A flexible model for the mean and variance functions, with application to medical cost data

Health Care Provider Clustering Using Fusion Penalty in Quasi-Likelihood

There has been growing research interest in developing methodology to evaluate the health care providers' performance with respect to a patient outcome. Random and fixed effects models are traditionally used for such a purpose. We propose a new method, using a fusion penalty to cluster health care providers based on quasi-likelihood. Without any priori knowledge of grouping information, our method provides a desirable data-driven approach for automatically clustering health care providers into different groups based on their performance. Further, the quasi-likelihood is more flexible and robust than the regular likelihood in that no distributional assumption is needed. An efficient alternating direction method of multipliers algorithm is developed to implement the proposed method. We show that the proposed method enjoys the oracle properties; namely, it performs as well as if the true group structure were known in advance. The consistency and asymptotic normality of the estimators are established. Simulation studies and analysis of the national kidney transplant registry data demonstrate the utility and validity of our method.

A flexible quasi-likelihood model for microbiome abundance count data

In this article, we present a flexible model for microbiome count data. We consider a quasi-likelihood framework, in which we do not make any assumptions on the distribution of the microbiome count except that its variance is an unknown but smooth function of the mean. By comparing our model to the negative binomial generalized linear model (GLM) and Poisson GLM in simulation studies, we show that our flexible quasi-likelihood method yields valid inferential results. Using a real microbiome study, we demonstrate the utility of our method by examining the relationship between adenomas and microbiota. We also provide an R package "fql" for the application of our method.

Estimating health care costs at scale in low- and middle-income countries: Mathematical notations and frameworks for the application of cost functions

Appropriate costing and economic modeling are major factors for the successful scale-up of health interventions. Various cost functions are currently being used to estimate costs of health interventions at scale in low- and middle-income countries (LMICs) potentially resulting in disparate cost projections. The aim of this study is to gain understanding of current methods used and provide guidance to inform the use of cost functions that is fit for purpose. We reviewed seven databases covering the economic and global health literature to identify studies reporting a quantitative analysis of costs informing the projected scale-up of a health intervention in LMICs between 2003 and 2019. Of the 8725 articles identified, 40 met the inclusion criteria. We classified studies according to the type of cost functions applied-accounting or econometric-and described the intended use of cost projections. Based on these findings, we developed new mathematical notations and cost function frameworks for the analysis of healthcare costs at scale in LMICs setting. These notations estimate variable returns to scale in cost projection methods, which is currently ignored in most studies. The frameworks help to balance simplicity versus accuracy and increase the overall transparency in reporting of methods.

A semi-parametric Bayesian model for semi-continuous longitudinal data

Semi-continuous data present challenges in both model fitting and interpretation. Parametric distributions may be inappropriate for extreme long right tails of the data. Mean effects of covariates, susceptible to extreme values, may fail to capture relevant information for most of the sample. We propose a two-component semi-parametric Bayesian mixture model, with the discrete component captured by a probability mass (typically at zero) and the continuous component of the density modeled by a mixture of B-spline densities that can be flexibly fit to any data distribution. The model includes random effects of subjects to allow for application to longitudinal data. We specify prior distributions on parameters and perform model inference using a Markov chain Monte Carlo (MCMC) Gibbs-sampling algorithm programmed in R. Statistical inference can be made for multiple quantiles of the covariate effects simultaneously providing a comprehensive view. Various MCMC sampling techniques are used to facilitate convergence. We demonstrate the performance and the interpretability of the model via simulations and analyses on the National Consortium on Alcohol and Neurodevelopment in Adolescence study (NCANDA) data on alcohol binge drinking.

Empirical likelihood-based inferences for median medical cost regression models with censored data

Recent studies show that medical cost data can be heavily censored and highly skewed, which leads to have more complex cost data analysis. In this paper, we propose influence function and empirical likelihood (EL)-based methods to construct confidence regions for regression parameters in median cost regression models with censored data. We further propose confidence intervals for the median cost with given covariates using the proposed EL-based confidence regions. Simulation studies are conducted to compare the proposed EL-based confidence regions with the existing normal approximation-based confidence regions in terms of coverage probabilities. The new EL-based methods are observed to have better finite sample performances than existing methods particularly when the censoring proportion is high. The new methods are also illustrated through a real data example.

Use of Claims Data for Cost and Cost-Effectiveness Research

Administrative claims data are big data generated from healthcare encounters. Claims data contain information on insurance payment as well as clinical diagnoses and procedure codes to ascertain medical conditions and treatments, making them valuable sources for economic evaluation research. This paper offers an introductory overview of the use of claims data for oncology-related cost-of-illness, cost comparison, and cost-effectiveness analyses. We reviewed analytical methods commonly employed in these analyses, such as the phase of care approach and net costing method for cost-of-illness studies, propensity score matching methods for cost comparison studies, and net benefit regression models for cost-effectiveness studies. We used published studies to explain each method and to discuss methodological challenges of conducting economic studies using claims data.

Simple Quasi-Bayes Approach for Modeling Mean Medical Costs

