2.4 Nested summary of diet

2.4.1 By cow:

It is of interest to note that each of the three diet types is given in each yard. Be aware of reverse causation - high protein milk does not necessarily influence diet type!

rm_covsum_nested(data = Milk, id = c("Cow"), covs = c("protein", "High_Protein", "Time", "Yard"), maincov = "Diet")
Full Sample (n=79) barley (n=25) barley+lupins (n=27) lupins (n=27) Unnested p-value Unnested Effect Size Unnested StatTest Nested p-value
protein 0.003 0.17 Kruskal Wallis, Eta sq 0.78
Mean (sd) 3.4 (0.2) 3.5 (0.2) 3.4 (0.1) 3.3 (0.2)
Median (Min,Max) 3.4 (2.9, 3.9) 3.6 (3.2, 3.9) 3.4 (3.2, 3.7) 3.3 (2.9, 3.7)
High_Protein 0.010 0.34 Chi Sq, Cramer’s V 0.80
N 63 (80) 15 (60) 23 (85) 25 (93)
Y 16 (20) 10 (40) 4 (15) 2 (7)
Time 0.97 0.001 Kruskal Wallis, Eta sq 0.99
Mean (sd) 9.0 (1.1) 9.1 (1.2) 9.0 (1.1) 9.0 (1.1)
Median (Min,Max) 10.0 (7.1, 10.4) 10.0 (7.1, 10.4) 10.0 (7.5, 10.1) 10.0 (7.5, 10.2)
Yard 0.66 0.19 Chi Sq, Cramer’s V 1.00
1 22 (28) 6 (24) 9 (33) 7 (26)
2 22 (28) 10 (40) 4 (15) 8 (30)
3 17 (22) 5 (20) 6 (22) 6 (22)
4 9 (11) 3 (12) 4 (15) 2 (7)
5 9 (11) 1 (4) 4 (15) 4 (15)

However, it appear that each cow was assigned to only one diet over all observation instances. Diet 1 = barley, 2 = barley+lupins and 3 = lupins.

Diet_t <- reshape2::dcast(Milk, Yard + Cow ~ Time, value.var="Diet_Type");
colnames(Diet_t)[-c(1:2)] <- paste("", colnames(Diet_t[, -c(1:2)]), "", sep = "");

options(knitr.kable.NA = '');
knitr::kable(Diet_t);
Yard Cow 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
1 B04 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
1 B08 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
1 B09 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
1 B13 1 1 1 1 1 1 1 1 1 1 1 1 1 1
1 B18 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
1 B21 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
1 BL01 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
1 BL05 2 2 2 2 2 2 2 2 2 2 2 2 2 2
1 BL08 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
1 BL14 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
1 BL16 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
1 BL20 2 2 2 2 2 2 2 2 2 2 2 2 2 2
1 BL21 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
1 BL22 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
1 BL25 2 2 2 2 2 2 2 2 2 2 2 2 2 2
1 L03 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
1 L14 3 3 3 3 3 3 3 3 3 3 3 3 3 3
1 L15 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
1 L17 3 3 3 3 3 3 3 3 3 3 3 3
1 L18 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
1 L22 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
1 L26 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
2 B02 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
2 B03 1 1 1 1 1 1 1 1 1 1 1 1 1 1
2 B07 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
2 B11 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
2 B12 1 1 1 1 1 1 1 1 1 1 1 1
2 B14 1 1 1 1 1 1 1 1 1 1 1 1 1 1
2 B15 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
2 B17 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
2 B22 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
2 B25 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
2 BL02 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 BL07 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 BL23 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 BL26 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 L01 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
2 L06 3 3 3 3 3 3 3 3 3 3 3 3 3 3
2 L10 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
2 L11 3 3 3 3 3 3 3 3 3 3 3 3 3 3
2 L13 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
2 L19 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
2 L25 3 3 3 3 3 3 3 3 3 3 3 3 3 3
2 L27 3 3 3 3 3 3 3 3 3 3 3 3 3 3
3 B01 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
3 B06 1 1 1 1 1 1 1 1 1 1 1 1 1 1
3 B19 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
3 B23 1 1 1 1 1 1 1 1 1 1 1 1 1 1
3 B24 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
3 BL03 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
3 BL06 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
3 BL09 2 2 2 2 2 2 2 2 2 2 2 2 2 2
3 BL10 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
3 BL12 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
3 BL17 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
3 L04 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
3 L05 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
3 L08 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
3 L09 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
3 L23 3 3 3 3 3 3 3 3 3 3 3 3 3 3
3 L24 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
4 B05 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
4 B10 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
4 B20 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
4 BL04 2 2 2 2 2 2 2 2 2 2 2 2 2 2
4 BL13 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
4 BL15 2 2 2 2 2 2 2 2 2 2 2 2 2 2
4 BL27 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
4 L02 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
4 L21 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
5 B16 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
5 BL11 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
5 BL18 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
5 BL19 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
5 BL24 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
5 L07 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
5 L12 3 3 3 3 3 3 3 3 3 3 3 3 3
5 L16 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
5 L20 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3