w=matrix(c(21,17,22,16,23,37,34,26,11),nrow=3)
rownames(w)=c("34 or younger","35-54","55 or older")
colnames(w)=c("worse","same","better")
chisq.test(w)
O=chisq.test(w)$observed
E=chisq.test(w)$expected
E
((O-E)^2)/E