<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 class="HOEnZb"><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></div>
<br>______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/<wbr>cgi-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-<wbr>guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div>