sexratio | R Documentation |
This function performs numerical analysis of a frequency dependent selection model of a hypothetical diploid sexually reproducing population in which sex is determined by the genotype at a biallelic genetic locus. Genotype AA are male, genotype aa are female, and genotype Aa might be male or female with probabilities that can be specified by the user. (Users may find, for instance, that setting sex.Aa=c(1,0)
will result in evolution towards an XY sex determination system; whereas sex.Aa=c(0,1)
will evolve towards a ZW system.)
sexratio(p0=0.01, time=40, show="p", pause=0, sex.Aa=c(0.5,0.5))
p0 |
Starting frequency for the A allele. Individuals with AA genotypes are male, while individuals with Aa genotypes are male or female with probability given by |
time |
Number of generations to run the analysis. |
show |
Two different options for plotting. |
pause |
Pause between generations. |
sex.Aa |
Probability that individuals with genotype Aa are male or female, respectively. |
The function creates one of four possible plots, depending on the value of show
. Numerical analysis of this model shows how frequency dependent selection should favor alleles that tend to produce the rarer sex in the population.
The function invisibly returns an object of class "sexratio"
that can be printed or re-plotted by the user. (See examples.)
Liam Revell liam.revell@umb.edu
freqdep
sexratio() sexratio(p0=0.001,show="sex-ratio") sexratio(p0=0.001,show="fitness") object<-sexratio(p0=0.001,sex.Aa=c(0.9,0.1), time=20) print(object) par(mfrow=c(2,1)) plot(object,lwd=4,type="s",show="sex-ratio") plot(object,lwd=4,type="s",show="genotypes") par(mfrow=c(1,1))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.