R2.test | R Documentation |
This function computes Ramos-Onsins and Rozas's test of neutrality for a set of DNA sequences.
R2.test(x, B = 1000, theta = 1, plot = TRUE, quiet = FALSE, ...)
x |
a DNA matrix (object of class |
B |
the number of replicates used for the simulation procedure. |
theta |
the value of the |
plot |
a logical value specifying whether to plot the results
( |
quiet |
a logical value specifying whether to not display the
progress of the simulations. The default is |
... |
further arguments passed to |
a list with two elements: R2
the value of the test statistic
R_2
, and P.val
the associated P-value. If
B = 0
a single value, the test statistic, is returned
The simulation procedure probably needs to be tested and improved. However the results make sense so far.
Emmanuel Paradis
Ramos-Onsins, R. and Rozas, R. (2002) Statistical properties of new neutrality tests against population growth. Molecular Biology and Evolution, 19, 2092–2100.
Sano, J. and Tachida, G. (2005) Gene genealogy and properties of test statistics of neutrality under population growth. Genetics, 169, 1687–1697.
read.dna
, dist.dna
data(woodmouse)
R2.test(woodmouse, quiet = TRUE)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.