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]]
______________________________________________
[email protected] 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.
