search for: modnl

...ining the greatest amount of variation appears to be da. I have used a log(da) to linearize the relationship, but an asymptotic relationship is more appropriate. I have used nls: B<-c(.007,1,3,2,2,2,2,2,2,2,2,2,.05,.001,.8) st<-list(ad=B[1],ahabdiv=B[2],abas=B[3:12],b=B[13],c=B[14],z=B[15]) modnls.a<-nls(sr~ad*log(d)+ahabdiv*habdiv+abas[bas]+(b/(c+(da^-z))), start=st,trace=T) I next used a random slope and intercept model using lmer from the package (lme4). modlme<-lmer(y~d+habdiv+log(da)+(log(da)|bas),method='ML') What I would like to do is use a similar model to the...