# 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)