Meu interesse é realizar um histograma da variável resposta com a densidade estimada conjuntamente com a densidade empírica de três modelos de regressão ajustados pelo pacote gamlss
, onde esses modelos possuem covariáveis. A ideia é obter um gráfico “próximo” ao exemplo mostrado abaixo.
Minha tentativa de realizar tal procedimento foi, no máximo, traçar a densidade empírica dessas distribuições, utilizando os seguintes códigos:
library(gamlss)
data(Orange)
mod1 = gamlss(circumference ~ age,
family=NO(mu.link = log), data = Orange)
mod2 = gamlss(circumference ~ age,
family=GA, data = Orange)
mod3 = gamlss(circumference ~ age,
family=WEI, data = Orange)
ymod1 <- predict(mod1)
ymod2 <- predict(mod2)
ymod3 <- predict(mod3)
dens1 <- density(exp(ymod1))
dens2 <- density(exp(ymod2))
dens3 <- density(exp(ymod3))
x11()
par(cex = 1.1)
hist(Orange$circumference, freq = FALSE,
xlab = "Circunferência", ylab = "Densidade", main=" ")
lines(dens1, col = "red", lwd = 2)
lines(dens2, col = "blue", lwd = 2)
lines(dens3, col = "green", lwd = 2)
legend(x = "topright", legend = c("Normal", "Gama", "Weibull"),
lty = c(1,1,1),
col = c("red", "blue", "green"))
No entanto, observe que, como mencionei, as densidades realizadas acima são empíricas, pois (a densidade de tais distribuições não são bimodais). Acredito que seja necessário utilizar os parâmetros estimados, mas não sei como fazer isso e/ou talvez expressar para uma única observação tal fato levando em consideração a estrutura do preditor linear e as especificidades algébricas em termos dos parâmetros das distribuições (Ex: parâmetro de forma da distribuição Weibull é dado por outra expressão de acordo com a reparametrização utilizada no R), o mesmo para a distribuição Gama.