Install necessary packages and load documents
require(agricolae)
## Loading required package: agricolae
Cabernet Sauvignon
setwd("C:/Users/19419/OneDrive - Washington State University (email.wsu.edu)/WSU/SMOKE/Stone Tree Bin Trials")
master <- read.csv("cabernet sauvignon gc.csv", header = T)
master$Location <- as.character(master$Location)
master$Wine.Treatment <- as.character(master$Wine.Treatment)
Guaiacol
By Location
guaiacol.aov <- aov(Guiacol~Location, data = master)
summary(guaiacol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 0.7 0.67 0.016 0.901
## Residuals 46 1979.5 43.03
By Wine Treatment
guaiacol.aov <- aov(Guiacol~Wine.Treatment, data = master)
summary(guaiacol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Wine.Treatment 3 1640.4 546.8 70.81 <2e-16 ***
## Residuals 44 339.8 7.7
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
guaiacol.tukey <- HSD.test(guaiacol.aov, c("Wine.Treatment"), group = T, console = T)
##
## Study: guaiacol.aov ~ c("Wine.Treatment")
##
## HSD Test for Guiacol
##
## Mean Square Error: 7.722331
##
## Wine.Treatment, means
##
## Guiacol std r Min Max
## bucket 8.34000 3.122525 12 4.55 12.72
## free 11.08833 1.236166 12 9.50 13.23
## press1 19.99250 3.708989 12 12.76 26.07
## press2 22.26833 2.419597 12 18.96 27.31
##
## Alpha: 0.05 ; DF Error: 44
## Critical Value of Studentized Range: 3.775958
##
## Minimun Significant Difference: 3.02908
##
## Treatments with the same letter are not significantly different.
##
## Guiacol groups
## press2 22.26833 a
## press1 19.99250 a
## free 11.08833 b
## bucket 8.34000 b
2-Way Interaction
guaiacol.aov <- aov(Guiacol~Location*Wine.Treatment, data = master)
summary(guaiacol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 0.7 0.7 0.092 0.764
## Wine.Treatment 3 1640.4 546.8 74.766 <2e-16 ***
## Location:Wine.Treatment 3 46.6 15.5 2.123 0.112
## Residuals 40 292.5 7.3
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Creosol
By Location
creosol.aov <- aov(Cresol~Location, data = master)
summary(creosol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 1.02 1.0150 1.123 0.295
## Residuals 46 41.56 0.9034
By Wine Treatment
creosol.aov <- aov(Cresol~Wine.Treatment, data = master)
summary(creosol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Wine.Treatment 3 31.13 10.38 39.89 1.3e-12 ***
## Residuals 44 11.44 0.26
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
creosol.tukey <- HSD.test(creosol.aov, c("Wine.Treatment"), group = T, console = T)
##
## Study: creosol.aov ~ c("Wine.Treatment")
##
## HSD Test for Cresol
##
## Mean Square Error: 0.2600989
##
## Wine.Treatment, means
##
## Cresol std r Min Max
## bucket 2.307500 0.9226259 12 1.04 3.65
## free 3.268333 0.3302845 12 2.89 4.13
## press1 1.306667 0.1769095 12 1.03 1.61
## press2 1.349167 0.2208438 12 0.96 1.63
##
## Alpha: 0.05 ; DF Error: 44
## Critical Value of Studentized Range: 3.775958
##
## Minimun Significant Difference: 0.5559117
##
## Treatments with the same letter are not significantly different.
##
## Cresol groups
## free 3.268333 a
## bucket 2.307500 b
## press2 1.349167 c
## press1 1.306667 c
2-Way Interaction
creosol.aov <- aov(Cresol~Location*Wine.Treatment, data = master)
summary(creosol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 1.015 1.015 4.133 0.0487 *
## Wine.Treatment 3 31.129 10.376 42.250 1.8e-12 ***
## Location:Wine.Treatment 3 0.606 0.202 0.822 0.4895
## Residuals 40 9.824 0.246
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
o-Cresol
By Location
ocresol.aov <- aov(oCresol~Location, data = master)
summary(ocresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 0.837 0.8374 1.341 0.253
## Residuals 46 28.724 0.6244
By Wine Treatment
ocresol.aov <- aov(oCresol~Wine.Treatment, data = master)
summary(ocresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Wine.Treatment 3 26.946 8.982 151.1 <2e-16 ***
## Residuals 44 2.616 0.059
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
ocresol.tukey <- HSD.test(ocresol.aov, c("Wine.Treatment"), group = T, console = T)
##
## Study: ocresol.aov ~ c("Wine.Treatment")
##
## HSD Test for oCresol
##
## Mean Square Error: 0.05945568
##
## Wine.Treatment, means
##
## oCresol std r Min Max
## bucket 2.908333 0.3547556 12 2.28 3.73
## free 3.521667 0.2126813 12 3.10 3.82
## press1 1.764167 0.1980568 12 1.48 2.04
## press2 1.797500 0.1658655 12 1.54 2.18
##
## Alpha: 0.05 ; DF Error: 44
## Critical Value of Studentized Range: 3.775958
##
## Minimun Significant Difference: 0.2657867
##
## Treatments with the same letter are not significantly different.
##
## oCresol groups
## free 3.521667 a
## bucket 2.908333 b
## press2 1.797500 c
## press1 1.764167 c
2-Way Interaction
ocresol.aov <- aov(oCresol~Location*Wine.Treatment, data = master)
summary(ocresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 0.837 0.837 19.461 7.55e-05 ***
## Wine.Treatment 3 26.946 8.982 208.736 < 2e-16 ***
## Location:Wine.Treatment 3 0.057 0.019 0.445 0.722
## Residuals 40 1.721 0.043
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
m-Cresol
By Location
mcresol.aov <- aov(mCresol~Location, data = master)
summary(mcresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 0.004 0.0037 0.008 0.927
## Residuals 46 19.925 0.4331
By Wine Treatment
mcresol.aov <- aov(mCresol~Wine.Treatment, data = master)
summary(mcresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Wine.Treatment 3 15.161 5.054 46.64 1.02e-13 ***
## Residuals 44 4.768 0.108
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mcresol.tukey <- HSD.test(mcresol.aov, c("Wine.Treatment"), group = T, console = T)
##
## Study: mcresol.aov ~ c("Wine.Treatment")
##
## HSD Test for mCresol
##
## Mean Square Error: 0.1083568
##
## Wine.Treatment, means
##
## mCresol std r Min Max
## bucket 2.703333 0.4716187 12 2.08 3.66
## free 3.120000 0.3155371 12 2.71 3.73
## press1 1.803333 0.2563851 12 1.39 2.22
## press2 1.851667 0.2137898 12 1.50 2.36
##
## Alpha: 0.05 ; DF Error: 44
## Critical Value of Studentized Range: 3.775958
##
## Minimun Significant Difference: 0.3588101
##
## Treatments with the same letter are not significantly different.
##
## mCresol groups
## free 3.120000 a
## bucket 2.703333 b
## press2 1.851667 c
## press1 1.803333 c
2-Way Interaction
mcresol.aov <- aov(mCresol~Location*Wine.Treatment, data = master)
summary(mcresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 0.004 0.004 0.035 0.853
## Wine.Treatment 3 15.161 5.054 47.908 2.62e-13 ***
## Location:Wine.Treatment 3 0.545 0.182 1.721 0.178
## Residuals 40 4.219 0.105
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
pCresol
By Location
pcresol.aov <- aov(pCresol~Location, data = master)
summary(pcresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 0.013 0.01333 0.046 0.831
## Residuals 46 13.266 0.28840
By Wine Treatment
pcresol.aov <- aov(pCresol~Wine.Treatment, data = master)
summary(pcresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Wine.Treatment 3 9.893 3.298 42.85 4.11e-13 ***
## Residuals 44 3.386 0.077
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
pcresol.tukey <- HSD.test(pcresol.aov, c("Wine.Treatment"), group = T, console = T)
##
## Study: pcresol.aov ~ c("Wine.Treatment")
##
## HSD Test for pCresol
##
## Mean Square Error: 0.07695871
##
## Wine.Treatment, means
##
## pCresol std r Min Max
## bucket 1.3108333 0.45299926 12 0.59 2.06
## free 1.1250000 0.26708868 12 0.79 1.69
## press1 0.3983333 0.07848953 12 0.27 0.51
## press2 0.2525000 0.15852301 12 0.08 0.59
##
## Alpha: 0.05 ; DF Error: 44
## Critical Value of Studentized Range: 3.775958
##
## Minimun Significant Difference: 0.3023886
##
## Treatments with the same letter are not significantly different.
##
## pCresol groups
## bucket 1.3108333 a
## free 1.1250000 a
## press1 0.3983333 b
## press2 0.2525000 b
2-Way Interaction
pcresol.aov <- aov(pCresol~Location*Wine.Treatment, data = master)
summary(pcresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 0.013 0.013 0.176 0.677
## Wine.Treatment 3 9.893 3.298 43.455 1.17e-12 ***
## Location:Wine.Treatment 3 0.337 0.112 1.481 0.234
## Residuals 40 3.036 0.076
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Eugenol
By Location
eugenol.aov <- aov(Eugenol~Location, data = master)
summary(eugenol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 0.864 0.8640 2.578 0.115
## Residuals 46 15.418 0.3352
By Wine Treatment
eugenol.aov <- aov(Eugenol~Wine.Treatment, data = master)
summary(eugenol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Wine.Treatment 3 13.430 4.477 69.05 <2e-16 ***
## Residuals 44 2.853 0.065
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
eugenol.tukey <- HSD.test(eugenol.aov, c("Wine.Treatment"), group = T, console = T)
##
## Study: eugenol.aov ~ c("Wine.Treatment")
##
## HSD Test for Eugenol
##
## Mean Square Error: 0.06483144
##
## Wine.Treatment, means
##
## Eugenol std r Min Max
## bucket 1.8600000 0.39691996 12 1.29 2.58
## free 1.4608333 0.25525240 12 1.16 1.94
## press1 0.5633333 0.07152029 12 0.40 0.65
## press2 0.7325000 0.17751440 12 0.46 1.12
##
## Alpha: 0.05 ; DF Error: 44
## Critical Value of Studentized Range: 3.775958
##
## Minimun Significant Difference: 0.2775424
##
## Treatments with the same letter are not significantly different.
##
## Eugenol groups
## bucket 1.8600000 a
## free 1.4608333 b
## press2 0.7325000 c
## press1 0.5633333 c
2-Way Interaction
eugenol.aov <- aov(Eugenol~Location*Wine.Treatment, data = master)
summary(eugenol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 0.864 0.864 18.83 9.45e-05 ***
## Wine.Treatment 3 13.430 4.477 97.54 < 2e-16 ***
## Location:Wine.Treatment 3 0.153 0.051 1.11 0.356
## Residuals 40 1.836 0.046
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Syrah
setwd("C:/Users/19419/OneDrive - Washington State University (email.wsu.edu)/WSU/SMOKE/Stone Tree Bin Trials")
master <- read.csv("syrah gc.csv", header = T)
master$Location <- as.character(master$Location)
master$Wine.Treatment <- as.character(master$Wine.Treatment)
Guaiacol
By Location
guaiacol.aov <- aov(Guiacol~Location, data = master)
summary(guaiacol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 1697 1697.4 66.47 1.76e-10 ***
## Residuals 46 1175 25.5
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
guaiacol.tukey <- HSD.test(guaiacol.aov, c("Location"), group = T, console = T)
##
## Study: guaiacol.aov ~ c("Location")
##
## HSD Test for Guiacol
##
## Mean Square Error: 25.53815
##
## Location, means
##
## Guiacol std r Min Max
## High 37.25125 4.289448 24 29.76 44.00
## Low 25.35792 5.716375 24 12.51 35.22
##
## Alpha: 0.05 ; DF Error: 46
## Critical Value of Studentized Range: 2.846664
##
## Minimun Significant Difference: 2.936469
##
## Treatments with the same letter are not significantly different.
##
## Guiacol groups
## High 37.25125 a
## Low 25.35792 b
By Wine Treatment
guaiacol.aov <- aov(Guiacol~Wine.Treatment, data = master)
summary(guaiacol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Wine.Treatment 3 418.4 139.47 2.501 0.0718 .
## Residuals 44 2453.8 55.77
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
2-Way Interaction
guaiacol.aov <- aov(Guiacol~Location*Wine.Treatment, data = master)
summary(guaiacol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 1697.4 1697.4 113.489 3.02e-13 ***
## Wine.Treatment 3 418.4 139.5 9.325 8.44e-05 ***
## Location:Wine.Treatment 3 158.1 52.7 3.523 0.0234 *
## Residuals 40 598.3 15.0
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
guaiacol.tukey <- HSD.test(guaiacol.aov, c("Wine.Treatment","Location"), group = T, console = T)
##
## Study: guaiacol.aov ~ c("Wine.Treatment", "Location")
##
## HSD Test for Guiacol
##
## Mean Square Error: 14.95667
##
## Wine.Treatment:Location, means
##
## Guiacol std r Min Max
## bucket:High 36.15667 5.536304 6 29.76 43.60
## bucket:Low 19.47333 4.614320 6 12.51 24.49
## free:High 36.04000 3.721204 6 30.89 39.82
## free:Low 21.94167 2.001444 6 19.72 24.40
## press1:High 38.67333 4.588711 6 34.06 44.00
## press1:Low 30.73667 3.431773 6 26.50 35.22
## press2:High 38.13500 3.538433 6 33.79 43.42
## press2:Low 29.28000 2.122216 6 27.43 32.47
##
## Alpha: 0.05 ; DF Error: 40
## Critical Value of Studentized Range: 4.520535
##
## Minimun Significant Difference: 7.137264
##
## Treatments with the same letter are not significantly different.
##
## Guiacol groups
## press1:High 38.67333 a
## press2:High 38.13500 a
## bucket:High 36.15667 ab
## free:High 36.04000 ab
## press1:Low 30.73667 b
## press2:Low 29.28000 b
## free:Low 21.94167 c
## bucket:Low 19.47333 c
Interaction Plot:
interaction.plot(master$Location,master$Wine.Treatment,master$Guiacol)
Creosol
By Location
creosol.aov <- aov(Cresol~Location, data = master)
summary(creosol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 149.35 149.35 79.19 1.45e-11 ***
## Residuals 46 86.76 1.89
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
creosol.tukey <- HSD.test(creosol.aov, c("Location"), group = T, console = T)
##
## Study: creosol.aov ~ c("Location")
##
## HSD Test for Cresol
##
## Mean Square Error: 1.885991
##
## Location, means
##
## Cresol std r Min Max
## High 6.195833 1.870598 24 4.03 11.24
## Low 2.667917 0.522344 24 1.51 3.62
##
## Alpha: 0.05 ; DF Error: 46
## Critical Value of Studentized Range: 2.846664
##
## Minimun Significant Difference: 0.7979955
##
## Treatments with the same letter are not significantly different.
##
## Cresol groups
## High 6.195833 a
## Low 2.667917 b
By Wine Treatment
creosol.aov <- aov(Cresol~Wine.Treatment, data = master)
summary(creosol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Wine.Treatment 3 25.92 8.639 1.808 0.16
## Residuals 44 210.19 4.777
2-Way Interaction
creosol.aov <- aov(Cresol~Location*Wine.Treatment, data = master)
summary(creosol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 149.35 149.35 149.341 4.4e-15 ***
## Wine.Treatment 3 25.92 8.64 8.638 0.000153 ***
## Location:Wine.Treatment 3 20.84 6.95 6.945 0.000717 ***
## Residuals 40 40.00 1.00
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
creosol.tukey <- HSD.test(creosol.aov, c("Location","Wine.Treatment"), group = T, console = T)
##
## Study: creosol.aov ~ c("Location", "Wine.Treatment")
##
## HSD Test for Cresol
##
## Mean Square Error: 1.000088
##
## Location:Wine.Treatment, means
##
## Cresol std r Min Max
## High:bucket 7.816667 2.2934050 6 5.50 11.24
## High:free 7.301667 1.0371387 6 5.94 8.35
## High:press1 4.778333 0.6615865 6 4.03 5.50
## High:press2 4.886667 0.4849192 6 4.29 5.43
## Low:bucket 2.490000 0.8171903 6 1.51 3.50
## Low:free 3.050000 0.2026820 6 2.70 3.21
## Low:press1 2.465000 0.2206128 6 2.20 2.75
## Low:press2 2.666667 0.4847130 6 2.34 3.62
##
## Alpha: 0.05 ; DF Error: 40
## Critical Value of Studentized Range: 4.520535
##
## Minimun Significant Difference: 1.845582
##
## Treatments with the same letter are not significantly different.
##
## Cresol groups
## High:bucket 7.816667 a
## High:free 7.301667 a
## High:press2 4.886667 b
## High:press1 4.778333 b
## Low:free 3.050000 bc
## Low:press2 2.666667 c
## Low:bucket 2.490000 c
## Low:press1 2.465000 c
Interaction Plot:
interaction.plot(master$Location,master$Wine.Treatment,master$Cresol)
o-Cresol
By Location
ocresol.aov <- aov(oCresol~Location, data = master)
summary(ocresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 44.56 44.56 58.42 9.98e-10 ***
## Residuals 46 35.09 0.76
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
ocresol.tukey <- HSD.test(ocresol.aov, c("Location"), group = T, console = T)
##
## Study: ocresol.aov ~ c("Location")
##
## HSD Test for oCresol
##
## Mean Square Error: 0.7628539
##
## Location, means
##
## oCresol std r Min Max
## High 3.992917 1.0968868 24 2.57 5.99
## Low 2.065833 0.5679323 24 1.31 3.09
##
## Alpha: 0.05 ; DF Error: 46
## Critical Value of Studentized Range: 2.846664
##
## Minimun Significant Difference: 0.5075178
##
## Treatments with the same letter are not significantly different.
##
## oCresol groups
## High 3.992917 a
## Low 2.065833 b
By Wine Treatment
ocresol.aov <- aov(oCresol~Wine.Treatment, data = master)
summary(ocresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Wine.Treatment 3 24.31 8.103 6.441 0.00103 **
## Residuals 44 55.35 1.258
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
ocresol.tukey <- HSD.test(ocresol.aov, c("Wine.Treatment"), group = T, console = T)
##
## Study: ocresol.aov ~ c("Wine.Treatment")
##
## HSD Test for oCresol
##
## Mean Square Error: 1.25789
##
## Wine.Treatment, means
##
## oCresol std r Min Max
## bucket 3.273333 1.2435604 12 1.94 5.58
## free 4.086667 1.3757995 12 2.57 5.99
## press1 2.407500 0.8808583 12 1.41 3.70
## press2 2.350000 0.9035385 12 1.31 3.73
##
## Alpha: 0.05 ; DF Error: 44
## Critical Value of Studentized Range: 3.775958
##
## Minimun Significant Difference: 1.222525
##
## Treatments with the same letter are not significantly different.
##
## oCresol groups
## free 4.086667 a
## bucket 3.273333 ab
## press1 2.407500 b
## press2 2.350000 b
2-Way Interaction
ocresol.aov <- aov(oCresol~Location*Wine.Treatment, data = master)
summary(ocresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 44.56 44.56 196.754 < 2e-16 ***
## Wine.Treatment 3 24.31 8.10 35.774 2.1e-11 ***
## Location:Wine.Treatment 3 1.72 0.57 2.537 0.0703 .
## Residuals 40 9.06 0.23
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
m-Cresol
By Location
mcresol.aov <- aov(mCresol~Location, data = master)
summary(mcresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 27.36 27.361 74.45 3.56e-11 ***
## Residuals 46 16.91 0.367
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mcresol.tukey <- HSD.test(mcresol.aov, c("Location"), group = T, console = T)
##
## Study: mcresol.aov ~ c("Location")
##
## HSD Test for mCresol
##
## Mean Square Error: 0.3674955
##
## Location, means
##
## mCresol std r Min Max
## High 2.812917 0.7868511 24 1.74 4.47
## Low 1.302917 0.3403768 24 0.80 1.92
##
## Alpha: 0.05 ; DF Error: 46
## Critical Value of Studentized Range: 2.846664
##
## Minimun Significant Difference: 0.3522546
##
## Treatments with the same letter are not significantly different.
##
## mCresol groups
## High 2.812917 a
## Low 1.302917 b
By Wine Treatment
mcresol.aov <- aov(mCresol~Wine.Treatment, data = master)
summary(mcresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Wine.Treatment 3 7.00 2.333 2.755 0.0536 .
## Residuals 44 37.27 0.847
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
2-Way Interaction
mcresol.aov <- aov(mCresol~Location*Wine.Treatment, data = master)
summary(mcresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 27.361 27.361 154.510 2.56e-15 ***
## Wine.Treatment 3 7.000 2.333 13.177 3.98e-06 ***
## Location:Wine.Treatment 3 2.821 0.940 5.311 0.00355 **
## Residuals 40 7.083 0.177
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mcresol.tukey <- HSD.test(mcresol.aov, c("Location","Wine.Treatment"), group = T, console = T)
##
## Study: mcresol.aov ~ c("Location", "Wine.Treatment")
##
## HSD Test for mCresol
##
## Mean Square Error: 0.1770833
##
## Location:Wine.Treatment, means
##
## mCresol std r Min Max
## High:bucket 3.461667 0.8547144 6 2.44 4.47
## High:free 3.351667 0.3877843 6 2.86 3.79
## High:press1 2.260000 0.4523273 6 1.74 2.88
## High:press2 2.178333 0.2430158 6 1.74 2.44
## Low:bucket 1.218333 0.2103727 6 0.80 1.38
## Low:free 1.703333 0.1910672 6 1.47 1.92
## Low:press1 1.143333 0.2859837 6 0.83 1.56
## Low:press2 1.146667 0.3309783 6 0.83 1.75
##
## Alpha: 0.05 ; DF Error: 40
## Critical Value of Studentized Range: 4.520535
##
## Minimun Significant Difference: 0.7766102
##
## Treatments with the same letter are not significantly different.
##
## mCresol groups
## High:bucket 3.461667 a
## High:free 3.351667 a
## High:press1 2.260000 b
## High:press2 2.178333 b
## Low:free 1.703333 bc
## Low:bucket 1.218333 c
## Low:press2 1.146667 c
## Low:press1 1.143333 c
Interaction Plot:
interaction.plot(master$Location,master$Wine.Treatment,master$mCresol)
pCresol
By Location
pcresol.aov <- aov(pCresol~Location, data = master)
summary(pcresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 1.200 1.2002 8.068 0.00669 **
## Residuals 46 6.843 0.1488
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
pcresol.tukey <- HSD.test(pcresol.aov, c("Location"), group = T, console = T)
##
## Study: pcresol.aov ~ c("Location")
##
## HSD Test for pCresol
##
## Mean Square Error: 0.1487644
##
## Location, means
##
## pCresol std r Min Max
## High 1.20000 0.4269711 24 0.58 2.04
## Low 0.88375 0.3394473 24 0.23 1.62
##
## Alpha: 0.05 ; DF Error: 46
## Critical Value of Studentized Range: 2.846664
##
## Minimun Significant Difference: 0.2241198
##
## Treatments with the same letter are not significantly different.
##
## pCresol groups
## High 1.20000 a
## Low 0.88375 b
By Wine Treatment
pcresol.aov <- aov(pCresol~Wine.Treatment, data = master)
summary(pcresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Wine.Treatment 3 4.277 1.4256 16.65 2.26e-07 ***
## Residuals 44 3.767 0.0856
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
pcresol.tukey <- HSD.test(pcresol.aov, c("Wine.Treatment"), group = T, console = T)
##
## Study: pcresol.aov ~ c("Wine.Treatment")
##
## HSD Test for pCresol
##
## Mean Square Error: 0.08560473
##
## Wine.Treatment, means
##
## pCresol std r Min Max
## bucket 1.2641667 0.3866278 12 0.73 2.04
## free 1.3875000 0.2772306 12 0.81 1.75
## press1 0.8725000 0.2671100 12 0.63 1.50
## press2 0.6433333 0.2115026 12 0.23 1.07
##
## Alpha: 0.05 ; DF Error: 44
## Critical Value of Studentized Range: 3.775958
##
## Minimun Significant Difference: 0.3189227
##
## Treatments with the same letter are not significantly different.
##
## pCresol groups
## free 1.3875000 a
## bucket 1.2641667 a
## press1 0.8725000 b
## press2 0.6433333 b
2-Way Interaction
pcresol.aov <- aov(pCresol~Location*Wine.Treatment, data = master)
summary(pcresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 1.200 1.2002 20.812 4.73e-05 ***
## Wine.Treatment 3 4.277 1.4256 24.720 3.26e-09 ***
## Location:Wine.Treatment 3 0.260 0.0866 1.501 0.229
## Residuals 40 2.307 0.0577
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Eugenol
By Location
eugenol.aov <- aov(Eugenol~Location, data = master)
summary(eugenol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 35.21 35.21 10.78 0.00196 **
## Residuals 46 150.19 3.26
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
eugenol.tukey <- HSD.test(eugenol.aov, c("Location"), group = T, console = T)
##
## Study: eugenol.aov ~ c("Location")
##
## HSD Test for Eugenol
##
## Mean Square Error: 3.264981
##
## Location, means
##
## Eugenol std r Min Max
## High 5.813750 2.150334 24 3.23 10.30
## Low 4.100833 1.380589 24 2.33 6.23
##
## Alpha: 0.05 ; DF Error: 46
## Critical Value of Studentized Range: 2.846664
##
## Minimun Significant Difference: 1.049956
##
## Treatments with the same letter are not significantly different.
##
## Eugenol groups
## High 5.813750 a
## Low 4.100833 b
By Wine Treatment
eugenol.aov <- aov(Eugenol~Wine.Treatment, data = master)
summary(eugenol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Wine.Treatment 3 117.93 39.31 25.64 9.54e-10 ***
## Residuals 44 67.46 1.53
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
eugenol.tukey <- HSD.test(eugenol.aov, c("Wine.Treatment"), group = T, console = T)
##
## Study: eugenol.aov ~ c("Wine.Treatment")
##
## HSD Test for Eugenol
##
## Mean Square Error: 1.533294
##
## Wine.Treatment, means
##
## Eugenol std r Min Max
## bucket 6.292500 2.0632592 12 3.66 10.30
## free 6.740000 1.0760534 12 5.24 8.05
## press1 3.454167 0.6236324 12 2.66 4.35
## press2 3.342500 0.5738724 12 2.33 4.12
##
## Alpha: 0.05 ; DF Error: 44
## Critical Value of Studentized Range: 3.775958
##
## Minimun Significant Difference: 1.349737
##
## Treatments with the same letter are not significantly different.
##
## Eugenol groups
## free 6.740000 a
## bucket 6.292500 a
## press1 3.454167 b
## press2 3.342500 b
2-Way Interaction
eugenol.aov <- aov(Eugenol~Location*Wine.Treatment, data = master)
summary(eugenol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 35.21 35.21 59.644 1.88e-09 ***
## Wine.Treatment 3 117.93 39.31 66.593 1.31e-15 ***
## Location:Wine.Treatment 3 8.64 2.88 4.881 0.00552 **
## Residuals 40 23.61 0.59
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
eugenol.tukey <- HSD.test(eugenol.aov, c("Location","Wine.Treatment"), group = T, console = T)
##
## Study: eugenol.aov ~ c("Location", "Wine.Treatment")
##
## HSD Test for Eugenol
##
## Mean Square Error: 0.5903154
##
## Location:Wine.Treatment, means
##
## Eugenol std r Min Max
## High:bucket 7.811667 1.6946081 6 6.10 10.30
## High:free 7.675000 0.4972625 6 6.95 8.05
## High:press1 3.950000 0.4449719 6 3.23 4.35
## High:press2 3.818333 0.2444109 6 3.45 4.12
## Low:bucket 4.773333 0.9771933 6 3.66 5.87
## Low:free 5.805000 0.4493885 6 5.24 6.23
## Low:press1 2.958333 0.2599551 6 2.66 3.31
## Low:press2 2.866667 0.3484059 6 2.33 3.21
##
## Alpha: 0.05 ; DF Error: 40
## Critical Value of Studentized Range: 4.520535
##
## Minimun Significant Difference: 1.417935
##
## Treatments with the same letter are not significantly different.
##
## Eugenol groups
## High:bucket 7.811667 a
## High:free 7.675000 a
## Low:free 5.805000 b
## Low:bucket 4.773333 bc
## High:press1 3.950000 cd
## High:press2 3.818333 cd
## Low:press1 2.958333 d
## Low:press2 2.866667 d
Interaction Plot:
interaction.plot(master$Location,master$Wine.Treatment,master$Eugenol)
Petit Verdot
setwd("C:/Users/19419/OneDrive - Washington State University (email.wsu.edu)/WSU/SMOKE/Stone Tree Bin Trials")
master <- read.csv("petite verdot gc.csv", header = T)
master$Location <- as.character(master$Location)
master$Wine.Treatment <- as.character(master$Wine.Treatment)
Guaiacol
By Location
guaiacol.aov <- aov(Guiacol~Location, data = master)
summary(guaiacol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 2686 2686.2 92.44 1.4e-12 ***
## Residuals 46 1337 29.1
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
guaiacol.tukey <- HSD.test(guaiacol.aov, c("Location"), group = T, console = T)
##
## Study: guaiacol.aov ~ c("Location")
##
## HSD Test for Guiacol
##
## Mean Square Error: 29.05913
##
## Location, means
##
## Guiacol std r Min Max
## High 48.05542 5.284538 24 35.01 56.58
## Low 33.09375 5.494717 24 22.82 42.05
##
## Alpha: 0.05 ; DF Error: 46
## Critical Value of Studentized Range: 2.846664
##
## Minimun Significant Difference: 3.132362
##
## Treatments with the same letter are not significantly different.
##
## Guiacol groups
## High 48.05542 a
## Low 33.09375 b
By Wine Treatment
guaiacol.aov <- aov(Guiacol~Wine.Treatment, data = master)
summary(guaiacol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Wine.Treatment 3 552 184.06 2.333 0.087 .
## Residuals 44 3471 78.88
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
2-Way Interaction
guaiacol.aov <- aov(Guiacol~Location*Wine.Treatment, data = master)
summary(guaiacol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 2686.2 2686.2 155.997 2.19e-15 ***
## Wine.Treatment 3 552.2 184.1 10.689 2.72e-05 ***
## Location:Wine.Treatment 3 95.8 31.9 1.854 0.153
## Residuals 40 688.8 17.2
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Creosol
By Location
creosol.aov <- aov(Cresol~Location, data = master)
summary(creosol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 278.0 277.97 45.25 2.34e-08 ***
## Residuals 46 282.6 6.14
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
creosol.tukey <- HSD.test(creosol.aov, c("Location"), group = T, console = T)
##
## Study: creosol.aov ~ c("Location")
##
## HSD Test for Cresol
##
## Mean Square Error: 6.142609
##
## Location, means
##
## Cresol std r Min Max
## High 13.241250 3.152189 24 8.83 18.23
## Low 8.428333 1.532618 24 6.65 12.03
##
## Alpha: 0.05 ; DF Error: 46
## Critical Value of Studentized Range: 2.846664
##
## Minimun Significant Difference: 1.440148
##
## Treatments with the same letter are not significantly different.
##
## Cresol groups
## High 13.241250 a
## Low 8.428333 b
By Wine Treatment
creosol.aov <- aov(Cresol~Wine.Treatment, data = master)
summary(creosol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Wine.Treatment 3 181.3 60.43 7.011 0.000591 ***
## Residuals 44 379.2 8.62
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
creosol.tukey <- HSD.test(creosol.aov, c("Wine.Treatment"), group = T, console = T)
##
## Study: creosol.aov ~ c("Wine.Treatment")
##
## HSD Test for Cresol
##
## Mean Square Error: 8.61927
##
## Wine.Treatment, means
##
## Cresol std r Min Max
## bucket 12.000000 3.795253 12 6.92 17.54
## free 13.423333 3.548232 12 8.97 18.23
## press1 8.998333 1.870133 12 6.65 12.43
## press2 8.917500 1.996443 12 6.72 11.80
##
## Alpha: 0.05 ; DF Error: 44
## Critical Value of Studentized Range: 3.775958
##
## Minimun Significant Difference: 3.200161
##
## Treatments with the same letter are not significantly different.
##
## Cresol groups
## free 13.423333 a
## bucket 12.000000 ab
## press1 8.998333 b
## press2 8.917500 b
2-Way Interaction
creosol.aov <- aov(Cresol~Location*Wine.Treatment, data = master)
summary(creosol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 277.97 277.97 152.330 3.21e-15 ***
## Wine.Treatment 3 181.28 60.43 33.115 6.32e-11 ***
## Location:Wine.Treatment 3 28.29 9.43 5.167 0.00411 **
## Residuals 40 72.99 1.82
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
creosol.tukey <- HSD.test(creosol.aov, c("Location","Wine.Treatment"), group = T, console = T)
##
## Study: creosol.aov ~ c("Location", "Wine.Treatment")
##
## HSD Test for Cresol
##
## Mean Square Error: 1.824785
##
## Location:Wine.Treatment, means
##
## Cresol std r Min Max
## High:bucket 15.070000 2.6327932 6 10.43 17.54
## High:free 16.690000 0.9578935 6 15.73 18.23
## High:press1 10.553333 1.2182063 6 9.25 12.43
## High:press2 10.651667 1.1708700 6 8.83 11.80
## Low:bucket 8.930000 1.4619576 6 6.92 10.63
## Low:free 10.156667 1.0815113 6 8.97 12.03
## Low:press1 7.443333 0.6379551 6 6.65 8.27
## Low:press2 7.183333 0.4244840 6 6.72 7.66
##
## Alpha: 0.05 ; DF Error: 40
## Critical Value of Studentized Range: 4.520535
##
## Minimun Significant Difference: 2.492988
##
## Treatments with the same letter are not significantly different.
##
## Cresol groups
## High:free 16.690000 a
## High:bucket 15.070000 a
## High:press2 10.651667 b
## High:press1 10.553333 b
## Low:free 10.156667 b
## Low:bucket 8.930000 bc
## Low:press1 7.443333 c
## Low:press2 7.183333 c
Interaction Plot:
interaction.plot(master$Location,master$Wine.Treatment,master$Cresol)
o-Cresol
By Location
ocresol.aov <- aov(oCresol~Location, data = master)
summary(ocresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 8.283 8.283 16.08 0.000221 ***
## Residuals 46 23.703 0.515
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
ocresol.tukey <- HSD.test(ocresol.aov, c("Location"), group = T, console = T)
##
## Study: ocresol.aov ~ c("Location")
##
## HSD Test for oCresol
##
## Mean Square Error: 0.5152752
##
## Location, means
##
## oCresol std r Min Max
## High 2.942083 0.8555699 24 1.61 4.15
## Low 2.111250 0.5463978 24 1.32 2.94
##
## Alpha: 0.05 ; DF Error: 46
## Critical Value of Studentized Range: 2.846664
##
## Minimun Significant Difference: 0.4171097
##
## Treatments with the same letter are not significantly different.
##
## oCresol groups
## High 2.942083 a
## Low 2.111250 b
By Wine Treatment
ocresol.aov <- aov(oCresol~Wine.Treatment, data = master)
summary(ocresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Wine.Treatment 3 18.24 6.079 19.46 3.53e-08 ***
## Residuals 44 13.75 0.312
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
ocresol.tukey <- HSD.test(ocresol.aov, c("Wine.Treatment"), group = T, console = T)
##
## Study: ocresol.aov ~ c("Wine.Treatment")
##
## HSD Test for oCresol
##
## Mean Square Error: 0.312458
##
## Wine.Treatment, means
##
## oCresol std r Min Max
## bucket 2.879167 0.5655160 12 2.18 3.67
## free 3.359167 0.6113838 12 2.56 4.15
## press1 1.948333 0.2665208 12 1.51 2.43
## press2 1.920000 0.6965630 12 1.32 3.69
##
## Alpha: 0.05 ; DF Error: 44
## Critical Value of Studentized Range: 3.775958
##
## Minimun Significant Difference: 0.6093016
##
## Treatments with the same letter are not significantly different.
##
## oCresol groups
## free 3.359167 a
## bucket 2.879167 a
## press1 1.948333 b
## press2 1.920000 b
2-Way Interaction
ocresol.aov <- aov(oCresol~Location*Wine.Treatment, data = master)
summary(ocresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 8.283 8.283 74.170 1.19e-10 ***
## Wine.Treatment 3 18.238 6.079 54.435 3.50e-14 ***
## Location:Wine.Treatment 3 0.998 0.333 2.977 0.0428 *
## Residuals 40 4.467 0.112
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
ocresol.tukey <- HSD.test(ocresol.aov, c("Location","Wine.Treatment"), group = T, console = T)
##
## Study: ocresol.aov ~ c("Location", "Wine.Treatment")
##
## HSD Test for oCresol
##
## Mean Square Error: 0.1116808
##
## Location:Wine.Treatment, means
##
## oCresol std r Min Max
## High:bucket 3.375000 0.2746452 6 2.98 3.67
## High:free 3.920000 0.2150349 6 3.61 4.15
## High:press1 2.128333 0.2142351 6 1.85 2.43
## High:press2 2.345000 0.7776310 6 1.61 3.69
## Low:bucket 2.383333 0.1952093 6 2.18 2.65
## Low:free 2.798333 0.1456594 6 2.56 2.94
## Low:press1 1.768333 0.1806008 6 1.51 1.99
## Low:press2 1.495000 0.1709678 6 1.32 1.75
##
## Alpha: 0.05 ; DF Error: 40
## Critical Value of Studentized Range: 4.520535
##
## Minimun Significant Difference: 0.6167421
##
## Treatments with the same letter are not significantly different.
##
## oCresol groups
## High:free 3.920000 a
## High:bucket 3.375000 ab
## Low:free 2.798333 bc
## Low:bucket 2.383333 cd
## High:press2 2.345000 cd
## High:press1 2.128333 d
## Low:press1 1.768333 de
## Low:press2 1.495000 e
Interaction Plot:
interaction.plot(master$Location,master$Wine.Treatment,master$oCresol)
m-Cresol
By Location
mcresol.aov <- aov(mCresol~Location, data = master)
summary(mcresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 8.96 8.961 6.168 0.0167 *
## Residuals 46 66.84 1.453
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mcresol.tukey <- HSD.test(mcresol.aov, c("Location"), group = T, console = T)
##
## Study: mcresol.aov ~ c("Location")
##
## HSD Test for mCresol
##
## Mean Square Error: 1.452937
##
## Location, means
##
## mCresol std r Min Max
## High 3.570833 1.5849096 24 -0.22 5.65
## Low 2.706667 0.6276434 24 1.83 3.79
##
## Alpha: 0.05 ; DF Error: 46
## Critical Value of Studentized Range: 2.846664
##
## Minimun Significant Difference: 0.7004128
##
## Treatments with the same letter are not significantly different.
##
## mCresol groups
## High 3.570833 a
## Low 2.706667 b
By Wine Treatment
mcresol.aov <- aov(mCresol~Wine.Treatment, data = master)
summary(mcresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Wine.Treatment 3 27.75 9.251 8.473 0.000149 ***
## Residuals 44 48.04 1.092
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mcresol.tukey <- HSD.test(mcresol.aov, c("Wine.Treatment"), group = T, console = T)
##
## Study: mcresol.aov ~ c("Wine.Treatment")
##
## HSD Test for mCresol
##
## Mean Square Error: 1.091869
##
## Wine.Treatment, means
##
## mCresol std r Min Max
## bucket 3.633333 0.9161712 12 2.39 5.30
## free 4.011667 0.7447615 12 2.66 5.25
## press1 2.024167 1.1390463 12 -0.22 3.54
## press2 2.885833 1.2946074 12 1.83 5.65
##
## Alpha: 0.05 ; DF Error: 44
## Critical Value of Studentized Range: 3.775958
##
## Minimun Significant Difference: 1.138995
##
## Treatments with the same letter are not significantly different.
##
## mCresol groups
## free 4.011667 a
## bucket 3.633333 a
## press2 2.885833 ab
## press1 2.024167 b
2-Way Interaction
mcresol.aov <- aov(mCresol~Location*Wine.Treatment, data = master)
summary(mcresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 8.96 8.961 11.270 0.00174 **
## Wine.Treatment 3 27.75 9.251 11.634 1.28e-05 ***
## Location:Wine.Treatment 3 7.27 2.424 3.049 0.03953 *
## Residuals 40 31.81 0.795
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mcresol.tukey <- HSD.test(mcresol.aov, c("Location","Wine.Treatment"), group = T, console = T)
##
## Study: mcresol.aov ~ c("Location", "Wine.Treatment")
##
## HSD Test for mCresol
##
## Mean Square Error: 0.7951875
##
## Location:Wine.Treatment, means
##
## mCresol std r Min Max
## High:bucket 4.335000 0.6943702 6 3.33 5.30
## High:free 4.558333 0.5788062 6 3.80 5.25
## High:press1 1.791667 1.6168539 6 -0.22 3.54
## High:press2 3.598333 1.5662620 6 1.83 5.65
## Low:bucket 2.931667 0.4276174 6 2.39 3.51
## Low:free 3.465000 0.4099146 6 2.66 3.79
## Low:press1 2.256667 0.3322449 6 1.83 2.65
## Low:press2 2.173333 0.1251666 6 2.01 2.34
##
## Alpha: 0.05 ; DF Error: 40
## Critical Value of Studentized Range: 4.520535
##
## Minimun Significant Difference: 1.645694
##
## Treatments with the same letter are not significantly different.
##
## mCresol groups
## High:free 4.558333 a
## High:bucket 4.335000 a
## High:press2 3.598333 ab
## Low:free 3.465000 ab
## Low:bucket 2.931667 abc
## Low:press1 2.256667 bc
## Low:press2 2.173333 bc
## High:press1 1.791667 c
Interaction Plot:
interaction.plot(master$Location,master$Wine.Treatment,master$mCresol)
pCresol
By Location
pcresol.aov <- aov(pCresol~Location, data = master)
summary(pcresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 0.064 0.0638 0.126 0.725
## Residuals 46 23.341 0.5074
By Wine Treatment
pcresol.aov <- aov(pCresol~Wine.Treatment, data = master)
summary(pcresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Wine.Treatment 3 11.29 3.763 13.66 1.96e-06 ***
## Residuals 44 12.12 0.275
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
pcresol.tukey <- HSD.test(pcresol.aov, c("Wine.Treatment"), group = T, console = T)
##
## Study: pcresol.aov ~ c("Wine.Treatment")
##
## HSD Test for pCresol
##
## Mean Square Error: 0.2753786
##
## Wine.Treatment, means
##
## pCresol std r Min Max
## bucket 1.863333 0.4630204 12 1.09 2.76
## free 1.905833 0.2714426 12 1.55 2.36
## press1 0.870000 0.3354508 12 0.24 1.54
## press2 0.965000 0.8372086 12 0.25 3.42
##
## Alpha: 0.05 ; DF Error: 44
## Critical Value of Studentized Range: 3.775958
##
## Minimun Significant Difference: 0.5720074
##
## Treatments with the same letter are not significantly different.
##
## pCresol groups
## free 1.905833 a
## bucket 1.863333 a
## press2 0.965000 b
## press1 0.870000 b
2-Way Interaction
pcresol.aov <- aov(pCresol~Location*Wine.Treatment, data = master)
summary(pcresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 0.064 0.064 0.229 0.635
## Wine.Treatment 3 11.288 3.763 13.481 3.19e-06 ***
## Location:Wine.Treatment 3 0.888 0.296 1.061 0.377
## Residuals 40 11.165 0.279
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Eugenol
By Location
eugenol.aov <- aov(Eugenol~Location, data = master)
summary(eugenol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 20.0 20.008 2.765 0.103
## Residuals 46 332.9 7.236
By Wine Treatment
eugenol.aov <- aov(Eugenol~Wine.Treatment, data = master)
summary(eugenol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Wine.Treatment 3 233.3 77.77 28.61 2.03e-10 ***
## Residuals 44 119.6 2.72
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
eugenol.tukey <- HSD.test(eugenol.aov, c("Wine.Treatment"), group = T, console = T)
##
## Study: eugenol.aov ~ c("Wine.Treatment")
##
## HSD Test for Eugenol
##
## Mean Square Error: 2.717746
##
## Wine.Treatment, means
##
## Eugenol std r Min Max
## bucket 10.622500 1.7553923 12 7.90 12.72
## free 10.608333 0.7285082 12 9.78 11.86
## press1 5.821667 1.0442730 12 4.15 8.31
## press2 6.673333 2.4836166 12 4.40 12.06
##
## Alpha: 0.05 ; DF Error: 44
## Critical Value of Studentized Range: 3.775958
##
## Minimun Significant Difference: 1.796971
##
## Treatments with the same letter are not significantly different.
##
## Eugenol groups
## bucket 10.622500 a
## free 10.608333 a
## press2 6.673333 b
## press1 5.821667 b
2-Way Interaction
eugenol.aov <- aov(Eugenol~Location*Wine.Treatment, data = master)
summary(eugenol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 20.01 20.01 9.421 0.00384 **
## Wine.Treatment 3 233.30 77.77 36.616 1.5e-11 ***
## Location:Wine.Treatment 3 14.62 4.87 2.295 0.09251 .
## Residuals 40 84.95 2.12
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Grenache
setwd("C:/Users/19419/OneDrive - Washington State University (email.wsu.edu)/WSU/SMOKE/Stone Tree Bin Trials")
master <- read.csv("grenache gc.csv", header = T)
master$Location <- as.character(master$Location)
master$Wine.Treatment <- as.character(master$Wine.Treatment)
Guaiacol
By Location
guaiacol.aov <- aov(Guiacol~Location, data = master)
summary(guaiacol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 2271 2271.4 45.52 2.69e-08 ***
## Residuals 44 2196 49.9
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
guaiacol.tukey <- HSD.test(guaiacol.aov, c("Location"), group = T, console = T)
##
## Study: guaiacol.aov ~ c("Location")
##
## HSD Test for Guiacol
##
## Mean Square Error: 49.89718
##
## Location, means
##
## Guiacol std r Min Max
## High 31.96917 8.054954 24 16.77 44.79
## Low 17.90182 5.786616 22 8.33 27.00
##
## Alpha: 0.05 ; DF Error: 44
## Critical Value of Studentized Range: 2.85016
##
## Groups according to probability of means differences and alpha level( 0.05 )
##
## Treatments with the same letter are not significantly different.
##
## Guiacol groups
## High 31.96917 a
## Low 17.90182 b
By Wine Treatment
guaiacol.aov <- aov(Guiacol~Wine.Treatment, data = master)
summary(guaiacol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Wine.Treatment 3 2053 684.5 11.91 8.9e-06 ***
## Residuals 42 2414 57.5
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
guaiacol.tukey <- HSD.test(guaiacol.aov, c("Wine.Treatment"), group = T, console = T)
##
## Study: guaiacol.aov ~ c("Wine.Treatment")
##
## HSD Test for Guiacol
##
## Mean Square Error: 57.46515
##
## Wine.Treatment, means
##
## Guiacol std r Min Max
## bucket 16.05000 6.360387 12 8.33 26.44
## free 23.49917 7.660778 12 13.97 32.58
## press1 34.67500 8.872523 10 21.86 44.79
## press2 28.31333 7.474070 12 19.16 37.97
##
## Alpha: 0.05 ; DF Error: 42
## Critical Value of Studentized Range: 3.782961
##
## Groups according to probability of means differences and alpha level( 0.05 )
##
## Treatments with the same letter are not significantly different.
##
## Guiacol groups
## press1 34.67500 a
## press2 28.31333 ab
## free 23.49917 bc
## bucket 16.05000 c
2-Way Interaction
guaiacol.aov <- aov(Guiacol~Location*Wine.Treatment, data = master)
summary(guaiacol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 2271.4 2271.4 257.403 < 2e-16 ***
## Wine.Treatment 3 1809.5 603.2 68.353 2.25e-15 ***
## Location:Wine.Treatment 3 50.6 16.9 1.912 0.144
## Residuals 38 335.3 8.8
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Creosol
By Location
creosol.aov <- aov(Cresol~Location, data = master)
summary(creosol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 184.90 184.90 169.8 <2e-16 ***
## Residuals 44 47.91 1.09
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
creosol.tukey <- HSD.test(creosol.aov, c("Location"), group = T, console = T)
##
## Study: creosol.aov ~ c("Location")
##
## HSD Test for Cresol
##
## Mean Square Error: 1.088836
##
## Location, means
##
## Cresol std r Min Max
## High 7.370417 1.2248034 24 5.22 9.74
## Low 3.356818 0.7989719 22 1.89 4.94
##
## Alpha: 0.05 ; DF Error: 44
## Critical Value of Studentized Range: 2.85016
##
## Groups according to probability of means differences and alpha level( 0.05 )
##
## Treatments with the same letter are not significantly different.
##
## Cresol groups
## High 7.370417 a
## Low 3.356818 b
By Wine Treatment
creosol.aov <- aov(Cresol~Wine.Treatment, data = master)
summary(creosol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Wine.Treatment 3 15.36 5.119 0.989 0.407
## Residuals 42 217.45 5.177
2-Way Interaction
creosol.aov <- aov(Cresol~Location*Wine.Treatment, data = master)
summary(creosol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 184.90 184.90 367.649 < 2e-16 ***
## Wine.Treatment 3 23.90 7.97 15.844 7.7e-07 ***
## Location:Wine.Treatment 3 4.89 1.63 3.243 0.0325 *
## Residuals 38 19.11 0.50
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
creosol.tukey <- HSD.test(creosol.aov, c("Location","Wine.Treatment"), group = T, console = T)
##
## Study: creosol.aov ~ c("Location", "Wine.Treatment")
##
## HSD Test for Cresol
##
## Mean Square Error: 0.5029333
##
## Location:Wine.Treatment, means
##
## Cresol std r Min Max
## High:bucket 7.903333 0.8788326 6 6.95 9.24
## High:free 7.478333 0.4932917 6 6.89 8.16
## High:press1 5.728333 0.4570959 6 5.22 6.21
## High:press2 8.371667 0.9394342 6 7.32 9.74
## Low:bucket 3.118333 0.9471096 6 1.89 4.22
## Low:free 3.235000 0.4332090 6 2.59 3.77
## Low:press1 2.805000 0.1537314 4 2.68 3.02
## Low:press2 4.085000 0.7850287 6 2.92 4.94
##
## Alpha: 0.05 ; DF Error: 38
## Critical Value of Studentized Range: 4.53321
##
## Groups according to probability of means differences and alpha level( 0.05 )
##
## Treatments with the same letter are not significantly different.
##
## Cresol groups
## High:press2 8.371667 a
## High:bucket 7.903333 a
## High:free 7.478333 a
## High:press1 5.728333 b
## Low:press2 4.085000 c
## Low:free 3.235000 c
## Low:bucket 3.118333 c
## Low:press1 2.805000 c
Interaction Plot:
interaction.plot(master$Location,master$Wine.Treatment,master$Cresol)
o-Cresol
By Location
ocresol.aov <- aov(oCresol~Location, data = master)
summary(ocresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 96.64 96.64 62.11 5.98e-10 ***
## Residuals 44 68.46 1.56
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
ocresol.tukey <- HSD.test(ocresol.aov, c("Location"), group = T, console = T)
##
## Study: ocresol.aov ~ c("Location")
##
## HSD Test for oCresol
##
## Mean Square Error: 1.555839
##
## Location, means
##
## oCresol std r Min Max
## High 8.267917 1.394729 24 5.16 9.72
## Low 5.366364 1.062694 22 2.72 6.45
##
## Alpha: 0.05 ; DF Error: 44
## Critical Value of Studentized Range: 2.85016
##
## Groups according to probability of means differences and alpha level( 0.05 )
##
## Treatments with the same letter are not significantly different.
##
## oCresol groups
## High 8.267917 a
## Low 5.366364 b
By Wine Treatment
ocresol.aov <- aov(oCresol~Wine.Treatment, data = master)
summary(ocresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Wine.Treatment 3 46.07 15.358 5.419 0.00304 **
## Residuals 42 119.02 2.834
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
ocresol.tukey <- HSD.test(ocresol.aov, c("Wine.Treatment"), group = T, console = T)
##
## Study: ocresol.aov ~ c("Wine.Treatment")
##
## HSD Test for oCresol
##
## Mean Square Error: 2.833796
##
## Wine.Treatment, means
##
## oCresol std r Min Max
## bucket 6.940000 1.636654 12 4.68 9.10
## free 7.721667 1.819565 12 5.69 9.72
## press1 5.063000 1.537914 10 2.72 6.73
## press2 7.493333 1.701573 12 5.56 9.67
##
## Alpha: 0.05 ; DF Error: 42
## Critical Value of Studentized Range: 3.782961
##
## Groups according to probability of means differences and alpha level( 0.05 )
##
## Treatments with the same letter are not significantly different.
##
## oCresol groups
## free 7.721667 a
## press2 7.493333 a
## bucket 6.940000 ab
## press1 5.063000 b
2-Way Interaction
ocresol.aov <- aov(oCresol~Location*Wine.Treatment, data = master)
summary(ocresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 96.64 96.64 370.636 < 2e-16 ***
## Wine.Treatment 3 57.64 19.21 73.696 6.73e-16 ***
## Location:Wine.Treatment 3 0.91 0.30 1.157 0.339
## Residuals 38 9.91 0.26
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
m-Cresol
By Location
mcresol.aov <- aov(mCresol~Location, data = master)
summary(mcresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 171.69 171.69 111.7 1.18e-13 ***
## Residuals 44 67.63 1.54
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mcresol.tukey <- HSD.test(mcresol.aov, c("Location"), group = T, console = T)
##
## Study: mcresol.aov ~ c("Location")
##
## HSD Test for mCresol
##
## Mean Square Error: 1.537086
##
## Location, means
##
## mCresol std r Min Max
## High 8.1425 1.5070220 24 4.84 10.19
## Low 4.2750 0.8562418 22 2.96 5.72
##
## Alpha: 0.05 ; DF Error: 44
## Critical Value of Studentized Range: 2.85016
##
## Groups according to probability of means differences and alpha level( 0.05 )
##
## Treatments with the same letter are not significantly different.
##
## mCresol groups
## High 8.1425 a
## Low 4.2750 b
By Wine Treatment
mcresol.aov <- aov(mCresol~Wine.Treatment, data = master)
summary(mcresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Wine.Treatment 3 23.38 7.794 1.516 0.224
## Residuals 42 215.94 5.141
2-Way Interaction
mcresol.aov <- aov(mCresol~Location*Wine.Treatment, data = master)
summary(mcresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 171.69 171.69 249.966 < 2e-16 ***
## Wine.Treatment 3 34.55 11.52 16.768 4.25e-07 ***
## Location:Wine.Treatment 3 6.98 2.33 3.388 0.0277 *
## Residuals 38 26.10 0.69
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
mcresol.tukey <- HSD.test(mcresol.aov, c("Location","Wine.Treatment"), group = T, console = T)
##
## Study: mcresol.aov ~ c("Location", "Wine.Treatment")
##
## HSD Test for mCresol
##
## Mean Square Error: 0.6868399
##
## Location:Wine.Treatment, means
##
## mCresol std r Min Max
## High:bucket 8.173333 0.9655189 6 6.65 9.56
## High:free 8.878333 1.1222730 6 7.42 9.94
## High:press1 6.120000 0.6962471 6 4.84 6.76
## High:press2 9.398333 0.5644614 6 8.62 10.19
## Low:bucket 4.496667 0.9297455 6 3.54 5.63
## Low:free 4.268333 0.4977315 6 3.64 4.82
## Low:press1 3.390000 0.5051072 4 2.96 4.00
## Low:press2 4.650000 0.9796122 6 3.51 5.72
##
## Alpha: 0.05 ; DF Error: 38
## Critical Value of Studentized Range: 4.53321
##
## Groups according to probability of means differences and alpha level( 0.05 )
##
## Treatments with the same letter are not significantly different.
##
## mCresol groups
## High:press2 9.398333 a
## High:free 8.878333 a
## High:bucket 8.173333 a
## High:press1 6.120000 b
## Low:press2 4.650000 bc
## Low:bucket 4.496667 c
## Low:free 4.268333 c
## Low:press1 3.390000 c
Interaction Plot:
interaction.plot(master$Location,master$Wine.Treatment,master$mCresol)
pCresol
By Location
pcresol.aov <- aov(pCresol~Location, data = master)
summary(pcresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 2.265 2.2653 5.497 0.0236 *
## Residuals 44 18.134 0.4121
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
pcresol.tukey <- HSD.test(pcresol.aov, c("Location"), group = T, console = T)
##
## Study: pcresol.aov ~ c("Location")
##
## HSD Test for pCresol
##
## Mean Square Error: 0.4121253
##
## Location, means
##
## pCresol std r Min Max
## High 1.983333 0.7428539 24 0.77 3.46
## Low 1.539091 0.5090318 22 0.72 2.62
##
## Alpha: 0.05 ; DF Error: 44
## Critical Value of Studentized Range: 2.85016
##
## Groups according to probability of means differences and alpha level( 0.05 )
##
## Treatments with the same letter are not significantly different.
##
## pCresol groups
## High 1.983333 a
## Low 1.539091 b
By Wine Treatment
pcresol.aov <- aov(pCresol~Wine.Treatment, data = master)
summary(pcresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Wine.Treatment 3 8.018 2.6727 9.067 9.52e-05 ***
## Residuals 42 12.381 0.2948
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
pcresol.tukey <- HSD.test(pcresol.aov, c("Wine.Treatment"), group = T, console = T)
##
## Study: pcresol.aov ~ c("Wine.Treatment")
##
## HSD Test for pCresol
##
## Mean Square Error: 0.2947754
##
## Wine.Treatment, means
##
## pCresol std r Min Max
## bucket 2.079167 0.8062873 12 0.72 3.46
## free 2.184167 0.4747336 12 1.19 2.74
## press1 1.099000 0.2458297 10 0.77 1.45
## press2 1.609167 0.4478729 12 1.19 2.27
##
## Alpha: 0.05 ; DF Error: 42
## Critical Value of Studentized Range: 3.782961
##
## Groups according to probability of means differences and alpha level( 0.05 )
##
## Treatments with the same letter are not significantly different.
##
## pCresol groups
## free 2.184167 a
## bucket 2.079167 a
## press2 1.609167 ab
## press1 1.099000 b
2-Way Interaction
pcresol.aov <- aov(pCresol~Location*Wine.Treatment, data = master)
summary(pcresol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 2.265 2.2653 12.353 0.00116 **
## Wine.Treatment 3 8.674 2.8915 15.768 8.1e-07 ***
## Location:Wine.Treatment 3 2.491 0.8302 4.527 0.00825 **
## Residuals 38 6.968 0.1834
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
pcresol.tukey <- HSD.test(pcresol.aov, c("Location","Wine.Treatment"), group = T, console = T)
##
## Study: pcresol.aov ~ c("Location", "Wine.Treatment")
##
## HSD Test for pCresol
##
## Mean Square Error: 0.1833787
##
## Location:Wine.Treatment, means
##
## pCresol std r Min Max
## High:bucket 2.681667 0.6126146 6 1.72 3.46
## High:free 2.263333 0.4494738 6 1.70 2.74
## High:press1 1.093333 0.2857038 6 0.77 1.45
## High:press2 1.895000 0.4902142 6 1.19 2.27
## Low:bucket 1.476667 0.4286102 6 0.72 1.82
## Low:free 2.105000 0.5279678 6 1.19 2.62
## Low:press1 1.107500 0.2123480 4 0.86 1.33
## Low:press2 1.323333 0.0700476 6 1.20 1.38
##
## Alpha: 0.05 ; DF Error: 38
## Critical Value of Studentized Range: 4.53321
##
## Groups according to probability of means differences and alpha level( 0.05 )
##
## Treatments with the same letter are not significantly different.
##
## pCresol groups
## High:bucket 2.681667 a
## High:free 2.263333 ab
## Low:free 2.105000 abc
## High:press2 1.895000 abcd
## Low:bucket 1.476667 bcde
## Low:press2 1.323333 cde
## Low:press1 1.107500 de
## High:press1 1.093333 e
interaction.plot(master$Location,master$Wine.Treatment,master$pCresol)
Eugenol
By Location
eugenol.aov <- aov(Eugenol~Location, data = master)
summary(eugenol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 0.685 0.6845 1.415 0.241
## Residuals 44 21.282 0.4837
By Wine Treatment
eugenol.aov <- aov(Eugenol~Wine.Treatment, data = master)
summary(eugenol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Wine.Treatment 3 16.95 5.649 47.26 1.6e-13 ***
## Residuals 42 5.02 0.120
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
eugenol.tukey <- HSD.test(eugenol.aov, c("Wine.Treatment"), group = T, console = T)
##
## Study: eugenol.aov ~ c("Wine.Treatment")
##
## HSD Test for Eugenol
##
## Mean Square Error: 0.1195185
##
## Wine.Treatment, means
##
## Eugenol std r Min Max
## bucket 2.580000 0.5352145 12 1.72 3.13
## free 1.305000 0.2233424 12 0.89 1.60
## press1 1.053000 0.2461278 10 0.77 1.44
## press2 2.106667 0.2654099 12 1.71 2.44
##
## Alpha: 0.05 ; DF Error: 42
## Critical Value of Studentized Range: 3.782961
##
## Groups according to probability of means differences and alpha level( 0.05 )
##
## Treatments with the same letter are not significantly different.
##
## Eugenol groups
## bucket 2.580000 a
## press2 2.106667 b
## free 1.305000 c
## press1 1.053000 c
2-Way Interaction
eugenol.aov <- aov(Eugenol~Location*Wine.Treatment, data = master)
summary(eugenol.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## Location 1 0.685 0.685 8.004 0.00741 **
## Wine.Treatment 3 17.363 5.788 67.676 2.64e-15 ***
## Location:Wine.Treatment 3 0.670 0.223 2.610 0.06549 .
## Residuals 38 3.250 0.086
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1