calculateCPProbabilities {EDISON} | R Documentation |
This function calculates the marginal changepoint probabilities from the changepoint samples taken during the MCMC simulation.
calculateCPProbabilities(network.samples)
network.samples |
List of network and changepoint samples collected
during the MCMC simulation by |
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)