基于回归分析的人口预测

统计系课程实验论文持。基于回归分析的人口数量预测学号：05姓名：李洋年级：2014级专业：统计学课程：回归分析指导教师：姜喜春完成日期：2016年6月19日iword格式支持编辑，如有帮助欢迎下载支持。持。目录错误!未定义书签。错误!未定义书签。第1章一元线性回归错误!未定义书签。1.1指标的选择错误!未定义书签。1.2样本确定错误!未定义书签。1.3—元回归分析错误!未定义书签。131绘制总人口与粮食产量的散点图错误!未定义书签。132设定理论模型错误!未定义书签。133回归诊断错误!未定义书签。第2章多元线性回归错误!未定义书签。2.1数据中心化标准化错误!未定义书签。2.2多元回归模型建立错误!未定义书签。2.3逐步回归法错误!未定义书签。2.4多重共线性错误!未定义书签。231多重共线性检测错误!未定义书签。2.4主成分分析错误!未定义书签。241主成分分析模型建立错误!未定义书签。第3章非线性模型错误!未定义书签。3.1曲线回归错误!未定义书签。3.1.1曲线拟合错误!未定义书签。3.2Logistic模型错误!未定义书签。结论错误!未定义书签。参考文献错误!未定义书签。iiword格式支持编辑，如有帮助欢迎下载支持。#t—a—前言最早的形式回归的方法是最小二乘法，这是在1805年出版的勒让德，和高斯在1809年。勒让德和高斯都采用的方法确定的问题，从天文观测，有关Sun的机构（主要是彗星，但后来也新发现的小行星）的轨道。1821年，高斯发表最小二乘法理论的进一步发展，在包括高斯-马尔可夫定理的一个版本。弗朗西斯•高尔顿在十九世纪的“回归”是杜撰来描述一种生物现象。这种现象是高度高大的祖先的后代往往倒退下来，对一个正常的平均水平（这种现象也被称为向均值回归）。对高尔顿,，回归只有这个生物意义，Udny圣诞节和皮尔逊但他的工作，后来扩展到更一般的统计范围内。在圣诞节和Pearson，工作的响应和解释变量的联合分布被假定为高斯。这个假设RA费舍尔在1922年和1925年，他的作品被削弱。费舍尔认为的响应变量的条件分布为高斯分布，但联合分布不一定要。在这方面，费舍尔的假设是高斯1821年制定的。在20世纪50年代和20世纪60年代，经济学家旧机电台计算器，计算回归。1970年以前，有时长达24小时接收从一个回归的结果。回归方法继续是一个活跃的研究领域。在最近的几十年中，新的方法已经制定了稳健回归，回归涉及的相关反应，如时间序列曲线和增长曲线，回归的预测或响应变量的曲线，图片，图表或其他复杂的数据对象，容纳不同的回归方法丢失的数据，非参数回归，贝叶斯方法进行回归，回归的预测变量的测量误差，预测变量的观测回归，回归和因果关系的推论与类型。第1章一元线性回归1.1指标的选择影响人口增长的主要因素经济因素，经济因素对人口自然增长的作用主要表现在它决定了人口的增殖条件和生存条件，通过改变人口的出生率和死亡率来影响人口的自然增率。一般情况下，当人口数量不能满足经济发展对劳动力的需求时，人口自身的再生产必将会刺激；当人口数量超越了经济发展所能提供的消费总数后，人口自身的再生产必将受到遏制。在现代生产力水平下，人口的自然增长率往往随着经济水平的提高而下降。经济因素对人口机械增长也有重要影响。通常情况下，经济发达或发展速度较快的地区，对人口具有一种吸引力和凝聚力，人口机械增长为正值；相反，经济落后或经济发展速度缓慢的地区，对人口会产生一种排斥力和离散力，人口机械增长一般为负值。与此同时粮食产量、出生率、死亡率，也是影响人口增长的因素。符号说明：用X、x、x、x表示粮食产量、GDP、出生率、死亡率。y表示总1234人口。样本确定通过查阅中国政府网，得到了1980年到2014年各因素的数据。表1-1样本数据年份粮食产量（万吨）GDP（亿元）出生率（%）死亡率（%）总人口（万人）201460702.61635910.212.377.16136782201360193.84588018.812.087.16136072201258957.9753412312.17.15135404201157120.85484123.511.937.14134735201054647.7140890311.97.11134091200953082.08345629.211.957.08133450200852870.92316751.712.147.06132802200750160.28268019.412.16.93132129200649804.23217656.612.096.81131448200548402.19185895.812.46.51130756200446946.95160714.412.296.42129988200343069.53136564.612.416.4129227200245705.7512100212.866.41128453200145263.67110270.413.386.43127627200046217.5299776.314.036.45126743199950838.5890187.714.646.46125786199851229.5384883.715.646.5124761199749417.179429.516.576.51123626199650453.571572.316.986.56122389199546661.861129.817.126.57121121199444510.148459.617.76.49119850199345648.835524.318.096.64118517199244265.827068.318.246.64117171199143529.321895.519.686.7115823199044624.318774.321.066.67114333198940754.917090.321.586.54112704198839408.115101.122.376.64111026198740297.712102.223.336.72109300198639151.210308.822.436.86107507198537910.89039.921.046.78105851198440730.57226.319.96.82104357198338727.55975.620.196.9103008198235450533322.286.61016541981325024898.120.916.36100072198032055.54551.618.216.34987051.3一元回归分析定义1.1回归分析(regressionanalysis)是确定两种或两种以上变量间相互依赖的定量关系的一种统计分析方法。如果在回归分析中，只包括一个自变量和一个因变量，且二者的关系可用一条直线近似表示，这种回归分析称为一元线性回归分析。绘制总人口与粮食产量的散点图图1-1粮食产量与总人口散点图设定理论模型根据图1-1随着粮食产量的增加，总人口的数量增加，且各样本点大致落在一条直线附近，故可以采用公式：y二卩。+卩iX+8(1-1)一元线性回归理论模型，对数据进行一元回归分析。1.3.3回归诊断表1T模型摘要模型RR平方調整後R平方標準偏斜度錯誤1.909a.826.8214985.99669a.預測值：(常數)，粮食产量(万吨)根据表1-1模型摘要表可以看到，r二0.909，说明以粮食产量为唯一因变量与总人口的拟合程度很高。表1-2系数非標準化係數標準化係數T顯著性

