# 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)
Run the code above in your browser using DataLab