Displaying 1 result from an estimated 1 matches for "growthfit".
2008 Jul 31
1
nls weights warning message
...9,170,163,197)
x<-c(1,1,1,1,1,2,2,3,3,3,3,3,3,4,4,4,5,5,5,5,6,6,6,6,6,6,7,7,7,7,7,7,7,7,7,8,8,8,
8,8,8,8,8,8,8,9,9,9,9,9,11,12)
Data<-data.frame(cbind(y,x))
n2<-length(x)
weight<-c(1,rep(0.1,n2-1))
Ta<-min(Data$x)
Tb<-max(Data$x)
G.st<-c(k=0.005, g1=50,g2=550)
#no weights
growthFit<-nls(y~((g1)*exp((log(g2/g1))*(1-exp(-k*(x-Ta)))
/(1-exp(-k*(Tb-Ta))))), data=Data, start=G.st, trace=T)
#with weights
growthFit.wts<-nls(y~((g1)*exp((log(g2/g1))*(1-exp(-k*(x-Ta)))
/(1-exp(-k*(Tb-Ta))))), data=Data, start=G.st, trace=T, weights=weight)
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~...