Quero comparar modelos para duas espécies de peixes. Consigo fazer usando nls mas
não gnls. Alguém tem alguma idéia de como fazer?
dados<-structure(list(Peixes = structure(c(1L, 1L,
1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L,
2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L),
.Label = c("1",
"2"), class = "factor"), Tempo =
c(0L, 0L, 0L, 30L, 30L, 30L,
60L, 60L, 60L, 90L, 90L, 90L, 240L, 240L, 240L, 270L,
270L, 270L,
0L, 0L, 0L, 30L, 30L, 30L, 60L, 60L, 60L, 90L, 90L,
90L, 240L,
240L, 240L,
270L, 270L, 270L), Peso = c(6.09, 6.09, 6.09, 65.73,
52, 53.13,
102.88, 115.8, 129.31, 161.06, 174.94, 213.63, 233.67,
291.93,
301.33, 258.47, 326.84, 337.04, 8.92, 8.92, 8.92, 78.8,
72.07, 85.8,
167.82, 146.13, 179.31, 255.75, 230.73, 207.06,
320.67, 280.2,
243.93, 376.55, 441.96, 349.84)), .Names = c("Peixes",
"Tempo",
"Peso"), row.names = c(NA, -36L), class = "data.frame")
library(nlme)
m1 <-gnls(Peso~a*(1-b*(exp(-c*Tempo)))^3 , weights=varPower(), data=dados, start=c(a=350,b=1, c=0.01))
m2<- nls(Peso~a*(1-b*(exp(-c*Tempo)))^3,
data=dados, start=c(a=300,b=10,
c=0.05))
m3<-gnls(Peso~a[Peixes]*(1-b[Peixes]*(exp(-c[Peixes]*Tempo)))^3,
, weights=varPower(), data=dados, start=list(a=c(300,300),b=c(10,10),
c=c(0.05,0.05)))
m4<-nls(Peso~a[Peixes]*(1-b[Peixes]*(exp(-c[Peixes]*Tempo)))^3,
data=dados, start=list(a=c(300,300),b=c(10,10),
c=c(0.05,0.05)))
anova(m1,m3)
anova(m2,m4)