Bias-corrected estimates for logistic regression models for complex surveys with application to the United States’ Nationwide Inpatient Sample

Abstract

For complex surveys with a binary outcome, logistic regression is widely used to model the outcome as a function of covariates. Complex survey sampling designs are typically stratified cluster samples, but consistent and asymptotically unbiased estimates of the logistic regression parameters can be obtained using weighted estimating equations (WEEs) under the naive assumption that subjects within a cluster are independent. Despite the relatively large samples typical of many complex surveys, with rare outcomes, many interaction terms, or analysis of subgroups, the logistic regression parameters estimates from WEE can be markedly biased, just as with independent samples. In this paper, we propose bias-corrected WEEs for complex survey data. The proposed method is motivated by a study of postoperative complications in laparoscopic cystectomy, using data from the 2009 United States’ Nationwide Inpatient Sample complex survey of hospitals.

Keywords

Binary responses bladder cancer population survey stratified cluster sampling weighted estimating equations

Get full access to this article

View all access options for this article.

References

Johnson CL, Dohrmann SM, Burt VL, et al. National Health and Nutrition Examination Survey: Sample design, 2011-2014. National Center for Health Statistics. Vital and Health Statistics Series 2, Number 168, 2014, http://www.cdc.gov/nchs/data/series/sr_02/sr02_162.pdf.

Albert

Anderson

. On the existence of maximum likelihood estimates in logistic regression models. Biometrika 1984; 71: 1–10.

Box

. Bias in nonlinear estimation. J R Stat Soc Series B (Methodol) 1971; 33: 171–201.

Firth

. Bias reduction of maximum likelihood estimates. Biometrika 1993; 80: 27–38.

Kosmidis

Firth

. Bias reduction in exponential family nonlinear models. Biometrika 2009; 96: 793–804.

Binder

. On the variance of asymptotically normal estimators from complex surveys. Int Stat Rev 1983; 51: 279–292.

Liang

Zeger

. Longitudinal data analysis using generalized linear models. Biometrika 1986; 73: 13–22.

Hevelone

Lipsitz

et al.

Comparative analysis of outcomes and costs following open radical cystectomy versus robot-assisted laparoscopic radical cystectomy: results from the US Nationwide Inpatient Sample. Eur Urol 2012; 61: 1239–1244.

HCUP, NIS. Database documentation. Healthcare Cost and Utilization Project (HCUP). SID Database Documentation. Rockville, MD: Agency for Healthcare Research and Quality, 2007.

10.

Kish

. Survey sampling, New York, NY: John Wiley, 1965.

11.

Shah BV, Barnwell BG and Bieler GS. SUDAAN, Software for the Statistical Analysis of Correlated Data: User’s Manual, Release 7.0. Research Triangle Park, NC: Research Triangle Institute, 1996.

12.

Heinze

Schemper

. Comparing the importance of prognostic factors in Cox and logistic regression using SAS. Comput Method Prog Biomed 2003; 71: 155–163.

13.

Wedderburn

RWM

. On the existence and uniqueness of the maximum likelihood estimates for certain generalized linear models. Biometrika 1976; 63: 27–32.

14.

Morel

Bokossa

Neerchal

. Small sample correction for the variance of GEE estimators. Biometric J 2003; 45: 395–409.

15.

Clogg

Rubin

Schenker

et al.

Multiple imputation of industry and occupation codes in census public-use samples using Bayesian logistic regression. J Am Statist Assoc 1991; 86: 68–78.

16.

Wang

Louis

. Matching conditional and marginal shapes in binary random intercept models using a bridge distribution function. Biometrika 2003; 90: 765–775.

17.

Parzen

et al.

A generalized linear mixed model for longitudinal binary data with a marginal logit link function. Ann Appl Stat 2011; 5: 449–467.

18.

Hauck

Donner

. Wald’s test as applied to hypotheses in logit analysis. J Am Statist Assoc 1977; 72: 851–853.

19.

Bull

Lewinger

Lee

. Confidence intervals for multinomial logistic regression in sparse data. Statist Med 2007; 26: 903–918.

20.

Agresti

Ryu

. Pseudo-score inference for parameters in discrete statistical models. Biometrika 2010; 97: 215–222.

21.

Little

. Calibrated Bayes, an alternative inferential paradigm for official statistics. J Off Statist 2012; 28: 309–334.

22.

Kosmidis I and Firth D. A generic algorithm for reducing bias in parametric estimation. Electron J Stat 2010; 4: 1097–1112.

23.

Kosmidis

Firth

. Multinomial logit bias reduction via the Poisson log-linear model. Biometrika 2011; 98: 755–759.