<div dir="ltr">Qual(is) resultado(s) você observa e por que você os considera errados?</div><div class="gmail_extra"><br><div class="gmail_quote">2015-10-31 12:17 GMT-02:00 ASANTOS <span dir="ltr"><<a href="mailto:alexandresantosbr@yahoo.com.br" target="_blank">alexandresantosbr@yahoo.com.br</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Prezados Listeiros,<br>
<br>
Consegui achar o problema, porém os resultados da comparação entre a distribuição observada e estimada estão errados, alguem poderia me dar uma luz? Segue novo CRM:<span class=""><br>
<br>
## Teste de distribuição de erros tabajara<br>
<br>
require(MASS)<br>
<br>
BD<-rpois(1000,10) ### Valores a serem testados<br>
<br>
Dist<-c("negative binomial","normal","Poisson")<br>
<br></span>
for(typeD in seq_along(Dist)){<br>
<br>
k <- fitdistr(BD,Dist[typeD])<br>
par <- k$estimate<br>
size <- par[1]#k<span class=""><br>
mu <- par[2]#Média<br>
SD<-sd(BD)<br>
N <- length(BD)<br>
<br></span>
if (Dist[typeD]=='negative binomial'){ ### Para binomial negativa<span class=""><br>
est <-N*dnbinom(BD,size=size,mu=mu) ## Estimativas<br>
fecdf <- ecdf(BD) ###ecdf- Função de distribuição cumulativa empírica<br>
knotsX <- knots(fecdf)<br>
emp <- fecdf(c(knotsX,Inf)) # Empírica<br></span>
print(Dist[typeD])<br>
print(chisq.test(table(emp),table(est),correct=TRUE)) ##Teste Chi - Quadrado para valiar a significância da distribuição<br>
}<br>
if (Dist[typeD]=='normal'){ ### Para normal<br>
est <-N*dnorm(BD,mean=size, sd=SD) ## Estimativas<span class=""><br>
fecdf <- ecdf(BD) ###ecdf- Função de distribuição cumulativa empírica<br>
knotsX <- knots(fecdf)<br>
emp <- fecdf(c(knotsX,Inf)) # Empírica<br></span>
print(Dist[typeD])<br>
print(chisq.test(table(emp),table(est),correct=TRUE)) ##Teste Chi - Quadrado para valiar a significância da distribuição<br>
}<br>
if (Dist[typeD]=='Poisson'){ ### Para Poisson<br>
est <-N*dpois(BD,lambda=size) ## Estimativas<span class=""><br>
fecdf <- ecdf(BD) ###ecdf- Função de distribuição cumulativa empírica<br>
knotsX <- knots(fecdf)<br>
emp <- fecdf(c(knotsX,Inf)) # Empírica<br></span>
print(Dist[typeD])<br>
print(chisq.test(table(emp),table(est),correct=TRUE))##Teste Chi - Quadrado para valiar a significância da distribuição<br>
}<br>
}<br>
#<br>
<br>
Obrigado<br><span class="im HOEnZb">
<br>
-- <br>
======================================================================<br>
Alexandre dos Santos<br>
Proteção Florestal<br>
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso<br>
Campus Cáceres<br>
Caixa Postal 244<br>
Avenida dos Ramires, s/n<br>
Bairro: Distrito Industrial<br>
Cáceres - MT CEP: 78.200-000<br>
Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO)<br>
<a href="mailto:e-mails%3Aalexandresantosbr@yahoo.com.br" target="_blank">e-mails:alexandresantosbr@yahoo.com.br</a><br>
<a href="mailto:alexandre.santos@cas.ifmt.edu.br" target="_blank">alexandre.santos@cas.ifmt.edu.br</a><br>
Lattes: <a href="http://lattes.cnpq.br/1360403201088680" rel="noreferrer" target="_blank">http://lattes.cnpq.br/1360403201088680</a><br>
======================================================================<br>
<br></span><div class="HOEnZb"><div class="h5">
Em 28/10/2015 15:25, ASANTOS escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Prezados listeiros,<br>
<br>
Estou tentando criar uma função para testar a distribuição de erros de normal, poisson e binomial negativa por enquanto, porém não sei o que estou fazendo de errado na minha função e solicito a ajuda de vocês para conseguir achar o erro, segue CRM:<br>
<br>
## Teste de distribuição de erros tabajara<br>
<br>
require(MASS)<br>
<br>
BD<-rpois(1000,10) ### Valores a serem testados<br>
<br>
Dist<-c("negative binomial","normal","Poisson")<br>
<br>
for(typeD in 1:Dist){<br>
k <- fitdistr(BD,typeD[Dist])<br>
par <- k$estimate<br>
size <- par[1]#Parametro<br>
mu <- par[2]#Média<br>
SD<-sd(BD)<br>
N <- length(BD)<br>
<br>
if (typeD[Dist]=='negative binomial'){ ### Para binomial negativa<br>
est <-N*dnbinom(BD,size=size,mu=mu) ## Estimativas<br>
fecdf <- ecdf(BD) ###ecdf- Função de distribuição cumulativa empírica<br>
knotsX <- knots(fecdf)<br>
emp <- fecdf(c(knotsX,Inf)) # Empírica<br>
chisq.test(table(emp),table(est),correct=TRUE) ##Teste Chi - Quadrado para valiar a significância da distribuição<br>
}<br>
if (typeD[Dist]=='normal'){ ### Para normal<br>
est <-N*dnorm(BD,mean=mu, sd=SD) ## Estimativas<br>
fecdf <- ecdf(BD) ###ecdf- Função de distribuição cumulativa empírica<br>
knotsX <- knots(fecdf)<br>
emp <- fecdf(c(knotsX,Inf)) # Empírica<br>
chisq.test(table(emp),table(est),correct=TRUE) ##Teste Chi - Quadrado para valiar a significância da distribuição<br>
}<br>
if (typeD[Dist]=='Poisson'){ ### Para Poisson<br>
est <-N*dpois(BD,lambda=mu) ## Estimativas<br>
fecdf <- ecdf(BD) ###ecdf- Função de distribuição cumulativa empírica<br>
knotsX <- knots(fecdf)<br>
emp <- fecdf(c(knotsX,Inf)) # Empírica<br>
chisq.test(table(emp),table(est),correct=TRUE) ##Teste Chi - Quadrado para valiar a significância da distribuição<br>
}<br>
}<br>
#<br>
<br>
Obrigado,<br>
<br>
</blockquote>
<br></div></div><div class="HOEnZb"><div class="h5">
_______________________________________________<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" rel="noreferrer" 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" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e fornea cdigo mnimo reproduzvel.</div></div></blockquote></div><br></div>