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.55 0.21 Chi Sq, Cramer’s V 1.00
1 11 (14) 2 (8) 4 (15) 5 (19)
2 14 (18) 3 (12) 8 (30) 3 (11)
3 13 (16) 6 (24) 4 (15) 3 (11)
4 22 (28) 7 (28) 6 (22) 9 (33)
5 19 (24) 7 (28) 5 (19) 7 (26)

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 B08 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
1 B23 1 1 1 1 1 1 1 1 1 1 1 1 1 1
1 BL05 2 2 2 2 2 2 2 2 2 2 2 2 2 2
1 BL12 2 2 2 2 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 BL24 2 2 2 2 2 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 L04 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
1 L05 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
1 L16 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
1 L27 3 3 3 3 3 3 3 3 3 3 3 3 3 3
2 B11 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
2 B18 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 BL04 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 BL06 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 BL11 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 BL14 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 BL16 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 BL18 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 BL20 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 L12 3 3 3 3 3 3 3 3 3 3 3 3 3
2 L18 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
2 L26 3 3 3 3 3 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 B03 1 1 1 1 1 1 1 1 1 1 1 1 1 1
3 B07 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
3 B09 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
3 B13 1 1 1 1 1 1 1 1 1 1 1 1 1 1
3 B20 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
3 BL09 2 2 2 2 2 2 2 2 2 2 2 2 2 2
3 BL13 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 BL19 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
3 L10 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
3 L21 3 3 3 3 3 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 B02 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
4 B04 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
4 B06 1 1 1 1 1 1 1 1 1 1 1 1 1 1
4 B14 1 1 1 1 1 1 1 1 1 1 1 1 1 1
4 B16 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
4 B17 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
4 B25 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
4 BL07 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
4 BL10 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 BL21 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
4 BL25 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 L06 3 3 3 3 3 3 3 3 3 3 3 3 3 3
4 L14 3 3 3 3 3 3 3 3 3 3 3 3 3 3
4 L15 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
4 L17 3 3 3 3 3 3 3 3 3 3 3 3
4 L19 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
4 L20 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
4 L22 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
4 L23 3 3 3 3 3 3 3 3 3 3 3 3 3 3
4 L25 3 3 3 3 3 3 3 3 3 3 3 3 3 3
5 B05 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
5 B10 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
5 B12 1 1 1 1 1 1 1 1 1 1 1 1
5 B15 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
5 B19 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
5 B21 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
5 B24 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
5 BL01 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
5 BL02 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
5 BL03 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
5 BL08 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
5 BL23 2 2 2 2 2 2 2 2 2 2 2 2 2 2
5 L01 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
5 L02 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
5 L07 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
5 L08 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
5 L09 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
5 L11 3 3 3 3 3 3 3 3 3 3 3 3 3 3
5 L13 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3