MPD dose dommer wheel and motion
MPDdoseDW <- read.csv("~/Desktop/experiment2statsD.csv")
MPDdoseDW
## Group.E2DmpdW Day.E2DmpdW Data.E2DmpdW Group.E2DmpdM Day.E2DmpdM
## 1 Saline 1 13.65313 Saline 1
## 2 Saline 1 35.68659 Saline 1
## 3 Saline 1 10.75247 Saline 1
## 4 2.5 1 0.00000 2.5 1
## 5 2.5 1 141.92391 2.5 1
## 6 2.5 1 73.71933 2.5 1
## 7 5 1 9.08067 5 1
## 8 5 1 152.17974 5 1
## 9 5 1 18.52850 5 1
## 10 10 1 93.70120 10 1
## 11 10 1 86.97943 10 1
## 12 10 1 1.65909 10 1
## 13 Saline 2 4.32358 Saline 2
## 14 Saline 2 36.91971 Saline 2
## 15 Saline 2 40.58958 Saline 2
## 16 2.5 2 0.00000 2.5 2
## 17 2.5 2 81.34606 2.5 2
## 18 2.5 2 0.07171 2.5 2
## 19 5 2 0.37602 5 2
## 20 5 2 68.39501 5 2
## 21 5 2 188.35905 5 2
## 22 10 2 108.22402 10 2
## 23 10 2 142.19247 10 2
## 24 10 2 34.17707 10 2
## 25 Saline 3 32.31135 Saline 3
## 26 Saline 3 26.03506 Saline 3
## 27 Saline 3 55.26250 Saline 3
## 28 2.5 3 0.00000 2.5 3
## 29 2.5 3 1.17120 2.5 3
## 30 2.5 3 5.12915 2.5 3
## 31 5 3 0.14861 5 3
## 32 5 3 45.11864 5 3
## 33 5 3 67.10686 5 3
## 34 10 3 65.83755 10 3
## 35 10 3 137.70446 10 3
## 36 10 3 17.11079 10 3
## 37 Saline 4 16.91048 Saline 4
## 38 Saline 4 33.04933 Saline 4
## 39 Saline 4 69.82811 Saline 4
## 40 2.5 4 0.00000 2.5 4
## 41 2.5 4 87.17062 2.5 4
## 42 2.5 4 26.14003 2.5 4
## 43 5 4 41.73102 5 4
## 44 5 4 45.38542 5 4
## 45 5 4 151.87475 5 4
## 46 10 4 94.57602 10 4
## 47 10 4 121.05845 10 4
## 48 10 4 1.61588 10 4
## Data.E2DmpdM
## 1 92.217
## 2 49.515
## 3 32.925
## 4 76.147
## 5 79.314
## 6 87.074
## 7 86.709
## 8 83.702
## 9 112.190
## 10 152.733
## 11 88.146
## 12 92.578
## 13 74.946
## 14 48.044
## 15 55.227
## 16 88.554
## 17 68.678
## 18 23.631
## 19 2.642
## 20 59.724
## 21 135.906
## 22 91.123
## 23 87.691
## 24 18.837
## 25 52.881
## 26 64.964
## 27 139.081
## 28 88.554
## 29 54.653
## 30 54.633
## 31 60.612
## 32 97.982
## 33 131.918
## 34 91.207
## 35 142.751
## 36 10.237
## 37 54.755
## 38 41.482
## 39 73.985
## 40 62.638
## 41 51.012
## 42 85.156
## 43 80.072
## 44 63.662
## 45 95.962
## 46 91.732
## 47 58.206
## 48 39.688
tapply(MPDdoseDW$Data.E2DmpdW, MPDdoseDW$Group.E2DmpdW, sd)
## 10 2.5 5 Saline
## 50.82 48.65 64.22 18.82
tapply(MPDdoseDW$Data.E2DmpdW, MPDdoseDW$Day.E2DmpdW, sd)
## 1 2 3 4
## 54.92 61.41 40.27 48.06
MPDdoseDW$Group.E2DmpdW <- as.factor(MPDdoseDW$Group.E2DmpdW)
MPDdoseDW$Day.E2DmpdW <- as.factor(MPDdoseDW$Day.E2DmpdW)
MPDdoseDW$Data.E2DmpdW <- as.numeric(MPDdoseDW$Data.E2DmpdW)
aov.MPDdoseDW = aov(Data.E2DmpdW ~ Group.E2DmpdW * Day.E2DmpdW, data = MPDdoseDW)
summary(aov.MPDdoseDW)
## Df Sum Sq Mean Sq F value Pr(>F)
## Group.E2DmpdW 3 17555 5852 2.10 0.12
## Day.E2DmpdW 3 3354 1118 0.40 0.75
## Group.E2DmpdW:Day.E2DmpdW 9 11057 1229 0.44 0.90
## Residuals 32 89298 2791
comparison <- aov(Data.E2DmpdW ~ Group.E2DmpdW * Day.E2DmpdW, data = MPDdoseDW)
TukeyHSD(comparison, "Group.E2DmpdW")
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = Data.E2DmpdW ~ Group.E2DmpdW * Day.E2DmpdW, data = MPDdoseDW)
##
## $Group.E2DmpdW
## diff lwr upr p adj
## 2.5-10 -40.680 -99.11 17.75 0.2539
## 5-10 -9.713 -68.14 48.72 0.9691
## Saline-10 -44.126 -102.56 14.30 0.1928
## 5-2.5 30.968 -27.46 89.40 0.4870
## Saline-2.5 -3.446 -61.88 54.98 0.9985
## Saline-5 -34.414 -92.84 24.02 0.3953
TukeyHSD(comparison, "Day.E2DmpdW")
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = Data.E2DmpdW ~ Group.E2DmpdW * Day.E2DmpdW, data = MPDdoseDW)
##
## $Day.E2DmpdW
## diff lwr upr p adj
## 2-1 5.593 -52.84 64.02 0.9938
## 3-1 -15.411 -73.84 43.02 0.8906
## 4-1 4.290 -54.14 62.72 0.9972
## 3-2 -21.003 -79.43 37.43 0.7651
## 4-2 -1.303 -59.73 57.13 0.9999
## 4-3 19.700 -38.73 78.13 0.7978
print(model.tables(aov.MPDdoseDW, "means"), digits = 9)
## Tables of means
## Grand mean
##
## 51.77
##
## Group.E2DmpdW
## Group.E2DmpdW
## 10 2.5 5 Saline
## 75.40 34.72 65.69 31.28
##
## Day.E2DmpdW
## Day.E2DmpdW
## 1 2 3 4
## 53.16 58.75 37.74 57.45
##
## Group.E2DmpdW:Day.E2DmpdW
## Day.E2DmpdW
## Group.E2DmpdW 1 2 3 4
## 10 60.78 94.86 73.55 72.42
## 2.5 71.88 27.14 2.10 37.77
## 5 59.93 85.71 37.46 79.66
## Saline 20.03 27.28 37.87 39.93
MPDdoseDM <- read.csv("~/Desktop/experiment2statsD.csv")
MPDdoseDM
## Group.E2DmpdW Day.E2DmpdW Data.E2DmpdW Group.E2DmpdM Day.E2DmpdM
## 1 Saline 1 13.65313 Saline 1
## 2 Saline 1 35.68659 Saline 1
## 3 Saline 1 10.75247 Saline 1
## 4 2.5 1 0.00000 2.5 1
## 5 2.5 1 141.92391 2.5 1
## 6 2.5 1 73.71933 2.5 1
## 7 5 1 9.08067 5 1
## 8 5 1 152.17974 5 1
## 9 5 1 18.52850 5 1
## 10 10 1 93.70120 10 1
## 11 10 1 86.97943 10 1
## 12 10 1 1.65909 10 1
## 13 Saline 2 4.32358 Saline 2
## 14 Saline 2 36.91971 Saline 2
## 15 Saline 2 40.58958 Saline 2
## 16 2.5 2 0.00000 2.5 2
## 17 2.5 2 81.34606 2.5 2
## 18 2.5 2 0.07171 2.5 2
## 19 5 2 0.37602 5 2
## 20 5 2 68.39501 5 2
## 21 5 2 188.35905 5 2
## 22 10 2 108.22402 10 2
## 23 10 2 142.19247 10 2
## 24 10 2 34.17707 10 2
## 25 Saline 3 32.31135 Saline 3
## 26 Saline 3 26.03506 Saline 3
## 27 Saline 3 55.26250 Saline 3
## 28 2.5 3 0.00000 2.5 3
## 29 2.5 3 1.17120 2.5 3
## 30 2.5 3 5.12915 2.5 3
## 31 5 3 0.14861 5 3
## 32 5 3 45.11864 5 3
## 33 5 3 67.10686 5 3
## 34 10 3 65.83755 10 3
## 35 10 3 137.70446 10 3
## 36 10 3 17.11079 10 3
## 37 Saline 4 16.91048 Saline 4
## 38 Saline 4 33.04933 Saline 4
## 39 Saline 4 69.82811 Saline 4
## 40 2.5 4 0.00000 2.5 4
## 41 2.5 4 87.17062 2.5 4
## 42 2.5 4 26.14003 2.5 4
## 43 5 4 41.73102 5 4
## 44 5 4 45.38542 5 4
## 45 5 4 151.87475 5 4
## 46 10 4 94.57602 10 4
## 47 10 4 121.05845 10 4
## 48 10 4 1.61588 10 4
## Data.E2DmpdM
## 1 92.217
## 2 49.515
## 3 32.925
## 4 76.147
## 5 79.314
## 6 87.074
## 7 86.709
## 8 83.702
## 9 112.190
## 10 152.733
## 11 88.146
## 12 92.578
## 13 74.946
## 14 48.044
## 15 55.227
## 16 88.554
## 17 68.678
## 18 23.631
## 19 2.642
## 20 59.724
## 21 135.906
## 22 91.123
## 23 87.691
## 24 18.837
## 25 52.881
## 26 64.964
## 27 139.081
## 28 88.554
## 29 54.653
## 30 54.633
## 31 60.612
## 32 97.982
## 33 131.918
## 34 91.207
## 35 142.751
## 36 10.237
## 37 54.755
## 38 41.482
## 39 73.985
## 40 62.638
## 41 51.012
## 42 85.156
## 43 80.072
## 44 63.662
## 45 95.962
## 46 91.732
## 47 58.206
## 48 39.688
tapply(MPDdoseDM$Data.E2DmpdM, MPDdoseDM$Group.E2DmpdM, sd)
## 10 2.5 5 Saline
## 43.28 19.86 36.07 28.39
tapply(MPDdoseDM$Data.E2DmpdM, MPDdoseDM$Day.E2DmpdM, sd)
## 1 2 3 4
## 29.42 36.91 40.52 18.85
MPDdoseDM$Group.E2DmpdM <- as.factor(MPDdoseDM$Group.E2DmpdM)
MPDdoseDM$Day.E2DmpdM <- as.factor(MPDdoseDM$Day.E2DmpdM)
MPDdoseDM$Data.E2DmpdM <- as.numeric(MPDdoseDM$Data.E2DmpdM)
aov.MPDdoseDM = aov(Data.E2DmpdM ~ Group.E2DmpdM * Day.E2DmpdM, data = MPDdoseDM)
summary(aov.MPDdoseDM)
## Df Sum Sq Mean Sq F value Pr(>F)
## Group.E2DmpdM 3 3100 1033 0.84 0.48
## Day.E2DmpdM 3 4748 1583 1.28 0.30
## Group.E2DmpdM:Day.E2DmpdM 9 3838 426 0.35 0.95
## Residuals 32 39538 1236
comparison <- aov(Data.E2DmpdM ~ Group.E2DmpdM * Day.E2DmpdM, data = MPDdoseDM)
TukeyHSD(comparison, "Group.E2DmpdM")
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = Data.E2DmpdM ~ Group.E2DmpdM * Day.E2DmpdM, data = MPDdoseDM)
##
## $Group.E2DmpdM
## diff lwr upr p adj
## 2.5-10 -12.074 -50.95 26.81 0.8343
## 5-10 3.846 -35.03 42.73 0.9931
## Saline-10 -15.409 -54.29 23.47 0.7077
## 5-2.5 15.920 -22.96 54.80 0.6865
## Saline-2.5 -3.335 -42.22 35.54 0.9955
## Saline-5 -19.255 -58.13 19.62 0.5439
TukeyHSD(comparison, "Day.E2DmpdM")
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = Data.E2DmpdM ~ Group.E2DmpdM * Day.E2DmpdM, data = MPDdoseDM)
##
## $Day.E2DmpdM
## diff lwr upr p adj
## 2-1 -23.187 -62.07 15.69 0.3844
## 3-1 -3.648 -42.53 35.23 0.9941
## 4-1 -19.575 -58.45 19.30 0.5303
## 3-2 19.539 -19.34 58.42 0.5319
## 4-2 3.612 -35.27 42.49 0.9943
## 4-3 -15.927 -54.81 22.95 0.6862
print(model.tables(aov.MPDdoseDM, "means"), digits = 9)
## Tables of means
## Grand mean
##
## 74.5
##
## Group.E2DmpdM
## Group.E2DmpdM
## 10 2.5 5 Saline
## 80.41 68.34 84.26 65.00
##
## Day.E2DmpdM
## Day.E2DmpdM
## 1 2 3 4
## 86.10 62.92 82.46 66.53
##
## Group.E2DmpdM:Day.E2DmpdM
## Day.E2DmpdM
## Group.E2DmpdM 1 2 3 4
## 10 111.15 65.88 81.40 63.21
## 2.5 80.84 60.29 65.95 66.27
## 5 94.20 66.09 96.84 79.90
## Saline 58.22 59.41 85.64 56.74