[R-br] Distribuição exponencial mista
Alisson Lucrécio
alisson.lucrecio em ifgoiano.edu.br
Segunda Setembro 19 10:54:25 BRT 2016
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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160919/764dfb8d/attachment.html>
Mais detalhes sobre a lista de discussão R-br