data("bodyfat_Aeolus")
fitPL <- PLreg(percentfat ~ 1, data = bodyfat_Aeolus,
family = "SN", zeta = 1.6)
fitPL
summary(fitPL)
coef(fitPL, model = "median")
vcov(fitPL)
logLik(fitPL)
AIC(fitPL)
Run the code above in your browser using DataLab