bodyfat {vasicekreg}

Percentage of body fat data set

Description

The body fat percentage of individuals assisted in a public hospital in Curitiba, Paraná, Brazil.

Usage

data(bodyfat, package = "vasicekreg")

Format

A data-frame with 298 observations and 9 columns:

• ARMS: arms fat percentage.

• LEGS: legs fat percentage.

• BODY: body fat percentage.

• ANDROID: android fat percentage.

• GYNECOID: ginecoid fat percentage.

• AGE: age of individuals.

• BMI: body mass index.

• SEX: 1 for female, 2 for male.

• IPAQ: 0 for IPAQ = sedentary, 1 for IPAQ = insufficiently active and 2 for IPAQ = active.

Author(s)

Josmar Mazucheli [e-mail]

Bruna Alves [e-mail]

Source

http://www.leg.ufpr.br/doku.php/publications:papercompanions:multquasibeta

References

Mazucheli, J., Leiva, V., Alves, B., and Menezes A. F. B., (2021). A new quantile regression for modeling bounded data under a unit Birnbaum-Saunders distribution with applications in medicine and politics. Symmetry, 13(4) 1–21.

Petterle, R. R., Bonat, W. H., Scarpin, C. T., Jonasson, T., and Borba, V. Z. C., (2020). Multivariate quasi-beta regression models for continuous bounded data. The International Journal of Biostatistics, 1–15, (preprint).

Examples

data(bodyfat, package = "vasicekreg")

bodyfat$BMI <- bodyfat$BMI / 100
bodyfat$SEX <- as.factor(bodyfat$SEX)
bodyfat$IPAQ<- as.factor(bodyfat$IPAQ)

library(gamlss)

# mean fit 
fitmean.logit <- gamlss(ARMS ~ AGE + BMI + SEX + IPAQ, 
data = bodyfat, family = VASIM(mu.link = "logit", sigma.link = "logit"))

fitmean.probit <- gamlss(ARMS ~ AGE + BMI + SEX + IPAQ, 
data = bodyfat, family = VASIM(mu.link = "probit", sigma.link = "logit"))

# quantile fit - tau = 0.5

tau <- 0.50
fitquant.logit <- gamlss(ARMS ~ AGE + BMI + SEX + IPAQ, data = bodyfat,
family = VASIQ(mu.link = "logit", sigma.link = "logit"))

fittaus <- lapply(c(0.10, 0.25, 0.50, 0.75, 0.90), function(Tau)
{
 tau <<- Tau; 
 gamlss(ARMS ~ AGE + BMI + SEX + IPAQ, data = bodyfat,
 family = VASIQ(mu.link = "logit", sigma.link = "logit"))
})

sapply(fittaus, summary, USE.NAMES = TRUE)