<div dir="ltr"><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><div class="gmail_extra"><br><div class="gmail_quote">2016-09-19 16:12 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">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><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On 19 September 2016 at 16:05, Cesar Rabak <span dir="ltr"><<a href="mailto:cesar.rabak@gmail.com" target="_blank">cesar.rabak@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">. . . e sua dúvida seria?</div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Sep 19, 2016 at 10:54 AM, Alisson Lucrécio via R-br <span dir="ltr"><<a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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>-- <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><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><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></font></span></blockquote></div><br></div>
</blockquote></div><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>
</div>
</div></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>