[R-br] Conflito de resultados + plotagem

[Fernando Mayer]

Rodrigo,

o excel é um programa de tabulação e visualização de dados, não um
pacote estatístico. Eu não confiaria nos resultados de um software que
não tem o propósito de análise estatística. Existe uma série de
problemas com resultados numéricos no excel, sugiro que vc de uma
olhada em

http://www.burns-stat.com/pages/Tutor/spreadsheet_addiction.html

principalmente na seção "Specific Problems with Excel". Além de tudo
isso, ele ainda é um programa de código fechado, ou seja, mesmo que vc
quisesse, não poderia olhar os cálculos que ele está fazendo - o que é
o oposto do R.

Para plotar a curva dos valores preditos sobreposta aos seus dados,
use a função predict() e veja o exemplo que tem em

?predict.nls

Usando os seus objetos, ficaria algo como:

plot(y ~ x, data = ds)
xx <- seq(min(ds$x), max(ds$x), length = 100) # quantos pontos forem necessarios
lines(xx, predict(m, newdata = list(x = xx, y = ds$y)))


[]s,

---
Fernando Mayer
URL: http://sites.google.com/site/fernandomayer
e-mail: fernandomayer [@] gmail.com



2011/5/19 Rodrigo Plei <rodrigo.plei em gmail.com>:
> 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
>
>
>
>
>
>
>
>
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>
>