Acidogenic Boluses and Welfare
Cainan Florentino, DVM
2026-05-26
PhD Dissertation - Chapter 5
1 - Descriptive Results
1.1 - Testing differences between groups - Baseline
| Variable | No-Bolus | Two-Bolus | P-value |
|---|---|---|---|
| Animals | |||
| Total Animals | 15 | 12 | |
| Udder Engorgement | |||
| Baseline Teat Distance (cm) | 6.64 (2.84) | 6.65 (3.09) | 0.995 |
| Baseline udder pressure at 20 cm (left) | 0.44 (0.26) | 0.48 (0.28) | 0.711 |
| Baseline udder pressure at 5 cm (left) | 0.28 (0.19) | 0.30 (0.26) | 0.792 |
| Baseline udder pressure at 20 cm (right) | 0.53 (0.61) | 0.46 (0.33) | 0.722 |
| Baseline udder pressure at 5 cm (right) | 0.32 (0.30) | 0.29 (0.20) | 0.719 |
| Acid-Base Balance | |||
| Baseline Urine pH | 8.23 (0.41) | 8.33 (0.27) | 0.490 |
| Baseline Blood pH | 7.48 (0.05) | 7.49 (0.06) | 0.707 |
| Production | |||
| Days in Milk at Dry-off | 305.60 (45.02) | 311.50 (25.53) | 0.690 |
| Milk Production Before Randomization (24h, kg) | 23.38 (8.74) | 25.42 (7.12) | 0.519 |
| Average Milk Production (week before Dry-off) | 26.27 (7.98) | 26.18 (5.85) | 0.976 |
| Average feed intake (week before Dry-off) | 39.48 (6.39) | 38.52 (5.39) | 0.712 |
| Behavior | |||
| Total Standing Time (day prior to Randomization) | 678.29 (137.29) | 698.15 (74.60) | 0.661 |
| Total Lying Time (day prior to Randomization) | 754.00 (137.56) | 729.33 (71.29) | 0.584 |
| Number of Bouts (day prior to Randomization) | 15.38 (5.94) | 12.00 (1.41) | 0.067 |
| Bout Duration (day prior to Randomization) | 54.29 (17.97) | 61.77 (10.74) | 0.224 |
| Note: Continuous variables are presented as mean (SD). P-values calculated using two-sample t-tests. | |||
2 - Teat Distance
2.1 - Teat Distance overtime
## $emmeans
## hours = 0:
## Treatment emmean SE df lower.CL upper.CL
## Control 6.74 0.704 20 5.27 8.21
## Treatment 5.83 0.867 20 4.02 7.64
##
## hours = 36:
## Treatment emmean SE df lower.CL upper.CL
## Control 7.81 0.691 20 6.37 9.25
## Treatment 6.20 0.856 20 4.41 7.98
##
## hours = 48:
## Treatment emmean SE df lower.CL upper.CL
## Control 9.96 0.686 20 8.53 11.39
## Treatment 7.81 0.856 20 6.03 9.60
##
## hours = 60:
## Treatment emmean SE df lower.CL upper.CL
## Control 9.16 0.686 20 7.73 10.59
## Treatment 7.65 0.856 20 5.86 9.43
##
## hours = 72:
## Treatment emmean SE df lower.CL upper.CL
## Control 9.65 0.686 20 8.22 11.08
## Treatment 8.06 0.856 20 6.27 9.85
##
## hours = 84:
## Treatment emmean SE df lower.CL upper.CL
## Control 9.07 0.687 20 7.63 10.50
## Treatment 7.81 0.857 20 6.03 9.60
##
## hours = 96:
## Treatment emmean SE df lower.CL upper.CL
## Control 8.21 0.687 20 6.77 9.64
## Treatment 7.60 0.857 20 5.81 9.39
##
## hours = 120:
## Treatment emmean SE df lower.CL upper.CL
## Control 8.20 0.687 20 6.77 9.64
## Treatment 7.10 0.857 20 5.31 8.89
##
## hours = 144:
## Treatment emmean SE df lower.CL upper.CL
## Control 7.63 0.687 20 6.20 9.07
## Treatment 7.26 0.857 20 5.47 9.04
##
## hours = 168:
## Treatment emmean SE df lower.CL upper.CL
## Control 7.70 0.687 20 6.27 9.14
## Treatment 6.42 0.857 20 4.64 8.21
##
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## hours = 0:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.906 1.12 20 0.811 0.4268
##
## hours = 36:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 1.613 1.10 20 1.467 0.1580
##
## hours = 48:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 2.152 1.10 20 1.962 0.0639
##
## hours = 60:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 1.513 1.10 20 1.379 0.1830
##
## hours = 72:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 1.589 1.10 20 1.449 0.1630
##
## hours = 84:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 1.254 1.10 20 1.142 0.2669
##
## hours = 96:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.604 1.10 20 0.550 0.5883
##
## hours = 120:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 1.103 1.10 20 1.005 0.3270
##
## hours = 144:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.376 1.10 20 0.343 0.7355
##
## hours = 168:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 1.281 1.10 20 1.167 0.2571
##
## Degrees-of-freedom method: containment## hours = 0:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.906 1.12 20 -1.424 3.24
##
## hours = 36:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 1.613 1.10 20 -0.681 3.91
##
## hours = 48:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 2.152 1.10 20 -0.136 4.44
##
## hours = 60:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 1.513 1.10 20 -0.775 3.80
##
## hours = 72:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 1.589 1.10 20 -0.699 3.88
##
## hours = 84:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 1.254 1.10 20 -1.036 3.54
##
## hours = 96:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.604 1.10 20 -1.686 2.89
##
## hours = 120:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 1.103 1.10 20 -1.187 3.39
##
## hours = 144:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.376 1.10 20 -1.914 2.67
##
## hours = 168:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 1.281 1.10 20 -1.010 3.57
##
## Degrees-of-freedom method: containment
## Confidence level used: 0.952.2 - Graph overtime
## Warning: `aes_()` was deprecated in ggplot2 3.0.0.
## ℹ Please use tidy evaluation idioms with `aes()`
## ℹ The deprecated feature was likely used in the emmeans package.
## Please report the issue at <https://github.com/rvlenth/emmeans/issues>.
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.
## Warning: Using `size` aesthetic for lines was deprecated in ggplot2 3.4.0.
## ℹ Please use `linewidth` instead.
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.
2.3 - Type III p-value
## Analysis of Deviance Table (Type III tests)
##
## Response: teatdistance
## Chisq Df Pr(>Chisq)
## (Intercept) 22.7719 1 1.824e-06 ***
## as.factor(Treatment) 0.6579 1 0.4173
## hours 89.1577 9 2.402e-15 ***
## weekbeforetreatmilk 1.6444 1 0.1997
## as.factor(Treatment):hours 10.7281 9 0.2948
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 12.4 - Testing normality of the model
2.5 - Average teat distance
## $emmeans
## Treatment emmean SE df lower.CL upper.CL
## Control 8.41 0.631 20 7.10 9.73
## Treatment 7.17 0.786 20 5.53 8.81
##
## Results are averaged over the levels of: hours
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## Control - Treatment 1.24 1.01 20 1.231 0.2327
##
## Results are averaged over the levels of: hours
## Degrees-of-freedom method: containment## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 1.24 1.01 20 -0.862 3.34
##
## Results are averaged over the levels of: hours
## Degrees-of-freedom method: containment
## Confidence level used: 0.953 - Pictures Model
3.1 - Udder Temperature overtime
## $emmeans
## hours = 36:
## Treatment emmean SE df lower.CL upper.CL
## Control 96.1 1.34 18 93.2 98.9
## Treatment 96.0 1.73 18 92.3 99.6
##
## hours = 48:
## Treatment emmean SE df lower.CL upper.CL
## Control 91.8 1.32 18 89.0 94.5
## Treatment 90.9 1.68 18 87.3 94.4
##
## hours = 60:
## Treatment emmean SE df lower.CL upper.CL
## Control 95.2 1.35 18 92.4 98.1
## Treatment 93.7 1.73 18 90.0 97.3
##
## hours = 72:
## Treatment emmean SE df lower.CL upper.CL
## Control 92.2 1.35 18 89.4 95.0
## Treatment 90.2 1.66 18 86.7 93.7
##
## hours = 84:
## Treatment emmean SE df lower.CL upper.CL
## Control 92.7 1.37 18 89.9 95.6
## Treatment 92.9 1.68 18 89.3 96.4
##
## hours = 96:
## Treatment emmean SE df lower.CL upper.CL
## Control 91.0 1.35 18 88.2 93.9
## Treatment 93.0 1.68 18 89.5 96.6
##
## hours = 120:
## Treatment emmean SE df lower.CL upper.CL
## Control 91.4 1.35 18 88.5 94.2
## Treatment 89.9 1.66 18 86.4 93.4
##
## hours = 144:
## Treatment emmean SE df lower.CL upper.CL
## Control 92.1 1.32 18 89.3 94.8
## Treatment 90.7 1.76 18 87.0 94.4
##
## hours = 168:
## Treatment emmean SE df lower.CL upper.CL
## Control 91.2 1.32 18 88.4 93.9
## Treatment 91.8 1.68 18 88.3 95.4
##
## Results are averaged over the levels of: Season
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## hours = 36:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.0805 1.68 18 0.048 0.9623
##
## hours = 48:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.9020 1.61 18 0.559 0.5831
##
## hours = 60:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 1.5665 1.70 18 0.922 0.3687
##
## hours = 72:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 1.9673 1.62 18 1.217 0.2394
##
## hours = 84:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -0.1207 1.67 18 -0.072 0.9433
##
## hours = 96:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -2.0113 1.65 18 -1.217 0.2394
##
## hours = 120:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 1.4475 1.61 18 0.896 0.3818
##
## hours = 144:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 1.3289 1.66 18 0.798 0.4351
##
## hours = 168:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -0.6772 1.61 18 -0.420 0.6794
##
## Results are averaged over the levels of: Season
## Degrees-of-freedom method: containment## hours = 36:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.0805 1.68 18 -3.45 3.61
##
## hours = 48:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.9020 1.61 18 -2.49 4.29
##
## hours = 60:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 1.5665 1.70 18 -2.00 5.14
##
## hours = 72:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 1.9673 1.62 18 -1.43 5.36
##
## hours = 84:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -0.1207 1.67 18 -3.64 3.40
##
## hours = 96:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -2.0113 1.65 18 -5.48 1.46
##
## hours = 120:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 1.4475 1.61 18 -1.94 4.84
##
## hours = 144:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 1.3289 1.66 18 -2.17 4.83
##
## hours = 168:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -0.6772 1.61 18 -4.06 2.71
##
## Results are averaged over the levels of: Season
## Degrees-of-freedom method: containment
## Confidence level used: 0.953.2 - Graph overtime
3.3 - Type III p-value
## Analysis of Deviance Table (Type III tests)
##
## Response: average
## Chisq Df Pr(>Chisq)
## (Intercept) 1093.9928 1 < 2.2e-16 ***
## as.factor(Treatment) 0.0023 1 0.9618
## hours 40.3938 8 2.705e-06 ***
## weekbeforetreatmilk 1.8371 1 0.1753
## Season 27.7507 2 9.419e-07 ***
## as.factor(Treatment):hours 7.6643 8 0.4669
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 13.4 - Testing normality of the model
3.5 - Average temperature
## $emmeans
## Treatment emmean SE df lower.CL upper.CL
## Control 92.7 1.09 18 90.4 95.0
## Treatment 92.2 1.37 18 89.3 95.1
##
## Results are averaged over the levels of: hours, Season
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.495 1.07 18 0.462 0.6493
##
## Results are averaged over the levels of: hours, Season
## Degrees-of-freedom method: containment4 - Pictures Model (C)
4.1 - Udder Temperature overtime
## $emmeans
## ampmh = 1.5 PM:
## Treatment emmean SE df lower.CL upper.CL
## Control 35.6 0.745 18 34.0 37.2
## Treatment 35.5 0.961 18 33.5 37.6
##
## ampmh = 2 AM:
## Treatment emmean SE df lower.CL upper.CL
## Control 33.2 0.734 18 31.7 34.7
## Treatment 32.7 0.934 18 30.7 34.7
##
## ampmh = 2.5 PM:
## Treatment emmean SE df lower.CL upper.CL
## Control 35.1 0.749 18 33.6 36.7
## Treatment 34.3 0.961 18 32.2 36.3
##
## ampmh = 3 AM:
## Treatment emmean SE df lower.CL upper.CL
## Control 33.4 0.749 18 31.9 35.0
## Treatment 32.4 0.925 18 30.4 34.3
##
## ampmh = 3.5 PM:
## Treatment emmean SE df lower.CL upper.CL
## Control 33.7 0.762 18 32.1 35.3
## Treatment 33.8 0.935 18 31.8 35.8
##
## ampmh = 4 AM:
## Treatment emmean SE df lower.CL upper.CL
## Control 32.8 0.749 18 31.2 34.4
## Treatment 33.9 0.935 18 31.9 35.9
##
## ampmh = 5 AM:
## Treatment emmean SE df lower.CL upper.CL
## Control 33.0 0.749 18 31.4 34.5
## Treatment 32.2 0.925 18 30.2 34.1
##
## ampmh = 6 AM:
## Treatment emmean SE df lower.CL upper.CL
## Control 33.4 0.734 18 31.8 34.9
## Treatment 32.6 0.978 18 30.6 34.7
##
## ampmh = 7 AM:
## Treatment emmean SE df lower.CL upper.CL
## Control 32.9 0.734 18 31.3 34.4
## Treatment 33.2 0.935 18 31.3 35.2
##
## Results are averaged over the levels of: Season
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## ampmh = 1.5 PM:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.0447 0.933 18 0.048 0.9623
##
## ampmh = 2 AM:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.5011 0.897 18 0.559 0.5831
##
## ampmh = 2.5 PM:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.8703 0.944 18 0.922 0.3687
##
## ampmh = 3 AM:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 1.0929 0.898 18 1.217 0.2394
##
## ampmh = 3.5 PM:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -0.0670 0.930 18 -0.072 0.9433
##
## ampmh = 4 AM:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -1.1174 0.918 18 -1.217 0.2394
##
## ampmh = 5 AM:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.8042 0.897 18 0.896 0.3818
##
## ampmh = 6 AM:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.7383 0.925 18 0.798 0.4351
##
## ampmh = 7 AM:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -0.3762 0.896 18 -0.420 0.6794
##
## Results are averaged over the levels of: Season
## Degrees-of-freedom method: containment## ampmh = 1.5 PM:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.0447 0.933 18 -1.916 2.006
##
## ampmh = 2 AM:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.5011 0.897 18 -1.383 2.385
##
## ampmh = 2.5 PM:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.8703 0.944 18 -1.112 2.853
##
## ampmh = 3 AM:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 1.0929 0.898 18 -0.794 2.980
##
## ampmh = 3.5 PM:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -0.0670 0.930 18 -2.020 1.886
##
## ampmh = 4 AM:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -1.1174 0.918 18 -3.047 0.812
##
## ampmh = 5 AM:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.8042 0.897 18 -1.080 2.689
##
## ampmh = 6 AM:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.7383 0.925 18 -1.205 2.681
##
## ampmh = 7 AM:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -0.3762 0.896 18 -2.258 1.506
##
## Results are averaged over the levels of: Season
## Degrees-of-freedom method: containment
## Confidence level used: 0.954.2 - Graph overtime
4.3 - Type III p-value
## Analysis of Deviance Table (Type III tests)
##
## Response: averagec
## Chisq Df Pr(>Chisq)
## (Intercept) 457.9192 1 < 2.2e-16 ***
## as.factor(Treatment) 0.0023 1 0.9618
## hours 40.3938 8 2.705e-06 ***
## weekbeforetreatmilk 1.8371 1 0.1753
## Season 27.7507 2 9.419e-07 ***
## as.factor(Treatment):hours 7.6643 8 0.4669
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 14.4 - Testing normality of the model
4.5 - Average temperature
## $emmeans
## Treatment emmean SE df lower.CL upper.CL
## Control 33.7 0.604 18 32.4 35
## Treatment 33.4 0.763 18 31.8 35
##
## Results are averaged over the levels of: hours, Season
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.275 0.595 18 0.462 0.6493
##
## Results are averaged over the levels of: hours, Season
## Degrees-of-freedom method: containment5 - Udder Pressure Higher in the Udder
5.1 - Pressure overtime
## $emmeans
## hours = 0:
## Treatment emmean SE df lower.CL upper.CL
## Control 0.621 0.169 20 0.2695 0.973
## Treatment 0.554 0.203 20 0.1300 0.978
##
## hours = 36:
## Treatment emmean SE df lower.CL upper.CL
## Control 0.931 0.165 20 0.5870 1.276
## Treatment 0.819 0.203 20 0.3959 1.242
##
## hours = 48:
## Treatment emmean SE df lower.CL upper.CL
## Control 1.332 0.163 20 0.9928 1.671
## Treatment 0.996 0.203 20 0.5730 1.419
##
## hours = 60:
## Treatment emmean SE df lower.CL upper.CL
## Control 1.739 0.163 20 1.3999 2.079
## Treatment 0.880 0.203 20 0.4573 1.304
##
## hours = 72:
## Treatment emmean SE df lower.CL upper.CL
## Control 1.293 0.163 20 0.9539 1.632
## Treatment 1.077 0.203 20 0.6537 1.500
##
## hours = 84:
## Treatment emmean SE df lower.CL upper.CL
## Control 1.551 0.163 20 1.2112 1.891
## Treatment 0.773 0.203 20 0.3496 1.197
##
## hours = 96:
## Treatment emmean SE df lower.CL upper.CL
## Control 0.958 0.163 20 0.6180 1.298
## Treatment 0.564 0.203 20 0.1404 0.988
##
## hours = 120:
## Treatment emmean SE df lower.CL upper.CL
## Control 0.827 0.163 20 0.4871 1.167
## Treatment 0.568 0.203 20 0.1443 0.992
##
## hours = 144:
## Treatment emmean SE df lower.CL upper.CL
## Control 0.551 0.163 20 0.2117 0.891
## Treatment 0.496 0.203 20 0.0726 0.920
##
## hours = 168:
## Treatment emmean SE df lower.CL upper.CL
## Control 0.582 0.163 20 0.2426 0.922
## Treatment 0.552 0.203 20 0.1283 0.976
##
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## hours = 0:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.0673 0.264 20 0.255 0.8014
##
## hours = 36:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.1124 0.262 20 0.430 0.6719
##
## hours = 48:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.3360 0.260 20 1.292 0.2111
##
## hours = 60:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.8587 0.260 20 3.303 0.0036
##
## hours = 72:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.2164 0.260 20 0.832 0.4151
##
## hours = 84:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.7776 0.260 20 2.987 0.0073
##
## hours = 96:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.3935 0.260 20 1.511 0.1463
##
## hours = 120:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.2588 0.260 20 0.994 0.3322
##
## hours = 144:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.0551 0.260 20 0.212 0.8345
##
## hours = 168:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.0303 0.260 20 0.116 0.9086
##
## Degrees-of-freedom method: containment## hours = 0:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.0673 0.264 20 -0.483 0.618
##
## hours = 36:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.1124 0.262 20 -0.433 0.658
##
## hours = 48:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.3360 0.260 20 -0.206 0.878
##
## hours = 60:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.8587 0.260 20 0.316 1.401
##
## hours = 72:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.2164 0.260 20 -0.326 0.759
##
## hours = 84:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.7776 0.260 20 0.234 1.321
##
## hours = 96:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.3935 0.260 20 -0.150 0.937
##
## hours = 120:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.2588 0.260 20 -0.284 0.802
##
## hours = 144:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.0551 0.260 20 -0.488 0.598
##
## hours = 168:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.0303 0.260 20 -0.513 0.573
##
## Degrees-of-freedom method: containment
## Confidence level used: 0.955.2 - Graph overtime
5.3 - Type III p-value
## Analysis of Deviance Table (Type III tests)
##
## Response: meantop
## Chisq Df Pr(>Chisq)
## (Intercept) 0.5982 1 0.4392554
## as.factor(Treatment) 0.0649 1 0.7988452
## hours 134.4199 9 < 2.2e-16 ***
## weekbeforetreatmilk 0.7174 1 0.3969936
## as.factor(Treatment):hours 29.6065 9 0.0005117 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 15.4 - Testing normality of the model
5.5 - Average pressure
## $emmeans
## Treatment emmean SE df lower.CL upper.CL
## Control 1.037 0.132 20 0.761 1.31
## Treatment 0.733 0.165 20 0.389 1.08
##
## Results are averaged over the levels of: hours
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.304 0.211 20 1.437 0.1662
##
## Results are averaged over the levels of: hours
## Degrees-of-freedom method: containment## $emmeans
## Treatment emmean SE df lower.CL upper.CL
## Control 1.037 0.132 20 0.761 1.31
## Treatment 0.733 0.165 20 0.389 1.08
##
## Results are averaged over the levels of: hours
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.304 0.211 20 -0.137 0.745
##
## Results are averaged over the levels of: hours
## Degrees-of-freedom method: containment
## Confidence level used: 0.956 - Udder Pressure Teat Cistern
6.1 - Pressure overtime
## $emmeans
## hours = 0:
## Treatment emmean SE df lower.CL upper.CL
## Control 0.377 0.144 20 0.0765 0.678
## Treatment 0.336 0.175 20 -0.0294 0.702
##
## hours = 36:
## Treatment emmean SE df lower.CL upper.CL
## Control 0.643 0.141 20 0.3484 0.938
## Treatment 0.468 0.174 20 0.1048 0.832
##
## hours = 48:
## Treatment emmean SE df lower.CL upper.CL
## Control 1.299 0.140 20 1.0066 1.591
## Treatment 0.882 0.175 20 0.5179 1.246
##
## hours = 60:
## Treatment emmean SE df lower.CL upper.CL
## Control 1.403 0.140 20 1.1113 1.695
## Treatment 0.912 0.175 20 0.5480 1.276
##
## hours = 72:
## Treatment emmean SE df lower.CL upper.CL
## Control 1.161 0.140 20 0.8700 1.453
## Treatment 0.918 0.174 20 0.5548 1.282
##
## hours = 84:
## Treatment emmean SE df lower.CL upper.CL
## Control 1.249 0.140 20 0.9563 1.542
## Treatment 0.768 0.175 20 0.4025 1.133
##
## hours = 96:
## Treatment emmean SE df lower.CL upper.CL
## Control 0.934 0.141 20 0.6412 1.227
## Treatment 0.577 0.175 20 0.2113 0.943
##
## hours = 120:
## Treatment emmean SE df lower.CL upper.CL
## Control 0.758 0.141 20 0.4653 1.052
## Treatment 0.570 0.175 20 0.2046 0.936
##
## hours = 144:
## Treatment emmean SE df lower.CL upper.CL
## Control 0.543 0.141 20 0.2497 0.836
## Treatment 0.497 0.175 20 0.1311 0.862
##
## hours = 168:
## Treatment emmean SE df lower.CL upper.CL
## Control 0.482 0.141 20 0.1889 0.775
## Treatment 0.525 0.175 20 0.1592 0.890
##
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## hours = 0:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.0413 0.227 20 0.182 0.8574
##
## hours = 36:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.1748 0.224 20 0.779 0.4451
##
## hours = 48:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.4165 0.224 20 1.861 0.0774
##
## hours = 60:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.4912 0.224 20 2.195 0.0401
##
## hours = 72:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.2431 0.223 20 1.088 0.2894
##
## hours = 84:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.4813 0.225 20 2.144 0.0445
##
## hours = 96:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.3574 0.225 20 1.591 0.1273
##
## hours = 120:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.1882 0.225 20 0.838 0.4122
##
## hours = 144:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.0461 0.225 20 0.205 0.8396
##
## hours = 168:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -0.0428 0.225 20 -0.190 0.8509
##
## Degrees-of-freedom method: containment## hours = 0:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.0413 0.227 20 -0.4320 0.515
##
## hours = 36:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.1748 0.224 20 -0.2932 0.643
##
## hours = 48:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.4165 0.224 20 -0.0502 0.883
##
## hours = 60:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.4912 0.224 20 0.0244 0.958
##
## hours = 72:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.2431 0.223 20 -0.2229 0.709
##
## hours = 84:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.4813 0.225 20 0.0130 0.950
##
## hours = 96:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.3574 0.225 20 -0.1112 0.826
##
## hours = 120:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.1882 0.225 20 -0.2805 0.657
##
## hours = 144:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.0461 0.225 20 -0.4225 0.515
##
## hours = 168:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -0.0428 0.225 20 -0.5114 0.426
##
## Degrees-of-freedom method: containment
## Confidence level used: 0.956.2 - Graph overtime
6.3 - Type III p-value
## Analysis of Deviance Table (Type III tests)
##
## Response: meanbot
## Chisq Df Pr(>Chisq)
## (Intercept) 0.1066 1 0.74407
## as.factor(Treatment) 0.0332 1 0.85552
## hours 152.9253 9 < 2e-16 ***
## weekbeforetreatmilk 0.6127 1 0.43378
## as.factor(Treatment):hours 17.4761 9 0.04176 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 16.4 - Testing normality of the model
6.5 - Average pressure
## $emmeans
## Treatment emmean SE df lower.CL upper.CL
## Control 0.877 0.119 20 0.628 1.126
## Treatment 0.659 0.149 20 0.348 0.969
##
## Results are averaged over the levels of: hours, rl
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.218 0.191 20 1.145 0.2656
##
## Results are averaged over the levels of: hours, rl
## Degrees-of-freedom method: containment## $emmeans
## Treatment emmean SE df lower.CL upper.CL
## Control 0.877 0.119 20 0.628 1.126
## Treatment 0.659 0.149 20 0.348 0.969
##
## Results are averaged over the levels of: hours, rl
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.218 0.191 20 -0.179 0.616
##
## Results are averaged over the levels of: hours, rl
## Degrees-of-freedom method: containment
## Confidence level used: 0.957 - Blood pH
7.1 - pH overtime
## $emmeans
## hours = 12:
## Treatment emmean SE df lower.CL upper.CL
## Control 7.486 0.0137 19 7.457 7.514
## Treatment 7.424 0.0137 19 7.395 7.452
##
## hours = 24:
## Treatment emmean SE df lower.CL upper.CL
## Control 7.458 0.0137 19 7.430 7.487
## Treatment 7.426 0.0137 19 7.398 7.455
##
## hours = 48:
## Treatment emmean SE df lower.CL upper.CL
## Control 7.455 0.0137 19 7.426 7.484
## Treatment 7.439 0.0137 19 7.410 7.467
##
## hours = 72:
## Treatment emmean SE df lower.CL upper.CL
## Control 7.446 0.0137 19 7.417 7.474
## Treatment 7.422 0.0137 19 7.394 7.451
##
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## hours = 12:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.0618 0.0194 19 3.186 0.0049
##
## hours = 24:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.0319 0.0194 19 1.645 0.1163
##
## hours = 48:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.0162 0.0194 19 0.835 0.4140
##
## hours = 72:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.0232 0.0194 19 1.196 0.2465
##
## Degrees-of-freedom method: containment## hours = 12:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.0618 0.0194 19 0.02122 0.1025
##
## hours = 24:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.0319 0.0194 19 -0.00869 0.0726
##
## hours = 48:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.0162 0.0194 19 -0.02441 0.0568
##
## hours = 72:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.0232 0.0194 19 -0.01741 0.0638
##
## Degrees-of-freedom method: containment
## Confidence level used: 0.957.2 - Graph overtime
7.3 - Type III p-value
## Analysis of Deviance Table (Type III tests)
##
## Response: bloodph
## Chisq Df Pr(>Chisq)
## (Intercept) 79.5194 1 < 2e-16 ***
## Treatment 10.1536 1 0.00144 **
## hours 4.6667 3 0.19789
## basebloodph 1.4164 1 0.23399
## Treatment:hours 3.1836 3 0.36418
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 17.4 - Testing normality of the model
7.5 - Average blood pH
## $emmeans
## Treatment emmean SE df lower.CL upper.CL
## Control 7.461 0.00734 19 7.446 7.477
## Treatment 7.428 0.00734 19 7.412 7.443
##
## Results are averaged over the levels of: hours
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.0337 0.0104 19 3.245 0.0043
##
## Results are averaged over the levels of: hours
## Degrees-of-freedom method: containment## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.0337 0.0104 19 0.012 0.0555
##
## Results are averaged over the levels of: hours
## Degrees-of-freedom method: containment
## Confidence level used: 0.958 - Urine pH
8.1 - pH overtime
## $emmeans
## hours = 12:
## Treatment emmean SE df lower.CL upper.CL
## Control 8.41 0.136 21 8.12 8.69
## Treatment 7.56 0.136 21 7.28 7.84
##
## hours = 24:
## Treatment emmean SE df lower.CL upper.CL
## Control 8.26 0.136 21 7.98 8.55
## Treatment 7.86 0.136 21 7.58 8.14
##
## hours = 48:
## Treatment emmean SE df lower.CL upper.CL
## Control 8.15 0.136 21 7.87 8.43
## Treatment 7.98 0.136 21 7.69 8.26
##
## hours = 72:
## Treatment emmean SE df lower.CL upper.CL
## Control 7.62 0.141 21 7.33 7.91
## Treatment 7.88 0.136 21 7.59 8.16
##
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## hours = 12:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.846 0.192 21 4.401 0.0002
##
## hours = 24:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.404 0.192 21 2.103 0.0477
##
## hours = 48:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.171 0.192 21 0.889 0.3843
##
## hours = 72:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -0.257 0.196 21 -1.309 0.2046
##
## Degrees-of-freedom method: containment## hours = 12:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.846 0.192 21 0.44609 1.245
##
## hours = 24:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.404 0.192 21 0.00442 0.804
##
## hours = 48:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.171 0.192 21 -0.22891 0.570
##
## hours = 72:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -0.257 0.196 21 -0.66544 0.151
##
## Degrees-of-freedom method: containment
## Confidence level used: 0.958.2 - Graph overtime
8.3 - Type III p-value
## Analysis of Deviance Table (Type III tests)
##
## Response: urineph
## Chisq Df Pr(>Chisq)
## (Intercept) 17.4383 1 2.968e-05 ***
## Treatment 19.3649 1 1.080e-05 ***
## hours 18.5941 3 0.0003317 ***
## baseurine 2.8152 1 0.0933746 .
## Treatment:hours 16.5803 3 0.0008620 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 18.4 - Testing normality of the model
8.5 - Average urine pH
## $emmeans
## Treatment emmean SE df lower.CL upper.CL
## Control 8.12 0.0944 21 7.92 8.31
## Treatment 7.82 0.0937 21 7.62 8.01
##
## Results are averaged over the levels of: hours
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.299 0.133 21 2.242 0.0359
##
## Results are averaged over the levels of: hours
## Degrees-of-freedom method: containment## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.299 0.133 21 0.0217 0.576
##
## Results are averaged over the levels of: hours
## Degrees-of-freedom method: containment
## Confidence level used: 0.959 - Lying Behavior
9.1 - Lying (min/d) overtime
## $emmeans
## daytotreatm = 0:
## Treatment emmean SE df lower.CL upper.CL
## Control 656 46.0 20 560 752
## Treatment 715 47.9 20 615 814
##
## daytotreatm = 1:
## Treatment emmean SE df lower.CL upper.CL
## Control 674 46.0 20 578 770
## Treatment 594 47.9 20 494 694
##
## daytotreatm = 2:
## Treatment emmean SE df lower.CL upper.CL
## Control 728 46.0 20 632 824
## Treatment 704 47.9 20 604 804
##
## daytotreatm = 3:
## Treatment emmean SE df lower.CL upper.CL
## Control 743 46.0 20 647 839
## Treatment 752 47.9 20 652 851
##
## daytotreatm = 4:
## Treatment emmean SE df lower.CL upper.CL
## Control 810 46.0 20 714 906
## Treatment 772 47.9 20 672 872
##
## daytotreatm = 5:
## Treatment emmean SE df lower.CL upper.CL
## Control 744 46.0 20 648 840
## Treatment 769 47.9 20 669 869
##
## daytotreatm = 6:
## Treatment emmean SE df lower.CL upper.CL
## Control 733 46.0 20 637 829
## Treatment 688 47.9 20 588 788
##
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## daytotreatm = 0:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -58.4 66.6 20 -0.877 0.3907
##
## daytotreatm = 1:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 80.3 66.6 20 1.205 0.2421
##
## daytotreatm = 2:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 23.7 66.6 20 0.356 0.7258
##
## daytotreatm = 3:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -8.2 66.6 20 -0.123 0.9032
##
## daytotreatm = 4:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 38.1 66.6 20 0.572 0.5738
##
## daytotreatm = 5:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -25.2 66.6 20 -0.379 0.7090
##
## daytotreatm = 6:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 45.0 66.6 20 0.676 0.5067
##
## Degrees-of-freedom method: containment## daytotreatm = 0:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -58.4 66.6 20 -197.3 80.5
##
## daytotreatm = 1:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 80.3 66.6 20 -58.6 219.1
##
## daytotreatm = 2:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 23.7 66.6 20 -115.2 162.6
##
## daytotreatm = 3:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -8.2 66.6 20 -147.1 130.7
##
## daytotreatm = 4:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 38.1 66.6 20 -100.8 177.0
##
## daytotreatm = 5:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -25.2 66.6 20 -164.1 113.7
##
## daytotreatm = 6:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 45.0 66.6 20 -93.9 183.9
##
## Degrees-of-freedom method: containment
## Confidence level used: 0.959.2 - Graph overtime
## Ignoring unknown labels:
## • linetype : "Treatment"
9.3 - Type III p-value
## Analysis of Deviance Table (Type III tests)
##
## Response: totalminlying
## Chisq Df Pr(>Chisq)
## (Intercept) 8.9215 1 0.002818 **
## as.factor(Treatment) 0.7697 1 0.380312
## as.factor(daytotreatm) 7.8295 6 0.250862
## baselying 0.0000 1 0.996175
## lastweekmilkaveragekg 0.7975 1 0.371833
## dim 0.1966 1 0.657443
## as.factor(Treatment):as.factor(daytotreatm) 4.9630 6 0.548565
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 19.4 - Testing normality of the model
9.5 - Average lying time
## $emmeans
## Treatment emmean SE df lower.CL upper.CL
## Control 725 26.6 16 669 781
## Treatment 701 32.0 16 633 769
##
## Results are averaged over the levels of: daytotreatm
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## Control - Treatment 24 42.1 16 0.571 0.5760
##
## Results are averaged over the levels of: daytotreatm
## Degrees-of-freedom method: containment## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 24 42.1 16 -65.2 113
##
## Results are averaged over the levels of: daytotreatm
## Degrees-of-freedom method: containment
## Confidence level used: 0.9510 - Standing Behavior
10.1 - Standing (min/d) overtime
## $emmeans
## daytotreatm = 0:
## Treatment emmean SE df lower.CL upper.CL
## Control 689 45.8 16 592 786
## Treatment 706 55.1 16 589 823
##
## daytotreatm = 1:
## Treatment emmean SE df lower.CL upper.CL
## Control 611 45.8 16 514 708
## Treatment 510 55.1 16 393 627
##
## daytotreatm = 2:
## Treatment emmean SE df lower.CL upper.CL
## Control 713 45.8 16 616 810
## Treatment 751 55.1 16 635 868
##
## daytotreatm = 3:
## Treatment emmean SE df lower.CL upper.CL
## Control 697 45.8 16 600 794
## Treatment 693 55.1 16 576 810
##
## daytotreatm = 4:
## Treatment emmean SE df lower.CL upper.CL
## Control 631 45.8 16 534 728
## Treatment 641 55.1 16 524 757
##
## daytotreatm = 5:
## Treatment emmean SE df lower.CL upper.CL
## Control 697 45.8 16 600 794
## Treatment 660 55.1 16 543 776
##
## daytotreatm = 6:
## Treatment emmean SE df lower.CL upper.CL
## Control 708 45.8 16 611 805
## Treatment 754 55.1 16 637 871
##
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## daytotreatm = 0:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -16.52 71.8 16 -0.230 0.8210
##
## daytotreatm = 1:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 101.25 71.8 16 1.410 0.1778
##
## daytotreatm = 2:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -38.46 71.8 16 -0.536 0.5997
##
## daytotreatm = 3:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 4.52 71.8 16 0.063 0.9505
##
## daytotreatm = 4:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -10.01 71.8 16 -0.139 0.8909
##
## daytotreatm = 5:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 37.68 71.8 16 0.525 0.6070
##
## daytotreatm = 6:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -45.87 71.8 16 -0.639 0.5321
##
## Degrees-of-freedom method: containment## daytotreatm = 0:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -16.52 71.8 16 -169 136
##
## daytotreatm = 1:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 101.25 71.8 16 -51 254
##
## daytotreatm = 2:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -38.46 71.8 16 -191 114
##
## daytotreatm = 3:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 4.52 71.8 16 -148 157
##
## daytotreatm = 4:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -10.01 71.8 16 -162 142
##
## daytotreatm = 5:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 37.68 71.8 16 -115 190
##
## daytotreatm = 6:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -45.87 71.8 16 -198 106
##
## Degrees-of-freedom method: containment
## Confidence level used: 0.9510.2 - Graph overtime
## Ignoring unknown labels:
## • linetype : "Treatment"
10.3 - Type III p-value
## Analysis of Deviance Table (Type III tests)
##
## Response: totalminstand
## Chisq Df Pr(>Chisq)
## (Intercept) 4.2248 1 0.03984 *
## as.factor(Treatment) 0.0529 1 0.81810
## as.factor(daytotreatm) 6.4162 6 0.37822
## basestand 0.3140 1 0.57526
## lastweekmilkaveragekg 0.6190 1 0.43143
## dim 0.7790 1 0.37746
## Intake.average.baseline 1.8258 1 0.17663
## as.factor(Treatment):as.factor(daytotreatm) 4.3388 6 0.63093
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 110.4 - Testing normality of the model
10.5 - Average standing time
## $emmeans
## Treatment emmean SE df lower.CL upper.CL
## Control 678 26.4 16 622 734
## Treatment 674 31.9 16 607 742
##
## Results are averaged over the levels of: daytotreatm
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 3.5 41.8 16 -85 92
##
## Results are averaged over the levels of: daytotreatm
## Degrees-of-freedom method: containment
## Confidence level used: 0.95## $emmeans
## Treatment emmean SE df lower.CL upper.CL
## Control 678 26.4 16 622 734
## Treatment 674 31.9 16 607 742
##
## Results are averaged over the levels of: daytotreatm
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## Control - Treatment 3.5 41.8 16 0.084 0.9344
##
## Results are averaged over the levels of: daytotreatm
## Degrees-of-freedom method: containment11 - Number of bouts
11.1 - Bouts (/d) overtime
## $emmeans
## daytotreatm = 0:
## Treatment emmean SE df lower.CL upper.CL
## Control 12.8 1.03 16 10.61 15.0
## Treatment 13.1 1.24 16 10.46 15.7
##
## daytotreatm = 1:
## Treatment emmean SE df lower.CL upper.CL
## Control 12.9 1.03 16 10.76 15.1
## Treatment 10.6 1.24 16 8.01 13.3
##
## daytotreatm = 2:
## Treatment emmean SE df lower.CL upper.CL
## Control 11.6 1.03 16 9.38 13.7
## Treatment 10.2 1.24 16 7.57 12.8
##
## daytotreatm = 3:
## Treatment emmean SE df lower.CL upper.CL
## Control 12.0 1.03 16 9.84 14.2
## Treatment 11.3 1.24 16 8.68 13.9
##
## daytotreatm = 4:
## Treatment emmean SE df lower.CL upper.CL
## Control 12.3 1.03 16 10.15 14.5
## Treatment 12.2 1.24 16 9.57 14.8
##
## daytotreatm = 5:
## Treatment emmean SE df lower.CL upper.CL
## Control 11.2 1.03 16 8.99 13.3
## Treatment 11.0 1.24 16 8.35 13.6
##
## daytotreatm = 6:
## Treatment emmean SE df lower.CL upper.CL
## Control 10.2 1.03 16 7.99 12.3
## Treatment 10.9 1.24 16 8.23 13.5
##
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## daytotreatm = 0:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -0.308 1.64 16 -0.187 0.8540
##
## daytotreatm = 1:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 2.291 1.64 16 1.393 0.1828
##
## daytotreatm = 2:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 1.350 1.64 16 0.821 0.4237
##
## daytotreatm = 3:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.701 1.64 16 0.426 0.6757
##
## daytotreatm = 4:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.120 1.64 16 0.073 0.9429
##
## daytotreatm = 5:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.188 1.64 16 0.114 0.9104
##
## daytotreatm = 6:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -0.701 1.64 16 -0.426 0.6758
##
## Degrees-of-freedom method: containment## daytotreatm = 0:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -0.308 1.64 16 -3.79 3.18
##
## daytotreatm = 1:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 2.291 1.64 16 -1.20 5.78
##
## daytotreatm = 2:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 1.350 1.64 16 -2.14 4.84
##
## daytotreatm = 3:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.701 1.64 16 -2.79 4.19
##
## daytotreatm = 4:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.120 1.64 16 -3.37 3.61
##
## daytotreatm = 5:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.188 1.64 16 -3.30 3.68
##
## daytotreatm = 6:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -0.701 1.64 16 -4.19 2.79
##
## Degrees-of-freedom method: containment
## Confidence level used: 0.9511.2 - Graph overtime
## Ignoring unknown labels:
## • linetype : "Treatment"
11.3 - Type III p-value
## Analysis of Deviance Table (Type III tests)
##
## Response: totalbouts
## Chisq Df Pr(>Chisq)
## (Intercept) 8.5373 1 0.0034793 **
## as.factor(Treatment) 0.0350 1 0.8516490
## as.factor(daytotreatm) 6.5307 6 0.3664281
## basebouts 11.3303 1 0.0007625 ***
## lastweekmilkaveragekg 2.4379 1 0.1184325
## dim 6.8459 1 0.0088846 **
## Intake.average.baseline 3.6361 1 0.0565379 .
## as.factor(Treatment):as.factor(daytotreatm) 3.2925 6 0.7713345
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 111.4 - Testing normality of the model
11.5 - Average bouts number
## $emmeans
## Treatment emmean SE df lower.CL upper.CL
## Control 11.8 0.640 16 10.46 13.2
## Treatment 11.4 0.784 16 9.72 13.0
##
## Results are averaged over the levels of: daytotreatm
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.431 1.06 16 0.406 0.6904
##
## Results are averaged over the levels of: daytotreatm
## Degrees-of-freedom method: containment## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.431 1.06 16 -1.82 2.69
##
## Results are averaged over the levels of: daytotreatm
## Degrees-of-freedom method: containment
## Confidence level used: 0.9512 - Bouts Length
12.1 - Bouts Length (/d) overtime
## $emmeans
## daytotreatm = 0:
## Treatment emmean SE df lower.CL upper.CL
## Control 53.3 5.50 21 41.9 64.7
## Treatment 58.5 5.73 21 46.6 70.4
##
## daytotreatm = 1:
## Treatment emmean SE df lower.CL upper.CL
## Control 55.2 5.50 21 43.8 66.6
## Treatment 57.2 5.73 21 45.3 69.1
##
## daytotreatm = 2:
## Treatment emmean SE df lower.CL upper.CL
## Control 62.1 5.50 21 50.6 73.5
## Treatment 76.5 5.73 21 64.6 88.4
##
## daytotreatm = 3:
## Treatment emmean SE df lower.CL upper.CL
## Control 61.1 5.50 21 49.6 72.5
## Treatment 71.4 5.73 21 59.5 83.3
##
## daytotreatm = 4:
## Treatment emmean SE df lower.CL upper.CL
## Control 64.2 5.50 21 52.8 75.7
## Treatment 74.0 5.73 21 62.1 85.9
##
## daytotreatm = 5:
## Treatment emmean SE df lower.CL upper.CL
## Control 64.5 5.50 21 53.0 75.9
## Treatment 80.6 5.73 21 68.7 92.5
##
## daytotreatm = 6:
## Treatment emmean SE df lower.CL upper.CL
## Control 68.9 5.50 21 57.5 80.4
## Treatment 73.1 5.73 21 61.2 85.0
##
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## daytotreatm = 0:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -5.16 7.95 21 -0.649 0.5232
##
## daytotreatm = 1:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -1.99 7.95 21 -0.251 0.8044
##
## daytotreatm = 2:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -14.46 7.95 21 -1.819 0.0832
##
## daytotreatm = 3:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -10.34 7.95 21 -1.301 0.2073
##
## daytotreatm = 4:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -9.74 7.95 21 -1.226 0.2339
##
## daytotreatm = 5:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -16.15 7.95 21 -2.032 0.0550
##
## daytotreatm = 6:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -4.15 7.95 21 -0.522 0.6074
##
## Degrees-of-freedom method: containment## daytotreatm = 0:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -5.16 7.95 21 -21.7 11.368
##
## daytotreatm = 1:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -1.99 7.95 21 -18.5 14.536
##
## daytotreatm = 2:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -14.46 7.95 21 -31.0 2.073
##
## daytotreatm = 3:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -10.34 7.95 21 -26.9 6.187
##
## daytotreatm = 4:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -9.74 7.95 21 -26.3 6.787
##
## daytotreatm = 5:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -16.15 7.95 21 -32.7 0.381
##
## daytotreatm = 6:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -4.15 7.95 21 -20.7 12.384
##
## Degrees-of-freedom method: containment
## Confidence level used: 0.9512.2 - Graph overtime
## Ignoring unknown labels:
## • linetype : "Treatment"
12.3 - Type III p-value
## Analysis of Deviance Table (Type III tests)
##
## Response: boutlength
## Chisq Df Pr(>Chisq)
## (Intercept) 0.3903 1 0.5322
## as.factor(Treatment) 0.4216 1 0.5161
## as.factor(daytotreatm) 9.4596 6 0.1493
## lastweekmilkaveragekg 0.0002 1 0.9881
## dim 1.3764 1 0.2407
## as.factor(Treatment):as.factor(daytotreatm) 6.4720 6 0.372412.4 - Testing normality of the model
12.5 - Average bouts length
## $emmeans
## Treatment emmean SE df lower.CL upper.CL
## Control 61.5 4.24 21 52.7 70.3
## Treatment 70.0 4.42 21 60.8 79.2
##
## Results are averaged over the levels of: daytotreatm
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -8.45 6.13 21 -21.2 4.3
##
## Results are averaged over the levels of: daytotreatm
## Degrees-of-freedom method: containment
## Confidence level used: 0.95## $emmeans
## Treatment emmean SE df lower.CL upper.CL
## Control 61.5 4.24 21 52.7 70.3
## Treatment 70.0 4.42 21 60.8 79.2
##
## Results are averaged over the levels of: daytotreatm
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## Control - Treatment -8.45 6.13 21 -1.379 0.1825
##
## Results are averaged over the levels of: daytotreatm
## Degrees-of-freedom method: containment13 - ‘As fed’ (pounds/d)
13.1 - Intake (lbs/d) overtime
## $emmeans
## Day.related.treat = 0:
## Treatment emmean SE df lower.CL upper.CL
## Control 78.6 3.39 19 71.5 85.7
## Treatment 79.2 4.23 19 70.4 88.1
##
## Day.related.treat = 1:
## Treatment emmean SE df lower.CL upper.CL
## Control 48.4 3.39 19 41.3 55.5
## Treatment 43.5 4.23 19 34.6 52.3
##
## Day.related.treat = 2:
## Treatment emmean SE df lower.CL upper.CL
## Control 38.8 3.39 19 31.7 45.9
## Treatment 43.1 4.23 19 34.3 52.0
##
## Day.related.treat = 3:
## Treatment emmean SE df lower.CL upper.CL
## Control 36.6 3.39 19 29.5 43.6
## Treatment 45.8 4.23 19 37.0 54.6
##
## Day.related.treat = 4:
## Treatment emmean SE df lower.CL upper.CL
## Control 42.8 3.39 19 35.7 49.9
## Treatment 50.7 4.23 19 41.8 59.5
##
## Day.related.treat = 5:
## Treatment emmean SE df lower.CL upper.CL
## Control 44.3 3.39 19 37.2 51.4
## Treatment 45.6 4.23 19 36.7 54.4
##
## Day.related.treat = 6:
## Treatment emmean SE df lower.CL upper.CL
## Control 47.6 3.39 19 40.5 54.6
## Treatment 46.6 4.23 19 37.7 55.4
##
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## Day.related.treat = 0:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -0.617 5.43 19 -0.114 0.9107
##
## Day.related.treat = 1:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 4.946 5.43 19 0.911 0.3735
##
## Day.related.treat = 2:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -4.363 5.43 19 -0.804 0.4314
##
## Day.related.treat = 3:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -9.244 5.43 19 -1.703 0.1048
##
## Day.related.treat = 4:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -7.919 5.43 19 -1.459 0.1609
##
## Day.related.treat = 5:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -1.307 5.43 19 -0.241 0.8122
##
## Day.related.treat = 6:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.978 5.43 19 0.180 0.8589
##
## Degrees-of-freedom method: containment## Day.related.treat = 0:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -0.617 5.43 19 -11.98 10.74
##
## Day.related.treat = 1:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 4.946 5.43 19 -6.41 16.31
##
## Day.related.treat = 2:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -4.363 5.43 19 -15.72 7.00
##
## Day.related.treat = 3:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -9.244 5.43 19 -20.60 2.12
##
## Day.related.treat = 4:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -7.919 5.43 19 -19.28 3.44
##
## Day.related.treat = 5:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -1.307 5.43 19 -12.67 10.05
##
## Day.related.treat = 6:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.978 5.43 19 -10.38 12.34
##
## Degrees-of-freedom method: containment
## Confidence level used: 0.9513.2 - Graph overtime
13.3 - Type III p-value
## Analysis of Deviance Table (Type III tests)
##
## Response: Intake
## Chisq Df Pr(>Chisq)
## (Intercept) 6.8869 1 0.0086828 **
## as.factor(Treatment) 0.0129 1 0.9094844
## Day.related.treat 114.5628 6 < 2.2e-16 ***
## Intake.average.baseline 11.5619 1 0.0006732 ***
## dimatdo 0.3336 1 0.5635424
## as.factor(Treatment):Day.related.treat 5.5003 6 0.4814243
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 113.4 - Testing normality of the model
13.5 - Average feed intake
## $emmeans
## Treatment emmean SE df lower.CL upper.CL
## Control 48.2 1.86 19 44.3 52.1
## Treatment 50.5 2.32 19 45.6 55.4
##
## Results are averaged over the levels of: Day.related.treat
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## Control - Treatment -2.25 2.99 19 -0.754 0.4603
##
## Results are averaged over the levels of: Day.related.treat
## Degrees-of-freedom method: containment## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -2.25 2.99 19 -8.51 4
##
## Results are averaged over the levels of: Day.related.treat
## Degrees-of-freedom method: containment
## Confidence level used: 0.9514 - ‘As fed’ (kg/d)
14.1 - Intake (kg/d) overtime
## $emmeans
## Day.related.treat = 0:
## Treatment emmean SE df lower.CL upper.CL
## Control 35.7 1.54 19 32.4 38.9
## Treatment 35.9 1.92 19 31.9 40.0
##
## Day.related.treat = 1:
## Treatment emmean SE df lower.CL upper.CL
## Control 22.0 1.54 19 18.7 25.2
## Treatment 19.7 1.92 19 15.7 23.7
##
## Day.related.treat = 2:
## Treatment emmean SE df lower.CL upper.CL
## Control 17.6 1.54 19 14.4 20.8
## Treatment 19.6 1.92 19 15.5 23.6
##
## Day.related.treat = 3:
## Treatment emmean SE df lower.CL upper.CL
## Control 16.6 1.54 19 13.4 19.8
## Treatment 20.8 1.92 19 16.8 24.8
##
## Day.related.treat = 4:
## Treatment emmean SE df lower.CL upper.CL
## Control 19.4 1.54 19 16.2 22.6
## Treatment 23.0 1.92 19 19.0 27.0
##
## Day.related.treat = 5:
## Treatment emmean SE df lower.CL upper.CL
## Control 20.1 1.54 19 16.9 23.3
## Treatment 20.7 1.92 19 16.7 24.7
##
## Day.related.treat = 6:
## Treatment emmean SE df lower.CL upper.CL
## Control 21.6 1.54 19 18.4 24.8
## Treatment 21.1 1.92 19 17.1 25.1
##
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## Day.related.treat = 0:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -0.280 2.46 19 -0.114 0.9107
##
## Day.related.treat = 1:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 2.243 2.46 19 0.911 0.3735
##
## Day.related.treat = 2:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -1.979 2.46 19 -0.804 0.4314
##
## Day.related.treat = 3:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -4.192 2.46 19 -1.703 0.1048
##
## Day.related.treat = 4:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -3.591 2.46 19 -1.459 0.1609
##
## Day.related.treat = 5:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -0.593 2.46 19 -0.241 0.8122
##
## Day.related.treat = 6:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 0.444 2.46 19 0.180 0.8589
##
## Degrees-of-freedom method: containment## Day.related.treat = 0:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -0.280 2.46 19 -5.43 4.872
##
## Day.related.treat = 1:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 2.243 2.46 19 -2.91 7.395
##
## Day.related.treat = 2:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -1.979 2.46 19 -7.13 3.173
##
## Day.related.treat = 3:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -4.192 2.46 19 -9.34 0.959
##
## Day.related.treat = 4:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -3.591 2.46 19 -8.74 1.560
##
## Day.related.treat = 5:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -0.593 2.46 19 -5.74 4.559
##
## Day.related.treat = 6:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 0.444 2.46 19 -4.71 5.595
##
## Degrees-of-freedom method: containment
## Confidence level used: 0.9514.2 - Graph overtime
14.3 - Type III p-value
## Analysis of Deviance Table (Type III tests)
##
## Response: Intakekg
## Chisq Df Pr(>Chisq)
## (Intercept) 6.8869 1 0.0086828 **
## as.factor(Treatment) 0.0129 1 0.9094844
## Day.related.treat 114.5628 6 < 2.2e-16 ***
## Intake.average.baseline 11.5619 1 0.0006732 ***
## dimatdo 0.3336 1 0.5635424
## as.factor(Treatment):Day.related.treat 5.5003 6 0.4814243
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 114.4 - Testing normality of the model
14.5 - Average feed intake
## $emmeans
## Treatment emmean SE df lower.CL upper.CL
## Control 21.9 0.844 19 20.1 23.6
## Treatment 22.9 1.050 19 20.7 25.1
##
## Results are averaged over the levels of: Day.related.treat
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## Control - Treatment -1.02 1.36 19 -0.754 0.4603
##
## Results are averaged over the levels of: Day.related.treat
## Degrees-of-freedom method: containment## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -1.02 1.36 19 -3.86 1.82
##
## Results are averaged over the levels of: Day.related.treat
## Degrees-of-freedom method: containment
## Confidence level used: 0.9515 - Milk Production
15.1 - Milk production Drop in Pounds
15.1.1 - Milk drop 24h after treatment
## $emmeans
## Treatment emmean SE df lower.CL upper.CL
## Control 2.40 3.44 23 -4.72 9.51
## Treatment -5.43 3.85 23 -13.39 2.53
##
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## Control - Treatment 7.83 5.17 23 1.514 0.1436## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 7.83 5.17 23 -2.87 18.5
##
## Confidence level used: 0.9515.1.2 - Bar Graph
## Treatment emmean SE df lower.CL upper.CL
## Control 2.40 3.44 23 -4.72 9.51
## Treatment -5.43 3.85 23 -13.39 2.53
##
## Confidence level used: 0.95
15.1.3 - Type III p-value
## Anova Table (Type III tests)
##
## Response: dropyield
## Sum Sq Df F value Pr(>F)
## (Intercept) 1.3 1 0.0075 0.9315
## Treatment 405.4 1 2.2923 0.1436
## dim 23.4 1 0.1324 0.7193
## lastweekmilkaverage 20.7 1 0.1169 0.7355
## Residuals 4067.2 2315.1.4 - Testing normality of the model
15.2 - Milk production overtime in Pounds
15.2.1 - Yield overtime
## $emmeans
## beforeafter = before:
## Treatment emmean SE df lower.CL upper.CL
## Control 51.4 2.56 23 46.1 56.7
## Treatment 56.2 2.86 23 50.3 62.2
##
## beforeafter = after:
## Treatment emmean SE df lower.CL upper.CL
## Control 53.7 2.56 23 48.4 59.0
## Treatment 50.9 2.86 23 45.0 56.8
##
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## beforeafter = before:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -4.84 3.84 23 -1.261 0.2199
##
## beforeafter = after:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 2.81 3.84 23 0.732 0.4718
##
## Degrees-of-freedom method: containment## beforeafter = before:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -4.84 3.84 23 -12.78 3.1
##
## beforeafter = after:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 2.81 3.84 23 -5.13 10.7
##
## Degrees-of-freedom method: containment
## Confidence level used: 0.9515.2.2 - Graph overtime
## beforeafter = before:
## Treatment emmean SE df lower.CL upper.CL
## Control 51.4 2.56 23 46.1 56.7
## Treatment 56.2 2.86 23 50.3 62.2
##
## beforeafter = after:
## Treatment emmean SE df lower.CL upper.CL
## Control 53.7 2.56 23 48.4 59.0
## Treatment 50.9 2.86 23 45.0 56.8
##
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
15.2.3 - Type III p-value
## Analysis of Deviance Table (Type III tests)
##
## Response: milk
## Chisq Df Pr(>Chisq)
## (Intercept) 0.6384 1 0.4243
## as.factor(Treatment) 1.5902 1 0.2073
## beforeafter 0.4872 1 0.4852
## dim 0.5009 1 0.4791
## lastweekmilkaverage 61.7629 1 3.874e-15 ***
## as.factor(Treatment):beforeafter 2.3553 1 0.1249
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 115.2.4 - Testing normality of the model
15.3 - Milk production Drop in Kg
15.3.1 - Milk drop 24h after treatment
## $emmeans
## Treatment emmean SE df lower.CL upper.CL
## Control 1.09 1.56 23 -2.14 4.31
## Treatment -2.46 1.74 23 -6.07 1.15
##
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## Control - Treatment 3.55 2.34 23 1.514 0.1436## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 3.55 2.34 23 -1.3 8.4
##
## Confidence level used: 0.9515.3.2 - Bar Graph
## Treatment emmean SE df lower.CL upper.CL
## Control 1.09 1.56 23 -2.14 4.31
## Treatment -2.46 1.74 23 -6.07 1.15
##
## Confidence level used: 0.95
15.3.3 - Type III p-value
## Anova Table (Type III tests)
##
## Response: dropyieldkg
## Sum Sq Df F value Pr(>F)
## (Intercept) 0.27 1 0.0075 0.9315
## Treatment 83.37 1 2.2923 0.1436
## dim 4.82 1 0.1324 0.7193
## lastweekmilkaveragekg 4.25 1 0.1169 0.7355
## Residuals 836.52 2315.3.4 - Testing normality of the model
15.4 - Milk production overtime in Kg
15.4.1 - Yield overtime
## $emmeans
## beforeafter = before:
## Treatment emmean SE df lower.CL upper.CL
## Control 23.3 1.16 23 20.9 25.7
## Treatment 25.5 1.30 23 22.8 28.2
##
## beforeafter = after:
## Treatment emmean SE df lower.CL upper.CL
## Control 24.4 1.16 23 22.0 26.8
## Treatment 23.1 1.30 23 20.4 25.8
##
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
##
## $contrasts
## beforeafter = before:
## contrast estimate SE df t.ratio p.value
## Control - Treatment -2.20 1.74 23 -1.261 0.2199
##
## beforeafter = after:
## contrast estimate SE df t.ratio p.value
## Control - Treatment 1.27 1.74 23 0.732 0.4718
##
## Degrees-of-freedom method: containment## beforeafter = before:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment -2.20 1.74 23 -5.80 1.41
##
## beforeafter = after:
## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 1.27 1.74 23 -2.33 4.88
##
## Degrees-of-freedom method: containment
## Confidence level used: 0.9515.4.2 - Graph overtime
## beforeafter = before:
## Treatment emmean SE df lower.CL upper.CL
## Control 23.3 1.16 23 20.9 25.7
## Treatment 25.5 1.30 23 22.8 28.2
##
## beforeafter = after:
## Treatment emmean SE df lower.CL upper.CL
## Control 24.4 1.16 23 22.0 26.8
## Treatment 23.1 1.30 23 20.4 25.8
##
## Degrees-of-freedom method: containment
## Confidence level used: 0.95
15.4.3 - Type III p-value
## Analysis of Deviance Table (Type III tests)
##
## Response: milkkg
## Chisq Df Pr(>Chisq)
## (Intercept) 0.6384 1 0.4243
## as.factor(Treatment) 1.5902 1 0.2073
## beforeafter 0.4872 1 0.4852
## dim 0.5009 1 0.4791
## lastweekmilkaveragekg 61.7629 1 3.874e-15 ***
## as.factor(Treatment):beforeafter 2.3553 1 0.1249
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 115.4.4 - Testing normality of the model
15.5 - Milk production percentage change
15.5.1 - Milk % drop 24h after treatment
## $emmeans
## Treatment emmean SE df lower.CL upper.CL
## Control 33.58 17.9 23 -3.51 70.7
## Treatment -9.89 20.1 23 -51.37 31.6
##
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## Control - Treatment 43.5 27 23 1.613 0.1204## contrast estimate SE df lower.CL upper.CL
## Control - Treatment 43.5 27 23 -12.3 99.2
##
## Confidence level used: 0.9515.5.2 - Bar Graph
## Treatment emmean SE df lower.CL upper.CL
## Control 33.58 17.9 23 -3.51 70.7
## Treatment -9.89 20.1 23 -51.37 31.6
##
## Confidence level used: 0.95
15.5.3 - Type III p-value
## Anova Table (Type III tests)
##
## Response: percentage_change
## Sum Sq Df F value Pr(>F)
## (Intercept) 5905 1 1.2285 0.27916
## Treatment 12503 1 2.6011 0.12043
## dim 50 1 0.0103 0.91992
## lastweekmilkaverage 35034 1 7.2882 0.01279 *
## Residuals 110560 23
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 115.5.4 - Testing normality of the model
