data <- read.table("c:\\temp\\chicks.txt",header=T)
attach(data)
names(data)
[1] "weight" "Time"
model <- nls(weight~SSfpl(Time, a, b, c, d))
xv <- seq(0,22,.2)
yv <- predict(model,list(Time=xv))
plot(weight~Time,pch=21,col="red",bg="green4")
lines(xv,yv,col="navy")
我尝试将自启动功能适合nls模型。但是我不知道如何创建公式和拟合模型。谢谢。
注意:上面的代码来自Michael J. Crawley的R书。