TY - JOUR
T1 - Augmenting beta regression for periodontal proportion data via the SAS NLMIXED procedure
AU - Lewis, Bradley R.
AU - Bandyopadhyay, Dipankar
AU - Desantis, Stacia M.
AU - John, Mike T.
N1 - Publisher Copyright:
© ISOSS Publications 2017.
PY - 2017
Y1 - 2017
N2 - Often in clinical dental research, clinical attachment level (CAL) is recorded at several sites throughout the mouth to assess the extent of periodontal disease (PD). One might be interested to quantify PD at the tooth-level via the proportion of diseased sites per tooth type (say, incisors, canines, pre-molars and molars) per subject. However, these studies might consist of relatively diseasefree and highly diseased subjects leading to the proportion responses distributed in the interval [0, 1]. While beta regression (BR) is often the model of choice to assess covariate effects for proportion data, the presence (and/or abundance) of zeros and/or ones makes it inapplicable here because the beta support is defined in the interval (0, 1). Avoiding ad hoc data transformation, we explore the potential of the augmented BR framework which augments the beta density with non-zero masses at zero and one while accounting for the clustering induced. Our classical estimation framework using maximum likelihood utilizes the potential of the SAS R Proc NLMIXED procedure. We explore our methodology via simulation studies and application to a real cross-sectional dataset on PD, and we assess the gain in model fit and parameter estimation over other ad hoc alternatives. This reveals newer insights into risk quantification on clustered proportion responses. Our methods can be implemented using standard SAS software routines. The augmented BR model results in a better fit to clustered periodontal proportion data over the standard beta model. We recommend using it as a parametric alternative for fitting proportion data, and avoid ad hoc data transformation.
AB - Often in clinical dental research, clinical attachment level (CAL) is recorded at several sites throughout the mouth to assess the extent of periodontal disease (PD). One might be interested to quantify PD at the tooth-level via the proportion of diseased sites per tooth type (say, incisors, canines, pre-molars and molars) per subject. However, these studies might consist of relatively diseasefree and highly diseased subjects leading to the proportion responses distributed in the interval [0, 1]. While beta regression (BR) is often the model of choice to assess covariate effects for proportion data, the presence (and/or abundance) of zeros and/or ones makes it inapplicable here because the beta support is defined in the interval (0, 1). Avoiding ad hoc data transformation, we explore the potential of the augmented BR framework which augments the beta density with non-zero masses at zero and one while accounting for the clustering induced. Our classical estimation framework using maximum likelihood utilizes the potential of the SAS R Proc NLMIXED procedure. We explore our methodology via simulation studies and application to a real cross-sectional dataset on PD, and we assess the gain in model fit and parameter estimation over other ad hoc alternatives. This reveals newer insights into risk quantification on clustered proportion responses. Our methods can be implemented using standard SAS software routines. The augmented BR model results in a better fit to clustered periodontal proportion data over the standard beta model. We recommend using it as a parametric alternative for fitting proportion data, and avoid ad hoc data transformation.
KW - Augmented distribution
KW - Beta regression
KW - Likelihood functions
KW - Periodontal disease
KW - Proc NLMIXED
KW - Proportion data
UR - http://www.scopus.com/inward/record.url?scp=85027339076&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85027339076&partnerID=8YFLogxK
M3 - Article
AN - SCOPUS:85027339076
SN - 1930-6792
VL - 12
SP - 49
EP - 66
JO - Journal of Applied Probability and Statistics
JF - Journal of Applied Probability and Statistics
IS - 1
ER -