You may try

model4 <- lme(logSSP ~ logM + K, random = list(species =
pdBlocked(list(pdSymm(~logM), pdSymm(~K)))), data=data1)

Alternatively another approach is

require(lme4)
model5 <- lmer(logSSP ~ logM + K + (logM|species) + (K|species), data=data1)

Since you didn't provide any data, I'm not so sure how these models pan
out. So let us know how it goes.

Cheers,
Gang


On Wed, Mar 14, 2012 at 4:41 PM, Bhogal, Harkiran <
harkiran.bhogal07@imperial.ac.uk> wrote:

> Please help!
>
>
>
> Got a few days left and I need to model a random effect of species on the
> body mass (logM) and temperature (K) slopes. This is what i've done so far
> that works:
>    model1<-lme(logSSP~logM + K,random=~1|species,data=data1)
>    model2<-lme(logSSP~logM + K,random=~K|species,data=data1)
>    model3<-lme(logSSP~logM + K,random=~logM|species,data=data1)
>
>
> The one I now want is:
>    model4<-lme(logSSP~logM + K,random=~logM|species,K|species,data=data1)
>    #I need the random effect of spp on both slopes of logM and K, but this
> code doesn't work so how do i change the code??????
>    :( Any help will be greatly appreciated
>
>        [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-mixed-models@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>

	[[alternative HTML version deleted]]


