Open Thinking2018 opened 5 years ago
I use coala package to simulate coalescent data
with fixed number of markers
The code is as follows m = 100 locus.n = 500 model <- coal_model(m, 0) + locus_averaged(locus.n,10)+ feat_mutation(1,fixed_number = T) + sumstat_seg_sites()
sumstats <- simulate(model,seed = 124+i) sites <- sumstats$seg_sites X <- sapply(sites,as.matrix)
When I package the above code into a function use `parSapply' to simulate more data matrix X, it shows that error message
ms simulation failed
I got a similar error. Perhaps you have solved it, but if not, adding the feature feat_pop_merge() may help.
I use coala package to simulate coalescent data
with fixed number of markers
The code is as follows m = 100 locus.n = 500 model <- coal_model(m, 0) + locus_averaged(locus.n,10)+ feat_mutation(1,fixed_number = T) + sumstat_seg_sites()
sumstats <- simulate(model,seed = 124+i) sites <- sumstats$seg_sites X <- sapply(sites,as.matrix)
When I package the above code into a function use `parSapply' to simulate more data matrix X, it shows that error message
ms simulation failed