jk is the expected hoof temperature, σ 2 e is the error variance, the b (k) i s are parameters and X ijk is the value of the ith factor (e.g. hoof, aspect or ambient temperature) for the jth observation on animal k. Betweenanimal variation was modelled by assuming that the parameters for each animal are drawn from higherorder distributions, such that,