[R-br] Conflito de resultados + plotagem

[Rodrigo Plei]

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 em gmail.com ; ocersm em lycos.com; rsmartins em usp.br

CV Lattes: http://lattes.cnpq.br/5350064124902777
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110519/8bfac450/attachment.html>