In: Statistics and Probability
unit anova (biostatistics) using R
on the table below Oxygen consumption you will find the rate of oxygen consumption of two species of limpets (Acmaea scabra & A. digitalis) in three different concentrations of salt water. The variable is given as μl O2/mg units per dry weight at 22°C. The researcher would like to know whether the rate of oxygen consumption differs between both species and the different concentrations of salt water.
Coding in Oxygen consumption:
Species 1 = Acmaea scabra
Species 2 = A. digitalis
Concentration 1 = 100%
Oxygen consumption:
consumption | species | concentration |
7.16 | 1 | 1 |
6.78 | 1 | 1 |
13.6 | 1 | 1 |
8.93 | 1 | 1 |
8.26 | 1 | 1 |
14 | 1 | 1 |
16.1 | 1 | 1 |
9.66 | 1 | 1 |
5.2 | 1 | 2 |
5.2 | 1 | 2 |
7.18 | 1 | 2 |
6.37 | 1 | 2 |
13.2 | 1 | 2 |
8.39 | 1 | 2 |
10.4 | 1 | 2 |
7.18 | 1 | 2 |
11.11 | 1 | 3 |
9.74 | 1 | 3 |
18.8 | 1 | 3 |
9.74 | 1 | 3 |
10.5 | 1 | 3 |
14.6 | 1 | 3 |
11.1 | 1 | 3 |
11.8 | 1 | 3 |
6.14 | 2 | 1 |
3.86 | 2 | 1 |
10.4 | 2 | 1 |
5.49 | 2 | 1 |
6.14 | 2 | 1 |
10 | 2 | 1 |
11.6 | 2 | 1 |
5.8 | 2 | 1 |
4.47 | 2 | 2 |
9.9 | 2 | 2 |
5.75 | 2 | 2 |
11.8 | 2 | 2 |
4.95 | 2 | 2 |
6.49 | 2 | 2 |
5.44 | 2 | 2 |
9.9 | 2 | 2 |
9.63 | 2 | 3 |
6.38 | 2 | 3 |
13.4 | 2 | 3 |
14.5 | 2 | 3 |
14.5 | 2 | 3 |
10.2 | 2 | 3 |
17.7 | 2 | 3 |
12.3 | 2 | 3 |
The r-codes are given below
data=read.csv(file.choose())
attach(data)
twanova<-aov(consumption~as.factor(species)*as.factor(concentration))
twanova
output |
|||||
Two-Way ANOVA | |||||
---|---|---|---|---|---|
df | Sum of Squares | Mean Squares | F | Pr(>F) | |
species | 1 | 16.638075 | 16.638075 | 1.74037878 | 0.19423806 |
concentration | 2 | 181.32101667 | 90.66050833 | 9.48328607 | 0.00039931 |
species: concentration | 2 | 23.9262 | 11.9631 | 1.25136624 | 0.29656157 |
Residuals | 42 | 401.5213 | 9.56003095 | NA | NA |
TukeyHSD(twanova)
For species |
||||
Tukey's HSD | ||||
---|---|---|---|---|
diff | lwr | upr | p adj | |
2-1 | -1.1775 | -2.97876567 | 0.62376567 | 0.19423806 |
For concentartion |
||||
Tukey's HSD | ||||
---|---|---|---|---|
diff | lwr | upr | p adj | |
2-1 | -1.38125 | -4.03708104 | 1.27458104 | 0.42339233 |
3-1 | 3.255 | 0.59916896 | 5.91083104 | 0.01306912 |
3-2 | 4.63625 | 1.98041896 | 7.29208104 | 0.00034582 |
From the we conclude that there is significant differen in the concetration only because the p-value associatiated with is 0.0003 which is significant ans species are insignificant because p-value > 5% level of significance.
In concentration the diffrence in solt of 3-a and 3-2 is significant, because p- value associated with it is < 0.05 level of significance