Skip to content
Prev 12199 / 20628 Next

MCMCglmm predict function output and interpretation

Hi Justine,

If you get the predictions on the link scale, and denote these as  
eta_1 and  and eta_635 for the first observation, then

Pr(A) = 1/(1+exp(eta_1)+exp(eta_635))
Pr(B) = exp(eta_1)/(1+exp(eta_1)+exp(eta_635))
Pr(C) = exp(eta_635)/(1+exp(eta_1)+exp(eta_635))

There is some code for doing this (and marginalising any random  
effects) in the CourseNotes (p97 after Eq. 5.7).

Cheers,

Jarrod


Quoting Justine <jsmith5 at ucsc.edu> on Mon, 28 Jul 2014 23:45:09 +0000 (UTC):