[R-br] dúvidas sobre utilização de disitribuição beta

[Augusto Ribas]

Mas um duvida que fiquei é a seguinte:

#
par(mfrow=c(4,4))

for(i in 1:4) {
  for(j in 1:4) {
    hist(rbeta(1000,i,j),main=
paste("shape1 =",i,"shape2 =",j),xlab="")
    função<-function(x) {dbeta(x,i,j)}
    curve(função(x)*100,0,1,add=T,lwd=2,col="red")
  }
}
#

Olhando o grafico acima, da pra comparar com como esta a distribuição dos
dados e pensar num chute bom para o start, que da pra inserir pelo
betareg.control().
Mas eu tenho dificuldade com matematica e não consegui fazer as estimativas
de média e phi que ele estima no ajuste3 voltar a ser o valores que usei no
shape1 e shape2 pra gerar os dados em rbeta().
O autor explica na parte 2 do artigo mas eu não entendi.



Em 19 de outubro de 2012 16:06, Wagner Bonat <wbonat em gmail.com> escreveu:

> Existe uma correção de borda para isto que funciona bem. O seu problema
> não parece ser com relação diretamente a isso, a mensagem diz respeito ao
> processo de maximização numerica, o betareg usa um esquema de busca de
> valores iniciais automatica quando esse procedimento não funciona ele solta
> a msg. Faz um ajuste preliminar e usa as estimativas deste ajuste novamente
> como chute inicial e verifique se os resultados estão sensíveis a isso. Vc
> pode disponibilizar seus dados ?
>
>
>
> --
> Wagner Hugo Bonat
> LEG - Laboratório de Estatística e Geoinformação
> UFPR - Universidade Federal do Paraná
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>



-- 
Grato
Augusto C. A. Ribas

Site Pessoal: http://augustoribas.heliohost.org
Lattes: http://lattes.cnpq.br/7355685961127056
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121019/d13b272a/attachment.html>