mtdt | R Documentation |
Transmission/disequilibrium test of a multiallelic marker
mtdt(x, n.sim = 0)
x |
the data table. |
n.sim |
the number of simulations. |
This function calculates transmission-disequilibrium statistics involving multiallelic marker. Inside the function are tril and triu used to obtain lower and upper triangular matrices.
It returned list contains the following components:
SE Spielman-Ewens Chi-square from the observed data.
ST Stuart or score Statistic from the observed data.
pSE the simulated p value.
sSE standard error of the simulated p value.
pST the simulated p value.
sST standard error of the simulated p value.
Mike Miller, Jing Hua Zhao
miller97gap
\insertRefsham95gap
\insertRefspielman96gap
\insertRefzhao99gap
bt
## Not run:
x <- matrix(c(0,0, 0, 2, 0,0, 0, 0, 0, 0, 0, 0,
0,0, 1, 3, 0,0, 0, 2, 3, 0, 0, 0,
2,3,26,35, 7,0, 2,10,11, 3, 4, 1,
2,3,22,26, 6,2, 4, 4,10, 2, 2, 0,
0,1, 7,10, 2,0, 0, 2, 2, 1, 1, 0,
0,0, 1, 4, 0,1, 0, 1, 0, 0, 0, 0,
0,2, 5, 4, 1,1, 0, 0, 0, 2, 0, 0,
0,0, 2, 6, 1,0, 2, 0, 2, 0, 0, 0,
0,3, 6,19, 6,0, 0, 2, 5, 3, 0, 0,
0,0, 3, 1, 1,0, 0, 0, 1, 0, 0, 0,
0,0, 0, 2, 0,0, 0, 0, 0, 0, 0, 0,
0,0, 1, 0, 0,0, 0, 0, 0, 0, 0, 0),nrow=12)
# See note to bt for the score test obtained by SAS
mtdt(x)
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.