calculateCPProbabilities {EDISON}

Calculate the changepoint probabilities.

Description

This function calculates the marginal changepoint probabilities from the changepoint samples taken during the MCMC simulation.

Usage

calculateCPProbabilities(network.samples)

Arguments

network.samples

List of network and changepoint samples collected during the MCMC simulation by EDISON.run and runDBN.

Value

Returns a matrix of dimension NumNodes by NumTimePoints, where each entry contains the marginal posterior probability of a changepoint for that node at that timepoint.

Author(s)

Frank Dondelinger

Examples

# Generate random gene network and simulate data from it
dataset = simulateNetwork()

# Run MCMC simulation to infer networks and changepoint locations
result = EDISON.run(dataset$sim_data, num.iter=500)

# Calculate posterior probabilities of changepoints
cps = calculateCPProbabilities(result)