icu=c(3,36,29,30)
died=c(0,6,5,33)
oth=c(84,1069,228,286)
IJ <- matrix(c(icu,died,oth),ncol=3)
IJ
fit <- chisq.test(IJ,correct=FALSE)
fit
fit$expected
fit$residual