Caro Colegas, Bom dia. Algum colega da poderia me ajudar na solução de uma distribução exponencial mista. # Script library(MASS) # Distribuição Gamma rg <- rgamma(100, shape = 1.5, rate = 1) fit_gamma <- fitdistr(rg,"gamma") par(bg="white", las=1, cex=1.1) plot(density(rg, bw=0.5, cut=0), las=1, lwd=2, xlim=c(0,5),col="steelblue") # Distribuição Weibull rw <- rweibull(100, scale=1, shape=1.5) par(bg="white", las=1, cex=1.1) plot(density(rw, bw=0.5, cut=0), las=1, lwd=2, xlim=c(0,5),col="steelblue") fit_weibull <- fitdistr(rw,"weibull") # Distribuição Exponential ex <- rexp(100, rate = 1) par(bg="white", las=1, cex=1.1) plot(density(ex, bw=0.5, cut=0), las=1, lwd=2, xlim=c(0,5),col="steelblue") fit_exp <- fitdistr(ex, "exponential") plot(density(rw, bw=0.5, cut=0), las=1, lwd=2, xlim=c(0,5),col="steelblue") lines(density(rg, bw=0.5, cut=0),col="green") lines(density(ex, bw=0.5, cut=0),col="red") # Distribuição Exponencial mista library(Renext) # f(x) = (alpha/beta1*exp(-x/beta1)) + ((1-alpha)/beta2*exp(-x/beta2)) mixexp <- rmixexp2(n = 100, prob1 = 0.5, rate1 = 1.0, rate2 = 3.0) fit_mixexp <- mom.mixexp2(mixexp) plot(density(mixexp, bw=0.5, cut=1), lwd=2, xlim=c(0, 5),col="steelblue") lines(density(ex, bw=0.5, cut=0),col="red") mod <- nls(mixexp ~ prob1 * dexp(x, rate1) + (1 - prob1) * dexp(x, rate2), start = list(prob1 = 1, rate1 = 1, rate2 = 3.0), control = nls.control(maxiter = 100)) Obrigado. -- Alisson Lucrecio da Costa