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)