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.

Effect Displays for Multinomial and Proportional-Odds Logit Models

Non Proportional Odds Models are Widely Dispensable -- Sparser Modeling based on Parametric and Additive Location-Shift Approaches

Discrete Exponential-Family Models for Multivariate Binary Outcomes

How Much Should We Trust Estimates from Multiplicative Interaction Models? Simple Tools to Improve Empirical Practice

Displaying Variation in Large Datasets: Plotting a Visual Summary of Effect Sizes

Multinomial Models with Linear Inequality Constraints: Overview and Improvements of Computational Methods for Bayesian Inference

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

On the "Poisson Trick" and its Extensions for Fitting Multinomial Regression Models

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

Power Estimation Using the Exemplary Dataset for Ordinal Data in the Presence of Nonproportional Odds

Testing the binomial fixed effects logit model, with an application to female labour supply

Ordinal Probit Functional Outcome Regression with Application to Computer-Use Behavior in Rhesus Monkeys

The Sorted Effects Method: Discovering Heterogeneous Effects Beyond Their Averages

Multivariate probit linear mixed models for multivariate longitudinal binary data

Individual-level probabilities and cluster-level proportions: Toward interpretable level 2 estimates in unconflated multilevel models for binary outcomes.

Logistic and multinomial-logit models: A brief review on their modifications and extensions

Generating contingency tables with fixed marginal probabilities and dependence structures described by loglinear models

Simple or complex statistical models: Non-traditional regression models with intuitive interpretations

A menu-driven software package of Bayesian nonparametric (and parametric) mixed models for regression analysis and density estimation

Mathematical Formalization and Applications to Data with Excess of Zeros and Ones of the Unit-Proportional Hazard Inflated Models

Copula-Based Random Effects Models for Clustered Data