I suppose I should've done some more experimenting before I posted.
Removing the "-1" from the random argument appears to specify the model as
desired.

Thanks anyway,

Chris




-----------------------------
Chris Gast
cmgast@gmail.com


On Sun, May 20, 2012 at 2:12 PM, Chris Gast <cmgast@gmail.com> wrote:

> Hello,
>
> I would like to specify the following model (in lmer() syntax) in
> MCMCglmm, but seem to be unable to do so. I receive the warning message
> "some fixed effects are not estimable and have been removed", although I do
> not believe this should be the case.
>
> Here's the lmer() model:
>
> fit = lmer(y ~ -1 + cyear:stratum + site + (-1+cyear|site),data=indat)
>
> And here's an attempt at fitting this in MCMCglmm():
>
> fitmcmc = MCMCglmm(y ~ -1 + cyear:stratum +
> site,random=~idh(-1+cyear):site,data=indat)
>
> The data structure is relatively simple: stratum and site are categorical
> (technically site is nested within stratum, so I'd give bonus points if
> someone can help me specify separate variance components for the
> site-specific slopes within each stratum), and "cyear" is the continuous
> predictor for which I would like stratum-specific slopes with random slopes
> for each site (again, within stratum, if possible).
>
>
> Thanks,
>
> Chris Gast
>
>
>
>
>
> -----------------------------
> Chris Gast
> cmgast@gmail.com
>

	[[alternative HTML version deleted]]

