# NOT RUN {
mysystem = createInSilicoSystem(G = 5, PC.p = 1, PC.PTM.p = 0.9, regcomplexes = "none", ploidy = 1)
mypop = createInSilicoPopulation(1, mysystem)
sim = simulateInSilicoSystem(mysystem, mypop, 100)
head(sim$Simulation)
mergedPTM = mergePTMAbundance(sim$Simulation)
head(mergedPTM)
# }
Run the code above in your browser using DataLab