[R-sig-ME] lme code

[Federico Calboli]

On 14 Mar 2012, at 20:41, Bhogal, Harkiran 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)

random ~ logM + species|K/species

????

please explain what you are trying to do in the random part, you code for model 4 does not just not make sense as R code, it just does not make sense period.

F


>    #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 at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models

--
Federico C. F. Calboli
Neuroepidemiology and Ageing Research
Imperial College, St. Mary's Campus
Norfolk Place, London W2 1PG

Tel +44 (0)20 75941602   Fax +44 (0)20 75943193

f.calboli [.a.t] imperial.ac.uk
f.calboli [.a.t] gmail.com