[R-sig-ME] Question about the predict function in MCMCglmm

[Jarrod Hadfield]

Hi,

Also

tcrossprod(beta_hat, X) + U

should be

(tcrossprod(beta_hat, X) + U)/sqrt(1 + k^2 * rowSums(mMCMC$VCV[,"units"]))

for pred2_hand

Cheers,

Jarrod



On 02/11/2016 16:07, Jarrod Hadfield wrote:
> tcrossprod(beta_hat, X) + U

-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-mixed-models/attachments/20161102/05af1b92/attachment-0001.pl>