Agu*_*cho 5 r curve-fitting lme4
我试图通过拟合非线性混合模型来绘制合成曲线.它应该像正常分布的曲线,但向右倾斜.我在这里和这里都遵循了以前的链接,但是当我使用我的数据时,我无法在不同的困难中实现它(见下文).
这是数据集 和代码
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...:~/
Run Code Online (Sandbox Code Playgroud)
之后,我想在原始数据上绘制曲线
with(vmax,plot(vmax ~ (TEM)))
x=vmax$TEM
lines(x, predict(nm1, newdata = data.frame(TEM = x, INDIVIDUO = "ACI5")))
Run Code Online (Sandbox Code Playgroud)
任何提示?
提前致谢
| 归档时间: |
|
| 查看次数: |
326 次 |
| 最近记录: |