MPD dose falafel wheel and motion
MPDdoseFW <- read.csv("~/Desktop/experiment2statsF.csv")
MPDdoseFW
## Group.E2FmpdW Day.E2FmpdW Data.E2FmpdW Group.E2FmpdM Day.E2FmpdM
## 1 Saline 1 66.9757 Saline 1
## 2 Saline 1 39.8084 Saline 1
## 3 Saline 1 71.8449 Saline 1
## 4 Saline 1 0.4504 Saline 1
## 5 2.5 1 120.3465 2.5 1
## 6 2.5 1 87.7458 2.5 1
## 7 2.5 1 135.0673 2.5 1
## 8 2.5 1 46.9731 2.5 1
## 9 2.5 1 96.2350 2.5 1
## 10 5 1 31.2685 5 1
## 11 5 1 68.8159 5 1
## 12 5 1 74.7980 5 1
## 13 5 1 113.3235 5 1
## 14 10 1 103.9607 10 1
## 15 10 1 63.7161 10 1
## 16 10 1 73.7223 10 1
## 17 10 1 95.2725 10 1
## 18 Saline 2 73.4696 Saline 2
## 19 Saline 2 77.5923 Saline 2
## 20 Saline 2 69.0582 Saline 2
## 21 Saline 2 1.3649 Saline 2
## 22 2.5 2 71.4790 2.5 2
## 23 2.5 2 55.0731 2.5 2
## 24 2.5 2 96.3686 2.5 2
## 25 2.5 2 60.0803 2.5 2
## 26 2.5 2 120.5474 2.5 2
## 27 5 2 63.1134 5 2
## 28 5 2 40.7544 5 2
## 29 5 2 58.2419 5 2
## 30 5 2 63.5639 5 2
## 31 10 2 84.7728 10 2
## 32 10 2 37.5777 10 2
## 33 10 2 64.3195 10 2
## 34 10 2 98.8629 10 2
## 35 Saline 3 54.9308 Saline 3
## 36 Saline 3 53.0260 Saline 3
## 37 Saline 3 56.6311 Saline 3
## 38 Saline 3 0.4504 Saline 3
## 39 2.5 3 41.7940 2.5 3
## 40 2.5 3 50.9133 2.5 3
## 41 2.5 3 117.2209 2.5 3
## 42 2.5 3 28.7446 2.5 3
## 43 2.5 3 102.9064 2.5 3
## 44 5 3 135.0030 5 3
## 45 5 3 55.5453 5 3
## 46 5 3 51.5457 5 3
## 47 5 3 40.2068 5 3
## 48 10 3 56.6067 10 3
## 49 10 3 37.1715 10 3
## 50 10 3 47.9529 10 3
## 51 10 3 59.5799 10 3
## Data.E2FmpdM
## 1 111.09
## 2 39.03
## 3 117.98
## 4 140.66
## 5 148.08
## 6 74.47
## 7 48.58
## 8 91.93
## 9 140.05
## 10 65.35
## 11 78.28
## 12 76.48
## 13 75.77
## 14 82.12
## 15 77.00
## 16 103.03
## 17 49.78
## 18 84.33
## 19 55.57
## 20 99.40
## 21 62.55
## 22 107.61
## 23 50.40
## 24 39.74
## 25 56.89
## 26 103.37
## 27 53.42
## 28 69.74
## 29 85.77
## 30 62.01
## 31 70.61
## 32 50.34
## 33 66.41
## 34 67.59
## 35 58.36
## 36 48.38
## 37 128.49
## 38 55.52
## 39 90.56
## 40 57.26
## 41 35.36
## 42 58.25
## 43 164.04
## 44 116.08
## 45 61.69
## 46 79.48
## 47 81.86
## 48 65.40
## 49 57.79
## 50 108.00
## 51 50.26
tapply(MPDdoseFW$Data.E2FmpdW, MPDdoseFW$Group.E2FmpdW, sd)
## 10 2.5 5 Saline
## 22.98 33.65 30.13 29.85
tapply(MPDdoseFW$Data.E2FmpdW, MPDdoseFW$Day.E2FmpdW, sd)
## 1 2 3
## 34.20 26.70 32.47
MPDdoseFW$Group.E2FmpdW <- as.factor(MPDdoseFW$Group.E2FmpdW)
MPDdoseFW$Day.E2FmpdW <- as.factor(MPDdoseFW$Day.E2FmpdW)
MPDdoseFW$Data.E2FmpdW <- as.numeric(MPDdoseFW$Data.E2FmpdW)
aov.MPDdoseFW = aov(Data.E2FmpdW ~ Group.E2FmpdW * Day.E2FmpdW, data = MPDdoseFW)
summary(aov.MPDdoseFW)
## Df Sum Sq Mean Sq F value Pr(>F)
## Group.E2FmpdW 3 8193 2731 2.96 0.044 *
## Day.E2FmpdW 2 2649 1325 1.44 0.250
## Group.E2FmpdW:Day.E2FmpdW 6 2825 471 0.51 0.797
## Residuals 39 35974 922
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
comparison <- aov(Data.E2FmpdW ~ Group.E2FmpdW * Day.E2FmpdW, data = MPDdoseFW)
TukeyHSD(comparison, "Group.E2FmpdW")
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = Data.E2FmpdW ~ Group.E2FmpdW * Day.E2FmpdW, data = MPDdoseFW)
##
## $Group.E2FmpdW
## diff lwr upr p adj
## 2.5-10 13.473 -18.09 45.037 0.6641
## 5-10 -2.278 -35.55 30.993 0.9978
## Saline-10 -21.493 -54.76 11.779 0.3207
## 5-2.5 -15.751 -47.32 15.813 0.5443
## Saline-2.5 -34.966 -66.53 -3.402 0.0249
## Saline-5 -19.215 -52.49 14.056 0.4186
TukeyHSD(comparison, "Day.E2FmpdW")
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = Data.E2FmpdW ~ Group.E2FmpdW * Day.E2FmpdW, data = MPDdoseFW)
##
## $Day.E2FmpdW
## diff lwr upr p adj
## 2-1 -9.064 -34.44 16.316 0.6621
## 3-1 -17.653 -43.03 7.727 0.2201
## 3-2 -8.589 -33.97 16.791 0.6903
print(model.tables(aov.MPDdoseFW, "means"), digits = 17)
## Tables of means
## Grand mean
##
## 67
##
## Group.E2FmpdW
## 10 2.5 5
## 68.626288375833312 82.099676402666631 66.348347914999977
## rep 12.000000000000000 15.000000000000000 12.000000000000000
## Saline
## 47.133547556999957
## rep 12.000000000000000
##
## Day.E2FmpdW
## 1 2 3
## 75.901436051294084 66.837628089411737 58.248780107176444
## rep 17.000000000000000 17.000000000000000 17.000000000000000
##
## Group.E2FmpdW:Day.E2FmpdW
## Day.E2FmpdW
## Group.E2FmpdW 1 2 3
## 10 84.17 71.38 50.33
## rep 4.00 4.00 4.00
## 2.5 97.27 80.71 68.32
## rep 5.00 5.00 5.00
## 5 72.05 56.42 70.58
## rep 4.00 4.00 4.00
## Saline 44.77 55.37 41.26
## rep 4.00 4.00 4.00
MPDdoseFM <- read.csv("~/Desktop/experiment2statsF.csv")
MPDdoseFM
## Group.E2FmpdW Day.E2FmpdW Data.E2FmpdW Group.E2FmpdM Day.E2FmpdM
## 1 Saline 1 66.9757 Saline 1
## 2 Saline 1 39.8084 Saline 1
## 3 Saline 1 71.8449 Saline 1
## 4 Saline 1 0.4504 Saline 1
## 5 2.5 1 120.3465 2.5 1
## 6 2.5 1 87.7458 2.5 1
## 7 2.5 1 135.0673 2.5 1
## 8 2.5 1 46.9731 2.5 1
## 9 2.5 1 96.2350 2.5 1
## 10 5 1 31.2685 5 1
## 11 5 1 68.8159 5 1
## 12 5 1 74.7980 5 1
## 13 5 1 113.3235 5 1
## 14 10 1 103.9607 10 1
## 15 10 1 63.7161 10 1
## 16 10 1 73.7223 10 1
## 17 10 1 95.2725 10 1
## 18 Saline 2 73.4696 Saline 2
## 19 Saline 2 77.5923 Saline 2
## 20 Saline 2 69.0582 Saline 2
## 21 Saline 2 1.3649 Saline 2
## 22 2.5 2 71.4790 2.5 2
## 23 2.5 2 55.0731 2.5 2
## 24 2.5 2 96.3686 2.5 2
## 25 2.5 2 60.0803 2.5 2
## 26 2.5 2 120.5474 2.5 2
## 27 5 2 63.1134 5 2
## 28 5 2 40.7544 5 2
## 29 5 2 58.2419 5 2
## 30 5 2 63.5639 5 2
## 31 10 2 84.7728 10 2
## 32 10 2 37.5777 10 2
## 33 10 2 64.3195 10 2
## 34 10 2 98.8629 10 2
## 35 Saline 3 54.9308 Saline 3
## 36 Saline 3 53.0260 Saline 3
## 37 Saline 3 56.6311 Saline 3
## 38 Saline 3 0.4504 Saline 3
## 39 2.5 3 41.7940 2.5 3
## 40 2.5 3 50.9133 2.5 3
## 41 2.5 3 117.2209 2.5 3
## 42 2.5 3 28.7446 2.5 3
## 43 2.5 3 102.9064 2.5 3
## 44 5 3 135.0030 5 3
## 45 5 3 55.5453 5 3
## 46 5 3 51.5457 5 3
## 47 5 3 40.2068 5 3
## 48 10 3 56.6067 10 3
## 49 10 3 37.1715 10 3
## 50 10 3 47.9529 10 3
## 51 10 3 59.5799 10 3
## Data.E2FmpdM
## 1 111.09
## 2 39.03
## 3 117.98
## 4 140.66
## 5 148.08
## 6 74.47
## 7 48.58
## 8 91.93
## 9 140.05
## 10 65.35
## 11 78.28
## 12 76.48
## 13 75.77
## 14 82.12
## 15 77.00
## 16 103.03
## 17 49.78
## 18 84.33
## 19 55.57
## 20 99.40
## 21 62.55
## 22 107.61
## 23 50.40
## 24 39.74
## 25 56.89
## 26 103.37
## 27 53.42
## 28 69.74
## 29 85.77
## 30 62.01
## 31 70.61
## 32 50.34
## 33 66.41
## 34 67.59
## 35 58.36
## 36 48.38
## 37 128.49
## 38 55.52
## 39 90.56
## 40 57.26
## 41 35.36
## 42 58.25
## 43 164.04
## 44 116.08
## 45 61.69
## 46 79.48
## 47 81.86
## 48 65.40
## 49 57.79
## 50 108.00
## 51 50.26
tapply(MPDdoseFM$Data.E2FmpdM, MPDdoseFM$Group.E2FmpdM, sd)
## 10 2.5 5 Saline
## 19.31 41.06 15.99 34.83
tapply(MPDdoseFM$Data.E2FmpdM, MPDdoseFM$Day.E2FmpdM, sd)
## 1 2 3
## 32.95 19.88 33.99
MPDdoseFM$Group.E2FmpdM <- as.factor(MPDdoseFM$Group.E2FmpdM)
MPDdoseFM$Day.E2FmpdM <- as.factor(MPDdoseFM$Day.E2FmpdM)
MPDdoseFM$Data.E2FmpdM <- as.numeric(MPDdoseFM$Data.E2FmpdM)
aov.MPDdoseFM = aov(Data.E2FmpdM ~ Group.E2FmpdM * Day.E2FmpdM, data = MPDdoseFM)
summary(aov.MPDdoseFM)
## Df Sum Sq Mean Sq F value Pr(>F)
## Group.E2FmpdM 3 1655 552 0.56 0.65
## Day.E2FmpdM 2 3331 1666 1.69 0.20
## Group.E2FmpdM:Day.E2FmpdM 6 1984 331 0.33 0.91
## Residuals 39 38541 988
comparison <- aov(Data.E2FmpdM ~ Group.E2FmpdM * Day.E2FmpdM, data = MPDdoseFM)
TukeyHSD(comparison, "Group.E2FmpdM")
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = Data.E2FmpdM ~ Group.E2FmpdM * Day.E2FmpdM, data = MPDdoseFM)
##
## $Group.E2FmpdM
## diff lwr upr p adj
## 2.5-10 13.7452 -18.93 46.42 0.6741
## 5-10 4.8008 -29.64 39.24 0.9819
## Saline-10 12.7510 -21.69 47.19 0.7539
## 5-2.5 -8.9444 -41.61 23.73 0.8825
## Saline-2.5 -0.9942 -33.66 31.68 0.9998
## Saline-5 7.9502 -26.49 42.39 0.9252
TukeyHSD(comparison, "Day.E2FmpdM")
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = Data.E2FmpdM ~ Group.E2FmpdM * Day.E2FmpdM, data = MPDdoseFM)
##
## $Day.E2FmpdM
## diff lwr upr p adj
## 2-1 -19.646 -45.91 6.624 0.1758
## 3-1 -11.937 -38.21 14.332 0.5155
## 3-2 7.709 -18.56 33.978 0.7562
print(model.tables(aov.MPDdoseFM, "means"), digits = 17)
## Tables of means
## Grand mean
##
## 78.87
##
## Group.E2FmpdM
## 10 2.5 5
## 70.695118094166673 84.440313270666678 75.495925015000012
## rep 12.000000000000000 15.000000000000000 12.000000000000000
## Saline
## 83.446075344166658
## rep 12.000000000000000
##
## Day.E2FmpdM
## 1 2 3
## 89.395136659999991 69.749579369999992 77.458232234705889
## rep 17.000000000000000 17.000000000000000 17.000000000000000
##
## Group.E2FmpdM:Day.E2FmpdM
## Day.E2FmpdM
## Group.E2FmpdM 1 2 3
## 10 77.98 63.74 70.37
## rep 4.00 4.00 4.00
## 2.5 100.62 71.60 81.09
## rep 5.00 5.00 5.00
## 5 73.97 67.74 84.78
## rep 4.00 4.00 4.00
## Saline 102.19 75.46 72.68
## rep 4.00 4.00 4.00