R/Chisqstat3.R

Defines functions Chisqstat3

Documented in Chisqstat3

Chisqstat3 <-
function(U, V, counts)

{

k = length(U)

n = sum(counts)

N = apply(counts, 2, sum)


if(N[1] == 0) 
{
chisq = t(U[1:(k-1)]) %*% solve(V[1:(k-1), 1:(k-1)]) %*% U[1:(k-1)] / n
df = k - 1
}
else
{
equal = 1

ratios = t(apply(counts, 1, "/", N))
r = 2

while(equal == 1 && r <= k)

{

if(any(ratios[r,] != ratios[1,]))

equal = 0

r = r + 1

}

if(equal == 1)

{

chisq = t(U[1:(k-1)]) %*% solve(V[1:(k-1), 1:(k-1)]) %*% U[1:(k-1)] / n
df = k - 1

}

else

{

chisq = t(U) %*% solve(V) %*% U / n

df = k

}
}

c(chisq, df)

}

Try the glrt package in your browser

Any scripts or data that you put into this service are public.

glrt documentation built on May 2, 2019, 8:31 a.m.