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.008 0.36 Chi Sq, Cramer’s V 0.041
1 19 (24) 5 (20) 12 (44) 2 (7)
2 11 (14) 2 (8) 4 (15) 5 (19)
3 17 (22) 10 (40) 2 (7) 5 (19)
4 11 (14) 4 (16) 1 (4) 6 (22)
5 21 (27) 4 (16) 8 (30) 9 (33)

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 B01 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
1 B14 1 1 1 1 1 1 1 1 1 1 1 1 1 1
1 B15 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
1 B20 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 BL02 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 BL07 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
1 BL09 2 2 2 2 2 2 2 2 2 2 2 2 2 2
1 BL11 2 2 2 2 2 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 BL17 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
1 BL18 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
1 BL19 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 BL22 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
1 BL27 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
1 L15 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 B07 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 BL01 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 BL04 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 BL08 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 BL10 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
2 L03 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 L16 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
2 L24 3 3 3 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 B03 1 1 1 1 1 1 1 1 1 1 1 1 1 1
3 B05 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 B08 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 B10 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
3 B11 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
3 B18 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
3 B21 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
3 B22 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 L09 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
3 L11 3 3 3 3 3 3 3 3 3 3 3 3 3 3
3 L13 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
3 L14 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
4 B02 1 1 1 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 B24 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 BL21 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
4 L01 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
4 L05 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
4 L10 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
4 L18 3 3 3 3 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 L22 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
5 B04 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
5 B13 1 1 1 1 1 1 1 1 1 1 1 1 1 1
5 B16 1 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 BL13 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
5 BL14 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
5 BL15 2 2 2 2 2 2 2 2 2 2 2 2 2 2
5 BL16 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 BL24 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
5 BL25 2 2 2 2 2 2 2 2 2 2 2 2 2 2
5 BL26 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
5 L02 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
5 L04 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 L12 3 3 3 3 3 3 3 3 3 3 3 3 3
5 L17 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
5 L23 3 3 3 3 3 3 3 3 3 3 3 3 3 3
5 L25 3 3 3 3 3 3 3 3 3 3 3 3 3 3