<div dir="ltr"><span style="font-size:12.8px">Jaildo Rocha,</span><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">Muito obrigado por sua ajuda.</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">Alisson</div></div><div class="gmail_extra"><br><div class="gmail_quote">2016-09-20 20:53 GMT-03:00 Alisson Lucrécio <span dir="ltr"><<a href="mailto:alisson.lucrecio@ifgoiano.edu.br" target="_blank">alisson.lucrecio@ifgoiano.edu.br</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Jaildo Rocha,<div><br></div><div>Muito obrigado por sua ajuda.</div><div><br></div><div>Alisson</div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="h5">2016-09-20 19:27 GMT-03:00 Jaildo Rocha <span dir="ltr"><<a href="mailto:j.vieirarocha@gmail.com" target="_blank">j.vieirarocha@gmail.com</a>></span>:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><div><div class="h5"><span><font face="Calibri" style="font-size:14px">Prezado Alisson,</font><div><font face="Calibri" style="font-size:14px"><br></font></div><div><font face="Calibri" style="font-size:14px">Para realizar as estimativas dos parâmetros via o package “Renext” utilize a função “EM.mixexp”. Apesar disso, é importante notar dois detalhes desta distribuição;</font></div><div><font face="Calibri" style="font-size:14px"><br></font></div><div><ul><li><font face="Calibri" style="font-size:14px">O primeiro é o caso do coeficiente de variação maior ou igual a um que tende a gerar “erros” (convergir para uma exponencial simples) quando a amostra apresenta tal coeficiente menor que um. </font></li></ul><div><font face="Calibri" style="font-size:14px"><br></font></div><ul><li><font face="Calibri" style="font-size:14px">Um segundo detalhe é que as estimativas dos parâmetros,no que diz respeito a convergência em relação a amostra, é bem mais “lenta” com o crescimento da amostra do que distribuições usuals como a gama, weibull, gumbel e outras  (o segundo detalhe é apenas um feeling particular).</font></li></ul></div><div><font face="Calibri" style="font-size:14px"><br></font></div><div><font face="Calibri" style="font-size:14px">Caso não queira usar  o “Renext”, outra possibilidade é utilizar o pacote “fitdistplus” para a realização das estimativas dos parâmetros, o único inconveniente é a necessidade de inserir um valor inicial.Por fim, segue os exemplos de estimativas utilizando os dois pacotes diferentes.</font></div><div><font face="Calibri" style="font-size:14px"><br></font></div><div><font face="Calibri" style="font-size:14px"><br></font></div><div><font face="Calibri" style="font-size:14px">#####</font></div><div><font face="Calibri" style="font-size:14px"><br></font></div><div><div><font face="Calibri" style="font-size:14px"># Via Renext</font></div><div><font face="Calibri" style="font-size:14px"><br></font></div><div><font face="Calibri" style="font-size:14px">library("Renext")</font></div><div><font face="Calibri" style="font-size:14px">mix <- rmixexp2(n = 1000, prob1 = 0.5, rate1 = 1.0, rate2 = 3.0)</font></div><div><font face="Calibri" style="font-size:14px">fit_mixexp <- EM.mixexp(x = mix)</font></div><div><font face="Calibri" style="font-size:14px"><br></font></div><div><font face="Calibri" style="font-size:14px"># Parâmetros</font></div><div><font face="Calibri" style="font-size:14px">fit_mixexp[[1]]</font></div><div><font face="Calibri" style="font-size:14px"><br></font></div><div><font face="Calibri" style="font-size:14px"><br></font></div><div><font face="Calibri" style="font-size:14px"># Via fitdistrplus</font></div><div><font face="Calibri" style="font-size:14px"><br></font></div><div><font face="Calibri" style="font-size:14px">library("Renext")</font></div><div><font face="Calibri" style="font-size:14px">library("fitdistrplus")</font></div><div><font face="Calibri" style="font-size:14px"><br></font></div><div><font face="Calibri" style="font-size:14px">fit_mixexp2 <- fitdist(mix, "mixexp2", start = list(prob1 = 0.8,rate1 = 1, rate2 = 2))</font></div><div><font face="Calibri" style="font-size:14px">fit_mixexp2[[1]]</font></div></div><div><div><font face="Calibri" style="font-size:14px"><br></font></div><div><font face="Calibri" style="font-size:14px"><br></font></div><div><font face="Calibri" style="font-size:14px">####</font></div><div><font face="Calibri" style="font-size:14px"><br></font></div><div><font face="Calibri" style="font-size:14px">Obs: No segundo exemplo, apesar de implícito, utilizei as funções pmixexp2, dmixexp2, qmixexp2 do pacote  “Renext” mas caso seja de sua preferência é possível você mesmo defini-las.  </font></div><div><font face="Calibri" style="font-size:14px"><br></font></div><div><font face="Calibri" style="font-size:14px">Att,</font></div><div><div><div style="word-wrap:break-word"><div style="word-wrap:break-word"><font face="Calibri" style="font-size:14px">-------</font></div><div style="word-wrap:break-word"><font face="Calibri" style="font-size:14px">Jaildo Rocha <br>MSc. pelo Programa Pós Graduação em Saneamento Meio Ambiente e Recursos  Hídricos –  UFMG<br>Área de concentração: Recursos Hídricos     <br></font></div></div></div><div><br></div></div></div><div><br></div></span></div></div><div><blockquote type="cite"><span class=""><span><div>Em 20 de set de 2016, à(s) 00:21, Alisson Lucrécio via R-br <<a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a>> escreveu:</div><br></span></span><div><span class=""><span><div dir="ltr" style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><div><span style="font-size:12.8px">Cesar,</span></div><div><span style="font-size:12.8px"><br></span></div><span style="font-size:12.8px">Preciso solucionar essa equação </span><span style="font-size:12.8px">f(x) = (alpha/beta1*exp(-x/beta1)) + ((1-alpha)/beta2*exp(-x/beta2)</span><span style="font-size:12.8px"><wbr>) para obter os paramentos da distribuição. Assim como solucionado pelo </span><span style="font-size:12.8px">fitdistr do pacote MASS para as ditribuções gamma, webull e exponencial. Observe que fiz uma tentativa com o pacote Renext, mas sem sucesso.</span><div style="font-size:12.8px"><span style="font-size:12.8px"><br></span></div><div style="font-size:12.8px"><span style="font-size:12.8px">Muito obrigado.</span></div></div></span></span><div class="gmail_extra" style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><br><div class="gmail_quote"><span class=""><span>2016-09-19 16:12 GMT-03:00 Alisson Lucrécio<span> </span><span dir="ltr"><<a href="mailto:alisson.lucrecio@ifgoiano.edu.br" target="_blank">alisson.lucrecio@ifg<wbr>oiano.edu.br</a>></span>:<br></span></span><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><span class=""><span><div dir="ltr">Como solucionar essa equação <span style="font-size:12.8px">f(x) = (alpha/beta1*exp(-x/beta1)) + ((1-alpha)/beta2*exp(-x/beta2)</span><span style="font-size:12.8px"><wbr>) para obter os paramentos da distribuição. Assim como é feito pelo </span><span style="font-size:12.8px">fitdistr do pacote MASS para as ditribuções gamma, webull e exponencial. Observe que fiz uma tentativa com o pacote Renext, mas sem sucesso.</span><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">Muito obrigado.<br></span><div><div><div><span style="font-size:12.8px"><br></span></div></div></div></div></div></span></span><div><div><div class="gmail_extra"><br><div class="gmail_quote"><span class=""><span>On 19 September 2016 at 16:05, Cesar Rabak<span> </span><span dir="ltr"><<a href="mailto:cesar.rabak@gmail.com" target="_blank">cesar.rabak@gmail.com</a>></span><span> </span><wbr>wrote:<br></span></span><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><span class=""><span><div dir="ltr">. . . e sua dúvida seria?</div></span></span><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="h5"><div><div>On Mon, Sep 19, 2016 at 10:54 AM, Alisson Lucrécio via R-br<span> </span><span dir="ltr"><<a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a><wbr>></span><span> </span>wrote:<br></div></div></div></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div><div class="h5"><div><div><div dir="ltr">Caro Colegas,<div><br></div><div>Bom dia.</div><div><br></div><div>Algum colega da poderia me ajudar na solução de uma distribução exponencial mista.</div><div><br></div><div># Script</div><div><br></div><div><div>library(MASS)</div><div><br></div><div># Distribuição Gamma</div><div><br></div><div>rg <- rgamma(100, shape = 1.5, rate = 1)</div><div><br></div><div>fit_gamma <- fitdistr(rg,"gamma")</div><div><br></div><div>par(bg="white", las=1, cex=1.1)</div><div>plot(density(rg, bw=0.5, cut=0), las=1, lwd=2,</div><div>     xlim=c(0,5),col="steelblue")</div><div><br></div><div># Distribuição Weibull</div><div><br></div><div>rw <- rweibull(100, scale=1, shape=1.5)</div><div><br></div><div>par(bg="white", las=1, cex=1.1)</div><div>plot(density(rw, bw=0.5, cut=0), las=1, lwd=2,</div><div>     xlim=c(0,5),col="steelblue")</div><div><br></div><div>fit_weibull <- fitdistr(rw,"weibull")</div><div><br></div><div># Distribuição Exponential</div><div><br></div><div>ex <- rexp(100, rate = 1)</div><div><br></div><div>par(bg="white", las=1, cex=1.1)</div><div>plot(density(ex, bw=0.5, cut=0), las=1, lwd=2,</div><div>     xlim=c(0,5),col="steelblue")</div><div><br></div><div>fit_exp <- fitdistr(ex, "exponential")</div><div><br></div><div>plot(density(rw, bw=0.5, cut=0), las=1, lwd=2,</div><div>     xlim=c(0,5),col="steelblue")</div><div>lines(density(rg, bw=0.5, cut=0),col="green")</div><div>lines(density(ex, bw=0.5, cut=0),col="red")</div><div><br></div><div># Distribuição Exponencial mista</div><div><br></div><div>library(Renext)</div><div><br></div><div># f(x) = (alpha/beta1*exp(-x/beta1)) + ((1-alpha)/beta2*exp(-x/beta2)<wbr>)</div><div><br></div><div>mixexp <- rmixexp2(n = 100, prob1 = 0.5, rate1 = 1.0, rate2 = 3.0)</div><div><br></div><div>fit_mixexp <- mom.mixexp2(mixexp)</div><div><br></div><div>plot(density(mixexp, bw=0.5, cut=1), lwd=2,</div><div>     xlim=c(0, 5),col="steelblue")</div><div>lines(density(ex, bw=0.5, cut=0),col="red")</div><div><br></div><div>mod <- nls(mixexp ~ prob1 * dexp(x, rate1) + (1 - prob1) * dexp(x, rate2),</div><div>           start = list(prob1 = 1, rate1 = 1, rate2 = 3.0), control = nls.control(maxiter = 100))</div></div><div><br></div><div>Obrigado.</div><span><font color="#888888"><span><font color="#888888"><div><div>--<span> </span><br><div><div dir="ltr"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;font-style:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;float:none;display:inline;background-color:rgb(255,255,255)">Alisson Lucrecio da Costa</span></div></div></div></div></font></span></font></span></div></div></div></div></div><span><font color="#888888"><br>______________________________<wbr>_________________<br>R-br mailing list<br><a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><span class=""><span><br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/cgi<wbr>-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-g<wbr>uia</a>) e forneça código mínimo reproduzível.<br></span></span></font></span></blockquote></div><br></div></blockquote></div><span class=""><span><br><br clear="all"><div><br></div>--<span> </span><br><div data-smartmail="gmail_signature"><div dir="ltr"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;font-style:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;float:none;background-color:rgb(255,255,255);display:inline!important">Alisson Lucrecio da Costa</span></div></div></span></span></div></div></div></blockquote></div><span class=""><span><br><br clear="all"><div><br></div>--<span> </span><br><div data-smartmail="gmail_signature"><div dir="ltr"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;font-style:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;float:none;background-color:rgb(255,255,255);display:inline!important">Alisson Lucrecio da Costa</span></div></div></span></span></div><span style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;float:none;display:inline!important">______________________________<wbr>_________________</span><br style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><span style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;float:none;display:inline!important">R-br mailing list</span><br style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><a href="mailto:R-br@listas.c3sl.ufpr.br" style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px" target="_blank">R-br@listas.c3sl.ufpr.br</a><span class=""><span><br style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px" target="_blank">https://listas.inf.ufpr.br/cgi<wbr>-bin/mailman/listinfo/r-br</a><br style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><span style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;float:none;display:inline!important">Leia o guia de postagem (</span><a href="http://www.leg.ufpr.br/r-br-guia" style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px" target="_blank">http://www.leg.ufpr.br/r-br-g<wbr>uia</a><span style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;float:none;display:inline!important">) e forneça código mínimo reproduzível.</span></span></span></div></blockquote></div><br></div></blockquote></div><span class=""><br><br clear="all"><div><br></div>-- <br><div data-smartmail="gmail_signature"><div dir="ltr"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;display:inline!important;float:none;background-color:rgb(255,255,255)">Alisson Lucrecio da Costa</span></div></div>
</span></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;display:inline!important;float:none;background-color:rgb(255,255,255)">Alisson Lucrecio da Costa</span></div></div>
</div>