Sal vs. IL falafel bout analysis
exp1FB <- read.csv("~/Desktop/experiment 1 bout analysis F.csv")
exp1FB
## Group.Fb Days.Fb Data.Fb Group.Fbl Days.Fbl Data.Fbl
## 1 Saline 1 86.78 Saline 1 59.68
## 2 Saline 1 99.01 Saline 1 48.03
## 3 Saline 1 69.70 Saline 1 60.06
## 4 Saline 1 91.30 Saline 1 72.18
## 5 Saline 1 76.53 Saline 1 71.87
## 6 Saline 1 73.86 Saline 1 73.74
## 7 Saline 1 104.76 Saline 1 95.59
## 8 Saline 1 102.56 Saline 1 45.29
## 9 Saline 1 82.09 Saline 1 84.02
## 10 Saline 2 95.04 Saline 2 73.58
## 11 Saline 2 74.26 Saline 2 79.86
## 12 Saline 2 75.76 Saline 2 85.09
## 13 Saline 2 82.61 Saline 2 78.95
## 14 Saline 2 96.94 Saline 2 97.09
## 15 Saline 2 73.86 Saline 2 59.97
## 16 Saline 2 123.81 Saline 2 93.14
## 17 Saline 2 111.11 Saline 2 56.92
## 18 Saline 2 115.67 Saline 2 91.88
## 19 Saline 3 86.78 Saline 3 60.25
## 20 Saline 3 74.26 Saline 3 76.97
## 21 Saline 3 103.03 Saline 3 96.10
## 22 Saline 3 73.91 Saline 3 53.38
## 23 Saline 3 91.84 Saline 3 97.86
## 24 Saline 3 85.23 Saline 3 75.36
## 25 Saline 3 90.48 Saline 3 73.53
## 26 Saline 3 102.56 Saline 3 66.71
## 27 Saline 3 123.13 Saline 3 121.77
## 28 IL-1B 1 76.92 IL-1B 1 62.38
## 29 IL-1B 1 48.91 IL-1B 1 14.63
## 30 IL-1B 1 37.50 IL-1B 1 28.65
## 31 IL-1B 1 33.02 IL-1B 1 17.33
## 32 IL-1B 1 18.02 IL-1B 1 16.94
## 33 IL-1B 1 47.06 IL-1B 1 23.48
## 34 IL-1B 1 47.62 IL-1B 1 23.27
## 35 IL-1B 1 51.20 IL-1B 1 18.37
## 36 IL-1B 1 79.17 IL-1B 1 29.82
## 37 IL-1B 2 91.35 IL-1B 2 60.40
## 38 IL-1B 2 48.91 IL-1B 2 20.33
## 39 IL-1B 2 46.88 IL-1B 2 33.85
## 40 IL-1B 2 47.17 IL-1B 2 43.73
## 41 IL-1B 2 58.56 IL-1B 2 51.55
## 42 IL-1B 2 64.71 IL-1B 2 49.90
## 43 IL-1B 2 71.43 IL-1B 2 42.75
## 44 IL-1B 2 54.22 IL-1B 2 52.42
## 45 IL-1B 2 83.33 IL-1B 2 52.28
## 46 IL-1B 3 134.62 IL-1B 3 121.78
## 47 IL-1B 3 114.13 IL-1B 3 116.26
## 48 IL-1B 3 115.62 IL-1B 3 118.23
## 49 IL-1B 3 127.36 IL-1B 3 115.51
## 50 IL-1B 3 99.10 IL-1B 3 95.73
## 51 IL-1B 3 100.00 IL-1B 3 105.68
## 52 IL-1B 3 104.76 IL-1B 3 162.88
## 53 IL-1B 3 78.31 IL-1B 3 61.38
## 54 IL-1B 3 112.50 IL-1B 3 122.11
tapply(exp1FB$Data.Fb, exp1FB$Group.Fb, sd)
## IL-1B Saline
## 31.36 15.68
tapply(exp1FB$Data.Fb, exp1FB$Days.Fb, sd)
## 1 2 3
## 25.52 23.58 17.93
exp1FB$Group.Fb <- as.factor(exp1FB$Group.Fb)
exp1FB$Days.Fb <- as.factor(exp1FB$Days.Fb)
exp1FB$Data.Fb <- as.numeric(exp1FB$Data.Fb)
aov.exp1FB = aov(Data.Fb ~ Group.Fb * Days.Fb, data = exp1FB)
summary(aov.exp1FB)
## Df Sum Sq Mean Sq F value Pr(>F)
## Group.Fb 1 4169 4169 14.8 0.00035 ***
## Days.Fb 2 10140 5070 18.0 1.5e-06 ***
## Group.Fb:Days.Fb 2 8299 4149 14.7 1.0e-05 ***
## Residuals 48 13517 282
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
comparison <- aov(Data.Fb ~ Group.Fb * Days.Fb, data = exp1FB)
TukeyHSD(comparison, "Group.Fb")
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = Data.Fb ~ Group.Fb * Days.Fb, data = exp1FB)
##
## $Group.Fb
## diff lwr upr p adj
## Saline-IL-1B 17.57 8.391 26.76 4e-04
TukeyHSD(comparison, "Days.Fb")
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = Data.Fb ~ Group.Fb * Days.Fb, data = exp1FB)
##
## $Days.Fb
## diff lwr upr p adj
## 2-1 10.53 -2.995 24.06 0.1546
## 3-1 32.87 19.339 46.39 0.0000
## 3-2 22.33 8.806 35.86 0.0006
print(model.tables(aov.exp1FB, "means"), digits = 17)
## Tables of means
## Grand mean
##
## 82.58
##
## Group.Fb
## Group.Fb
## IL-1B Saline
## 73.79 91.37
##
## Days.Fb
## Days.Fb
## 1 2 3
## 68.11 78.64 100.98
##
## Group.Fb:Days.Fb
## Days.Fb
## Group.Fb 1 2 3
## IL-1B 48.82 62.95 109.60
## Saline 87.40 94.34 92.36
exp1FBl <- read.csv("~/Desktop/experiment 1 bout analysis F.csv")
exp1FBl
## Group.Fb Days.Fb Data.Fb Group.Fbl Days.Fbl Data.Fbl
## 1 Saline 1 86.78 Saline 1 59.68
## 2 Saline 1 99.01 Saline 1 48.03
## 3 Saline 1 69.70 Saline 1 60.06
## 4 Saline 1 91.30 Saline 1 72.18
## 5 Saline 1 76.53 Saline 1 71.87
## 6 Saline 1 73.86 Saline 1 73.74
## 7 Saline 1 104.76 Saline 1 95.59
## 8 Saline 1 102.56 Saline 1 45.29
## 9 Saline 1 82.09 Saline 1 84.02
## 10 Saline 2 95.04 Saline 2 73.58
## 11 Saline 2 74.26 Saline 2 79.86
## 12 Saline 2 75.76 Saline 2 85.09
## 13 Saline 2 82.61 Saline 2 78.95
## 14 Saline 2 96.94 Saline 2 97.09
## 15 Saline 2 73.86 Saline 2 59.97
## 16 Saline 2 123.81 Saline 2 93.14
## 17 Saline 2 111.11 Saline 2 56.92
## 18 Saline 2 115.67 Saline 2 91.88
## 19 Saline 3 86.78 Saline 3 60.25
## 20 Saline 3 74.26 Saline 3 76.97
## 21 Saline 3 103.03 Saline 3 96.10
## 22 Saline 3 73.91 Saline 3 53.38
## 23 Saline 3 91.84 Saline 3 97.86
## 24 Saline 3 85.23 Saline 3 75.36
## 25 Saline 3 90.48 Saline 3 73.53
## 26 Saline 3 102.56 Saline 3 66.71
## 27 Saline 3 123.13 Saline 3 121.77
## 28 IL-1B 1 76.92 IL-1B 1 62.38
## 29 IL-1B 1 48.91 IL-1B 1 14.63
## 30 IL-1B 1 37.50 IL-1B 1 28.65
## 31 IL-1B 1 33.02 IL-1B 1 17.33
## 32 IL-1B 1 18.02 IL-1B 1 16.94
## 33 IL-1B 1 47.06 IL-1B 1 23.48
## 34 IL-1B 1 47.62 IL-1B 1 23.27
## 35 IL-1B 1 51.20 IL-1B 1 18.37
## 36 IL-1B 1 79.17 IL-1B 1 29.82
## 37 IL-1B 2 91.35 IL-1B 2 60.40
## 38 IL-1B 2 48.91 IL-1B 2 20.33
## 39 IL-1B 2 46.88 IL-1B 2 33.85
## 40 IL-1B 2 47.17 IL-1B 2 43.73
## 41 IL-1B 2 58.56 IL-1B 2 51.55
## 42 IL-1B 2 64.71 IL-1B 2 49.90
## 43 IL-1B 2 71.43 IL-1B 2 42.75
## 44 IL-1B 2 54.22 IL-1B 2 52.42
## 45 IL-1B 2 83.33 IL-1B 2 52.28
## 46 IL-1B 3 134.62 IL-1B 3 121.78
## 47 IL-1B 3 114.13 IL-1B 3 116.26
## 48 IL-1B 3 115.62 IL-1B 3 118.23
## 49 IL-1B 3 127.36 IL-1B 3 115.51
## 50 IL-1B 3 99.10 IL-1B 3 95.73
## 51 IL-1B 3 100.00 IL-1B 3 105.68
## 52 IL-1B 3 104.76 IL-1B 3 162.88
## 53 IL-1B 3 78.31 IL-1B 3 61.38
## 54 IL-1B 3 112.50 IL-1B 3 122.11
tapply(exp1FBl$Data.Fbl, exp1FBl$Group.Fbl, sd)
## IL-1B Saline
## 42.22 17.85
tapply(exp1FBl$Data.Fbl, exp1FBl$Days.Fbl, sd)
## 1 2 3
## 26.20 21.79 29.00
exp1FBl$Group.Fbl <- as.factor(exp1FBl$Group.Fbl)
exp1FBl$Days.Fbl <- as.factor(exp1FBl$Days.Fbl)
exp1FBl$Data.Fbl <- as.numeric(exp1FBl$Data.Fbl)
aov.exp1FBl = aov(Data.Fbl ~ Group.Fbl * Days.Fbl, data = exp1FBl)
summary(aov.exp1FBl)
## Df Sum Sq Mean Sq F value Pr(>F)
## Group.Fbl 1 2777 2777 8.35 0.0058 **
## Days.Fbl 2 23376 11688 35.13 4e-10 ***
## Group.Fbl:Days.Fbl 2 15294 7647 22.99 1e-07 ***
## Residuals 48 15970 333
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
comparison <- aov(Data.Fbl ~ Group.Fbl * Days.Fbl, data = exp1FBl)
TukeyHSD(comparison, "Group.Fbl")
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = Data.Fbl ~ Group.Fbl * Days.Fbl, data = exp1FBl)
##
## $Group.Fbl
## diff lwr upr p adj
## Saline-IL-1B 14.34 4.361 24.32 0.0058
TukeyHSD(comparison, "Days.Fbl")
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = Data.Fbl ~ Group.Fbl * Days.Fbl, data = exp1FBl)
##
## $Days.Fbl
## diff lwr upr p adj
## 2-1 15.46 0.76 30.17 0.0373
## 3-1 49.79 35.08 64.49 0.0000
## 3-2 34.32 19.62 49.03 0.0000
print(model.tables(aov.exp1FBl, "means"), digits = 17)
## Tables of means
## Grand mean
##
## 68.71
##
## Group.Fbl
## Group.Fbl
## IL-1B Saline
## 61.54 75.88
##
## Days.Fbl
## Days.Fbl
## 1 2 3
## 46.96 62.43 96.75
##
## Group.Fbl:Days.Fbl
## Days.Fbl
## Group.Fbl 1 2 3
## IL-1B 26.10 45.24 113.28
## Saline 67.83 79.61 80.22