discoal_sim: discoal_sim function
In deponent-verb/popgen.stats: Popgen summary statistics calculator

Description Usage Arguments Value Examples

View source: R/discoal_sim.R

discoal_sim function

1 2	discoal_sim(mu, recomb_rate, Ne, nSites, samplesize, s = 0, discoal_path, fix_generation, seed, sweep)

`mu`	mutation rate per base per generation. Warning! Check this again later in the discoal manual!
`recomb_rate`	The recombination rate per base per generation
`Ne`	The effective population size
`nSites`	The number of bases to simulate per genome
`samplesize`	Number of samples to take from the population
`s`	selection coefficient for the selected mutation
`discoal_path`	path to your discoal program
`fix_generation`	number of generations ago when the selected mutation was fixed
`seed`	vector of 2 numbers used for the simulations
`sweep`	the kind of selective sweep. Input "hard" or "neutral"

an object of class sim_obj. Here are the features. cmd is the command. Seeds: the seeds used in the discoal simulation. num_seg: number of segregating sites in the sampled population. pos: vector of the positions of every seg site (infinite sites model) sweep: the kind of selective sweep modelled. s: the selection coefficient

1	sim<-discoal_sim(mu=mu,recomb_rate=recomb_rate,Ne=Ne,nSites=nSites,samplesize=samplesize,s=s,discoal_path=discoal_path,fix_generation=fix,sweep=sweep)