Dear list I have a conceptual question. I run a MCMCglmm model (model 1) and extracted the posterior distribution of the random effects. This is a matrix of 75 (number of individuals) x 5000 (number of iteractions of the model). Now I'd like to use that posterior distribution (75x5000) as a explanatory variable X (fixed effect) in another mixed-effect model (model 2). I have thought of two options: - Run model 2, 5000 times, one with each of the values of the distribution of variable X, and then get a distribution of responses - Run model two, e.g. 1000 times randomly picking one value for the distribution of variable X on each run, and then get a distribution of responses Could you please give me some advise on how to proceed, and if my thoughts are correct? I'd like to know how to do it both with MCMCglmm (Bayesian) and lme or lme4 (frequentist). Thanks! David
use a posterior distribution from one model as fixed effect in another model
1 message · David Villegas Ríos