Abstract:BACKGROUND:Observational studies incorporating real-world data from multiple institutions facilitate study of rare outcomes or exposures and improve generalizability of results. Due to privacy concerns surrounding patient-level data sharing across institutions, methods for performing regression analyses distributively are desirable. Meta-analysis of institution-specific estimates is commonly used, but has been shown to produce biased estimates in certain settings. While distributed regression methods are increasingly available, methods for analyzing count outcomes are currently limited. Count data in practice are commonly subject to overdispersion, exhibiting greater variability than expected under a given statistical model.OBJECTIVE:We propose a novel computational method, a one-shot distributed algorithm for quasi-Poisson regression (ODAP), to distributively model count outcomes while accounting for overdispersion.METHODS:ODAP incorporates a surrogate likelihood approach to perform distributed quasi-Poisson regression without requiring patient-level data sharing, only requiring sharing of aggregate data from each participating institution. ODAP requires at most three rounds of non-iterative communication among institutions to generate coefficient estimates and corresponding standard errors. In simulations, we evaluate ODAP under several data scenarios possible in multi-site analyses, comparing ODAP and meta-analysis estimates in terms of error relative to pooled regression estimates, considered the gold standard. In a proof-of-concept real-world data analysis, we similarly compare ODAP and meta-analysis in terms of relative error to pooled estimatation using data from the OneFlorida Clinical Research Consortium, modeling length of stay in COVID-19 patients as a function of various patient characteristics. In a second proof-of-concept analysis, using the same outcome and covariates, we incorporate data from the UnitedHealth Group Clinical Discovery Database together with the OneFlorida data in a distributed analysis to compare estimates produced by ODAP and meta-analysis.RESULTS:In simulations, ODAP exhibited negligible error relative to pooled regression estimates across all settings explored. Meta-analysis estimates, while largely unbiased, were increasingly variable as heterogeneity in the outcome increased across institutions. When baseline expected count was 0.2, relative error for meta-analysis was above 5% in 25% of iterations (250/1000), while the largest relative error for ODAP in any iteration was 3.59%. In our proof-of-concept analysis using only OneFlorida data, ODAP estimates were closer to pooled regression estimates than those produced by meta-analysis for all 15 covariates. In our distributed analysis incorporating data from both OneFlorida and the UnitedHealth Group Clinical Discovery Database, ODAP and meta-analysis estimates were largely similar, while some differences in estimates (as large as 13.8%) could be indicative of bias in meta-analytic estimates.CONCLUSIONS:ODAP performs privacy-preserving, communication-efficient distributed quasi-Poisson regression to analyze count outcomes using data stored within multiple institutions. Our method produces estimates nearly matching pooled regression estimates and sometimes more accurate than meta-analysis estimates, most notably in settings with relatively low counts and high outcome heterogeneity across institutions.

Distributed Estimation for Large-Scale Cox Regression with Poisson Subsampling

Communication-Efficient Distributed Estimation and Inference for Cox's Model

Learning from Local to Global - an Efficient Distributed Algorithm for Modeling Time-to-event Data

Distributed Bootstrap Simultaneous Inference for High-Dimensional Quantile Regression

Optimal subsampling for the Cox proportional hazards model with massive survival data

Approximating Partial Likelihood Estimators via Optimal Subsampling

An optimal subsampling design for large-scale Cox model with censored data

Optimal Cox Regression Subsampling Procedure with Rare Events

Distributed Cox proportional hazards regression using summary-level information.

Multicenter Privacy-Preserving Cox Analysis Based on Homomorphic Encryption

Privacy enhanced collaborative inference in the Cox proportional hazards model for distributed data

Optimal Distributed Subsampling for Maximum Quasi-Likelihood Estimators with Massive Data

Combining Cox Regressions Across a Heterogeneous Distributed Research Network Facing Small and Zero Counts

Communication-efficient Distributed Statistical Inference on Zero-Inflated Poisson Models

Distributed Quasi-Poisson Regression Algorithm for Modeling Multi-Site Count Outcomes in Distributed Data Networks

Bayesian Cox Regression for Large-scale Inference with Applications to Electronic Health Records

Sampling-Based Estimation For Massive Survival Data With Additive Hazards Model

On the Maximum Penalized Full Likelihood Approach for Cox Model with Extreme Value for Heavily Censored Survival Data

A uniformly ergodic Gibbs sampler for Bayesian survival analysis

VERTICOX: Vertically Distributed Cox Proportional Hazards Model Using the Alternating Direction Method of Multipliers

Distributed Subsampling for Multiplicative Regression