A=c(483,1101) B=c(477,1121) EX <- matrix(c(A,B),ncol=2) EX fit <- chisq.test(EX,correct=FALSE) fit fit$expected fit$residual