Abstract:BackgroundIn biomedical research, response variables are often encountered which have bounded support on the open unit interval - (0,1). Traditionally, researchers have attempted to estimate covariate effects on these types of response data using linear regression. Alternative modelling strategies may include: beta regression, variable-dispersion beta regression, and fractional logit regression models. This study employs a Monte Carlo simulation design to compare the statistical properties of the linear regression model to that of the more novel beta regression, variable-dispersion beta regression, and fractional logit regression models.MethodsIn the Monte Carlo experiment we assume a simple two sample design. We assume observations are realizations of independent draws from their respective probability models. The randomly simulated draws from the various probability models are chosen to emulate average proportion/percentage/rate differences of pre-specified magnitudes. Following simulation of the experimental data we estimate average proportion/percentage/rate differences. We compare the estimators in terms of bias, variance, type-1 error and power. Estimates of Monte Carlo error associated with these quantities are provided.ResultsIf response data are beta distributed with constant dispersion parameters across the two samples, then all models are unbiased and have reasonable type-1 error rates and power profiles. If the response data in the two samples have different dispersion parameters, then the simple beta regression model is biased. When the sample size is small (N0 = N1 = 25) linear regression has superior type-1 error rates compared to the other models. Small sample type-1 error rates can be improved in beta regression models using bias correction/reduction methods. In the power experiments, variable-dispersion beta regression and fractional logit regression models have slightly elevated power compared to linear regression models. Similar results were observed if the response data are generated from a discrete multinomial distribution with support on (0,1).ConclusionsThe linear regression model, the variable-dispersion beta regression model and the fractional logit regression model all perform well across the simulation experiments under consideration. When employing beta regression to estimate covariate effects on (0,1) response data, researchers should ensure their dispersion sub-model is properly specified, else inferential errors could arise.

Regression analysis of proportional data using simplex distribution

Simplexreg: an R Package for Regression Analysis of Proportional Data Using the Simplex Distribution

Modified Poisson Regression Model for Data of Prospective Studies with Common Outcomes

Marginal Models for Longitudinal Continuous Proportional Data

A new regression model for rates and proportions data with applications

Time series regression models for zero-inflated proportions

PROPORTIONAL HAZARDS MODELS WITH CONTINUOUS MARKS

A Monte Carlo simulation study comparing linear regression, beta regression, variable-dispersion beta regression and fractional logit regression at recovering average difference measures in a two sample design

Regression analysis of bivariate current status data under the proportional hazards model

A Simulation Study for Propensity Score in Dealing with Col- linearity Data

A bootstrap semiparametric homogeneity test for the distributions of multigroup proportional data, with applications to analysis of quality of life outcomes in clinical trials

Robust Transformation Mixed‐effects Models for Longitudinal Continuous Proportional Data

A Partially Linear Regression Model for Data from an Outcome-Dependent Sampling Design

Applications of Proportional Odds Model to Multivariate Interval-Censored Data

A comparison of methods for the analysis of binomial proportion data in behavioral research

Extended arcsine distribution to proportional data: Properties and applications

Proportional Hazards Model with Varying Coefficients for Length-Biased Data

Double Probability Integral Transform Residuals for Regression Models with Discrete Outcomes

Logistic regression analysis with standardized markers

Proportional rates models for multivariate panel count data

Using a monotonic density ratio model to increase the power of the goodness-of-fit test for logistic regression models with case-control data