在原始拟合数据上绘制非线性混合模型

时间:2016-10-16 23:32:28

标签: r curve-fitting lme4

我试图通过拟合非线性混合模型来绘制合成曲线。它应该像正常分布的曲线,但向右倾斜。我在这里和这里都遵循了以前的链接,但是当我使用我的数据时,我不能在不同的困难中实现它(见下文)。

这是dataset 和代码

s=read.csv("GRVMAX tadpoles.csv") 
t=s[s$SPP== levels(s$SPP)[1],]
head(t)
vmax=t[t$PERFOR=="VMAX",]
colnames(vmax)[6]="vmax"
vmax$TEM=as.numeric(as.character(vmax$TEM));
require(lme4)
start =c(TEM=25)
is.numeric(start)
nm1 <- nlmer ( vmax ~ deriv(TEM)~TEM|INDIVIDUO,nlpars=start, nAGQ  =0,data= vmax)# this gives an error suggesting nlpars is not numeric, despite start is numeric...:~/

之后,我想在原始数据上绘制曲线

 with(vmax,plot(vmax ~ (TEM)))
 x=vmax$TEM
 lines(x, predict(nm1, newdata = data.frame(TEM = x, INDIVIDUO = "ACI5")))

任何提示?

提前致谢

0 个答案:

没有答案