Muito obrigado Paulo, fico feliz com a sua ajuda.<br>Abraços.<br><br><div class="gmail_quote">Em 21 de novembro de 2012 14:35, Paulo Justiniano <span dir="ltr"><<a href="mailto:paulojus@leg.ufpr.br" target="_blank">paulojus@leg.ufpr.br</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Para comparação viisual voce poode adaptar do exemplo a seguir:<br>
<br>
hist(x, prob=T)<br>
lines(density(x))<br>
rug(x)<br>
fit.n <- fitdistr(x,"normal")<br>
with(fit.n, curve(dnorm(x, m=estimate["mean"], sd=estimate["sd"]), from=150, to=250, add=T, col=4))<br>
fit.t <- fitdistr(x,"t")<br>
mydt <- function(x, m, s, df) dt((x-m)/s, df)/s<br>
with(fit.t, curve(mydt(x, m=estimate["m"], s=estimate["s"], df=estimate["df"]), from=150, to=250, col=2, add=T))<br>
<br>
<br>
Para comparacao para indicar melhor ajuste use a varossimilahnca retornada pelos ajustes (penalizando modelos com mais parametros, por ex. via AIC)<div class="HOEnZb"><div class="h5"><br>
<br>
<br>
<br>
<br>
On Wed, 21 Nov 2012, Matheus Monteiro wrote:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
obrigado Paulo, qual comando voce acha melhor para fazer tal comparação?<br>
<br>
Em 21 de novembro de 2012 13:20, Paulo Justiniano <<a href="mailto:paulojus@leg.ufpr.br" target="_blank">paulojus@leg.ufpr.br</a>> escreveu:<br>
fora dos parentesis estao as estimativas e dentro os erros padrao.<br>
<br>
Portanto voce deve pegar as estimativas e usar para tracar curvas de cada distribuição<br>
<br>
<br>
On Wed, 21 Nov 2012, Matheus Monteiro wrote:<br>
<br>
Obtive as densidades aproximadas da 't' e da 'normal', gostaria de compara-las através de graficos dessas densidades<br>
estimadas.<br>
Qual seria o melhor comando para isso? Uso o que está entre parênteses ou o outro? Obrigado.<br>
<br>
t <br>
<br>
> fitdistr(x[,1],"t")<br>
<br>
m s df <br>
<br>
-0.001874414 0.009122738 3.522064232 <br>
<br>
( 0.001726835) ( 0.001862692) ( 2.234067399)<br>
<br>
normal<br>
<br>
fitdistr(x[,1],"normal")<br>
<br>
mean sd<br>
<br>
0.0001040816 0.0128201843<br>
<br>
(0.0018314549) (0.0012950342)<br>
<br>
<br>
x são dados simulados. <br>
<br>
<br>
<br>
______________________________<u></u>_________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/<u></u>cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-<u></u>guia</a>) e forneça código mínimo reproduzível.<br>
<br>
<br>
<br>
</blockquote>
</div></div><br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br>