# NOT RUN {
mysystem = createInSilicoSystem(G = 5, PC.p = 1, PC.TC.p = 1, ploidy = 1)
mysystem = addComplex(mysystem, c(1, 2))
mypop = createInSilicoPopulation(1, mysystem)
sim = simulateInSilicoSystem(mysystem, mypop, 100)
head(sim$Simulation)
mergedComplex = mergeComplexesAbundance(sim$Simulation)
head(mergedComplex)
# }
Run the code above in your browser using DataLab