D.11 Jersey: Milk Yield

D.11.1: Full Period: 1998-2023

Model Summary

A. parametric coefficientsEstimateStd. Errort-valuep-value
(Intercept)13.76210.712319.3200\(< 0.0001\)
parityprimiparous-2.01960.0125-162.2099\(< 0.0001\)
year1999-0.57880.7821-0.74010.4592
year20000.17060.79080.21580.8292
year20010.68720.76110.90290.3666
year20021.13450.76211.48860.1366
year20031.62740.76882.11690.0343
year20042.13810.74072.88680.0039
year20052.32940.72833.19860.0014
year20062.61620.72763.59570.0003
year20072.43840.72393.36830.0008
year20082.67800.72313.70370.0002
year20092.93020.72814.0244\(< 0.0001\)
year20103.15020.73384.2930\(< 0.0001\)
year20113.35870.72754.6169\(< 0.0001\)
year20123.62060.72105.0218\(< 0.0001\)
year20133.69430.71855.1417\(< 0.0001\)
year20144.30000.72225.9540\(< 0.0001\)
year20154.41020.72096.1173\(< 0.0001\)
year20164.63090.71996.4324\(< 0.0001\)
year20175.01410.72376.9289\(< 0.0001\)
year20185.33530.72257.3848\(< 0.0001\)
year20195.68060.71837.9088\(< 0.0001\)
year20205.90190.72238.1711\(< 0.0001\)
year20215.95780.72258.2465\(< 0.0001\)
year20225.92430.72228.2033\(< 0.0001\)
year20236.02940.71568.4255\(< 0.0001\)

B. smooth termsedfRef.dfF-valuep-value
s(thi_mean_t0_3d):paritymultiparous8.06368.06361133.7258\(< 0.0001\)
s(thi_mean_t0_3d):parityprimiparous8.57118.571144.5719\(< 0.0001\)
s(days_in_milk_t):paritymultiparous14.554714.554757761.6129\(< 0.0001\)
s(days_in_milk_t):parityprimiparous13.843613.84369080.5596\(< 0.0001\)
Table D.61: Jersey: Milk Yield - 1998-2023 - GAMM model summary without random effect terms.

Smooth Term Fixed EffectEst.SEzp
s(thi_mean_t0_3d):multiFx10.30760.08663.55\(< 0.0004\)
s(thi_mean_t0_3d):primiFx10.14740.14591.010.3124
s(days_in_milk_):multiFx14.91060.397812.34\(< 1e-34\)
s(days_in_milk_):primiFx14.12490.53237.75\(< 1e-14\)

Variance ComponentEstimated \( \sigma \)
\( \sigma_\alpha\)2.3991
\( \sigma_\iota\)0.9500
\( \sigma_\phi\)2.9764
s(thi_mean_t0_3d):multi1.5575
s(days_in_milk_):primi6.9995
s(days_in_milk_):multi8.4095
s(thi_mean_t0_3d):primi3.9762
Residual2.9355
Table D.62: Jersey: Milk Yield - 1998-2023 - Mixed Model Summary - Smooth Terms and Random Effects.

Model Diagnostics

Figure D.61: Jersey: Milk Yield - 1998-2023 - Diagnostic Plot
Figure D.61: Jersey: Milk Yield - 1998-2023 - Diagnostic Plot

THI Effect and Lactation Curve

Figure D.62: Jersey: Milk Yield - 1998 - 2023 - THI Effect and Lactation Curve

D.11.2: Split Period: Until 2010 - After 2010

D.11.2.1: Split Period: 1998 - 2010

Model Summary
A. parametric coefficientsEstimateStd. Errort-valuep-value
(Intercept)13.84640.615822.4846\(< 0.0001\)
parityprimiparous-1.59620.0195-81.7385\(< 0.0001\)
year1999-0.46700.6855-0.68120.4958
year20000.29870.67510.44240.6582
year20010.95510.65241.46390.1432
year20021.54170.65312.36040.0183
year20032.04550.64553.16880.0015
year20042.61150.63854.0902\(< 0.0001\)
year20052.86750.63314.5293\(< 0.0001\)
year20063.22720.62975.1248\(< 0.0001\)
year20073.12390.62914.9658\(< 0.0001\)
year20083.47960.62955.5276\(< 0.0001\)
year20093.88850.62876.1845\(< 0.0001\)
year20104.24740.62616.7842\(< 0.0001\)

