Description Usage Arguments Details Value References See Also Examples
This function computes the MQLS proposed by Thormton. see the details T. Thornton et al. 2004.
1 | Mqls.test1(data, lstk, fam.id, fl, k)
|
data |
Data frame. |
lstk |
List of the family kinship matrices |
fam.id |
Family identifying |
fl |
The formula |
k |
The disease prevalence. By default, it's equal to zero. A second statistic is evaluated while optimizing the parameter k by Davies method. see D. Davies et al. (1987) |
The fonction implantes the MQLS published by T. Thornton et al (2004). The function uses the Davies method D. Davies et al. (1987) to evaluate a second statistic while optimizing the parameter k. A Software in C++ was also developed for the test. see https://galton.uchicago.edu/~mcpeek/software/MQLS/index.html.
Return a list containing
Rstt |
The result matrix containing the statistic of the test, the number of freedom degree, the p-value computed with the parameter k and the corrected p-value computed by the Davies's method. |
Rstd |
The result matrix containing the data description, the number of stratums, the number of freedom degree, the empiric allelic frequency, the estimated of the allelic frequency computed by BLUE method and the allelic frequency variance. |
sg |
The number of kinship matrix that are singulars. |
fMqls |
The function of the test statistic which depends of the nuisance parameter k. |
Timothy Thornton and Mary sara McPeek Case-Control Association Testing with Related Individuals: A More Powerful Quasi-likelihood Score Test American journal of human genetics vol 81 August 2007.
Robert B. Davies. Hypothesis testing when a nuisance parameter is present only under the alternative Biometrika (1987), 74, 1, pp.33-43.
To See Also as Mqls.test2
1 2 3 4 5 | data(DataGPh)
fl=formula(X2~Y2)
lst_mkinshp2=kinshAprior("fam.id","subject.id","father.id","mother.id","sex",DataGPh)
lR<-Mqls.test1(DataGPh,lst_mkinshp2,"fam.id",fl,0.008)
lR
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.