1、江苏省昆山市人口统计学指标变化趋势及预测分析 胡文斌 秦威 史建国 仝岚 邱和泉 周杰 金亦徐 张婷 罗晓明 昆山市疾病预防控制中心慢性非传染性疾病预防控制科 昆山高新区江浦社区服务中心 摘 要： 目的 分析江苏省昆山市 19812014 年户籍人口变动情况以及时间序列特点。方法 历年年均人口数、出生率、死亡率、人口自然增长率及机械增长率资料来源于 19812014 年昆山市统计年鉴。使用统计描述、时间趋势与自回归模型结合的方法对出生率、死亡率进行时间序列建模, 其模型的识别、估计以及预测用 SAS 软件实现。结果 江苏省昆山市 19812014 年男女合计平均观察人年数由 1981 年的 5
2、27 337 人增加到 2014 年的 761 336 人, 上升趋势差异有统计学意义 (APC=1.1%, 95%CI:1.0% 1.2%) 。2000 年以后粗出生率明显上升, 粗死亡率在 19812014 年期间明显下降 (APC=-0.3%, 95%CI:-0.5% -0.1%) 。时间趋势与自回归模型对出生率和死亡率进行预测, 结果显示, 未来年份人口出生率急剧上升, 而死亡率则缓慢下降。结论 江苏省昆山市年均人口数在人口自然增长和机械增长双重作用下逐年上升, 同时也呈现出粗死亡率下降和老龄化进程加速的趋势。关键词： 出生率; 死亡率; 人口自然增长率; 年均人口数; 时间序列; 收
3、稿日期：2016-11-16基金：昆山市社会发展科技专项 (KS1655) Demographic Trends and Prediction in Kunshan City, Jiangsu ProvinceHU Wenbin ZHANG Ting LUO Xiaoming Department of Chronic Non-communicable Disease Preventive and Control, Kunshan Center for Disease Control and Prevention; Abstract： Objectives To explore the dem
4、ographic trend and the time serie characteristic in Kunshan city, Jiangsu province, 1981-2014.Methods Average annual population, birth rate, death rate, population mechanical growth rate, and natural population growth rate were derived from Kunshan Statistical Yearbook.Statistical description and st
5、epwise autoregressive method were used to modeling the demographic change indicators, and the identify, estimate, and prediction were conducted in SAS software.Results Annual average population was 527 337 in 1981 and 761 336 in2014 in Kunshan city, the upward trend was statistically significant (AP
6、C=1.1%, 95%CI:1.0% 1.2%) .Crude birth rate was increasing gradually after the year of 2000, and crude death rate was decreased significantly during1981 to 2014 (APC=-0.3%, 95%CI:-0.5% -0.1%) .The stepwise autoregressive method of time series model demonstrated that the birth rate showed upward trend
7、 in the following years, but the death rate showed downward trend.Conclusions The average population in Kunshan city, Jiangsu province, have been driven upwards by the natural population growth and mechanical growth.Moreover, the decline of death rate and the acceleration of the aging process were c
8、oexist.Keyword： Birth rate; Death rate; Natural population growth rate; Mechanical growth rate; Time series; Received： 2016-11-16人口变动是社会学与经济学关注的课题, 可影响到社会各个方面1。人口迁移对居民总体健康指标影响不容忽视2-3。人口机械增长升高可提高平均期望寿命、降低全死因死亡率、明显改变年龄构成等。因此, 新时期深入研究人口变动是客观看待居民健康指标的关键。江苏省昆山市地处上海和苏州之间, 特殊的区位优势使得社会经济快速发展、人口持续变动4-5。既往研究提
9、示, 人口老龄化对昆山市人群死因影响较大6-9, 但户籍人口变动情况对其影响如何、年龄结构如何、人口发展趋势如何, 目前尚未有系统研究结果。为此, 本研究采用时间序列分析方法, 分析昆山市 19812014 年户籍人口变动情况, 不仅可对同类研究结果提供相互比较依据, 也对今后与户籍人口相关的疾病监测具有重要参考意义。1 资料与方法1.1 资料来源江苏省昆山市 19812014 年年均人口数、出生数、出生率、死亡数、死亡率、自然增长人数、人口自然增长率、迁入人口和迁出人口、人口机械增长率等数据来源于昆山市统计局历年发布的统计年鉴10。户籍人口自然增长率是指在一年内户籍人口自然增加数 (出生人数
10、减死亡人数) 与该时期内平均人口数 (或期中人口数) 之比, 用千分率表示。人口机械增长率= (迁入人口数-迁出人口数) /该时期内平均人口数 (或期中人口数) , 用千分率表示。19812010 年人口年龄构成来自于昆山市疾病预防控制中心每年 10%人口抽样调查;20112014 年人口年龄结构来源于昆山市公安局。人口老龄化指标为60 岁的比例与65 岁的比例。1.2 统计学分析对昆山市 19812014 年出生率、死亡率及年均人口数的时间序列预测模型拟合中, 使用时间趋势与自回归模型结合的方法 (stepwise autoregressive method) 11-12, 其一般的数学模型
11、如下:x 1=b0+b1t+b2t=ut, 其中 ut=a1ut-1+a2ut-2+aput-p= t。模型评价指标采用决定系数 (R) 、施瓦兹-贝叶斯判别准则 (SBC) 、赤池信息准则 (AIC) 及相关系数 (CORR) 等。该模型中, 趋势模型适合于获取长期行为, 而自回归模型则更加适合获取短期波动11-12。其模型的识别、估计, 以及预测用 SAS 软件给出的模块“FORECAST procedure”实现13。采用对数线性模型估计年度变化百分比 (annual percentage change, APC) 及其 95%CI, 来评价出生率、死亡率及年均人口数的时间趋势14。采用
12、 SAS 9.3 进行 APC 检验、时间序列建模及绘图。2 结果2.1 昆山市 19812014 年人口统计学主要指标变化趋势江苏省昆山市 19812014 年男女合计平均观察人年数由 1981 年的 527 337 人年增加到 2014 年的 761 336 人年, 上升趋势差异有统计学意义 (APC=1.1%, 95%CI:1.0%1.2%) 。男性人群中平均观察人年数由 1981 年的 265 251 人年增加到 2014 年的 378 477 人年, 上升趋势差异有统计学意义 (APC=1.0%, 95%CI:0.9%1.1%) ;女性人群中由 1981 年的 262 086 人年增
13、加到 2014 年的382 859 人年, 上升趋势差异有统计学意义 (APC=1.1%, 95%CI:1.0%1.2%) 。见表 1。人口老龄化指标中, 60 岁人口的比例由 1981 年的 8.2%上升到 2014 年的20.6%, 上升趋势明显 (APC=2.6%, 95%CI:2.3%2.8%) ;65 岁人口的比例由1981 年的 5.1%上升到 2014 年的 14.0% (APC=3.0%, 95%CI:2.6%3.4%) 。见图1。粗出生率在 19812014 年期间, 呈现先下降后上升的变化。在 2001 年出生率下降至最低点, 为 6.25, 之后的年份逐渐上升, 至 20
14、14 年粗出生率为13.09。死亡数则在特定水平略有上升, 上升趋势差异有统计学意义 (APC=0.8%, 95%CI:0.7%0.9%) ;而粗死亡率在 19812014 年期间明显下降, 趋势差异有统计学意义 (APC=-0.3%, 95%CI:-0.5%-0.1%) 。江苏省昆山市人口变动情况中, 从 1990 年开始人口机械增长逐年增加, 2000 年以后随着迁入人口的大幅增加, 人口机械增长也急剧上升, 至 2014 年人口机械增长率为14.85, 机械增长人数为达 11 306 人。见表 1。2.2 江苏省昆山市 19812014 年出生率和死亡率时间序列分析江苏省昆山市出生率时间
15、趋势与自回归模型建模中, R=0.92、AIC=-13.02、SBC=-9.24、CORR=0.96, 表明模型拟合较好。预测结果显示, 20152024 年昆山市出生率分别为13.57、14.78、15.77、16.82、17.54、18.45、19.41、20.59、21.75和 22.95。见图 2。死亡率时间趋势与自回归模型建模中, R=0.60、AIC=-101.07、SBC=-96.49、CORR=0.77, 表示模型拟合尚可。预测结果显示, 20152024 年昆山市死亡率分别为5.35、5.25、5.14、5.03、4.92、4.80、4.67、4.54、4.41和 4.27。
16、见图 2。3 讨论随着社会经济发展及城市化进程, 人口变动持续动态发展1。分析过去 34 年户籍人口变动情况, 为昆山市户籍人口相关的疾病监测提供人口学特征变化依据, 也为深入研究健康期望寿命、科学评价昆山市居民疾病负担奠定基础。本研究显示, 江苏省昆山市人口出生率急剧上升、死亡率缓慢下降, 人口自然增长率持续上升;人口机械增长率也在高水平波动, 导致年均人口数逐年上升, 同时伴随着老龄化进程加速。本研究显示, 人口出生率在 19812000 年明显下降, 而 2000 年之后则明显上升。可能与国家人口宏观政策有关, 与赖平耀等1研究结果一致。人口死亡率相对平稳, 但总体上呈现了平稳偏下降态势
17、。分开来看, 人口死亡数呈现上升趋势 (APC=0.8%, 95%CI:0.7%0.9%) , 由于总人口数上升的速率远高于死亡数, 而导致死亡率明显下降 (APC=-0.3%, 95%CI:-0.5%-0.1%) , 也可能在未来年份死亡率会进一步下降。在国家全面二胎政策的背景下, 有理由相信未来年份出生率会明显上升;这也与本文时间序列预测结果一致。李辉等15对宁波市人口变动情况分析显示, 20002010 年期间, 宁波市人口数量急剧上升, 流动人口、老年人口、014 岁儿童的比例在 20002010 年期间发生了巨大变化, 这些特殊人群在基本公共服务等方面亟需加强。总体来看, 江苏省昆山
18、市年均人口数上升速率在 2003 年左右出现拐点, 考虑可能与 20032008 年人口机械增长高峰有关。但同时也注意到, 死亡率逐年下降、人口老龄化进程在加速。人口老龄化程度一直是社会各个部门关注的重要指标, 除了对经济产生影响2, 人口老龄化与疾病也息息相关。表 1 江苏省昆山市 19812014 年人口统计学指标变化趋势分析 下载原表 图 1 江苏省昆山市 19812014 年户籍人口年龄构成趋势图 下载原图图 2 江苏省昆山市出生率和死亡率时间序列预测 下载原图本研究显示, 昆山市历年60 岁人口的比例由 1981 年的 8.2%上升到 2014 年的20.6%。在昆山市平均期望寿命居
19、历史新高的背景下16, 人口学特征对糖尿病相关的死亡6和对肿瘤死亡7,9的影响不容忽视。既往研究提示, 昆山市高血压17-18和糖尿病19等慢性病患病率随着年龄增加而逐渐上升;肿瘤粗死亡率明显上升20, 是影响昆山居民寿命的第二位疾病21, 表现尤为突出的是肺癌发病与死亡率处于明显上升趋势9,22-23, 肺癌和胃癌等依旧是造成巨大疾病负担的主要癌种24-25。本研究依据时间序列分析方法对昆山市出生率、死亡率等进行分析, 旨在探索过去昆山市人口变动情况及将来的发展。时间趋势与自回归模型结合的方法, 在 SAS 9.3 软件给出的模块“FORECAST procedure”可以实现13, 在较早
20、的研究中也鲜有应用11-12。本研究中对出生率和死亡率时间序列分析中, 使用了非线性 (3 次方) 趋势, 较为稳健地拟合了时间序列预测模型。本研究初步呈现了昆山市 19812014 年户籍人口变动情况, 但研究中也有几点缺陷:一是数据来源的局限性, 不能对死亡率和出生率进行年龄标化。二是数据来源的差异, 在一定程度上数据同质性较差, 难以精确呈现户籍人口变动。如出生数与出生率、死亡数与死亡率、迁入和迁出人口数来源于历年昆山市统计年鉴;而人口年龄构成来源于疾控中心抽样调查。三是在说明老龄化现象时昆山市 19812010 人口年龄构成可能与实际人口年龄结构存在差距。参考文献1赖平耀, 钟甫宁.中
21、国的人口变动:1978-2012J.南京农业大学学报:社会科学版, 2014, 14 (1) :99-110.1赖平耀, 钟甫宁.中国的人口变动:1978-2012J.南京农业大学学报:社会科学版, 2014, 14 (1) :99-110. 2王美凤.我国老年人口年龄结构变动及其对经济增长的影响J.人口与社会, 2015 (3) :19-26.2王美凤.我国老年人口年龄结构变动及其对经济增长的影响J.人口与社会, 2015 (3) :19-26. 3袁小平, 梁海艳.中国人口年龄结构变动对出生率的影响研究J.西北人口, 2014, 35 (6) :49-53.3袁小平, 梁海艳.中国人口年龄
22、结构变动对出生率的影响研究J.西北人口, 2014, 35 (6) :49-53. 4滕琪.改革开放以来昆山市人口增长、分布及迁移研究D.南京:南京师范大学, 2013.4滕琪.改革开放以来昆山市人口增长、分布及迁移研究D.南京:南京师范大学, 2013. 5姚鑫, 杨桂山, 孙洪波, 等.昆山市人口指标的数学模型J.长江流域资源与环境, 2010, 19 (3) :244-248.5姚鑫, 杨桂山, 孙洪波, 等.昆山市人口指标的数学模型J.长江流域资源与环境, 2010, 19 (3) :244-248. 6沙汝明.江苏省昆山市人口老龄化及危险因素改变对糖尿病死亡率的影响J.中华流行病学杂
23、志, 2005, 26 (10) :810-810.6沙汝明.江苏省昆山市人口老龄化及危险因素改变对糖尿病死亡率的影响J.中华流行病学杂志, 2005, 26 (10) :810-810. 7沙汝明.昆山市人口老龄化及危险因素改变对恶性肿瘤死亡率的影响研究J.现代预防医学, 2005, 32 (10) :1343-1344.7沙汝明.昆山市人口老龄化及危险因素改变对恶性肿瘤死亡率的影响研究J.现代预防医学, 2005, 32 (10) :1343-1344. 8张瑾, 沙汝明.昆山市人口老龄化及危险因素改变对前 10 位死因及其死亡率的影响研究J.安徽预防医学杂志, 2003, 9 (4) :
24、217-219.8张瑾, 沙汝明.昆山市人口老龄化及危险因素改变对前 10 位死因及其死亡率的影响研究J.安徽预防医学杂志, 2003, 9 (4) :217-219. 9胡文斌, 张婷, 秦威, 等.江苏省昆山市 19932014 年肺癌死亡趋势及差别分解分析J.中国公共卫生, 2015, 31 (11) :1460-1463.9胡文斌, 张婷, 秦威, 等.江苏省昆山市 19932014 年肺癌死亡趋势及差别分解分析J.中国公共卫生, 2015, 31 (11) :1460-1463. 10昆山市统计局.昆山统计年鉴M/OL.2016-01-16.http:/ 11Barry F.Fitt
25、ing autoregressive trend stationary models with finite samplesJ.Int J Forecast, 1999, 15 (1) :11-25.11Barry F.Fitting autoregressive trend stationary models with finite samplesJ.Int J Forecast, 1999, 15 (1) :11-25. 12Pflaumer P.Forecasting U.S.population totals with the Box-Jenkins approachJ.Int J F
26、orecast, 1992, 8 (3) :329-338.12Pflaumer P.Forecasting U.S.population totals with the Box-Jenkins approachJ.Int J Forecast, 1992, 8 (3) :329-338. 13Brocklebank JC, Dickey DA.SAS for Forecasting Time Series, Second EditionM.Cary NC:SAS Institute Inc, USA, 2003.13Brocklebank JC, Dickey DA.SAS for Fore
27、casting Time Series, Second EditionM.Cary NC:SAS Institute Inc, USA, 2003. 14Fay MP, Tiwari RC, Feuer EJ, et al.Estimating average annual percent change for disease rates without assuming constant changeJ.Biometrics, 2006, 62 (3) :847-854.14Fay MP, Tiwari RC, Feuer EJ, et al.Estimating average annua
28、l percent change for disease rates without assuming constant changeJ.Biometrics, 2006, 62 (3) :847-854. 15李辉, 崔军, 王永, 等.2000 2010 年宁波市人口发展趋势分析J.中国卫生统计, 2014, 31 (5) :876-877, 881.15李辉, 崔军, 王永, 等.20002010 年宁波市人口发展趋势分析J.中国卫生统计, 2014, 31 (5) :876-877, 881. 16胡文斌, 张婷, 秦威, 等.江苏省昆山市 19812014 年期望寿命时间序列分析J.
29、中国全科医学, 2016, 19 (8) :971-975.16胡文斌, 张婷, 秦威, 等.江苏省昆山市 19812014 年期望寿命时间序列分析J.中国全科医学, 2016, 19 (8) :971-975. 17胡文斌, 张婷, 史建国, 等.BMI 与高血压关联强度剂量-反应关系分析J.中国卫生统计, 2015, 32 (6) :971-974.17胡文斌, 张婷, 史建国, 等.BMI 与高血压关联强度剂量-反应关系分析J.中国卫生统计, 2015, 32 (6) :971-974. 18胡文斌, 张婷, 史建国, 等.男性吸烟与高血压病的剂量-反应关系J.中华心血管病杂志, 201
30、4, 42 (9) :773-777.18胡文斌, 张婷, 史建国, 等.男性吸烟与高血压病的剂量-反应关系J.中华心血管病杂志, 2014, 42 (9) :773-777. 19胡文斌, 张婷, 秦威, 等.昆山市 2012 年 2 型糖尿病患病现况及其危险因素J.江苏预防医学, 2015, 26 (3) :87-89.19胡文斌, 张婷, 秦威, 等.昆山市 2012 年 2 型糖尿病患病现况及其危险因素J.江苏预防医学, 2015, 26 (3) :87-89. 20胡文斌, 张婷, 秦威, 等.江苏省昆山市 1981-2013 年恶性肿瘤死亡趋势分析J.中国肿瘤, 2015, 24
31、(9) :752-757.20胡文斌, 张婷, 秦威, 等.江苏省昆山市 1981-2013 年恶性肿瘤死亡趋势分析J.中国肿瘤, 2015, 24 (9) :752-757. 21朱士新, 郑建中, 陆秋明, 等.昆山市居民 19821999 年潜在寿命损失年数分析J.疾病控制杂志, 2001, 5 (2) :121-122.21朱士新, 郑建中, 陆秋明, 等.昆山市居民 19821999 年潜在寿命损失年数分析J.疾病控制杂志, 2001, 5 (2) :121-122. 22秦威, 许寒冰.1981-2005 年昆山市居民肺癌死亡率趋势J.中国慢性病预防与控制, 2009, 17 (1
32、) :37-38.22秦威, 许寒冰.1981-2005 年昆山市居民肺癌死亡率趋势J.中国慢性病预防与控制, 2009, 17 (1) :37-38. 23胡文斌, 张婷, 史建国, 等.江苏省昆山市 20062012 年肺癌发病趋势分析J.中国肿瘤, 2014, 23 (10) :810-814.23胡文斌, 张婷, 史建国, 等.江苏省昆山市 20062012 年肺癌发病趋势分析J.中国肿瘤, 2014, 23 (10) :810-814. 24胡文斌, 张婷, 史建国, 等.江苏省昆山市 2010 年肿瘤死亡率与疾病负担分析J.江苏预防医学, 2015, 26 (5) :42-44.24胡文斌, 张婷, 史建国, 等.江苏省昆山市 2010 年肿瘤死亡率与疾病负担分析J.江苏预防医学, 2015, 26 (5) :42-44. 25胡文斌, 张婷, 秦威, 等.江苏省昆山市 2011 年恶性肿瘤死亡率与疾病负担分析J.中国肿瘤, 2015, 24 (11) :905-909.25胡文斌, 张婷, 秦威, 等.江苏省昆山市 2011 年恶性肿瘤死亡率与疾病负担分析J.中国肿瘤, 2015, 24 (11) :905-909.
