# NOT RUN {
data(bodyfat, package = "ugomquantreg")
library(gamlss)
tau <- 0.50
fit.logit <- gamlss(ARMS ~ AGE + I(BMI / 100) + as.factor(SEX) + as.factor(IPAQ),
data = bodyfat, family = UGOM(mu.link = "logit", sigma.link = "log"))
tau <- 0.50;
fit.probit <- gamlss(ARMS ~ AGE + I(BMI / 100) + as.factor(SEX) + as.factor(IPAQ),
data = bodyfat, family = UGOM(mu.link = "probit", sigma.link = "log"))
# }
Run the code above in your browser using DataLab