Displaying 2 results from an estimated 2 matches for "revordershrunk".
2001 Sep 07
3
fitting models with gnls
...mmands, got an error message:
> growth.gnls <- gnls(lt ~ Linf*(1-exp(-K*(age-t0))),
+ data=growth.dat, params= Linf +K + t0 ~ 1, start=list(Linf=500,K=0.2,t0=0),
+ control = list(returnObject = T), corr = corAR1(form=~fish|age))
Error in "[<-.factor"(*tmp*, , value = grpShrunk[revOrderShrunk]) :
Argument "i" is missing, with no default
If I use nls function I can estimate the parameters:
> growth.nls <- nls(lt ~ Linf*(1-exp(-K*(age-t0))),
+ data=growth.dat, start=list(Linf=500,K=0.2,t0=0))
> summary(growth.nls)
Formula: lt ~ Linf * (1 - exp(-K * (age - t...
2001 Jun 01
1
nls works but not gnls
This works fine:
fit42<-nls(Vfs~SSlogis(Months,Asym.Int+Asym.Group*Groupdum,xmid,scal),
data=df,
start=c(Asym.Int=22,Asym.Group=5,xmid=2,scal=6),
na.action=na.omit)
But this, identical except using gnls, doesn't converge:
fit43<-gnls(Vfs~SSlogis(Months,Asym.Int+Asym.Group*Groupdum,xmid,scal),
data=df,
start=c(Asym.Int=22,Asym.Group=5,xmid=2,scal=6),
na.action=na.omit)
Error in gnls(Vfs