nMelN: integer, number of Apis mellifera mellifera North individuals to simulate
nMelS: integer, number of Apis mellifera mellifera South individuals to simulate
nCar: integer, number of Apis mellifera carnica individuals to simulate
nLig: integer, number of Apis mellifera ligustica individuals to simulate
Ne: integer, effective size of the simulated population. Currently set to 170,000, according to Wallberg et al., 2014. Would discourage you to change it since it is linked to the parameters of the demographic model we use for the simulation. However, there might be some edge cases when using a different Ne is necessary, but proceed with caution.
ploidy: integer, the ploidy of the individuals
nChr: integer, number of chromosomes to simulate
nSegSites: integer, number of segregating sites to keep per chromosome
nBp: integer, base pair length of chromosome
genLen: numeric, genetic length of chromosome in Morgans
mutRate: numeric, per base pair mutation rate
recRate: numeric, per base pair recombination rate
nThreads: integer, if OpenMP is available, this will allow for simulating chromosomes in parallel. If NULL, the number of threads is automatically detected
Returns
MapPop-class
References
Wallberg, A., Bunikis, I., Pettersson, O.V. et al. A hybrid de novo genome assembly of the honeybee, Apis mellifera, with chromosome-length scaffolds. 2019, BMC Genomics 20:275. tools:::Rd_expr_doi("10.1186/s12864-019-5642-0")
Beye M, Gattermeier I, Hasselmann M, et al. Exceptionally high levels of recombination across the honey bee genome. 2006, Genome Res 16(11):1339-1344. tools:::Rd_expr_doi("10.1101/gr.5680406")
Wallberg, A., Han, F., Wellhagen, G. et al. A worldwide survey of genome sequence variation provides insight into the evolutionary history of the honeybee Apis mellifera. 2014, Nat Genet 46:1081–1088. tools:::Rd_expr_doi("10.1038/ng.3077")
Yang S, Wang L, Huang J, Zhang X, Yuan Y, Chen JQ, Hurst LD, Tian D. Parent-progeny sequencing indicates higher mutation rates in heterozygotes. 2015, Nature 523(7561):463-7. tools:::Rd_expr_doi("10.1038/nature14649") .
See Also
Due to the computational time and resources required to run this function, we do not include an example here, but we demonstrate its use in the Honeybee biology vignette.