AbstractSeveral statistical issues associated with health care costs, such as heteroscedasticity and severe skewness, make it challenging to estimate or predict medical costs. When the interest is modeling the mean cost, it is desirable to make no assumption on the density function or higher order moments. Another challenge in developing cost prediction models is the presence of many covariates, making it necessary to apply variable selection methods to achieve a balance of prediction accuracy and model simplicity. We propose Spike-or-Slab priors for Bayesian variable selection based on asymptotic normal estimates of the full model parameters that are consistent as long as the assumption on the mean cost is satisfied. In addition, the scope of model searching can be reduced by ranking the Z-statistics. This method possesses four advantages simultaneously: robust (due to avoiding assumptions on the density function or higher order moments), parsimonious (feature of variable selection), informative (due to its Bayesian flavor, which can compare posterior probabilities of candidate models) and efficient (by reducing model searching scope with the use of Z-ranking). We apply this method to the Medical Expenditure Panel Survey dataset.

A flexible quantile regression model for medical costs with application to Medical Expenditure Panel Survey Study

Medical costs are often skewed to the right and heteroscedastic, having a sophisticated relation with covariates. Mean function regression models with low-dimensional covariates have been extensively considered in the literature. However, it is important to develop a robust alternative to find the underlying relationship between medical costs and high-dimensional covariates. In this paper, we propose a new quantile regression model to analyze medical costs. We also consider variable selection, using an adaptive lasso penalized variable selection method to identify significant factors of the covariates. Simulation studies are conducted to illustrate the performance of the estimation method. We apply our method to the analysis of the Medical Expenditure Panel Survey dataset.

Semiparametric Estimation of Longitudinal Medical Cost Trajectory

Estimating the average monthly medical costs from disease diagnosis to a terminal event such as death for an incident cohort of patients is a topic of immense interest to researchers in health policy and health economics because patterns of average monthly costs over time reveal how medical costs vary across phases of care. The statistical challenges to estimating monthly medical costs longitudinally are multifold; the longitudinal cost trajectory (formed by plotting the average monthly costs from diagnosis to the terminal event) is likely to be nonlinear, with its shape depending on the time of the terminal event, which can be subject to right censoring. The goal of this paper is to tackle this statistically challenging topic by estimating the conditional mean cost at any month t given the time of the terminal event s. The longitudinal cost trajectories with different terminal event times form a bivariate surface of t and s, under the constraint t ≤ s. We propose to estimate this surface using bivariate penalized splines in an Expectation-Maximization algorithm that treats the censored terminal event times as missing data. We evaluate the proposed model and estimation method in simulations and apply the method to the medical cost data of an incident cohort of stage IV breast cancer patients from the Surveillance, Epidemiology and End Results-Medicare Linked Database.

A marginalized two-part model with heterogeneous variance for semicontinuous data

Semicontinuous data, characterized by a point mass at zero followed by a positive, continuous distribution, arise frequently in medical research. These data are typically analyzed using two-part mixtures that separately model the probability of incurring a positive outcome and the distribution of positive values among those who incur them. In such a conditional specification, however, standard two-part models do not provide a marginal interpretation of covariate effects on the overall population. We have previously proposed a marginalized two-part model that yields more interpretable effect estimates by parameterizing the model in terms of the marginal mean. In the original formulation, a constant variance was assumed for the positive values. We now extend this model to a more general framework by allowing non-constant variance to be explicitly modeled as a function of covariates, and incorporate this variance into two flexible distributional assumptions, log-skew-normal and generalized gamma, both of which take the log-normal distribution as a special case. Using simulation studies, we compare the performance of each of these models with respect to bias, coverage, and efficiency. We illustrate the proposed modeling framework by evaluating the effect of a behavioral weight loss intervention on health care expenditures in the Veterans Affairs health system.

The Portion of Health Care Costs Associated With Lifestyle-Related Modifiable Health Risks Based on a Sample of 223,461 Employees in Seven Industries: The UM-HMRC Study

OBJECTIVE

This study estimates the percent of health care costs associated with employees' modifiable health risks.

METHODS

Cross-sectional multivariate analysis of 223,461 employees from seven industries who completed a health risk assessment during 2007 to 2012.

RESULTS

Modifiable health risks were associated with 26.0% of health care costs ($761/person) among employees with no self-reported medical conditions and 25.4% among employees with a medical condition ($2598/person). The prevalence and relative costs of each of the 10 risks were different for those without and with medical conditions, but high body mass index was the most prevalent risk for both groups (41.0% and 63.9%) and also contributed the largest percentage of excess costs (7.2% and 7.3%).

CONCLUSIONS

This study, coupled with past work, gives an employer a sense of the magnitude that might be saved if modifiable health risks could be eliminated.

A flexible model for correlated medical costs, with application to medical expenditure panel survey data

We propose a flexible model for correlated medical cost data with several appealing features. First, the mean function is partially linear. Second, the distributional form for the response is not specified. Third, the covariance structure of correlated medical costs has a semiparametric form. We use extended generalized estimating equations to simultaneously estimate all parameters of interest. B-splines are used to estimate unknown functions, and a modification to Akaike information criterion is proposed for selecting knots in spline bases. We apply the model to correlated medical costs in the Medical Expenditure Panel Survey dataset. Simulation studies are conducted to assess the performance of our method.