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)