<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><font face="Calibri" style="font-size: 14px;" class="">Prezado Alisson,</font><div class=""><font face="Calibri" style="font-size: 14px;" class=""><br class=""></font></div><div class=""><font face="Calibri" style="font-size: 14px;" class="">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 class=""><font face="Calibri" style="font-size: 14px;" class=""><br class=""></font></div><div class=""><ul class="MailOutline"><li class=""><font face="Calibri" style="font-size: 14px;" class="">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 class=""><font face="Calibri" style="font-size: 14px;" class=""><br class=""></font></div><ul class="MailOutline"><li class=""><font face="Calibri" style="font-size: 14px;" class="">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 class=""><font face="Calibri" style="font-size: 14px;" class=""><br class=""></font></div><div class=""><font face="Calibri" style="font-size: 14px;" class="">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 class=""><font face="Calibri" style="font-size: 14px;" class=""><br class=""></font></div><div class=""><font face="Calibri" style="font-size: 14px;" class=""><br class=""></font></div><div class=""><font face="Calibri" style="font-size: 14px;" class="">#####</font></div><div class=""><font face="Calibri" style="font-size: 14px;" class=""><br class=""></font></div><div class=""><div class=""><font face="Calibri" style="font-size: 14px;" class=""># Via Renext</font></div><div class=""><font face="Calibri" style="font-size: 14px;" class=""><br class=""></font></div><div class=""><font face="Calibri" style="font-size: 14px;" class="">library("Renext")</font></div><div class=""><font face="Calibri" style="font-size: 14px;" class="">mix <- rmixexp2(n = 1000, prob1 = 0.5, rate1 = 1.0, rate2 = 3.0)</font></div><div class=""><font face="Calibri" style="font-size: 14px;" class="">fit_mixexp <- EM.mixexp(x = mix)</font></div><div class=""><font face="Calibri" style="font-size: 14px;" class=""><br class=""></font></div><div class=""><font face="Calibri" style="font-size: 14px;" class=""># Parâmetros</font></div><div class=""><font face="Calibri" style="font-size: 14px;" class="">fit_mixexp[[1]]</font></div><div class=""><font face="Calibri" style="font-size: 14px;" class=""><br class=""></font></div><div class=""><font face="Calibri" style="font-size: 14px;" class=""><br class=""></font></div><div class=""><font face="Calibri" style="font-size: 14px;" class=""># Via fitdistrplus</font></div><div class=""><font face="Calibri" style="font-size: 14px;" class=""><br class=""></font></div><div class=""><font face="Calibri" style="font-size: 14px;" class="">library("Renext")</font></div><div class=""><font face="Calibri" style="font-size: 14px;" class="">library("fitdistrplus")</font></div><div class=""><font face="Calibri" style="font-size: 14px;" class=""><br class=""></font></div><div class=""><font face="Calibri" style="font-size: 14px;" class="">fit_mixexp2 <- fitdist(mix, "mixexp2", start = list(prob1 = 0.8,rate1 = 1, rate2 = 2))</font></div><div class=""><font face="Calibri" style="font-size: 14px;" class="">fit_mixexp2[[1]]</font></div></div><div class=""><div class=""><font face="Calibri" style="font-size: 14px;" class=""><br class=""></font></div><div class=""><font face="Calibri" style="font-size: 14px;" class=""><br class=""></font></div><div class=""><font face="Calibri" style="font-size: 14px;" class="">####</font></div><div class=""><font face="Calibri" style="font-size: 14px;" class=""><br class=""></font></div><div class=""><font face="Calibri" style="font-size: 14px;" class="">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 class=""><font face="Calibri" style="font-size: 14px;" class=""><br class=""></font></div><div class=""><font face="Calibri" style="font-size: 14px;" class="">Att,</font></div><div class=""><div class=""><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><font face="Calibri" style="font-size: 14px;" class="">-------</font></div><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><font face="Calibri" style="font-size: 14px;" class="">Jaildo Rocha <br class="">MSc. pelo Programa Pós Graduação em Saneamento Meio Ambiente e Recursos Hídricos – UFMG<br class="">Área de concentração: Recursos Hídricos <br class=""></font></div></div></div><div class=""><br class="webkit-block-placeholder"></div></div></div><div class=""><br class=""></div><div><blockquote type="cite" class=""><div class="">Em 20 de set de 2016, à(s) 00:21, Alisson Lucrécio via R-br <<a href="mailto:r-br@listas.c3sl.ufpr.br" class="">r-br@listas.c3sl.ufpr.br</a>> escreveu:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><div class=""><span style="font-size: 12.8px;" class="">Cesar,</span></div><div class=""><span style="font-size: 12.8px;" class=""><br class=""></span></div><span style="font-size: 12.8px;" class="">Preciso solucionar essa equação </span><span style="font-size: 12.8px;" class="">f(x) = (alpha/beta1*exp(-x/beta1)) + ((1-alpha)/beta2*exp(-x/beta2)</span><span style="font-size: 12.8px;" class=""><wbr class="">) para obter os paramentos da distribuição. Assim como solucionado pelo </span><span style="font-size: 12.8px;" class="">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;" class=""><span style="font-size: 12.8px;" class=""><br class=""></span></div><div style="font-size: 12.8px;" class=""><span style="font-size: 12.8px;" class="">Muito obrigado.</span></div></div><div class="gmail_extra" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;"><br class=""><div class="gmail_quote">2016-09-19 16:12 GMT-03:00 Alisson Lucrécio<span class="Apple-converted-space"> </span><span dir="ltr" class=""><<a href="mailto:alisson.lucrecio@ifgoiano.edu.br" target="_blank" class="">alisson.lucrecio@ifgoiano.edu.br</a>></span>:<br class=""><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 dir="ltr" class="">Como solucionar essa equação <span style="font-size: 12.8px;" class="">f(x) = (alpha/beta1*exp(-x/beta1)) + ((1-alpha)/beta2*exp(-x/beta2)</span><span style="font-size: 12.8px;" class=""><wbr class="">) para obter os paramentos da distribuição. Assim como é feito pelo </span><span style="font-size: 12.8px;" class="">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 class=""><span style="font-size: 12.8px;" class=""><br class=""></span></div><div class=""><span style="font-size: 12.8px;" class="">Muito obrigado.<br class=""></span><div class=""><div class=""><div class=""><span style="font-size: 12.8px;" class=""><br class=""></span></div></div></div></div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br class=""><div class="gmail_quote">On 19 September 2016 at 16:05, Cesar Rabak<span class="Apple-converted-space"> </span><span dir="ltr" class=""><<a href="mailto:cesar.rabak@gmail.com" target="_blank" class="">cesar.rabak@gmail.com</a>></span><span class="Apple-converted-space"> </span>wrote:<br class=""><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 dir="ltr" class="">. . . e sua dúvida seria?</div><div class="gmail_extra"><br class=""><div class="gmail_quote">On Mon, Sep 19, 2016 at 10:54 AM, Alisson Lucrécio via R-br<span class="Apple-converted-space"> </span><span dir="ltr" class=""><<a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank" class="">r-br@listas.c3sl.ufpr.br</a>></span><span class="Apple-converted-space"> </span>wrote:<br class=""><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 dir="ltr" class="">Caro Colegas,<div class=""><br class=""></div><div class="">Bom dia.</div><div class=""><br class=""></div><div class="">Algum colega da poderia me ajudar na solução de uma distribução exponencial mista.</div><div class=""><br class=""></div><div class=""># Script</div><div class=""><br class=""></div><div class=""><div class="">library(MASS)</div><div class=""><br class=""></div><div class=""># Distribuição Gamma</div><div class=""><br class=""></div><div class="">rg <- rgamma(100, shape = 1.5, rate = 1)</div><div class=""><br class=""></div><div class="">fit_gamma <- fitdistr(rg,"gamma")</div><div class=""><br class=""></div><div class="">par(bg="white", las=1, cex=1.1)</div><div class="">plot(density(rg, bw=0.5, cut=0), las=1, lwd=2,</div><div class=""> xlim=c(0,5),col="steelblue")</div><div class=""><br class=""></div><div class=""># Distribuição Weibull</div><div class=""><br class=""></div><div class="">rw <- rweibull(100, scale=1, shape=1.5)</div><div class=""><br class=""></div><div class="">par(bg="white", las=1, cex=1.1)</div><div class="">plot(density(rw, bw=0.5, cut=0), las=1, lwd=2,</div><div class=""> xlim=c(0,5),col="steelblue")</div><div class=""><br class=""></div><div class="">fit_weibull <- fitdistr(rw,"weibull")</div><div class=""><br class=""></div><div class=""># Distribuição Exponential</div><div class=""><br class=""></div><div class="">ex <- rexp(100, rate = 1)</div><div class=""><br class=""></div><div class="">par(bg="white", las=1, cex=1.1)</div><div class="">plot(density(ex, bw=0.5, cut=0), las=1, lwd=2,</div><div class=""> xlim=c(0,5),col="steelblue")</div><div class=""><br class=""></div><div class="">fit_exp <- fitdistr(ex, "exponential")</div><div class=""><br class=""></div><div class="">plot(density(rw, bw=0.5, cut=0), las=1, lwd=2,</div><div class=""> xlim=c(0,5),col="steelblue")</div><div class="">lines(density(rg, bw=0.5, cut=0),col="green")</div><div class="">lines(density(ex, bw=0.5, cut=0),col="red")</div><div class=""><br class=""></div><div class=""># Distribuição Exponencial mista</div><div class=""><br class=""></div><div class="">library(Renext)</div><div class=""><br class=""></div><div class=""># f(x) = (alpha/beta1*exp(-x/beta1)) + ((1-alpha)/beta2*exp(-x/beta2)<wbr class="">)</div><div class=""><br class=""></div><div class="">mixexp <- rmixexp2(n = 100, prob1 = 0.5, rate1 = 1.0, rate2 = 3.0)</div><div class=""><br class=""></div><div class="">fit_mixexp <- mom.mixexp2(mixexp)</div><div class=""><br class=""></div><div class="">plot(density(mixexp, bw=0.5, cut=1), lwd=2,</div><div class=""> xlim=c(0, 5),col="steelblue")</div><div class="">lines(density(ex, bw=0.5, cut=0),col="red")</div><div class=""><br class=""></div><div class="">mod <- nls(mixexp ~ prob1 * dexp(x, rate1) + (1 - prob1) * dexp(x, rate2),</div><div class=""> start = list(prob1 = 1, rate1 = 1, rate2 = 3.0), control = nls.control(maxiter = 100))</div></div><div class=""><br class=""></div><div class="">Obrigado.</div><span class=""><font color="#888888" class=""><span class=""><font color="#888888" class=""><div class=""><div class="">--<span class="Apple-converted-space"> </span><br class=""><div class=""><div dir="ltr" class=""><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);" class="">Alisson Lucrecio da Costa</span></div></div></div></div></font></span></font></span></div><span class=""><font color="#888888" class=""><br class="">______________________________<wbr class="">_________________<br class="">R-br mailing list<br class=""><a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank" class="">R-br@listas.c3sl.ufpr.br</a><br class=""><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" target="_blank" class="">https://listas.inf.ufpr.br/cgi<wbr class="">-bin/mailman/listinfo/r-br</a><br class="">Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="noreferrer" target="_blank" class="">http://www.leg.ufpr.br/r-br-g<wbr class="">uia</a>) e forneça código mínimo reproduzível.<br class=""></font></span></blockquote></div><br class=""></div></blockquote></div><br class=""><br clear="all" class=""><div class=""><br class=""></div>--<span class="Apple-converted-space"> </span><br class=""><div data-smartmail="gmail_signature" class=""><div dir="ltr" class=""><span style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 13px; font-style: normal; font-variant-ligatures: normal; font-variant-position: normal; font-variant-caps: normal; font-variant-numeric: normal; font-variant-alternates: normal; font-variant-east-asian: 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;" class="">Alisson Lucrecio da Costa</span></div></div></div></div></div></blockquote></div><br class=""><br clear="all" class=""><div class=""><br class=""></div>--<span class="Apple-converted-space"> </span><br class=""><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr" class=""><span style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 13px; font-style: normal; font-variant-ligatures: normal; font-variant-position: normal; font-variant-caps: normal; font-variant-numeric: normal; font-variant-alternates: normal; font-variant-east-asian: 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;" class="">Alisson Lucrecio da Costa</span></div></div></div><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">_______________________________________________</span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">R-br mailing list</span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><a href="mailto:R-br@listas.c3sl.ufpr.br" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class="">R-br@listas.c3sl.ufpr.br</a><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class="">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">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-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class="">http://www.leg.ufpr.br/r-br-guia</a><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">) e forneça código mínimo reproduzível.</span></div></blockquote></div><br class=""></body></html>