[R-br] Comparar modelos com gnls

[Jonas Brito]

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)
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120329/53465e0d/attachment.html>