rm(list = ls())
laz=matrix(c(
1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,
0,0,0,0,0,0,1,1,1,1,1,1,0,0,0,0,
0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,
1,0,0,0,0,0,1,0,0,0,0,0,1,0,0,0,
0,1,1,0,0,0,0,1,1,0,0,0,0,1,0,0,
0,0,0,1,0,0,0,0,0,1,1,0,0,0,1,0,
0,0,0,0,1,1,0,0,0,0,0,1,0,0,0,1),16,7)
laz

# fit indep/homogen model to table
fit <- chisq.test(laz,correct=FALSE)
fit
fit$expected
fit$residual
t(fit$residual)%*%fit$residual
t(laz)%*%laz
XtX=t(laz)%*%laz
fitt = chisq.test(XtX[1:3,4:7],correct=F)
fitt$residual