Progressively simpler models are subsequently defined by removing terms from previous models beginning with interaction effects (e.g. 1 vs. 2). Nested models were tested evaluated using likelihood ratio tests to select the most parsimonious model. The reported p-values correspond to results of the likelihood ratio tests wherein a value greater than the nominal level of significance indicates that removal of the model term does not increase model deviance significantly. The final, parsimonious model that was selected is indicated in bold characters (Env = environment, EBV = EBV group). In the case of GHR, removal of the EBVxEnv term (model 5) resulted in a marginally non-significant likelihood ratio test. However, the interaction term was retained in the final model (model 4 AIC = 185.02), given that its exclusion resulted in a greater penalised likelihood score (model 5 AIC = 186.11).