[R-br] Comparar modelos com gnls

[Walmes Zeviani]

A gnls() não entende a[Peixes], b[Peixes], enfim, separar parâmetros por
variáveis indicadoras não é como na nls(). Você deve fazer isso por meio do
argumento params= da gnls(). No seu CMR, o m3 ficaria assim

m3 <- gnls(Peso~a*(1-b*(exp(-c*Tempo)))^3,
           params=list(a+b+c~Peixes),
           weights=varPower(), data=dados,
           start=list(a=c(300,0), b=c(10,0), c=c(0.05,0)))

Veja a documentação da função para maiores informações.

À disposição.
Walmes.

==========================================================================
Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
Departamento de Estatística - Universidade Federal do Paraná
fone: (+55) 41 3361 3573
VoIP: (3361 3600) 1053 1173
e-mail: walmes em ufpr.br
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120329/6ece520b/attachment.html>