RStudio_distreg.vis_unititled_9B4DE965

library("distreg.vis")
library("gamlss")
library("bamlss")
library("dplyr", warn.conflicts = FALSE)
### Fit a Lognormal model
wage <- ISLR::Wage
model <- gamlss(wage ~ year + ps(age) + education, 
                ~ year + ps(age), data = wage,
                family = "LOGNO")
ndata <- wage %>%
  select(age, year, education) %>%
  sample_n(5)
plot_moments(model, int_var = "age", pred_data = ndata)

## Inequality function
ineq <- function(par) {
  2 * pnorm((par[["sigma"]] / 2) * sqrt(2)) - 1
}

## Pred_params
pred <- preds(model, ndata)
apply(pred, 1, FUN = ineq)