# summary(data[1:3]) .Table <- xtabs(~row name+column name, data=data) .Table # Chi-squared test .Test <- chisq.test(.Table, correct=TRUE) .Test res <- chisq.test(.Table, correct=TRUE) # Expected count res$expected # Standardized Residual res$residuals # Adjusted Residual res$stdres # Residual analysis (Haberman 1973) pnorm(abs(res$stdres), lower.tail = FALSE) * 2 # Bonferroni correction p.value.matrix <- pnorm(abs(res$stdres), lower.tail=FALSE)*2 matrix(p.adjust(p.value.matrix,"bonferronih),row number,column number) # Holm correction p.value.matrix <- pnorm(abs(res$stdres), lower.tail=FALSE)*2 matrix(p.adjust(p.value.matrix,"holm"),row number,column number) # BH correction p.value.matrix <- pnorm(abs(res$stdres), lower.tail=FALSE)*2 matrix(p.adjust(p.value.matrix,"BH"),row number,column number)