Prezados,<div><br></div><div>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:</div>
<div><br></div><div>m <- nls(y~a*x^b, data = ds, start = c(a=.1, b=2))</div><div><br></div><div>Obtive as seguintes constantes: a = 0.0003967 , b = 2.3933797 e o R2 = 0.9681</div><div><br></div><div>Ocorre que o mesmo ajuste no Excel 2007 resulta em a = 0.0002 , b = 2.6097 e R2 = 0.9819</div>
<div><br></div><div>O mesmo ajuste no Statistica resulta em constantes e R2 idênticas aos encontrados no R.</div><div><br></div><div>Não deveria ter dado o mesmo resultado nos três programas?</div><div><br></div><div>Não sei em qual resultado confiar.</div>
<div><br></div><div>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).</div><div><br></div><div>Agradeço desde logo.</div><div><br></div>
<div>Rodrigo</div><div><br>-- <br>=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8--><br><br>Rodrigo Silvestre Martins, PhD<br>Bolsista FAPESP Pós-Doutorado<br><br>Laboratório de Ecossistemas Pesqueiros (LabPesq)<br>
Universidade de São Paulo, Instituto Oceanográfico<br>Praça do Oceanográfico, 191. Cidade Universitaria (sala 107-A/B)<br>Butantã - São Paulo/SP, Brasil<br>05508-900<br>Tel: +55 11 3091 6549<br>Email: <a href="mailto:rodrigo.plei@gmail.com" target="_blank">rodrigo.plei@gmail.com</a> ; <a href="mailto:ocersm@lycos.com" target="_blank">ocersm@lycos.com</a>; <a href="mailto:rsmartins@usp.br" target="_blank">rsmartins@usp.br</a><br>
<br>CV Lattes: <a href="http://lattes.cnpq.br/5350064124902777" target="_blank">http://lattes.cnpq.br/5350064124902777</a><br><br><br><br><br><br><br><br><br><br>
</div>