Dear list,

I am fitting a model with lmer, using the following code:

lmer(Y ~ b_0 + b_1*X1 + b_2*X2 + b_3*X1:X2 + (1  | groups) + (0 + X1  |
groups), data=mydata)

Where X2 is measured at the 2nd level (named groups in the code above).

Then, I calculate the marginal fixed effects of X1:  b_1 + b_3*X2. But I
would also like to calculate the HPD interval and the pMCMC for such
marginal effect, similar to what the pvals.fnc function does for all the
parameters of the lmer model. What would be the best way to do so, once
pvals.fnc, mcmcsamp and HPDintervals only work with complete fitted models,
not with particular singular values?

Thanks for your time.

Tom

	[[alternative HTML version deleted]]

