<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))</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><div><div>-- <br><div class="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-ligatures:normal;font-variant-caps: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></div>