data(obesity.levels)
# short-hand variable names.
r = obesity.levels[, "NObeyesdad"]
s = obesity.levels[, c("Gender", "Age")]
p = obesity.levels[, setdiff(names(obesity.levels), c("NObeyesdad", "Gender", "Age"))]
if (FALSE) {
# the lambda = 0.1 is very helpful in making model estimation succeed.
m = fgrrm(response = r, sensitive = s, predictors = p, ,
family = "multinomial", unfairness = 0.05, lambda = 0.1)
summary(m)
}
Run the code above in your browser using DataLab