B標準錯誤Beta1（常數）53054.4045500.0139.646.000粮食产量（万吨）1.468.117.90912.513.000a.應變數\：人□（万人）根据表1-2系数表，将系数带入公式（1-1）可得出回归模型公式：用一元回归模型对2014年的总人口进行预测，的到的预测值y142165.8355第2章多元线性回归定义2.1在回归分析中，如果有两个或两个以上的自变量，就称为多元回归。事实上，一种现象常常是与多个因素相联系的，由多个自变量的最优组合共同来预测或估计因变量，比只用一个自变量进行预测或估计更有效，更符合实际。数据中心化标准化数据中心化和标准化在回归分析中的意义是取消由于量纲不同、自身变异或者数值相差较大所引起的误差。图2-1标准化结果多元回归模型建立多元回归模型公式：对数据进行多元回归分析，结果如图2.2表2-1模型摘要模型RR平方調整後R平方標準偏斜度錯誤1.972a.944.9362967.56830a.預測值：（常數），死亡率（%）,出生率（％）,粮食产量（万吨），GDP（亿元）表2-2變異數分析a模型平方和df平均值平方F顯著性1迴歸43.984410.996126.295.000b殘差7.616308806461.587總計91.60034應變數：人□（万人）預測值：（常數），死亡率（％）,出生率（%）,粮食产量（万吨），GDP（亿元）表2-3係數a模型非標準,化係數標準化係數T顯著性B標準錯2口誤Beta1（常數）124574.66721501.7295.794.000

