This contains the code from Jiang et al (2021), in package form so it is easier for me to compare in simulations. You can see the results of these comparisons in Gerard (2022).
The original repo that contains the code from Jiang et al (2021) is located at https://github.com/LiboJiang/DoubleReduction.
The main functions are
main_multi()
: Fit main_p2()
across many loci.main_p2()
: Calculate maximum likelihood esimates of double
reduction and allele frequency in tetraploids using all genotypes.main_p3()
: Calculate maximum likelihood esimates of double
reduction and allele frequency in tetraploids using just
heterozygous/homozygous.sim_p()
: Simulate genotypes at equilibrium in tetraploids.
However, there is a mistake in Jiang et al (2021) and these
genotypes are not actually simulated at equilibrium. See
Gerard (2022) for details.sim_p3()
: Simulate heterozygous/homozygous at equilibrium in
tetraploids. However, there is a mistake in Jiang et al (2021) and
these are not actually simulated at equilibrium. See Gerard (2022)
for details.David Gerard’s work was supported by the National Science Foundation under Grant No. 2132247. The opinions, findings, and conclusions or recommendations expressed are those of the author and do not necessarily reflect the views of the National Science Foundation.
The work of Jiang et al (2021) was supported by grant 201404102 from the State Administration of Forestry of China and grant 32070601 from the National Natural Science Foundation of China.
Please note that the phwelike project is released with a Contributor Code of Conduct. By contributing to this project, you agree to abide by its terms.
Gerard D (2022). “Double reduction estimation and equilibrium tests in natural autopolyploid populations.” Biometrics In press. doi:10.1111/biom.13722.
Jiang, Libo, Xiangyu Ren, and Rongling Wu. 2021. “Computational Characterization of Double Reduction in Autotetraploid Natural Populations.” The Plant Journal 105 (6): 1703–9. doi:10.1111/tpj.15126
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.