Mas um duvida que fiquei é a seguinte:<br><br>#<br>
par(mfrow=c(4,4))<br><br>for(i in 1:4) {<br>  for(j in 1:4) {<br>    hist(rbeta(1000,i,j),main=<div id=":1vw">paste("shape1 =",i,"shape2 =",j),xlab="")<br>    função<-function(x) {dbeta(x,i,j)}<br>
    curve(função(x)*100,0,1,add=T,lwd=2,col="red")<br>
  }<br>}<br>#<br><br>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().<br>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().<br>
O autor explica na parte 2 do artigo mas eu não entendi.<br><br></div><br><br><div class="gmail_quote">Em 19 de outubro de 2012 16:06, Wagner Bonat <span dir="ltr"><<a href="mailto:wbonat@gmail.com" target="_blank">wbonat@gmail.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div>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 ?</div>
<span class="HOEnZb"><font color="#888888">

<div> </div><br clear="all"><br>-- <br>Wagner Hugo Bonat<br>LEG - Laboratório de Estatística e Geoinformação<br>UFPR - Universidade Federal do Paraná<br>
</font></span><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><br clear="all"><br>-- <br><div>Grato<br>
Augusto C. A. Ribas</div>
<div> </div>
<div>Site Pessoal: <a href="http://augustoribas.heliohost.org" target="_blank">http://augustoribas.heliohost.org</a></div>
<div>Lattes: <a href="http://lattes.cnpq.br/7355685961127056" target="_blank">http://lattes.cnpq.br/7355685961127056</a><br></div><br>