Description Usage Format Author(s) Source References Examples
The body fat percentage of individuals assisted in a public hospital in Curitiba, Paraná, Brazil.
1 |
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.
Josmar Mazucheli jmazucheli@gmail.com
Bruna Alves pg402900@uem.br
http://www.leg.ufpr.br/doku.php/publications:papercompanions:multquasibeta
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).
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 | 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)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.