E=c(21,3,7) G=c(11,2,1) F=c(4,2,1) IJ <- matrix(c(E,G,F),ncol=3) IJ fit <- chisq.test(IJ,correct=FALSE) fit fit$expected fit$residual