MPD dose falafel wheel and motion hourly post hoc
MPD1hourlyFW <- read.csv("~/Desktop/exp 2 posthochourFalafelWheel.csv")
MPD1hourlyFW
## Group.1hFw Day.1hFw Data.1hFw Group.1.5hFw Day.1.5hFw Data.1.5hFw
## 1 Saline 1 327.869 Saline 1 15.556
## 2 Saline 1 5.734 Saline 1 1.646
## 3 Saline 1 14.985 Saline 1 40.424
## 4 Saline 1 3.165 Saline 1 1.096
## 5 Saline 2 49.180 Saline 2 2.963
## 6 Saline 2 0.000 Saline 2 16.462
## 7 Saline 2 28.900 Saline 2 31.303
## 8 Saline 2 0.000 Saline 2 0.000
## 9 Saline 3 8.197 Saline 3 3.333
## 10 Saline 3 0.000 Saline 3 0.000
## 11 Saline 3 21.541 Saline 3 18.705
## 12 Saline 3 0.000 Saline 3 0.000
## 13 2.5 mg/kg 1 1030.303 2.5 mg/kg 1 117.978
## 14 2.5 mg/kg 1 319.018 2.5 mg/kg 1 145.198
## 15 2.5 mg/kg 1 139.501 2.5 mg/kg 1 195.688
## 16 2.5 mg/kg 1 52.174 2.5 mg/kg 1 21.016
## 17 2.5 mg/kg 1 116.390 2.5 mg/kg 1 132.719
## 18 2.5 mg/kg 2 151.515 2.5 mg/kg 2 44.944
## 19 2.5 mg/kg 2 0.000 2.5 mg/kg 2 34.463
## 20 2.5 mg/kg 2 162.261 2.5 mg/kg 2 101.575
## 21 2.5 mg/kg 2 65.217 2.5 mg/kg 2 56.042
## 22 2.5 mg/kg 2 84.917 2.5 mg/kg 2 145.968
## 23 2.5 mg/kg 3 333.333 2.5 mg/kg 3 42.135
## 24 2.5 mg/kg 3 144.172 2.5 mg/kg 3 59.322
## 25 2.5 mg/kg 3 300.294 2.5 mg/kg 3 228.027
## 26 2.5 mg/kg 3 73.913 2.5 mg/kg 3 56.918
## 27 2.5 mg/kg 3 62.055 2.5 mg/kg 3 24.078
## 28 5.0 mg/kg 1 51.724 5.0 mg/kg 1 25.862
## 29 5.0 mg/kg 1 790.598 5.0 mg/kg 1 154.613
## 30 5.0 mg/kg 1 1230.964 5.0 mg/kg 1 229.462
## 31 5.0 mg/kg 1 3390.323 5.0 mg/kg 1 1405.544
## 32 5.0 mg/kg 2 206.897 5.0 mg/kg 2 103.448
## 33 5.0 mg/kg 2 380.342 5.0 mg/kg 2 73.982
## 34 5.0 mg/kg 2 253.807 5.0 mg/kg 2 79.792
## 35 5.0 mg/kg 2 816.129 5.0 mg/kg 2 348.049
## 36 5.0 mg/kg 3 68.966 5.0 mg/kg 3 86.207
## 37 5.0 mg/kg 3 540.598 5.0 mg/kg 3 105.154
## 38 5.0 mg/kg 3 456.853 5.0 mg/kg 3 92.540
## 39 5.0 mg/kg 3 835.484 5.0 mg/kg 3 286.448
## 40 10.0 mg/kg 1 1652.597 10.0 mg/kg 1 1035.206
## 41 10.0 mg/kg 1 1389.716 10.0 mg/kg 1 474.755
## 42 10.0 mg/kg 1 3305.654 10.0 mg/kg 1 663.499
## 43 10.0 mg/kg 1 3943.396 10.0 mg/kg 1 536.537
## 44 10.0 mg/kg 2 865.260 10.0 mg/kg 2 519.383
## 45 10.0 mg/kg 2 598.782 10.0 mg/kg 2 221.364
## 46 10.0 mg/kg 2 1351.590 10.0 mg/kg 2 250.655
## 47 10.0 mg/kg 2 3379.717 10.0 mg/kg 2 453.056
## 48 10.0 mg/kg 3 308.442 10.0 mg/kg 3 276.899
## 49 10.0 mg/kg 3 928.281 10.0 mg/kg 3 316.880
## 50 10.0 mg/kg 3 1349.823 10.0 mg/kg 3 250.655
## 51 10.0 mg/kg 3 1653.302 10.0 mg/kg 3 317.759
## Group.2hFw Day.2hFw Data.2hFw
## 1 Saline 1 7.7036
## 2 Saline 1 1.0274
## 3 Saline 1 43.0528
## 4 Saline 1 0.9158
## 5 Saline 2 29.7139
## 6 Saline 2 58.0479
## 7 Saline 2 45.2976
## 8 Saline 2 0.0000
## 9 Saline 3 1.6508
## 10 Saline 3 0.0000
## 11 Saline 3 27.5124
## 12 Saline 3 0.0000
## 13 2.5 mg/kg 1 96.7742
## 14 2.5 mg/kg 1 144.1013
## 15 2.5 mg/kg 1 203.9683
## 16 2.5 mg/kg 1 37.8049
## 17 2.5 mg/kg 1 93.6873
## 18 2.5 mg/kg 2 140.5530
## 19 2.5 mg/kg 2 25.7324
## 20 2.5 mg/kg 2 133.7302
## 21 2.5 mg/kg 2 39.6341
## 22 2.5 mg/kg 2 109.1345
## 23 2.5 mg/kg 3 213.1336
## 24 2.5 mg/kg 3 41.5677
## 25 2.5 mg/kg 3 272.6190
## 26 2.5 mg/kg 3 87.8049
## 27 2.5 mg/kg 3 60.8967
## 28 5.0 mg/kg 1 16.6667
## 29 5.0 mg/kg 1 87.0379
## 30 5.0 mg/kg 1 176.6618
## 31 5.0 mg/kg 1 759.4235
## 32 5.0 mg/kg 2 66.6667
## 33 5.0 mg/kg 2 42.1151
## 34 5.0 mg/kg 2 81.6474
## 35 5.0 mg/kg 2 192.3503
## 36 5.0 mg/kg 3 61.1111
## 37 5.0 mg/kg 3 59.1951
## 38 5.0 mg/kg 3 75.1445
## 39 5.0 mg/kg 3 183.4812
## 40 10.0 mg/kg 1 852.1739
## 41 10.0 mg/kg 1 506.0541
## 42 10.0 mg/kg 1 355.6772
## 43 10.0 mg/kg 1 456.4450
## 44 10.0 mg/kg 2 380.5797
## 45 10.0 mg/kg 2 243.9459
## 46 10.0 mg/kg 2 130.8140
## 47 10.0 mg/kg 2 389.4938
## 48 10.0 mg/kg 3 267.8261
## 49 10.0 mg/kg 3 318.1980
## 50 10.0 mg/kg 3 152.0178
## 51 10.0 mg/kg 3 238.0812
tapply(MPD1hourlyFW$Data.1hFw, MPD1hourlyFW$Group.1hFw, sd)
## 10.0 mg/kg 2.5 mg/kg 5.0 mg/kg Saline
## 1175.88 250.29 903.27 92.42
tapply(MPD1hourlyFW$Data.1hFw, MPD1hourlyFW$Day.1hFw, sd)
## 1 2 3
## 1311.5 836.4 499.7
MPD1hourlyFW$Group.1hFw <- as.factor(MPD1hourlyFW$Group.1hFw)
MPD1hourlyFW$Day.1hFw <- as.factor(MPD1hourlyFW$Day.1hFw)
MPD1hourlyFW$Data.1hFw <- as.numeric(MPD1hourlyFW$Data.1hFw)
aov.MPD1hourlyFW = aov(Data.1hFw ~ Group.1hFw * Day.1hFw, data = MPD1hourlyFW)
summary(aov.MPD1hourlyFW)
## Df Sum Sq Mean Sq F value Pr(>F)
## Group.1hFw 3 21545254 7181751 15.60 8.1e-07 ***
## Day.1hFw 2 3991010 1995505 4.33 0.02 *
## Group.1hFw:Day.1hFw 6 3207445 534574 1.16 0.35
## Residuals 39 17957150 460440
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
comparison <- aov(Data.1hFw ~ Group.1hFw * Day.1hFw, data = MPD1hourlyFW)
TukeyHSD(comparison, "Group.1hFw")
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = Data.1hFw ~ Group.1hFw * Day.1hFw, data = MPD1hourlyFW)
##
## $Group.1hFw
## diff lwr upr p adj
## 2.5 mg/kg-10.0 mg/kg -1524.9 -2230.1 -819.68 0.0000
## 5.0 mg/kg-10.0 mg/kg -975.3 -1718.7 -231.98 0.0059
## Saline-10.0 mg/kg -1688.9 -2432.3 -945.57 0.0000
## 5.0 mg/kg-2.5 mg/kg 549.6 -155.6 1254.75 0.1739
## Saline-2.5 mg/kg -164.0 -869.2 541.16 0.9237
## Saline-5.0 mg/kg -713.6 -1456.9 29.75 0.0639
TukeyHSD(comparison, "Day.1hFw")
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = Data.1hFw ~ Group.1hFw * Day.1hFw, data = MPD1hourlyFW)
##
## $Day.1hFw
## diff lwr upr p adj
## 2-1 -551.15 -1118 15.88 0.0583
## 3-1 -628.17 -1195 -61.13 0.0270
## 3-2 -77.02 -644 490.02 0.9415
print(model.tables(aov.MPD1hourlyFW, "means"), digits = 9)
## Tables of means
## Grand mean
##
## 651.8
##
## Group.1hFw
## 10.0 mg/kg 2.5 mg/kg 5.0 mg/kg Saline
## 1727.21336 202.337526 751.890377 38.2976014
## rep 12.00000 15.000000 12.000000 12.0000000
##
## Day.1hFw
## 1 2 3
## 1044.9477 493.794991 416.779602
## rep 17.0000 17.000000 17.000000
##
## Group.1hFw:Day.1hFw
## Day.1hFw
## Group.1hFw 1 2 3
## 10.0 mg/kg 2572.841 1548.837 1059.962
## rep 4.000 4.000 4.000
## 2.5 mg/kg 331.477 92.782 182.753
## rep 5.000 5.000 5.000
## 5.0 mg/kg 1365.902 414.294 475.475
## rep 4.000 4.000 4.000
## Saline 87.938 19.520 7.435
## rep 4.000 4.000 4.000
MPD1.5hourlyFW <- read.csv("~/Desktop/exp 2 posthochourFalafelWheel.csv")
MPD1.5hourlyFW
## Group.1hFw Day.1hFw Data.1hFw Group.1.5hFw Day.1.5hFw Data.1.5hFw
## 1 Saline 1 327.869 Saline 1 15.556
## 2 Saline 1 5.734 Saline 1 1.646
## 3 Saline 1 14.985 Saline 1 40.424
## 4 Saline 1 3.165 Saline 1 1.096
## 5 Saline 2 49.180 Saline 2 2.963
## 6 Saline 2 0.000 Saline 2 16.462
## 7 Saline 2 28.900 Saline 2 31.303
## 8 Saline 2 0.000 Saline 2 0.000
## 9 Saline 3 8.197 Saline 3 3.333
## 10 Saline 3 0.000 Saline 3 0.000
## 11 Saline 3 21.541 Saline 3 18.705
## 12 Saline 3 0.000 Saline 3 0.000
## 13 2.5 mg/kg 1 1030.303 2.5 mg/kg 1 117.978
## 14 2.5 mg/kg 1 319.018 2.5 mg/kg 1 145.198
## 15 2.5 mg/kg 1 139.501 2.5 mg/kg 1 195.688
## 16 2.5 mg/kg 1 52.174 2.5 mg/kg 1 21.016
## 17 2.5 mg/kg 1 116.390 2.5 mg/kg 1 132.719
## 18 2.5 mg/kg 2 151.515 2.5 mg/kg 2 44.944
## 19 2.5 mg/kg 2 0.000 2.5 mg/kg 2 34.463
## 20 2.5 mg/kg 2 162.261 2.5 mg/kg 2 101.575
## 21 2.5 mg/kg 2 65.217 2.5 mg/kg 2 56.042
## 22 2.5 mg/kg 2 84.917 2.5 mg/kg 2 145.968
## 23 2.5 mg/kg 3 333.333 2.5 mg/kg 3 42.135
## 24 2.5 mg/kg 3 144.172 2.5 mg/kg 3 59.322
## 25 2.5 mg/kg 3 300.294 2.5 mg/kg 3 228.027
## 26 2.5 mg/kg 3 73.913 2.5 mg/kg 3 56.918
## 27 2.5 mg/kg 3 62.055 2.5 mg/kg 3 24.078
## 28 5.0 mg/kg 1 51.724 5.0 mg/kg 1 25.862
## 29 5.0 mg/kg 1 790.598 5.0 mg/kg 1 154.613
## 30 5.0 mg/kg 1 1230.964 5.0 mg/kg 1 229.462
## 31 5.0 mg/kg 1 3390.323 5.0 mg/kg 1 1405.544
## 32 5.0 mg/kg 2 206.897 5.0 mg/kg 2 103.448
## 33 5.0 mg/kg 2 380.342 5.0 mg/kg 2 73.982
## 34 5.0 mg/kg 2 253.807 5.0 mg/kg 2 79.792
## 35 5.0 mg/kg 2 816.129 5.0 mg/kg 2 348.049
## 36 5.0 mg/kg 3 68.966 5.0 mg/kg 3 86.207
## 37 5.0 mg/kg 3 540.598 5.0 mg/kg 3 105.154
## 38 5.0 mg/kg 3 456.853 5.0 mg/kg 3 92.540
## 39 5.0 mg/kg 3 835.484 5.0 mg/kg 3 286.448
## 40 10.0 mg/kg 1 1652.597 10.0 mg/kg 1 1035.206
## 41 10.0 mg/kg 1 1389.716 10.0 mg/kg 1 474.755
## 42 10.0 mg/kg 1 3305.654 10.0 mg/kg 1 663.499
## 43 10.0 mg/kg 1 3943.396 10.0 mg/kg 1 536.537
## 44 10.0 mg/kg 2 865.260 10.0 mg/kg 2 519.383
## 45 10.0 mg/kg 2 598.782 10.0 mg/kg 2 221.364
## 46 10.0 mg/kg 2 1351.590 10.0 mg/kg 2 250.655
## 47 10.0 mg/kg 2 3379.717 10.0 mg/kg 2 453.056
## 48 10.0 mg/kg 3 308.442 10.0 mg/kg 3 276.899
## 49 10.0 mg/kg 3 928.281 10.0 mg/kg 3 316.880
## 50 10.0 mg/kg 3 1349.823 10.0 mg/kg 3 250.655
## 51 10.0 mg/kg 3 1653.302 10.0 mg/kg 3 317.759
## Group.2hFw Day.2hFw Data.2hFw
## 1 Saline 1 7.7036
## 2 Saline 1 1.0274
## 3 Saline 1 43.0528
## 4 Saline 1 0.9158
## 5 Saline 2 29.7139
## 6 Saline 2 58.0479
## 7 Saline 2 45.2976
## 8 Saline 2 0.0000
## 9 Saline 3 1.6508
## 10 Saline 3 0.0000
## 11 Saline 3 27.5124
## 12 Saline 3 0.0000
## 13 2.5 mg/kg 1 96.7742
## 14 2.5 mg/kg 1 144.1013
## 15 2.5 mg/kg 1 203.9683
## 16 2.5 mg/kg 1 37.8049
## 17 2.5 mg/kg 1 93.6873
## 18 2.5 mg/kg 2 140.5530
## 19 2.5 mg/kg 2 25.7324
## 20 2.5 mg/kg 2 133.7302
## 21 2.5 mg/kg 2 39.6341
## 22 2.5 mg/kg 2 109.1345
## 23 2.5 mg/kg 3 213.1336
## 24 2.5 mg/kg 3 41.5677
## 25 2.5 mg/kg 3 272.6190
## 26 2.5 mg/kg 3 87.8049
## 27 2.5 mg/kg 3 60.8967
## 28 5.0 mg/kg 1 16.6667
## 29 5.0 mg/kg 1 87.0379
## 30 5.0 mg/kg 1 176.6618
## 31 5.0 mg/kg 1 759.4235
## 32 5.0 mg/kg 2 66.6667
## 33 5.0 mg/kg 2 42.1151
## 34 5.0 mg/kg 2 81.6474
## 35 5.0 mg/kg 2 192.3503
## 36 5.0 mg/kg 3 61.1111
## 37 5.0 mg/kg 3 59.1951
## 38 5.0 mg/kg 3 75.1445
## 39 5.0 mg/kg 3 183.4812
## 40 10.0 mg/kg 1 852.1739
## 41 10.0 mg/kg 1 506.0541
## 42 10.0 mg/kg 1 355.6772
## 43 10.0 mg/kg 1 456.4450
## 44 10.0 mg/kg 2 380.5797
## 45 10.0 mg/kg 2 243.9459
## 46 10.0 mg/kg 2 130.8140
## 47 10.0 mg/kg 2 389.4938
## 48 10.0 mg/kg 3 267.8261
## 49 10.0 mg/kg 3 318.1980
## 50 10.0 mg/kg 3 152.0178
## 51 10.0 mg/kg 3 238.0812
tapply(MPD1.5hourlyFW$Data.1.5hFw, MPD1.5hourlyFW$Group.1.5hFw, sd)
## 10.0 mg/kg 2.5 mg/kg 5.0 mg/kg Saline
## 233.03 64.71 376.59 13.68
tapply(MPD1.5hourlyFW$Data.1.5hFw, MPD1.5hourlyFW$Day.1.5hFw, sd)
## 1 2 3
## 401.8 159.5 121.4
MPD1.5hourlyFW$Group.1.5hFw <- as.factor(MPD1.5hourlyFW$Group.1.5hFw)
MPD1.5hourlyFW$Day.1.5hFw <- as.factor(MPD1.5hourlyFW$Day.1.5hFw)
MPD1.5hourlyFW$Data.1.5hFw <- as.numeric(MPD1.5hourlyFW$Data.1.5hFw)
aov.MPD1.5hourlyFW = aov(Data.1.5hFw ~ Group.1.5hFw * Day.1.5hFw, data = MPD1.5hourlyFW)
summary(aov.MPD1.5hourlyFW)
## Df Sum Sq Mean Sq F value Pr(>F)
## Group.1.5hFw 3 1333773 444591 10.70 2.9e-05 ***
## Day.1.5hFw 2 326570 163285 3.93 0.028 *
## Group.1.5hFw:Day.1.5hFw 6 271000 45167 1.09 0.387
## Residuals 39 1620447 41550
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
comparison <- aov(Data.1.5hFw ~ Group.1.5hFw * Day.1.5hFw, data = MPD1.5hourlyFW)
TukeyHSD(comparison, "Group.1.5hFw")
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = Data.1.5hFw ~ Group.1.5hFw * Day.1.5hFw, data = MPD1.5hourlyFW)
##
## $Group.1.5hFw
## diff lwr upr p adj
## 2.5 mg/kg-10.0 mg/kg -349.32 -561.16 -137.5 0.0004
## 5.0 mg/kg-10.0 mg/kg -193.80 -417.10 29.5 0.1089
## Saline-10.0 mg/kg -432.10 -655.40 -208.8 0.0000
## 5.0 mg/kg-2.5 mg/kg 155.52 -56.32 367.4 0.2168
## Saline-2.5 mg/kg -82.78 -294.62 129.1 0.7222
## Saline-5.0 mg/kg -238.30 -461.60 -15.0 0.0326
TukeyHSD(comparison, "Day.1.5hFw")
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = Data.1.5hFw ~ Group.1.5hFw * Day.1.5hFw, data = MPD1.5hourlyFW)
##
## $Day.1.5hFw
## diff lwr upr p adj
## 2-1 -159.61 -329.9 10.728 0.0702
## 3-1 -178.34 -348.7 -8.001 0.0384
## 3-2 -18.73 -189.1 151.608 0.9613
print(model.tables(aov.MPD1.5hourlyFW, "means"), digits = 9)
## Tables of means
## Grand mean
##
## 193
##
## Group.1.5hFw
## 10.0 mg/kg 2.5 mg/kg 5.0 mg/kg Saline
## 443.054016 93.7379936 249.258451 10.9573318
## rep 12.000000 15.0000000 12.000000 12.0000000
##
## Day.1.5hFw
## 1 2 3
## 305.694017 146.085299 127.356418
## rep 17.000000 17.000000 17.000000
##
## Group.1.5hFw:Day.1.5hFw
## Day.1.5hFw
## Group.1.5hFw 1 2 3
## 10.0 mg/kg 677.499 361.115 290.548
## rep 4.000 4.000 4.000
## 2.5 mg/kg 122.520 76.598 82.096
## rep 5.000 5.000 5.000
## 5.0 mg/kg 453.870 151.318 142.587
## rep 4.000 4.000 4.000
## Saline 14.680 12.682 5.509
## rep 4.000 4.000 4.000
MPD2hourlyFW <- read.csv("~/Desktop/exp 2 posthochourFalafelWheel.csv")
MPD2hourlyFW
## Group.1hFw Day.1hFw Data.1hFw Group.1.5hFw Day.1.5hFw Data.1.5hFw
## 1 Saline 1 327.869 Saline 1 15.556
## 2 Saline 1 5.734 Saline 1 1.646
## 3 Saline 1 14.985 Saline 1 40.424
## 4 Saline 1 3.165 Saline 1 1.096
## 5 Saline 2 49.180 Saline 2 2.963
## 6 Saline 2 0.000 Saline 2 16.462
## 7 Saline 2 28.900 Saline 2 31.303
## 8 Saline 2 0.000 Saline 2 0.000
## 9 Saline 3 8.197 Saline 3 3.333
## 10 Saline 3 0.000 Saline 3 0.000
## 11 Saline 3 21.541 Saline 3 18.705
## 12 Saline 3 0.000 Saline 3 0.000
## 13 2.5 mg/kg 1 1030.303 2.5 mg/kg 1 117.978
## 14 2.5 mg/kg 1 319.018 2.5 mg/kg 1 145.198
## 15 2.5 mg/kg 1 139.501 2.5 mg/kg 1 195.688
## 16 2.5 mg/kg 1 52.174 2.5 mg/kg 1 21.016
## 17 2.5 mg/kg 1 116.390 2.5 mg/kg 1 132.719
## 18 2.5 mg/kg 2 151.515 2.5 mg/kg 2 44.944
## 19 2.5 mg/kg 2 0.000 2.5 mg/kg 2 34.463
## 20 2.5 mg/kg 2 162.261 2.5 mg/kg 2 101.575
## 21 2.5 mg/kg 2 65.217 2.5 mg/kg 2 56.042
## 22 2.5 mg/kg 2 84.917 2.5 mg/kg 2 145.968
## 23 2.5 mg/kg 3 333.333 2.5 mg/kg 3 42.135
## 24 2.5 mg/kg 3 144.172 2.5 mg/kg 3 59.322
## 25 2.5 mg/kg 3 300.294 2.5 mg/kg 3 228.027
## 26 2.5 mg/kg 3 73.913 2.5 mg/kg 3 56.918
## 27 2.5 mg/kg 3 62.055 2.5 mg/kg 3 24.078
## 28 5.0 mg/kg 1 51.724 5.0 mg/kg 1 25.862
## 29 5.0 mg/kg 1 790.598 5.0 mg/kg 1 154.613
## 30 5.0 mg/kg 1 1230.964 5.0 mg/kg 1 229.462
## 31 5.0 mg/kg 1 3390.323 5.0 mg/kg 1 1405.544
## 32 5.0 mg/kg 2 206.897 5.0 mg/kg 2 103.448
## 33 5.0 mg/kg 2 380.342 5.0 mg/kg 2 73.982
## 34 5.0 mg/kg 2 253.807 5.0 mg/kg 2 79.792
## 35 5.0 mg/kg 2 816.129 5.0 mg/kg 2 348.049
## 36 5.0 mg/kg 3 68.966 5.0 mg/kg 3 86.207
## 37 5.0 mg/kg 3 540.598 5.0 mg/kg 3 105.154
## 38 5.0 mg/kg 3 456.853 5.0 mg/kg 3 92.540
## 39 5.0 mg/kg 3 835.484 5.0 mg/kg 3 286.448
## 40 10.0 mg/kg 1 1652.597 10.0 mg/kg 1 1035.206
## 41 10.0 mg/kg 1 1389.716 10.0 mg/kg 1 474.755
## 42 10.0 mg/kg 1 3305.654 10.0 mg/kg 1 663.499
## 43 10.0 mg/kg 1 3943.396 10.0 mg/kg 1 536.537
## 44 10.0 mg/kg 2 865.260 10.0 mg/kg 2 519.383
## 45 10.0 mg/kg 2 598.782 10.0 mg/kg 2 221.364
## 46 10.0 mg/kg 2 1351.590 10.0 mg/kg 2 250.655
## 47 10.0 mg/kg 2 3379.717 10.0 mg/kg 2 453.056
## 48 10.0 mg/kg 3 308.442 10.0 mg/kg 3 276.899
## 49 10.0 mg/kg 3 928.281 10.0 mg/kg 3 316.880
## 50 10.0 mg/kg 3 1349.823 10.0 mg/kg 3 250.655
## 51 10.0 mg/kg 3 1653.302 10.0 mg/kg 3 317.759
## Group.2hFw Day.2hFw Data.2hFw
## 1 Saline 1 7.7036
## 2 Saline 1 1.0274
## 3 Saline 1 43.0528
## 4 Saline 1 0.9158
## 5 Saline 2 29.7139
## 6 Saline 2 58.0479
## 7 Saline 2 45.2976
## 8 Saline 2 0.0000
## 9 Saline 3 1.6508
## 10 Saline 3 0.0000
## 11 Saline 3 27.5124
## 12 Saline 3 0.0000
## 13 2.5 mg/kg 1 96.7742
## 14 2.5 mg/kg 1 144.1013
## 15 2.5 mg/kg 1 203.9683
## 16 2.5 mg/kg 1 37.8049
## 17 2.5 mg/kg 1 93.6873
## 18 2.5 mg/kg 2 140.5530
## 19 2.5 mg/kg 2 25.7324
## 20 2.5 mg/kg 2 133.7302
## 21 2.5 mg/kg 2 39.6341
## 22 2.5 mg/kg 2 109.1345
## 23 2.5 mg/kg 3 213.1336
## 24 2.5 mg/kg 3 41.5677
## 25 2.5 mg/kg 3 272.6190
## 26 2.5 mg/kg 3 87.8049
## 27 2.5 mg/kg 3 60.8967
## 28 5.0 mg/kg 1 16.6667
## 29 5.0 mg/kg 1 87.0379
## 30 5.0 mg/kg 1 176.6618
## 31 5.0 mg/kg 1 759.4235
## 32 5.0 mg/kg 2 66.6667
## 33 5.0 mg/kg 2 42.1151
## 34 5.0 mg/kg 2 81.6474
## 35 5.0 mg/kg 2 192.3503
## 36 5.0 mg/kg 3 61.1111
## 37 5.0 mg/kg 3 59.1951
## 38 5.0 mg/kg 3 75.1445
## 39 5.0 mg/kg 3 183.4812
## 40 10.0 mg/kg 1 852.1739
## 41 10.0 mg/kg 1 506.0541
## 42 10.0 mg/kg 1 355.6772
## 43 10.0 mg/kg 1 456.4450
## 44 10.0 mg/kg 2 380.5797
## 45 10.0 mg/kg 2 243.9459
## 46 10.0 mg/kg 2 130.8140
## 47 10.0 mg/kg 2 389.4938
## 48 10.0 mg/kg 3 267.8261
## 49 10.0 mg/kg 3 318.1980
## 50 10.0 mg/kg 3 152.0178
## 51 10.0 mg/kg 3 238.0812
tapply(MPD2hourlyFW$Data.2hFw, MPD2hourlyFW$Group.2hFw, sd)
## 10.0 mg/kg 2.5 mg/kg 5.0 mg/kg Saline
## 192.70 72.55 200.40 21.59
tapply(MPD2hourlyFW$Data.2hFw, MPD2hourlyFW$Day.2hFw, sd)
## 1 2 3
## 268.1 117.1 107.0
MPD2hourlyFW$Group.2hFw <- as.factor(MPD2hourlyFW$Group.2hFw)
MPD2hourlyFW$Day.2hFw <- as.factor(MPD2hourlyFW$Day.2hFw)
MPD2hourlyFW$Data.2hFw <- as.numeric(MPD2hourlyFW$Data.2hFw)
aov.MPD2hourlyFW = aov(Data.2hFw ~ Group.2hFw * Day.2hFw, data = MPD2hourlyFW)
summary(aov.MPD2hourlyFW)
## Df Sum Sq Mean Sq F value Pr(>F)
## Group.2hFw 3 744148 248049 15.09 1.1e-06 ***
## Day.2hFw 2 120764 60382 3.67 0.035 *
## Group.2hFw:Day.2hFw 6 167138 27856 1.69 0.148
## Residuals 39 641159 16440
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
comparison <- aov(Data.2hFw ~ Group.2hFw * Day.2hFw, data = MPD2hourlyFW)
TukeyHSD(comparison, "Group.2hFw")
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = Data.2hFw ~ Group.2hFw * Day.2hFw, data = MPD2hourlyFW)
##
## $Group.2hFw
## diff lwr upr p adj
## 2.5 mg/kg-10.0 mg/kg -244.20 -377.45 -110.947 0.0001
## 5.0 mg/kg-10.0 mg/kg -207.48 -347.94 -67.023 0.0017
## Saline-10.0 mg/kg -339.70 -480.16 -199.238 0.0000
## 5.0 mg/kg-2.5 mg/kg 36.72 -96.54 169.968 0.8806
## Saline-2.5 mg/kg -95.50 -228.75 37.753 0.2352
## Saline-5.0 mg/kg -132.21 -272.68 8.246 0.0714
TukeyHSD(comparison, "Day.2hFw")
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = Data.2hFw ~ Group.2hFw * Day.2hFw, data = MPD2hourlyFW)
##
## $Day.2hFw
## diff lwr upr p adj
## 2-1 -101.748 -208.9 5.397 0.0657
## 3-1 -104.643 -211.8 2.502 0.0568
## 3-2 -2.895 -110.0 104.250 0.9976
print(model.tables(aov.MPD2hourlyFW, "means"), digits = 9)
## Tables of means
## Grand mean
##
## 157
##
## Group.2hFw
## 10.0 mg/kg 2.5 mg/kg 5.0 mg/kg Saline
## 357.608885 113.409474 150.125112 17.9101757
## rep 12.000000 15.000000 12.000000 12.0000000
##
## Day.2hFw
## 1 2 3
## 225.833859 124.085673 121.190596
## rep 17.000000 17.000000 17.000000
##
## Group.2hFw:Day.2hFw
## Day.2hFw
## Group.2hFw 1 2 3
## 10.0 mg/kg 542.588 286.208 244.031
## rep 4.000 4.000 4.000
## 2.5 mg/kg 115.267 89.757 135.204
## rep 5.000 5.000 5.000
## 5.0 mg/kg 259.947 95.695 94.733
## rep 4.000 4.000 4.000
## Saline 13.175 33.265 7.291
## rep 4.000 4.000 4.000
MPDhourly1FM <- read.csv("~/Desktop/exp 2 posthochourFalafelMotion.csv")
MPDhourly1FM
## Group.1hFm Day.1hFm Data.1hFm Group.1.5hFm Day.1.5hFm Data.1.5hFm
## 1 Saline 1 70.423 Saline 1 62.75
## 2 Saline 1 30.749 Saline 1 31.09
## 3 Saline 1 72.770 Saline 1 98.46
## 4 Saline 1 185.185 Saline 1 84.62
## 5 Saline 2 32.864 Saline 2 53.92
## 6 Saline 2 5.348 Saline 2 26.65
## 7 Saline 2 56.338 Saline 2 60.00
## 8 Saline 2 37.037 Saline 2 30.77
## 9 Saline 3 18.779 Saline 3 38.24
## 10 Saline 3 0.000 Saline 3 12.69
## 11 Saline 3 56.338 Saline 3 61.54
## 12 Saline 3 55.556 Saline 3 38.46
## 13 2.5 mg/kg 1 274.809 2.5 mg/kg 1 131.52
## 14 2.5 mg/kg 1 49.559 2.5 mg/kg 1 56.54
## 15 2.5 mg/kg 1 74.703 2.5 mg/kg 1 57.06
## 16 2.5 mg/kg 1 374.150 2.5 mg/kg 1 198.74
## 17 2.5 mg/kg 1 147.601 2.5 mg/kg 1 155.17
## 18 2.5 mg/kg 2 137.405 2.5 mg/kg 2 149.27
## 19 2.5 mg/kg 2 13.216 2.5 mg/kg 2 47.00
## 20 2.5 mg/kg 2 61.121 2.5 mg/kg 2 52.05
## 21 2.5 mg/kg 2 23.810 2.5 mg/kg 2 53.63
## 22 2.5 mg/kg 2 16.605 2.5 mg/kg 2 65.27
## 23 2.5 mg/kg 3 374.046 2.5 mg/kg 3 195.20
## 24 2.5 mg/kg 3 60.573 2.5 mg/kg 3 48.37
## 25 2.5 mg/kg 3 21.222 2.5 mg/kg 3 26.53
## 26 2.5 mg/kg 3 95.238 2.5 mg/kg 3 72.56
## 27 2.5 mg/kg 3 103.321 2.5 mg/kg 3 76.35
## 28 5.0 mg/kg 1 312.844 5.0 mg/kg 1 246.68
## 29 5.0 mg/kg 1 574.419 5.0 mg/kg 1 176.06
## 30 5.0 mg/kg 1 144.467 5.0 mg/kg 1 94.32
## 31 5.0 mg/kg 1 215.328 5.0 mg/kg 1 165.45
## 32 5.0 mg/kg 2 179.817 5.0 mg/kg 2 132.69
## 33 5.0 mg/kg 2 469.767 5.0 mg/kg 2 145.97
## 34 5.0 mg/kg 2 119.877 5.0 mg/kg 2 81.25
## 35 5.0 mg/kg 2 215.328 5.0 mg/kg 2 120.00
## 36 5.0 mg/kg 3 321.101 5.0 mg/kg 3 254.52
## 37 5.0 mg/kg 3 444.186 5.0 mg/kg 3 126.12
## 38 5.0 mg/kg 3 119.877 5.0 mg/kg 3 97.73
## 39 5.0 mg/kg 3 379.562 5.0 mg/kg 3 250.91
## 40 10.0 mg/kg 1 217.277 10.0 mg/kg 1 234.38
## 41 10.0 mg/kg 1 403.955 10.0 mg/kg 1 250.73
## 42 10.0 mg/kg 1 307.047 10.0 mg/kg 1 211.61
## 43 10.0 mg/kg 1 15.432 10.0 mg/kg 1 16.42
## 44 10.0 mg/kg 2 340.314 10.0 mg/kg 2 262.50
## 45 10.0 mg/kg 2 220.339 10.0 mg/kg 2 151.60
## 46 10.0 mg/kg 2 384.228 10.0 mg/kg 2 180.40
## 47 10.0 mg/kg 2 30.864 10.0 mg/kg 2 73.13
## 48 10.0 mg/kg 3 86.387 10.0 mg/kg 3 92.19
## 49 10.0 mg/kg 3 316.384 10.0 mg/kg 3 183.67
## 50 10.0 mg/kg 3 803.691 10.0 mg/kg 3 324.59
## 51 10.0 mg/kg 3 21.605 10.0 mg/kg 3 55.22
## Group.2hFm Day.2hFm Data.2hFm
## 1 Saline 1 48.343
## 2 Saline 1 24.000
## 3 Saline 1 98.780
## 4 Saline 1 102.151
## 5 Saline 2 58.702
## 6 Saline 2 44.444
## 7 Saline 2 70.732
## 8 Saline 2 21.505
## 9 Saline 3 27.624
## 10 Saline 3 8.889
## 11 Saline 3 74.390
## 12 Saline 3 26.882
## 13 2.5 mg/kg 1 83.236
## 14 2.5 mg/kg 1 57.143
## 15 2.5 mg/kg 1 50.933
## 16 2.5 mg/kg 1 169.625
## 17 2.5 mg/kg 1 155.914
## 18 2.5 mg/kg 2 136.205
## 19 2.5 mg/kg 2 47.429
## 20 2.5 mg/kg 2 45.101
## 21 2.5 mg/kg 2 57.199
## 22 2.5 mg/kg 2 79.749
## 23 2.5 mg/kg 3 181.024
## 24 2.5 mg/kg 3 43.429
## 25 2.5 mg/kg 3 40.047
## 26 2.5 mg/kg 3 98.619
## 27 2.5 mg/kg 3 129.032
## 28 5.0 mg/kg 1 185.554
## 29 5.0 mg/kg 1 107.675
## 30 5.0 mg/kg 1 76.692
## 31 5.0 mg/kg 1 135.910
## 32 5.0 mg/kg 2 112.457
## 33 5.0 mg/kg 2 124.069
## 34 5.0 mg/kg 2 81.834
## 35 5.0 mg/kg 2 100.998
## 36 5.0 mg/kg 3 198.097
## 37 5.0 mg/kg 3 74.888
## 38 5.0 mg/kg 3 89.117
## 39 5.0 mg/kg 3 223.192
## 40 10.0 mg/kg 1 225.537
## 41 10.0 mg/kg 1 221.963
## 42 10.0 mg/kg 1 158.750
## 43 10.0 mg/kg 1 41.037
## 44 10.0 mg/kg 2 200.477
## 45 10.0 mg/kg 2 136.682
## 46 10.0 mg/kg 2 130.417
## 47 10.0 mg/kg 2 90.713
## 48 10.0 mg/kg 3 97.852
## 49 10.0 mg/kg 3 171.729
## 50 10.0 mg/kg 3 232.500
## 51 10.0 mg/kg 3 69.114
tapply(MPDhourly1FM$Data.1hFm, MPDhourly1FM$Group.1hFm, sd)
## 10.0 mg/kg 2.5 mg/kg 5.0 mg/kg Saline
## 222.84 122.51 150.40 48.23
tapply(MPDhourly1FM$Data.1hFm, MPDhourly1FM$Day.1hFm, sd)
## 1 2 3
## 154.6 144.2 216.7
MPDhourly1FM$Group.1hFm <- as.factor(MPDhourly1FM$Group.1hFm)
MPDhourly1FM$Day.1hFm <- as.factor(MPDhourly1FM$Day.1hFm)
MPDhourly1FM$Data.1hFm <- as.numeric(MPDhourly1FM$Data.1hFm)
aov.MPDhourly1FM = aov(Data.1hFm ~ Group.1hFm * Day.1hFm, data = MPDhourly1FM)
summary(aov.MPDhourly1FM)
## Df Sum Sq Mean Sq F value Pr(>F)
## Group.1hFm 3 478030 159343 6.53 0.0011 **
## Day.1hFm 2 42699 21350 0.87 0.4252
## Group.1hFm:Day.1hFm 6 35709 5952 0.24 0.9590
## Residuals 39 952346 24419
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
comparison <- aov(Data.1hFm ~ Group.1hFm * Day.1hFm, data = MPDhourly1FM)
TukeyHSD(comparison, "Group.1hFm")
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = Data.1hFm ~ Group.1hFm * Day.1hFm, data = MPDhourly1FM)
##
## $Group.1hFm
## diff lwr upr p adj
## 2.5 mg/kg-10.0 mg/kg -140.47 -302.870 21.93 0.1107
## 5.0 mg/kg-10.0 mg/kg 29.09 -142.099 200.27 0.9681
## Saline-10.0 mg/kg -210.51 -381.698 -39.33 0.0107
## 5.0 mg/kg-2.5 mg/kg 169.56 7.154 331.96 0.0378
## Saline-2.5 mg/kg -70.04 -232.445 92.36 0.6568
## Saline-5.0 mg/kg -239.60 -410.785 -68.41 0.0030
TukeyHSD(comparison, "Day.1hFm")
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = Data.1hFm ~ Group.1hFm * Day.1hFm, data = MPDhourly1FM)
##
## $Day.1hFm
## diff lwr upr p adj
## 2-1 -66.26 -196.84 64.32 0.4394
## 3-1 -11.34 -141.93 119.24 0.9756
## 3-2 54.92 -75.67 185.50 0.5660
print(model.tables(aov.MPDhourly1FM, "means"), digits = 9)
## Tables of means
## Grand mean
##
## 178.3
##
## Group.1hFm
## 10.0 mg/kg 2.5 mg/kg 5.0 mg/kg Saline
## 262.293725 121.825224 291.381154 51.7821476
## rep 12.000000 15.000000 12.000000 12.0000000
##
## Day.1hFm
## 1 2 3
## 204.159922 137.898656 192.815696
## rep 17.000000 17.000000 17.000000
##
## Group.1hFm:Day.1hFm
## Day.1hFm
## Group.1hFm 1 2 3
## 10.0 mg/kg 235.93 243.94 307.02
## rep 4.00 4.00 4.00
## 2.5 mg/kg 184.16 50.43 130.88
## rep 5.00 5.00 5.00
## 5.0 mg/kg 311.76 246.20 316.18
## rep 4.00 4.00 4.00
## Saline 89.78 32.90 32.67
## rep 4.00 4.00 4.00
MPDhourly1.5FM <- read.csv("~/Desktop/exp 2 posthochourFalafelMotion.csv")
MPDhourly1.5FM
## Group.1hFm Day.1hFm Data.1hFm Group.1.5hFm Day.1.5hFm Data.1.5hFm
## 1 Saline 1 70.423 Saline 1 62.75
## 2 Saline 1 30.749 Saline 1 31.09
## 3 Saline 1 72.770 Saline 1 98.46
## 4 Saline 1 185.185 Saline 1 84.62
## 5 Saline 2 32.864 Saline 2 53.92
## 6 Saline 2 5.348 Saline 2 26.65
## 7 Saline 2 56.338 Saline 2 60.00
## 8 Saline 2 37.037 Saline 2 30.77
## 9 Saline 3 18.779 Saline 3 38.24
## 10 Saline 3 0.000 Saline 3 12.69
## 11 Saline 3 56.338 Saline 3 61.54
## 12 Saline 3 55.556 Saline 3 38.46
## 13 2.5 mg/kg 1 274.809 2.5 mg/kg 1 131.52
## 14 2.5 mg/kg 1 49.559 2.5 mg/kg 1 56.54
## 15 2.5 mg/kg 1 74.703 2.5 mg/kg 1 57.06
## 16 2.5 mg/kg 1 374.150 2.5 mg/kg 1 198.74
## 17 2.5 mg/kg 1 147.601 2.5 mg/kg 1 155.17
## 18 2.5 mg/kg 2 137.405 2.5 mg/kg 2 149.27
## 19 2.5 mg/kg 2 13.216 2.5 mg/kg 2 47.00
## 20 2.5 mg/kg 2 61.121 2.5 mg/kg 2 52.05
## 21 2.5 mg/kg 2 23.810 2.5 mg/kg 2 53.63
## 22 2.5 mg/kg 2 16.605 2.5 mg/kg 2 65.27
## 23 2.5 mg/kg 3 374.046 2.5 mg/kg 3 195.20
## 24 2.5 mg/kg 3 60.573 2.5 mg/kg 3 48.37
## 25 2.5 mg/kg 3 21.222 2.5 mg/kg 3 26.53
## 26 2.5 mg/kg 3 95.238 2.5 mg/kg 3 72.56
## 27 2.5 mg/kg 3 103.321 2.5 mg/kg 3 76.35
## 28 5.0 mg/kg 1 312.844 5.0 mg/kg 1 246.68
## 29 5.0 mg/kg 1 574.419 5.0 mg/kg 1 176.06
## 30 5.0 mg/kg 1 144.467 5.0 mg/kg 1 94.32
## 31 5.0 mg/kg 1 215.328 5.0 mg/kg 1 165.45
## 32 5.0 mg/kg 2 179.817 5.0 mg/kg 2 132.69
## 33 5.0 mg/kg 2 469.767 5.0 mg/kg 2 145.97
## 34 5.0 mg/kg 2 119.877 5.0 mg/kg 2 81.25
## 35 5.0 mg/kg 2 215.328 5.0 mg/kg 2 120.00
## 36 5.0 mg/kg 3 321.101 5.0 mg/kg 3 254.52
## 37 5.0 mg/kg 3 444.186 5.0 mg/kg 3 126.12
## 38 5.0 mg/kg 3 119.877 5.0 mg/kg 3 97.73
## 39 5.0 mg/kg 3 379.562 5.0 mg/kg 3 250.91
## 40 10.0 mg/kg 1 217.277 10.0 mg/kg 1 234.38
## 41 10.0 mg/kg 1 403.955 10.0 mg/kg 1 250.73
## 42 10.0 mg/kg 1 307.047 10.0 mg/kg 1 211.61
## 43 10.0 mg/kg 1 15.432 10.0 mg/kg 1 16.42
## 44 10.0 mg/kg 2 340.314 10.0 mg/kg 2 262.50
## 45 10.0 mg/kg 2 220.339 10.0 mg/kg 2 151.60
## 46 10.0 mg/kg 2 384.228 10.0 mg/kg 2 180.40
## 47 10.0 mg/kg 2 30.864 10.0 mg/kg 2 73.13
## 48 10.0 mg/kg 3 86.387 10.0 mg/kg 3 92.19
## 49 10.0 mg/kg 3 316.384 10.0 mg/kg 3 183.67
## 50 10.0 mg/kg 3 803.691 10.0 mg/kg 3 324.59
## 51 10.0 mg/kg 3 21.605 10.0 mg/kg 3 55.22
## Group.2hFm Day.2hFm Data.2hFm
## 1 Saline 1 48.343
## 2 Saline 1 24.000
## 3 Saline 1 98.780
## 4 Saline 1 102.151
## 5 Saline 2 58.702
## 6 Saline 2 44.444
## 7 Saline 2 70.732
## 8 Saline 2 21.505
## 9 Saline 3 27.624
## 10 Saline 3 8.889
## 11 Saline 3 74.390
## 12 Saline 3 26.882
## 13 2.5 mg/kg 1 83.236
## 14 2.5 mg/kg 1 57.143
## 15 2.5 mg/kg 1 50.933
## 16 2.5 mg/kg 1 169.625
## 17 2.5 mg/kg 1 155.914
## 18 2.5 mg/kg 2 136.205
## 19 2.5 mg/kg 2 47.429
## 20 2.5 mg/kg 2 45.101
## 21 2.5 mg/kg 2 57.199
## 22 2.5 mg/kg 2 79.749
## 23 2.5 mg/kg 3 181.024
## 24 2.5 mg/kg 3 43.429
## 25 2.5 mg/kg 3 40.047
## 26 2.5 mg/kg 3 98.619
## 27 2.5 mg/kg 3 129.032
## 28 5.0 mg/kg 1 185.554
## 29 5.0 mg/kg 1 107.675
## 30 5.0 mg/kg 1 76.692
## 31 5.0 mg/kg 1 135.910
## 32 5.0 mg/kg 2 112.457
## 33 5.0 mg/kg 2 124.069
## 34 5.0 mg/kg 2 81.834
## 35 5.0 mg/kg 2 100.998
## 36 5.0 mg/kg 3 198.097
## 37 5.0 mg/kg 3 74.888
## 38 5.0 mg/kg 3 89.117
## 39 5.0 mg/kg 3 223.192
## 40 10.0 mg/kg 1 225.537
## 41 10.0 mg/kg 1 221.963
## 42 10.0 mg/kg 1 158.750
## 43 10.0 mg/kg 1 41.037
## 44 10.0 mg/kg 2 200.477
## 45 10.0 mg/kg 2 136.682
## 46 10.0 mg/kg 2 130.417
## 47 10.0 mg/kg 2 90.713
## 48 10.0 mg/kg 3 97.852
## 49 10.0 mg/kg 3 171.729
## 50 10.0 mg/kg 3 232.500
## 51 10.0 mg/kg 3 69.114
tapply(MPDhourly1.5FM$Data.1.5hFm, MPDhourly1.5FM$Group.1.5hFm, sd)
## 10.0 mg/kg 2.5 mg/kg 5.0 mg/kg Saline
## 94.23 57.30 62.54 25.03
tapply(MPDhourly1.5FM$Data.1.5hFm, MPDhourly1.5FM$Day.1.5hFm, sd)
## 1 2 3
## 77.67 63.72 92.93
MPDhourly1.5FM$Group.1.5hFm <- as.factor(MPDhourly1.5FM$Group.1.5hFm)
MPDhourly1.5FM$Day.1.5hFm <- as.factor(MPDhourly1.5FM$Day.1hFm)
MPDhourly1.5FM$Data.1.5hFm <- as.numeric(MPDhourly1.5FM$Data.1hFm)
aov.MPDhourly1.5FM = aov(Data.1.5hFm ~ Group.1.5hFm * Day.1.5hFm, data = MPDhourly1.5FM)
summary(aov.MPDhourly1.5FM)
## Df Sum Sq Mean Sq F value Pr(>F)
## Group.1.5hFm 3 478030 159343 6.53 0.0011 **
## Day.1.5hFm 2 42699 21350 0.87 0.4252
## Group.1.5hFm:Day.1.5hFm 6 35709 5952 0.24 0.9590
## Residuals 39 952346 24419
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
comparison <- aov(Data.1.5hFm ~ Group.1.5hFm * Day.1.5hFm, data = MPDhourly1.5FM)
TukeyHSD(comparison, "Group.1.5hFm")
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = Data.1.5hFm ~ Group.1.5hFm * Day.1.5hFm, data = MPDhourly1.5FM)
##
## $Group.1.5hFm
## diff lwr upr p adj
## 2.5 mg/kg-10.0 mg/kg -140.47 -302.870 21.93 0.1107
## 5.0 mg/kg-10.0 mg/kg 29.09 -142.099 200.27 0.9681
## Saline-10.0 mg/kg -210.51 -381.698 -39.33 0.0107
## 5.0 mg/kg-2.5 mg/kg 169.56 7.154 331.96 0.0378
## Saline-2.5 mg/kg -70.04 -232.445 92.36 0.6568
## Saline-5.0 mg/kg -239.60 -410.785 -68.41 0.0030
TukeyHSD(comparison, "Day.1.5hFm")
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = Data.1.5hFm ~ Group.1.5hFm * Day.1.5hFm, data = MPDhourly1.5FM)
##
## $Day.1.5hFm
## diff lwr upr p adj
## 2-1 -66.26 -196.84 64.32 0.4394
## 3-1 -11.34 -141.93 119.24 0.9756
## 3-2 54.92 -75.67 185.50 0.5660
print(model.tables(aov.MPDhourly1.5FM, "means"), digits = 9)
## Tables of means
## Grand mean
##
## 178.3
##
## Group.1.5hFm
## 10.0 mg/kg 2.5 mg/kg 5.0 mg/kg Saline
## 262.293725 121.825224 291.381154 51.7821476
## rep 12.000000 15.000000 12.000000 12.0000000
##
## Day.1.5hFm
## 1 2 3
## 204.159922 137.898656 192.815696
## rep 17.000000 17.000000 17.000000
##
## Group.1.5hFm:Day.1.5hFm
## Day.1.5hFm
## Group.1.5hFm 1 2 3
## 10.0 mg/kg 235.93 243.94 307.02
## rep 4.00 4.00 4.00
## 2.5 mg/kg 184.16 50.43 130.88
## rep 5.00 5.00 5.00
## 5.0 mg/kg 311.76 246.20 316.18
## rep 4.00 4.00 4.00
## Saline 89.78 32.90 32.67
## rep 4.00 4.00 4.00
MPDhourly2FM <- read.csv("~/Desktop/exp 2 posthochourFalafelMotion.csv")
MPDhourly2FM
## Group.1hFm Day.1hFm Data.1hFm Group.1.5hFm Day.1.5hFm Data.1.5hFm
## 1 Saline 1 70.423 Saline 1 62.75
## 2 Saline 1 30.749 Saline 1 31.09
## 3 Saline 1 72.770 Saline 1 98.46
## 4 Saline 1 185.185 Saline 1 84.62
## 5 Saline 2 32.864 Saline 2 53.92
## 6 Saline 2 5.348 Saline 2 26.65
## 7 Saline 2 56.338 Saline 2 60.00
## 8 Saline 2 37.037 Saline 2 30.77
## 9 Saline 3 18.779 Saline 3 38.24
## 10 Saline 3 0.000 Saline 3 12.69
## 11 Saline 3 56.338 Saline 3 61.54
## 12 Saline 3 55.556 Saline 3 38.46
## 13 2.5 mg/kg 1 274.809 2.5 mg/kg 1 131.52
## 14 2.5 mg/kg 1 49.559 2.5 mg/kg 1 56.54
## 15 2.5 mg/kg 1 74.703 2.5 mg/kg 1 57.06
## 16 2.5 mg/kg 1 374.150 2.5 mg/kg 1 198.74
## 17 2.5 mg/kg 1 147.601 2.5 mg/kg 1 155.17
## 18 2.5 mg/kg 2 137.405 2.5 mg/kg 2 149.27
## 19 2.5 mg/kg 2 13.216 2.5 mg/kg 2 47.00
## 20 2.5 mg/kg 2 61.121 2.5 mg/kg 2 52.05
## 21 2.5 mg/kg 2 23.810 2.5 mg/kg 2 53.63
## 22 2.5 mg/kg 2 16.605 2.5 mg/kg 2 65.27
## 23 2.5 mg/kg 3 374.046 2.5 mg/kg 3 195.20
## 24 2.5 mg/kg 3 60.573 2.5 mg/kg 3 48.37
## 25 2.5 mg/kg 3 21.222 2.5 mg/kg 3 26.53
## 26 2.5 mg/kg 3 95.238 2.5 mg/kg 3 72.56
## 27 2.5 mg/kg 3 103.321 2.5 mg/kg 3 76.35
## 28 5.0 mg/kg 1 312.844 5.0 mg/kg 1 246.68
## 29 5.0 mg/kg 1 574.419 5.0 mg/kg 1 176.06
## 30 5.0 mg/kg 1 144.467 5.0 mg/kg 1 94.32
## 31 5.0 mg/kg 1 215.328 5.0 mg/kg 1 165.45
## 32 5.0 mg/kg 2 179.817 5.0 mg/kg 2 132.69
## 33 5.0 mg/kg 2 469.767 5.0 mg/kg 2 145.97
## 34 5.0 mg/kg 2 119.877 5.0 mg/kg 2 81.25
## 35 5.0 mg/kg 2 215.328 5.0 mg/kg 2 120.00
## 36 5.0 mg/kg 3 321.101 5.0 mg/kg 3 254.52
## 37 5.0 mg/kg 3 444.186 5.0 mg/kg 3 126.12
## 38 5.0 mg/kg 3 119.877 5.0 mg/kg 3 97.73
## 39 5.0 mg/kg 3 379.562 5.0 mg/kg 3 250.91
## 40 10.0 mg/kg 1 217.277 10.0 mg/kg 1 234.38
## 41 10.0 mg/kg 1 403.955 10.0 mg/kg 1 250.73
## 42 10.0 mg/kg 1 307.047 10.0 mg/kg 1 211.61
## 43 10.0 mg/kg 1 15.432 10.0 mg/kg 1 16.42
## 44 10.0 mg/kg 2 340.314 10.0 mg/kg 2 262.50
## 45 10.0 mg/kg 2 220.339 10.0 mg/kg 2 151.60
## 46 10.0 mg/kg 2 384.228 10.0 mg/kg 2 180.40
## 47 10.0 mg/kg 2 30.864 10.0 mg/kg 2 73.13
## 48 10.0 mg/kg 3 86.387 10.0 mg/kg 3 92.19
## 49 10.0 mg/kg 3 316.384 10.0 mg/kg 3 183.67
## 50 10.0 mg/kg 3 803.691 10.0 mg/kg 3 324.59
## 51 10.0 mg/kg 3 21.605 10.0 mg/kg 3 55.22
## Group.2hFm Day.2hFm Data.2hFm
## 1 Saline 1 48.343
## 2 Saline 1 24.000
## 3 Saline 1 98.780
## 4 Saline 1 102.151
## 5 Saline 2 58.702
## 6 Saline 2 44.444
## 7 Saline 2 70.732
## 8 Saline 2 21.505
## 9 Saline 3 27.624
## 10 Saline 3 8.889
## 11 Saline 3 74.390
## 12 Saline 3 26.882
## 13 2.5 mg/kg 1 83.236
## 14 2.5 mg/kg 1 57.143
## 15 2.5 mg/kg 1 50.933
## 16 2.5 mg/kg 1 169.625
## 17 2.5 mg/kg 1 155.914
## 18 2.5 mg/kg 2 136.205
## 19 2.5 mg/kg 2 47.429
## 20 2.5 mg/kg 2 45.101
## 21 2.5 mg/kg 2 57.199
## 22 2.5 mg/kg 2 79.749
## 23 2.5 mg/kg 3 181.024
## 24 2.5 mg/kg 3 43.429
## 25 2.5 mg/kg 3 40.047
## 26 2.5 mg/kg 3 98.619
## 27 2.5 mg/kg 3 129.032
## 28 5.0 mg/kg 1 185.554
## 29 5.0 mg/kg 1 107.675
## 30 5.0 mg/kg 1 76.692
## 31 5.0 mg/kg 1 135.910
## 32 5.0 mg/kg 2 112.457
## 33 5.0 mg/kg 2 124.069
## 34 5.0 mg/kg 2 81.834
## 35 5.0 mg/kg 2 100.998
## 36 5.0 mg/kg 3 198.097
## 37 5.0 mg/kg 3 74.888
## 38 5.0 mg/kg 3 89.117
## 39 5.0 mg/kg 3 223.192
## 40 10.0 mg/kg 1 225.537
## 41 10.0 mg/kg 1 221.963
## 42 10.0 mg/kg 1 158.750
## 43 10.0 mg/kg 1 41.037
## 44 10.0 mg/kg 2 200.477
## 45 10.0 mg/kg 2 136.682
## 46 10.0 mg/kg 2 130.417
## 47 10.0 mg/kg 2 90.713
## 48 10.0 mg/kg 3 97.852
## 49 10.0 mg/kg 3 171.729
## 50 10.0 mg/kg 3 232.500
## 51 10.0 mg/kg 3 69.114
tapply(MPDhourly2FM$Data.2hFm, MPDhourly2FM$Group.2hFm, sd)
## 10.0 mg/kg 2.5 mg/kg 5.0 mg/kg Saline
## 64.63 49.99 50.26 30.81
tapply(MPDhourly2FM$Data.2hFm, MPDhourly2FM$Day.2hFm, sd)
## 1 2 3
## 63.37 45.50 71.81
MPDhourly2FM$Group.2hFm <- as.factor(MPDhourly2FM$Group.2hFm)
MPDhourly2FM$Day.2hFm <- as.factor(MPDhourly2FM$Day.2hFm)
MPDhourly2FM$Data.2hFm <- as.numeric(MPDhourly2FM$Data.2hFm)
aov.MPDhourly2FM = aov(Data.2hFm ~ Group.1hFm * Day.2hFm, data = MPDhourly2FM)
summary(aov.MPDhourly2FM)
## Df Sum Sq Mean Sq F value Pr(>F)
## Group.1hFm 3 65606 21869 7.78 0.00034 ***
## Day.2hFm 2 4894 2447 0.87 0.42672
## Group.1hFm:Day.2hFm 6 4648 775 0.28 0.94505
## Residuals 39 109632 2811
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
comparison <- aov(Data.2hFm ~ Group.2hFm * Day.2hFm, data = MPDhourly2FM)
TukeyHSD(comparison, "Group.2hFm")
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = Data.2hFm ~ Group.2hFm * Day.2hFm, data = MPDhourly2FM)
##
## $Group.2hFm
## diff lwr upr p adj
## 2.5 mg/kg-10.0 mg/kg -56.42 -111.52 -1.317 0.0430
## 5.0 mg/kg-10.0 mg/kg -22.19 -80.27 35.891 0.7358
## Saline-10.0 mg/kg -97.53 -155.61 -39.445 0.0003
## 5.0 mg/kg-2.5 mg/kg 34.23 -20.87 89.329 0.3546
## Saline-2.5 mg/kg -41.11 -96.21 13.993 0.2048
## Saline-5.0 mg/kg -75.34 -133.42 -17.255 0.0066
TukeyHSD(comparison, "Day.2hFm")
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = Data.2hFm ~ Group.2hFm * Day.2hFm, data = MPDhourly2FM)
##
## $Day.2hFm
## diff lwr upr p adj
## 2-1 -23.796 -68.10 20.51 0.3991
## 3-1 -9.224 -53.53 35.08 0.8683
## 3-2 14.571 -29.73 58.88 0.7045
print(model.tables(aov.MPDhourly2FM, "means"), digits = 9)
## Tables of means
## Grand mean
##
## 103.3
##
## Group.1hFm
## 10.0 mg/kg 2.5 mg/kg 5.0 mg/kg Saline
## 148.064231 91.6457104 125.873537 50.5368263
## rep 12.000000 15.0000000 12.000000 12.0000000
##
## Day.2hFm
## 1 2 3
## 114.308398 90.512447 105.083907
## rep 17.000000 17.000000 17.000000
##
## Group.1hFm:Day.2hFm
## Day.2hFm
## Group.1hFm 1 2 3
## 10.0 mg/kg 161.82 139.57 142.80
## rep 4.00 4.00 4.00
## 2.5 mg/kg 103.37 73.14 98.43
## rep 5.00 5.00 5.00
## 5.0 mg/kg 126.46 104.84 146.32
## rep 4.00 4.00 4.00
## Saline 68.32 48.85 34.45
## rep 4.00 4.00 4.00