2 x 3 table

classic Classic list List threaded Threaded
1 message Options
Reply | Threaded
Open this post in threaded view
|

2 x 3 table

Jim Silverton
I am trying to find the likelihood ratio for a 2 x3 table (actually SNP
data). I do not want the log likelihood. Is this the correct R code for the
following table? Row totals are fixed at 100 each (case control situation).
aa     Aa     AA   Total
10     20    70     100
20     40    50     100
30     60    120


So  the row totals are fixed. Is this the correct code3?

Likelihood.ratio =
dmulitnom(100, c(0.1, 0.2, 0.7))*dmultinom(100, c(0.2, 0.4, 0.5))
divided by
dmulitnom(100, c(0.33, 0.33, 0.34))*dmultinom(100, c(0.33, 0.33, 0.34))

What about if the null came from the HWE proportions how will this work for
the case control SNP situation?

--
Thanks,
Jim.

        [[alternative HTML version deleted]]

______________________________________________
[hidden email] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.