Estou interessado em realizar um histograma da variável resposta com a densidade estimada de dois modelos de regressão ajustado pelo pacote gamlss
.
A ideia consiste em se obter um gráfico "próximo" ao exemplo que está sendo apresentado abaixo.
Porém, como já foi descrito, o meu interesse é realizar o traço da densidade estimada por meio dos modelos de regressão ajustado com diferentes distribuições para a variável resposta. Neste fórum, encontra-se inúmeros exemplos considerando simplesmente o traçado para checar a aderência da distribuição e/ou exemplos quando o modelo de regressão não contém covariáveis.
Abaixo encontra-se a rotina computacional de minha tentativa.
library(gamlss)
head(Orange)
mod1 = gamlss(circumference ~ age,
family=NO, data = Orange)
mod2 = gamlss(circumference ~ age,
family=GA, data = Orange)
x11()
hist(Orange$circumference)
lines(predict(mod1))
lines(predict(mod2))
Todavia, não estou obtendo sucesso na construção de tal gráfico.