Prezados,
Utilizei a função "nls" do R para estimar as constantes de uma função potencial ajustada aos meus dados de comprimento e peso de lulas utilizando o método dos mínimos quadrados:
m <- nls(y~a*x^b, data = ds, start = c(a=.1, b=2))
Obtive as seguintes constantes: a = 0.0003967 , b = 2.3933797 e o R2 = 0.9681
Ocorre que o mesmo ajuste no Excel 2007 resulta em a = 0.0002 , b = 2.6097 e R2 = 0.9819
O mesmo ajuste no Statistica resulta em constantes e R2 idênticas aos encontrados no R.
Não deveria ter dado o mesmo resultado nos três programas?
Não sei em qual resultado confiar.
Também gostaria de saber como faço para plotar a curva dessa função potencial sobreposta a minha dispersão de dados (peso vs. comprimeto).
Agradeço desde logo.
Rodrigo
--
=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->
Rodrigo Silvestre Martins, PhD
Bolsista FAPESP Pós-Doutorado
Laboratório de Ecossistemas Pesqueiros (LabPesq)
Universidade de São Paulo, Instituto Oceanográfico
Praça do Oceanográfico, 191. Cidade Universitaria (sala 107-A/B)
Butantã - São Paulo/SP, Brasil
05508-900
Tel: +55 11 3091 6549
Email:
rodrigo.plei@gmail.com ;
ocersm@lycos.com;
rsmartins@usp.br
CV Lattes:
http://lattes.cnpq.br/5350064124902777