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.

Beta-Binomial Regression and Its Application

Variance Analyzed under E-Bayesian Estimation and Hierarchical Bayesian Estimation for Diagnosing Failure Rates

Logistic Regression Model for a Bivariate Binomial Distribution with Applications in Baseball Data Analysis

The Tilted Beta-Binomial Distribution in Overdispersed Data: Maximum Likelihood and Bayesian Estimation

Application of Beta-Distribution and Combined Uniform and Binomial Methods in Longitudinal Modeling of Bounded Outcome Score Data

Robust Inference under the Beta Regression Model with Application to Health Care Studies

Application of Binomial Distribution Function with Discrete Random Variables

A robust regression model for bounded count health data

Thinking Inside the Bounds: Improved Error Distributions for Indifference Point Data Analysis and Simulation Via Beta Regression using Common Discounting Functions

Discrete Beta and Shifted Beta-Binomial models for rating and ranking data

On Beta-Product Convolutions

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

Robust beta regression through the logit transformation

Testing Multivariate Normality Based on Beta-Representative Points

A generalized correlated binomial distribution with application in multiple testing problems

Modeling microbial abundances and dysbiosis with beta-binomial regression

Equivalence Test in Multi-dimensional Space with Applications in A/B Testing

Beta regression for double‐bounded response with correlated high‐dimensional covariates

A new probabilistic model: Theory, simulation and applications to sports and failure times data

Likelihood analysis for a class of beta mixed models

Logistic-beta processes for modeling dependent random probabilities with beta marginals