B. smooth termsedfRef.dfF-valuep-value
s(thi_mean_t0_3d):paritymultiparous7.46877.4687367.5611\(< 0.0001\)
s(thi_mean_t0_3d):parityprimiparous5.92915.929131.4722\(< 0.0001\)
s(days_in_milk_t):paritymultiparous13.960013.960017839.8816\(< 0.0001\)
s(days_in_milk_t):parityprimiparous12.882212.88223856.3426\(< 0.0001\)
Table D.63: Jersey: Milk Yield - 1998-2010 - GAMM model summary without random effect terms.

Smooth Term Fixed EffectEst.SEzp
s(thi_mean_t0_3d):multiFx10.38330.14332.680.0075
s(thi_mean_t0_3d):primiFx10.08680.16310.530.5943
s(days_in_milk_):multiFx14.21880.67496.25\(< 1e-09\)
s(days_in_milk_):primiFx13.15640.74554.23\(< 1e-04\)

Variance ComponentEstimated \( \sigma \)
\( \sigma_\alpha\)0.9524
\( \sigma_\iota\)2.2355
\( \sigma_\phi\)2.6202
s(thi_mean_t0_3d):multi1.8172
s(days_in_milk_):primi7.1500
s(days_in_milk_):multi9.1468
s(thi_mean_t0_3d):primi1.3871
Residual2.6539
Table D.64: Jersey: Milk Yield - 1998-2010 - Mixed Model Summary - Smooth Terms and Random Effects.

Model Diagnostics
Figure D.63: Jersey: Milk Yield - 1998-2010 - Diagnostic Plot
Figure D.63: Jersey: Milk Yield - 1998-2010 - Diagnostic Plot
THI Effect and Lactation Curve
Figure D.64: Jersey: Milk Yield - 1998 - 2010 - THI Effect and Lactation Curve

D.11.2.2: Split Period: 2010 - 2023

Model Summary
A. parametric coefficientsEstimateStd. Errort-valuep-value
(Intercept)16.93740.1584106.9310\(< 0.0001\)
parityprimiparous-2.08850.0163-128.1973\(< 0.0001\)
year20120.20600.17961.14710.2513
year20130.28420.21081.34810.1776
year20140.87490.20544.2601\(< 0.0001\)
year20150.99390.20454.8606\(< 0.0001\)
year20161.22440.19326.3379\(< 0.0001\)
year20171.61260.19928.0952\(< 0.0001\)
year20181.93560.189510.2118\(< 0.0001\)
year20192.28700.184012.4266\(< 0.0001\)
year20202.50840.200112.5361\(< 0.0001\)
year20212.56360.199412.8587\(< 0.0001\)
year20222.52370.200812.5669\(< 0.0001\)
year20232.61250.179114.5866\(< 0.0001\)

B. smooth termsedfRef.dfF-valuep-value
s(thi_mean_t0_3d):paritymultiparous8.10138.1013856.5674\(< 0.0001\)
s(thi_mean_t0_3d):parityprimiparous5.90655.906553.5215\(< 0.0001\)
s(days_in_milk_t):paritymultiparous14.512914.512942118.2414\(< 0.0001\)
s(days_in_milk_t):parityprimiparous13.645213.64525835.7977\(< 0.0001\)
Table D.65: Jersey: Milk Yield - 2011-2023 - GAMM model summary without random effect terms.

Smooth Term Fixed EffectEst.SEzp
s(thi_mean_t0_3d):multiFx1-0.31120.1041-2.990.0028
s(thi_mean_t0_3d):primiFx1-0.16580.1121-1.480.1391
s(days_in_milk_):multiFx14.76960.452410.54\(< 1e-25\)
s(days_in_milk_):primiFx13.92230.60906.44\(< 1e-09\)

Variance ComponentEstimated \( \sigma \)
\( \sigma_\alpha\)2.3907
\( \sigma_\iota\)0.9756
\( \sigma_\phi\)3.1309
s(thi_mean_t0_3d):multi1.9309
s(days_in_milk_):primi7.4722
s(days_in_milk_):multi9.2408
s(thi_mean_t0_3d):primi1.0857
Residual2.9808
Table D.66: Jersey: Milk Yield - 2011-2023 - Mixed Model Summary - Smooth Terms and Random Effects.

Model Diagnostics
Figure D.65: Jersey: Milk Yield - 2011-2023 - Diagnostic Plot
Figure D.65: Jersey: Milk Yield - 2011-2023 - Diagnostic Plot
THI Effect and Lactation Curve
Figure D.66: Jersey: Milk Yield - 2010 - 2023 - THI Effect and Lactation Curve