> music.ls <- aov(prod ~ row + col + trt)
> anova(music.ls)
Analysis of Variance Table

Response: prod
          Df Sum Sq Mean Sq F value    Pr(>F)    
row        4   82.0    20.5  1.3057 0.3225681    
col        4  477.2   119.3  7.5987 0.0027266 ** 
trt        4  664.4   166.1 10.5796 0.0006597 ***
Residuals 12  188.4    15.7                      
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 
> 
> TukeyHSD(music.ls, "trt")
  Tukey multiple comparisons of means
    95% family-wise confidence level

Fit: aov(formula = prod ~ row + col + trt)

$trt
    diff         lwr        upr     p adj
2-1 15.2   7.2123177 23.1876823 0.0004403
3-1  8.2   0.2123177 16.1876823 0.0432844
4-1 12.4   4.4123177 20.3876823 0.0025340
5-1 10.2   2.2123177 18.1876823 0.0110146
3-2 -7.0 -14.9876823  0.9876823 0.0966985
4-2 -2.8 -10.7876823  5.1876823 0.7947419
5-2 -5.0 -12.9876823  2.9876823 0.3242324
4-3  4.2  -3.7876823 12.1876823 0.4820666
5-3  2.0  -5.9876823  9.9876823 0.9262383
5-4 -2.2 -10.1876823  5.7876823 0.8996466