calculateCPProbabilities: Calculate the changepoint probabilities.
In EDISON: Network Reconstruction and Changepoint Detection

Description Usage Arguments Value Author(s) Examples

View source: R/calculateCPProbabilities.R

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

1	calculateCPProbabilities(network.samples)

network.samples

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

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.

Frank Dondelinger

# 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)