### Independence test for a 4 x 2 matrix
data(Anderson)
fisher.test(Anderson)
Anderson = Anderson[(c("Heimlich", "Bloom", "Dougal", "Cobblestone")),]
PT = pairwiseNominalMatrix(Anderson,
fisher = TRUE,
gtest = FALSE,
chisq = FALSE)$Adjusted
PT
library(multcompView)
multcompLetters(PT)
Run the code above in your browser using DataLab