Displaying 1 result from an estimated 1 matches for "dcum".
Did you mean:
dcu
2005 Nov 08
1
Can someone Help in nls() package
...248.2,266,
277,291.4,392,461.2,470,418,410.8,395.6,365.2,189.4,43.4,33.2,32,29,26,26,25.6,24.6,
24,23.4,18.4,17.6,15.8,13.6,10.6,8.4,6.4,5.4,4.2,3.6,2.7))
times<-c(0:60)
f <- function(t,xx,pars){
Cur<- xx[1]
Cum<- xx[2]
a <- pars[1]
b <- pars[2]
dCur <- a*Cur-b*Cur*Cum^2
dCum <- Cur
list(c(dCur,dCum))
}
require(nls)
fit <- nls(Cur ~ lsoda(c(Cur=1,Cum=1), times, f, c(a=a, b=b))[,2],
data=dat,
start=list(a=2.5, b=.01),
trace=T
)
print(coef(fit))
It runs the first iterations and shows the output, after the intial
evelua...