粮食产量（万吨）1.045.152.6476.893.000GDP（亿元）-.005.008-.078-.635.531出生率（％）-1460.798262.207-.496-5.571.000死亡率（％）-4060.2723649.700-.091-1.112.275a.應變數\：人□（万人）根据表2-1可知R=0.972，拟合度高，所以能用该模型进行预测，同时模型的检验P值sigv0.05，说明该模型显著可信；根据表2-3得到多远线性回归方程：同时，根据方程得到14年总人口预测值y二138182.2403，与真实值相差e=1400.25403.2.3逐步回归法定义2.2逐步回归法：逐步回归的基本思想是有进有出。具体做法是将变量一个一个的引入，每引入一个变量后，对已引入变量进行逐个检验，当原引入变量因后引入变量变得不显著时，对其进行剔除，保证最后所得的回归子集是最优回归子集。根据图2.2可知，变量GDP与死亡率未通过显著性检验，因此采用逐步回归的方法对方程进行优化。表2-4逐步回归係數a模型非標準化係數標準化係數T顯著性B標準錯誤Beta1（常數）165146.3883543.31346.608.000出牛率（％）-2687.949210.086-.912-12.794.0002（常數）108407.8058355.07212.975.000出生率（％）-1534.270211.155-.521-7.266.000粮食产量（万吨）.816.116.5057.052.0003（常數）134678.87814308.5429.412.000出生率（％）-1366.202213.615-.464-6.396.000粮食产量（万吨）1.012.141.6267.184.000死亡率（％）-5681.5572580.897-.127-2.201.035a.應變數\：人□（万人）根据图表2-4,经过逐步回归，剔除了GDP这一变量得到一个拟合度更好的方程:同时得到2014年总人口的预测值y二138506.21662与真实值相差e二172.21662对比一般多元回归方法的出的模型预测值与运用了逐步回归后的预测值，发现不用逐步回归的预测值更接近真实值；其原因主要有三方面：我国的人口普查为10年一次，其数据主要也是通过预测的出，所以数据本身存在较大的误差；在实际问题中，自变量之间通常存在相关性，当相关程度严重时称为多重共线性。自变量之间的多重共线性会对回归产生极大的影响；人口预测不适合采用线性多元回归的方法进行预测。2.4多重共线性定义2.3多重共线性（Multicollinearity）是指线性回归模型中的解释变量之间由于存在精确相关关系或高度相关关系而使模型估计失真或难以估计准确。多重共线性产生原因：（1）经济变量相关的共同趋势；（2）滞后变量的引入；（3）样本资料的限制。影响：（1）完全共线性下参数估计量不存在；（2）近似共线性下OLS估计量非有效，多重共线性使参数估计值的方差增大，1/（1-⑵为方差膨胀因子（VarianceInflationFactor,VIF）如果方差膨胀因子值越大，说明共线性越强。相反因为，容许度是方差膨胀因子的倒数，所以，容许度越小，共线性越强。可以这样记忆：容许度代表容许，也就是许可，如果，值越小，代表在数值上越不容许，就是越小，越不要。而共线性是一个负面指标，在分析中都是不希望它出现，将共线性和容许度联系在一起，容许度越小，越不要，实际情况越不好，共线性这个“坏蛋”越强。进一步，方差膨胀因子因为是容许度倒数，所以反过来；（3）参数估计量经济含义不合理；（4）变量的显著性检验失去意义，可能将重要的解释变量排除在模型之外；（5）模型的预测功能失效。变大的方差容易使区间预测的“区间”变大，使预测失去意义。2.3.1多重共线性检测用SPSS22.0对数据进行多重共线性检测表2-5共線性診斷a模型維度特徵值條件指數變異數比例（常數）出牛率（％）粮食产量（万吨）死亡率（％）111.9721.000.01.012.0288.447.99.99212.9291.000.00.00.002.0686.557.00.18.063.00333.7981.00.81.94313.9271.000.00.00.00.002.0707.506.00.17.03.003.00336.436.15.79.74.034.00186.828.85.04.24.97a.應變數：人□（万人）检测结果表2-5显示：（1）特征值有3个接近0，表明存在多重共线性；（2）条件指数有3个＞10，表明可能存在多重共线性；（3）在相关系数矩阵中，死亡率数值接近1可能存在多重共线性表2-6排除共线性变量后的係數a模型非標準化係數標準化係數T顯著性共線性統計資料B標準錯誤Beta允差VIF1（常數）165146.3883543.31346.608.000出牛率（％）-2687.949210.086-.912-12.794.0001.0001.0002（常數）108407.8058355.07212.975.000出生率（％）-1534.270211.155-.521-7.266.000.4002.502粮食产量（万吨）.816.116.5057.052.000.4002.5023（常數）134678.87814308.5429.412.000出生率（％）-1366.202213.615-.464-6.396.000.3492.868粮食产量（万吨）1.012.141.6267.184.000.2414.146死亡率（％）-5681.5572580.897-.127-2.201.035.5531.808a.應變數\：人□（万人）根据表2-6,剔除GDP这一变量后，其余变量的VIF全部小于10,排除变量间存在多重共线性的原因，不需要再次剔除变量。得到剔除共线性后的多元线性回归方程：y二134678.878+1.012x-1366.202x-5681.557x1342.4主成分分析定义2.4主成分分析（PrincipalComponentAnalysis,PCA），是一种统计方法。通过正交变换将一组可能存在相关性的变量转换为一组线性不相关的变量，转换后的这组变量叫主成分。在实际课题中,为了全面分析问题,往往提出很多与此有关的变量（或因素）,因为每个变量都在不同程度上反映这个课题的某些信息。主成分分析首先是由K.皮尔森对非随机变量引入的，尔后H.霍特林将此方法推广到随机向量的情形。信息的大小通常用离差平方和或方差来衡量。主成分分析的主要原理是寻找一个适当的线性变换：（1）将彼此相关的变量转变为彼此独立的新变量；（2）方差较大的几个新变量就能综合反应原多个变量所包含的主要信息；（3）新变量各自带有独特的专业含义。住成分分析的作用是：（1）减少指标变量的个数；（2）决多重相关性问题主成分分析模型建立运用SPSS22.0对数据进行主成分分析表2-7各变量所占权重元件起始特徵值擷取平方和載入總計變異的%累加%總計變異的%累加%13.04076.00476.0043.04076.00476.0042.72718.17694.1793.1523.79397.9724.0812.028100.000擷取方法：主體元件分析。根据主成分分析结果可以看出前两个变量所占比重最多，二者的和所占比例为94.18%，所以可以采用前两个变量建立回归模型。表2-8回归係數a模型非標準化係數標準化係數T顯著性B標準錯誤Beta1（常數）53075.7119733.1915.453.000粮食产量（万吨）1.468.235.9086.252.000GDP（亿元）2.494E-5.009.000.003.998a.應變數\：人□（万人）根据表2-8,得到回归方程：根据回归方程预测出2014年人口y142180.52727万人。因为：1.主成分分析中，我们首先应保证所提取的前几个主成分的累计贡献率达到一个较高的水平（即变量降维后的信息量须保持在一个较高水平上），其次对这些被提取的主成分必须都能够给出符合实际背景和意义的解释（否则主成分将空有信息量而无实际含义）。2.主成分的解释其含义一般多少带有点模糊性，不像原始变量的含义那么清楚、确切，这是变量降维过程中不得不付出的代价。因此，提取的主成分个数m通常应明显小于原始变量个数p（除非p本身较小），否则维数降低的“利”可能抵不过主成分含义不如原始变量清楚的“弊”。所以预测的结果出现了更大的误差第3章非线性模型曲线回归定义3.1非线性回归是在掌握大量观察数据的基础上，利用数理统计方法建立因变量与自变量之间的回归关系函数表达式（称回归方程式）。回归分析中，当研究的因果关系只涉及因变量和一个自变量时，叫做一元回归分析；当研究的因果关系涉及因变量和两个或两个以上自变量时，叫做多元回归分析。确定两个变数间数量变化的某种特定的规则或规律；估计表示该种曲线关系特点的一些重要参数，如回归参数、极大值、极小值和渐近值等；为生产预测或试验控制进行内插，或在论据充足时作出理论上的外推。3.1.1曲线拟合对国内总人口的拟合，选取总人口指标为因变量，单位为万人，拟合总人口关于时间t的趋势曲线。以1980年为基准年，取值为t=1,2014年t=35。绘制总人口与变量t的散点图，如图3-1所示。图3-1总人口对t的散点图从散点图可以看到，总人口大致符合三次函数形式，当人口的增长速度大致相同时，其趋势线就是三次函数形式。图3-2拟合曲线图通过观察图3-2，发现三次曲线模型拟合度最好，其次为线性和复合模型，故根据公式：建立三次曲线模型。表3-1模型摘要RR平方調整後R平方標準偏斜度仝±±兰口市錯n誤1.000.999.999365.498自變數為时间。表3-2變異數分析平方和df平均值平方F顯著性迴歸35.390345.13011749.73.0002殘差4141256.2131133588.9100總計91.60034自變數為时间。表3-3三次曲线系数非標準化係數標準化係數T顯著性B標準錯誤Beta时间1849.09465.5011.60928.230.000时间-17.6834.196-.571-4.214.000**2时间-.070.077-.076-.908.371**3（常96103.58276.15348.00.000數）874根据图表3-1,得到三次函数模型的R二1.000，说明拟合程度非常好；同时根据表3-2可知，方程整体都通过了显著性检验；根据表3-3可知，不但方程整体通过了显著性检验，每个不同次幂的t也经过了显著性检验，所以建立的三次方程式完全符合实际情况的。根据第三个表建立出三次函数方程：同时得出2014年预测值y二136174.94841万人，预测值与真实值相差e二607.0516万人。预测值与真实值相差程度远小于运用线性多元回归方法、逐步回归法和主成分分析后的预测值与真实值之间的差异，进一步证明建立的三次函数模型符合实际情况。Logistic模型Logistic模型增长公式为：其中p为时刻的人口总数，p为人口极限规模，e为自然对数的底，t为时刻长tm度，a、b为待定参数。Logistic模型考虑到人口总数增长的有限性，提出了人口总数增长的规律即随着人口总数的增长，人口增长率逐渐下降，但对于在短期内如30-50年内人口增长可能呈上升趋势如人口生育率上升、死亡率下降等原因而导致人口呈上升趋势。Logistic模型在应用中对时间长，人口数据变化大，因此误差较大且不稳定。而小城镇人口的变化就存在人口数据变化较大的特点，所以Logistic模型对小城镇人口的预测并不适合。用对人口进行Logistic曲线拟合，运用命令：x=[20142013201220112010200920082007200620052004200320022001200019991998199719961995199419931992199119901989198819871986198519841983198219811980]y=[13678213607213540413473513409113345013280213212913144813075612998812922712845312762712674312578612476112362612238912112111985011851711717111582311433311270411102610930010750710585110435710300810165410007298705]x=x';y=y';st_=[20001276271];(在x,y内任意取的数)ft_=fittype('a/(1+b*exp(-k*(x-1980)))',...'dependent',{'y'},'independent',{'x'},...'coefficients',{'a','b','k'});cf_=fit(x,y,ft_,'Startpoint',st_)最后运行出的结果：cf_=Generalmodel:cf_(x)=a/(1+b*exp(-k*(x-1980)))Coefficients(with95%confidencebounds):a=1.217e+005(1.144e+005,1.29e+005)b=1.176e+005(-4.741e+009,4.742e+009)k=13.21(-4.032e+004,4.035e+004)根据运行结果，得到Logistic拟合方程：其中a为人口上线的估计值，因为在1995年人口的实际值就已经超过了预测上限，所以该模型不是和用于预测未来的人口。造成模型不准确的原因主要是数据过少，且人口数据多数来自抽样调查，数据本身存在一定误差。结论通过不同的模型建立方法对获得的人口数据建立了多个不同的预测模型；经过对

比发现三次曲线模型是最符合实际运用的；导致其他模型不适用的原因主要为：（1）目前我国的人口出生率低、死亡率低、自然增长率低．这表明我国将进入“低、低、低”现代人口再生产类型的行列；即我国今后的人口数量趋势大致会成为一条水平线。（2）自20