Displaying 1 result from an estimated 1 matches for "meansqresid".
2003 Apr 22
4
fisher exact vs. simulated chi-square
...on's Chi-squared test
data: ta
X-squared = 9.6976, df = 4, p-value = 0.04584
Warning message:
Chi-squared approximation may be incorrect in: chisq.test(ta)
# simulate values by hand, based on r2dtable example
> expected <- outer(rowSums(ta), colSums(ta), "*") / sum(ta)
> meanSqResid <- function(x) mean((x - expected) ^ 2 / expected)
> sum(sapply(r2dtable(100000, rowSums(ta), colSums(ta)), meanSqResid)
>= meanSqResid(ta))/ 100000
[1] 0.03939
# is similar to
> sum(sapply(r2dtable(100000, rowSums(ta), colSums(ta)),
function(x) { chisq.test(x)$sta...