Agradeço a todos os comentários e ao Augusto por ajudar a explicar melhor todo o problema! <br><br>Wagner, também sou bastante limitado em matemática por isso do tópico. <br><br>Aqui estão os dados que estamos testando e o modelo. Ajuda sobre a melhor forma de análisa-los é muito bem vinda!<br>
Grato <br>Nicolay<br><br><br><div id=":4td">_____________________________________________________________________________<br><br><br>#dados<br><br>index<-c(0.749962502, 0.629968502, 0.990454495, 0.97559371, 0.959552189, <br>
0.992212907, 0.991207709, 0.977094002, 0.960647357, 0.965756886, <br>
0.710047142, 0.887653978, 0, 0.840074171, 0.991626994, 0.926527555, <br>0.855316235, 0.93723019, 0.779908921, 0.971128419, 0.84293306, <br>0.69501627, 0.927165576, 0.685086853, 0.97300735, 0, 0.985987221, <br>0, 0.856096633, 0, 0.749962502, 0.830983382, 0.868087428, 0, <br>
0, 0.692213048, 0.746629335, 0.984632505, 0.969449446, 0.990320893, <br>0, 0, 0.793822786, 0.925777711, 0.832891421, 0.835937039, 0.083261143, <br>0.994466643, 0.990346881, 0.995345641, 0.995227688, 0.942101655, <br>0.987427629, 0.990576035, 0.926467349, 0.897074488, 0.534408541, <br>
0.638856946, 0.996976761, 0.901154942, 0.464509646, 0.822662948, <br>0.963468058, 0.615255249, 0.884813406, 0.907454627, 0.950820473, <br>0.44401526, 0.952543568, 0.749516361, 0.691836864, 0.944174245, <br>0.983889943)<br>
<br><br>periodicity<-structure(c(2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, <br>2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, <br>1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, <br>
2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, <br>3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 1L, 1L, 2L), .Label = c("ephemerous", <br>"permanent", "permanent_instable"), class = "factor")<br>
<br>size<-c(1000, 200, 200, 2000, 100, 2000, 2000, 200, 200, 500, 100, <br>200, 25, 200, 50, 200, 100, 100, 100, 500, 100, 500, 100, 500, <br>500, 25, 200, 25, 100, 100, 25, 25, 200, 25, 25, 100, 200, 200, <br>500, 1000, 25, 50, 500, 500, 500, 100, 50, 500, 1000, 2000, 200, <br>
2000, 2000, 500, 100, 50, 2000, 2000, 2000, 2000, 2000, 2000, <br>2000, 2000, 2000, 2000, 2000, 2000, 500, 50, 50, 50, 200)<br><br>#modelo<br><br>library(betareg)<br>index[index=='1']<-0.9999<br>index[findex=='0']<-0.0001<br>
<br>evenn_logit<-betareg(index~log(size)+periodicity)<br>coef(evenn_logit)<br>summary(evenn_logit)</div><br><br><div class="gmail_quote">2012/10/19 Augusto Ribas <span dir="ltr"><<a href="mailto:ribas.aca@gmail.com" target="_blank">ribas.aca@gmail.com</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Mas um duvida que fiquei é a seguinte:<div class="im"><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><div><div class="im">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></div>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><div class="h5"><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><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></div></div><div class="im">_______________________________________________<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/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></div></blockquote></div><br><br clear="all"><div class="HOEnZb">
<div class="h5"><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>
</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><br clear="all"><br>-- <br>Nicolay Leme da Cunha<br>
<br>Biólogo, Mestre, Doutorando em Ecologia e Conservação<br>Universidade Federal de Mato Grosso do Sul, 79070-900<br>Campo Grande, MS, Brasil<br>E-mail: <a href="mailto:nicolaycunha@gmail.com" target="_blank">nicolaycunha@gmail.com</a><div>
<a href="http://lattes.cnpq.br/5916316648872099" target="_blank">lattes.cnpq.br/5916316648872099</a> </div><br>