Dataset <- 
  read.table("/Users/minaka/TeachingMaterials/R-statistics/Win/R/data/Box1_R.txt",
             header=TRUE, sep="", na.strings="NA", dec=".", strip.white=TRUE)
attach(Dataset)
pairwise.t.test(DATA, TRT, p.adj="none")

## 
##  Pairwise comparisons using t tests with pooled SD 
## 
## data:  DATA and TRT 
## 
##     AZO     Con     DB      DDT     DK      DM1    
## Con 0.00124 -       -       -       -       -      
## DB  0.14214 0.03874 -       -       -       -      
## DDT 0.06508 1.2e-05 0.00228 -       -       -      
## DK  0.05269 0.10842 0.60284 0.00065 -       -      
## DM1 0.99547 0.00125 0.14358 0.06435 0.05330 -      
## DM2 0.01962 3.3e-06 0.00057 0.56811 0.00016 0.01937
## 
## P value adjustment method: none

pairwise.t.test(DATA, TRT, p.adj="bonferroni")

## 
##  Pairwise comparisons using t tests with pooled SD 
## 
## data:  DATA and TRT 
## 
##     AZO     Con     DB      DDT     DK      DM1    
## Con 0.02596 -       -       -       -       -      
## DB  1.00000 0.81348 -       -       -       -      
## DDT 1.00000 0.00026 0.04786 -       -       -      
## DK  1.00000 1.00000 1.00000 0.01364 -       -      
## DM1 1.00000 0.02632 1.00000 1.00000 1.00000 -      
## DM2 0.41194 7e-05   0.01206 1.00000 0.00341 0.40684
## 
## P value adjustment method: bonferroni

pairwise.t.test(DATA, TRT, p.adj="holm")

## 
##  Pairwise comparisons using t tests with pooled SD 
## 
## data:  DATA and TRT 
## 
##     AZO     Con     DB      DDT     DK      DM1    
## Con 0.01978 -       -       -       -       -      
## DB  0.71072 0.42611 -       -       -       -      
## DDT 0.52695 0.00025 0.03191 -       -       -      
## DK  0.52695 0.65051 1.00000 0.01104 -       -      
## DM1 1.00000 0.01978 0.71072 0.52695 0.52695 -      
## DM2 0.25185 7e-05   0.01033 1.00000 0.00309 0.25185
## 
## P value adjustment method: holm