row=c("1-10","11-19","20") rownames=row colnames = c("AfAm","Asian","Hisp","White") z = read.table("smokepreg.txt",sep=",",row.names=rownames, col.names=colnames) t = chisq.test(z) t$observed t$expected t$observed - t$expected t$residuals t$residuals^2 drows=c("Y","N") dcols = c("oldest","middle","youngest","only") z2 = read.table("delinq.txt",sep=",",row.names=drows, col.names=dcols) t2 = chisq.test(z2) t2$observed t2$expected t2$observed - t2$expected t2$residuals t2$residuals^2