

FOLLOWUS
1. 1Department of Epidemiology and Statistics, Institute of Basic Medical Sciences, Chinese Academy of Medical Sciences & School of Basic Medicine, Peking Union Medical College, Beijing 100005, China;
2. 2Department of physiopathology, Institute of Basic Medical Sciences, Chinese Academy of Medical Sciences & School of Basic Medicine, Peking Union Medical College, Beijing 100005, China;
*Corresponding author Tel: 86-10-69156408, E-mail: hansm1@vip.sina.com
Received:30 March 2017,
Published Online:30 December 2017,
Published:2017-12
Scan QR Code
Xu Tao, Zhu Guangjin, Han Shaomei. Study of Zero-Inflated Regression Models in a Large-Scale Population Survey of Sub-Health Status and Its Influencing Factors△[J]. Chinese medical sciences journal, 2017, 32(4): 218-225.
Xu Tao, Zhu Guangjin, Han Shaomei. Study of Zero-Inflated Regression Models in a Large-Scale Population Survey of Sub-Health Status and Its Influencing Factors△[J]. Chinese medical sciences journal, 2017, 32(4): 218-225. DOI: 10.24920/J1001-9294.2017.054.
目的
2
亚健康近年来引起了越来越多的医务人员和普通公众的注意。把亚健康症状数作为计数资料进行处理能够获得更完整和准确的分析结果。本研究旨在探讨计数资料模型的拟合优度以确定最佳的亚健康研究模型。
方法
2
研究样本来源于一项2007年到2014年在中国四省两自治区进行的一项关于生理和心理的大规模人群调查研究。用SAS软件拟合四个计数资料模型:Poisson回归、负二项回归(NB)、零频数过多的Poisson回归(ZIP)和零频数过多的Poisson回归负二项回归模型(ZINB)。亚健康症状数是主要结局变量。用α系数和
O
检验验证是否过度离散
用Vuong检验验证零频数是否过多。用预测概率曲线和似然比检验评价模型的拟合优度。
结果
2
在78
307名受试者中
38.53%的受试者没有任何亚健康症状。亚健康症状数的平均水平为2.98±3.72。
O
检验统计量为720.995 (P
<
0.001)。α系数为0.618 (95% 可信区间: 0.600-0.636)。Vuong检验的统计量
Z
=45.487。ZINB模型有最大的似然比(-167519)、最小的赤池信息准则系数(335112) 和贝叶斯信息准则系数 (335455)。ZINB模型对大多数症状数的预测概率都优于其他模型。ZINB模型的logit部分显示年龄、性别、职业、吸烟、饮酒、民族、肥胖是发生亚健康症状的危险因素。负二项部分显示性别、职业、吸烟、饮酒、民族、婚姻状况和肥胖对于亚健康的程度有影响。
结论
2
各种拟合优度检验和预测概率曲线的综合结果表明ZINB模型是探索亚健康症状数的影响因素的最佳模型。
Objective
2
Sub-health status has progressively gained more attention from both medical professionals and the publics. Treating the number of sub-health symptoms as count data rather than dichotomous data helps to completely and accurately analyze findings in sub-healthy population. This study aims to compare the goodness of fit for count outcome models to identify the optimum model for sub-health study.
Methods
2
The sample of the study derived from a large-scale population survey on physiological and psychological constants from 2007 to 2011 in 4 provinces and 2 autonomous regions in China. We constructed four count outcome models using SAS: Poisson model
negative binomial (NB) model
zero-inflated Poisson (ZIP) model and zero-inflated negative binomial (ZINB) model. The number of sub-health symptoms was used as the main outcome measure. The alpha dispersion parameter and
O
test were used to identify over-dispersed data
and Vuong test was used to evaluate the excessive zero count. The goodness of fit of regression models were determined by predictive probability curves and statistics of likelihood ratio test.
Results
2
Of all 78 307 respondents
38.53% reported no sub-health symptoms. The mean number of sub-health symptoms was 2.98
and the standard deviation was 3.72. The statistic
O
in over-dispersion test was 720.995 (
P
<
0.001); the estimated alpha was 0.618 (95%
CI
: 0.600-0.636) comparing ZINB model and ZIP model; Vuong test statistic Z was 45.487. These results indicated over-dispersion of the data and excessive zero counts in this sub-health study. ZINB model had the largest log likelihood (-167 519)
the smallest Akaike’s Information Criterion coefficient (335 112) and the smallest Bayesian information criterion coefficient (335455)
indicating its best goodness of fit. The predictive probabilities for most counts in ZINB model fitted the observed counts best. The logit section of ZINB model analysis showed that age
sex
occupation
smoking
alcohol drinking
ethnicity and obesity were determinants for presence of sub-health symptoms; the binomial negative section of ZINB model analysis showed that sex
occupation
smoking
alcohol drinking
ethnicity
marital status and obesity had significant effect on the severity of sub-health.
Conclusions
2
All tests for goodness of fit and the predictive probability curve produced the same finding that ZINB model was the optimum model for exploring the influencing factors of sub-health symptoms.
Medicine CAOC. The TCM clinical guidelines of suboptimal health status. Beijing : China Press Traditional Chinese Medicine , 2006 .
Wang W , Yan Y. Suboptimal health: a new health dimension for translational medicine. Clin Transl Med 2012 ; 1 : 28 . doi: 10.1186/2001-1326-1-28 https://doi.org/10.1186/2001-1326-1-28 .
Wang DX , Zhou HL . Health, disease and subhealth . Medicine and Society 2007 ; 20 : 5 - 8 . Chinese.
Chen QS , Wang SY , Jing CX , Dong XM , Chi GB , Zhu L. Evaluation on diagnostic criterion of sub-health with Delphi method . Chinese Journal of Public Health . 2003 ; 19 : 1467 - 8 . Chinese. doi: 10.3321/j.issn:1001-0580.2003. 12.030 https://doi.org/10.3321/j.issn:1001-0580.2003.12.030 .
Xu T , Han SM , Liu JT , Zhu GJ , Mao M. Comparison of sub-health status between Tibetan people and Han people Chin Med J 2009 ; 89 : 2671 - 4 . Chinese. doi: 10.3760/cma.j.issn.0376-2491.2009.38.003 https://doi.org/10.3760/cma.j.issn.0376-2491.2009.38.003 .
Xu T , Liu JT , Han SM , Zhu GJ , Mao M. Analyses for risk factors of sub-health status with logistic model and binomial model-survey in Tibetan people. J Clin Rehabilitative Tis Engineer Res 2009 ; 13 : 6597 - 600 . Chinese. doi: 10.3969/j.issn.1673-8225.2009.33.043 https://doi.org/10.3969/j.issn.1673-8225.2009.33.043 .
Zhu L , Wang SY , Fan CX , Xiao YJ , Ou CF. Logistic regression analysis of risk factors on subhealth of young teachers in institution of higher learning. Chin J Public Health . 2003 ; 19 : 595 - 6 . Chinese. doi: 10.3321/j.issn: 1001-0580.2003.05.048 https://doi.org/10.3321/j.issn:1001-0580.2003.05.048 .
Wang Y , Chen Q. A cross-sectional study on sub-health of students in a university in Guangzhou . Chin General Practice 2005 ; 8 : 738 - 40 . Chinese. doi: 10.3969/j.issn.1007-9572.2005.09.019 https://doi.org/10.3969/j.issn.1007-9572.2005.09.019 .
Pang H , Wang P. Research on sub-health condition and influencing factors of middle school teacher in Xinjiang . China Sport Science and Technology 2008 ; 44 : 47 - 50 . Chinese. doi: 10.3969/j.issn.1002-9826.2008.05.008 https://doi.org/10.3969/j.issn.1002-9826.2008.05.008 .
Nie XL , Xue Q , Lai MH , Chen J , Zhao XS , Luo R . Subhealth and its influential factors among civil servants in taxation department . Chinese Journal of Public Health 2010 ; 26 : 634 - 5 . Chinese.
Hall DB. Zero-inflated Poisson and binomial regression with random effects: a case study. Biometric 2000 ; 56 : 1030 - 9 . doi: 10.1111/j.0006-341X.2000.01030.x https://doi.org/10.1111/j.0006-341X.2000.01030.x .
Agresti A. Generalized Linear Models. In: Agresti A, editors. An introduction to categorical data analysis. New York: Chichester Wiley, 1996 . p 74 - 83 .
Cheung YB. Zero-inflated models for regression analysis of count data: a study of growth and development. Stat Med 2002 ; 21 : 1461 - 9 . doi: 10.1002/sim.1088 https://doi.org/10.1002/sim.1088 .
Solinas G , Campus G , Maida C , Sotgiu G , Cagetti MG , Lesaffre E. et al. What statistical method should be used to evaluate risk factors associated with dmfs index? Evidence from the National Pathfinder Survey of 4-year-old Italian children. Community Dent Oral Epidemiol 2009 ; 37 : 539 - 46 . doi: 10.1111/j.1600-0528.2009.00500.x https://doi.org/10.1111/j.1600-0528.2009.00500.x .
Kipnis V , Midthune D , Buckman DW , Dodd KW , Guenther PM , Krebs-Smith SM , et al. Modeling data with excess zeros and measurement error: application to evaluating relationships between episodically consumed foods and health outcomes. Biometrics 2009 ; 65 : 1003 - 10 . doi: 10. 1111/j.1541-0420.2009.01223.x https://doi.org/10.1111/j.1541-0420.2009.01223.x .
Wang H , Heitjan DF. Modeling heaping in self-reported cigarette counts. Stat Med 2008 ; 27 : 3789 - 804 . doi: 10.1002/sim.3281 https://doi.org/10.1002/sim.3281 .
Sheu ML , Hu TW , Keeler TE , Ong M , Sung H Y. The effect of a major cigarette price change on smoking behavior in California: a zero-inflated negative binomial model. Health Economics 2004 ; 13 : 781 - 91 . doi: 10.1002/hec.849 https:/10.1002/doi.org//hec.849 .
Denwood MJ , Stear MJ , Matthews L , Reid SW , Toft N , Innocent GT. The distribution of the pathogenic nematode Nematodirus battus in lambs is zero-inflated. Parasitology 2008 ; 135 : 1225 - 35 . doi: 10.1017/S003118-2008004708 https://doi.org/10.1017/S003118-2008004708 .
Lewsey JD , Thomson WM. The utility of the zero-inflated Poisson and zero-inflated negative binomial models: a case study of cross-sectional and longitudinal DMF data examining the effect of socio-economic status. Community Dent Oral Epidemiol 2004 ; 32 : 183 - 9 . doi: 10.1111/j. 1600-0528.2004.00155.x https://doi.org/10.1111/j.1600-0528.2004.00155.x .
Javali SB , Pandit PV. Using zero inflated models to analyze dental caries with many zeroes. Indian J Dent Res 2010 ; 21 : 480 - 5 . doi: 10.4103/0970-9290.74210 https://doi.org/10.4103/0970-9290.74210 .
Lambert D. Zero-Inflated Poisson regression, with an application to defects in manufacturing. Technometrics 1992 ; 34 : 1 - 14 . doi: 10.2307/1269547 https://doi.org/10.2307/1269547 .
Joseph M , Hilbe. Problems with zero counts . In: Joseph M, Hilbe. Negative Binomial Regression. London: Cambridge Univ Pr 2007 : 173 - 77 .
Bohning D , Dietz E , Schlattmann P. Zero-inflated count models and their applications in public health and social science . In: Rost J, Langeheine R eds. Applications of latent trait and latent class models in the Social Sciences. 1st ed. Germany: Munster Waxmann 1997 : 333 - 44 .
Vuong QH. Likelihood ratio tests for model selection and non-nested hypothesis. Econometrica 1989 ; 57 : 307 - 33 . doi: 10.2307/1912557 https://doi.org/10.2307/1912557 .
Rose CE , Martin SW , Wannemuehler KA , Plikaytis BD. On the use of zero-inflated and hurdle models for modeling vaccine adverse event count data. J Biopharm Stat 2006 ; 16 : 463 - 81 . doi: 10.1080/10543400600719384 https://doi.org/10.1080/10543400600719384 .
Akram K , Pedersen-Bjergaard U , Carstensen B , Borch-Johnsen K , Thorsteinsson B. Frequency and risk factors of severe hypoglycaemia in insulin-treated Type 2 diabetes: a cross-sectional survey. Diabet Med 2006 ; 23 : 750 - 6 . doi: 10.1111/j.1464-5491.2006.01880.x https://doi.org/10.1111/j.1464-5491.2006.01880.x .
Zaninotto P , Falaschetti E. Comparison of methods for modelling a count outcome with excess zeros: application to Activities of Daily Living (ADL-s). J Epidemiol Community Health 2011 ; 65 : 205 - 10 . doi: 10.1136/jech. 2008. 079640 https://doi.org/10.1136/jech. 2008. 079640 .
Turner AN , Miller WC , Padian NS , Kaufman JS , Behets FM , Chipato T , et al. Unprotected sex following HIV testing among women in Uganda and Zimbabwe: short- and long-term comparisons with pre-test behaviour. Int J Epidemiol 2009 ; 38 : 997 - 1007 . doi: 10.1093/ije/dyp171 https://doi.org/10.1093/ije/dyp171 .
Carrel M , Voss P , Streatfield PK , Yunus M , Emch M. Protection from annual flooding is correlated with increased cholera prevalence in Bangladesh: a zero-inflated regression analysis. Environ Health 2010 ; 9 : 13 . doi: 10.1186/1476-069X-9-13 https://doi.org/10.1186/1476-069X-9-13 .
Lacruz ME , Emeny RT , Haefner S , Zimmermann AK , Linkohr B , Holle R , et al. Relation between depressed mood, somatic comorbidities and health service utilisation in older adults: results from the KORA-Age study. Age Ageing 2012 ; 41 : 183 - 90 . doi: 10.1093/ageing/afr162 https://doi.org/10.1093/ageing/afr162 .
Meszaros ZS , Dimmock JA , Ploutz-Snyder RJ , Abdul-Malak Y , Leontiera L , Canfield K , et al. Predictors of smoking severity in patients with schizophrenia and alcohol use disorders. Am J Addict 2011 ; 20 : 462 - 7 . doi: 10.1111/ j.1521-0391.2011.00150.x https://doi.org/10.1111/j.1521-0391.2011.00150.x .
Bandiera FC , Arheart KL , Caban-Martinez AJ , Fleming LE , McCollister K, Dietz NA, et al. Secondhand smoke exposure and depressive symptoms. Psychosom Med 2010 ; 72 : 68 - 72 . doi: 10.1097/PSY.0b013e3181c6c8b5 https://doi.org/10.1097/PSY.0b013e3181c6c8b5 .
Slymen DJ , Ayala GX , Arredondo EM , Elder JP. A demonstration of modeling count data with an application to physical activity. Epidemiol Perspect Innov 2006 ; 3 : 3 . doi: 10.1186/1742-5573-3-3 https://doi.org/10.1186/1742-5573-3-3 .
Karazsia BT , van Dulmen MH . Regression models for count data: illustrations using longitudinal predictors of childhood injury. J Pediatr Psychol 2008 ; 33 : 1076 - 84 . doi: 10.1093/jpepsy/jsn055 https://doi.org/10.1093/jpepsy/jsn055 .
Bergemann TL , Huang Z. A new method to account for missing data in case-parent triad studies. Hum Hered 2009 ; 68 : 268 - 77 . doi: 10.1159/000228924 https://doi.org/10.1159/000228924 .
Ceppi M , Biasotti B , Fenech M , Bonassi S. Human population studies with the exfoliated buccal micronucleus assay: statistical and epidemiological issues. Mutat Res 2010 ; 705 : 11 - 9 . doi: 10.1016/j.mrrev.2009.11.001 https://doi.org/10.1016/j.mrrev.2009.11.001 .
Zhang HJ , Min J , Wang P , et al. A statistical method to extra zero in the data of field survey. Acta Universities Medicinalis Nanjing 2007 , 27 : 634 - 636 . doi: 10.3969/j.issn.1007-4368.2007.06.035 https://doi.org/10.3969/j.issn.1007-4368.2007.06.035 .
Marioni RE , Matthews FE , Brayne C. The association between late-life cognitive test scores and retrospective informant interview data. Int Psychogeriatr 2011 ; 23 : 274 - 9 . doi: 10.1017/S1041610210001201 https://doi.org/10.1017/S1041610210001201 .
Publicity Resources
Related Articles
Related Author
Related Institution
京公网安备11010802024621