多元线性回归分析医学论文

实验三 多元回归模型【实验目的】掌握建立多元回归模型和比较、筛选模型的方法。【实验内容】建立我国国有独立核算工业企业生产函数。根据生产函数理论，生产函数的基本形式为： 。其中，L、K分别为生产过程中投入的劳动与资金，时间变量 反映技术进步的影响。表3-1列出了我国1978-1994年期间国有独立核算工业企业的有关统计资料；其中产出Y为工业总产值（可比价），L、K分别为年末职工人数和固定资产净值（可比价）。表3-1 我国国有独立核算工业企业统计资料年份 时间 工业总产值Y（亿元） 职工人数L（万人） 固定资产K（亿元）1978 1 3139 2 3208 3 3334 4 3488 5 3582 6 3632 7 3669 8 3815 9 3955 10 4086 11 4229 12 4273 13 4364 14 4472 15 4521 16 4498 17 4545 资料来源：根据《中国统计年鉴－1995》和《中国工业经济年鉴-1995》计算整理【实验步骤】一、建立多元线性回归模型一建立包括时间变量的三元线性回归模型；在命令窗口依次键入以下命令即可：⒈建立工作文件： CREATE A 78 94⒉输入统计资料： DATA Y L K⒊生成时间变量 ： GENR T=@TREND(77)⒋建立回归模型： LS Y C T L K则生产函数的估计结果及有关信息如图3-1所示。 图3-1 我国国有独立核算工业企业生产函数的估计结果因此，我国国有独立工业企业的生产函数为： （模型1） ＝() () () () 模型的计算结果表明，我国国有独立核算工业企业的劳动力边际产出为，资金的边际产出为，技术进步的影响使工业总产值平均每年递增亿元。回归系数的符号和数值是较为合理的。 ，说明模型有很高的拟合优度，F检验也是高度显著的，说明职工人数L、资金K和时间变量 对工业总产值的总影响是显著的。从图3-1看出，解释变量资金K的 统计量值为，表明资金对企业产出的影响是显著的。但是，模型中其他变量（包括常数项）的 统计量值都较小，未通过检验。因此，需要对以上三元线性回归模型做适当的调整，按照统计检验程序，一般应先剔除 统计量最小的变量（即时间变量）而重新建立模型。二建立剔除时间变量的二元线性回归模型； 命令：LS Y C L K则生产函数的估计结果及有关信息如图3-2所示。 图3-2 剔除时间变量后的估计结果因此，我国国有独立工业企业的生产函数为： （模型2） ＝() () () 从图3-2的结果看出，回归系数的符号和数值也是合理的。劳动力边际产出为，资金的边际产出为，表明这段时期劳动力投入的增加对我国国有独立核算工业企业的产出的影响最为明显。模型2的拟合优度较模型1并无多大变化，F检验也是高度显著的。这里，解释变量、常数项的 检验值都比较大，显著性概率都小于，因此模型2较模型1更为合理。三建立非线性回归模型——C-D生产函数。C-D生产函数为： ，对于此类非线性函数，可以采用以下两种方式建立模型。方式1：转化成线性模型进行估计；在模型两端同时取对数，得： 在EViews软件的命令窗口中依次键入以下命令：GENR LNY=log（Y）GENR LNL=log（L）GENR LNK=log（K）LS LNY C LNL LNK则估计结果如图3-3所示。 图3-3 线性变换后的C-D生产函数估计结果即可得到C-D生产函数的估计式为： （模型3） ＝ () () () 即： 从模型3中看出，资本与劳动的产出弹性都是在0到1之间，模型的经济意义合理，而且拟合优度较模型2还略有提高，解释变量都通过了显著性检验。方式2：迭代估计非线性模型，迭代过程中可以作如下控制：⑴在工作文件窗口中双击序列C，输入参数的初始值；⑵在方程描述框中点击Options，输入精度控制值。控制过程：①参数初值：0，0，0；迭代精度：10－3；则生产函数的估计结果如图3-4所示。 图3-4 生产函数估计结果此时，函数表达式为： （模型4） ＝()(－)() 可以看出，模型4中劳动力弹性 ＝，资金的产出弹性 ＝，很显然模型的经济意义不合理，因此，该模型不能用来描述经济变量间的关系。而且模型的拟合优度也有所下降，解释变量L的显著性检验也未通过，所以应舍弃该模型。②参数初值：0，0，0；迭代精度：10－5； 图3-5 生产函数估计结果从图3-5看出，将收敛的误差精度改为10－5后，迭代100次后仍报告不收敛，说明在使用迭代估计法时参数的初始值与误差精度或迭代次数设置不当，会直接影响模型的估计结果。③参数初值：0，0，0；迭代精度：10－5，迭代次数1000； 图3-6 生产函数估计结果此时，迭代953次后收敛，函数表达式为： （模型5） ＝()()() 从模型5中看出，资本与劳动的产出弹性都是在0到1之间，模型的经济意义合理， ，具有很高的拟合优度，解释变量都通过了显著性检验。将模型5与通过方式1所估计的模型3比较，可见两者是相当接近的。④参数初值：1，1，1；迭代精度：10－5，迭代次数100； 图3-7 生产函数估计结果此时，迭代14次后收敛，估计结果与模型5相同。比较方式2的不同控制过程可见，迭代估计过程的收敛性及收敛速度与参数初始值的选取密切相关。若选取的初始值与参数真值比较接近，则收敛速度快；反之，则收敛速度慢甚至发散。因此，估计模型时最好依据参数的经济意义和有关先验信息，设定好参数的初始值。二、比较、选择最佳模型估计过程中，对每个模型检验以下内容，以便选择出一个最佳模型：一回归系数的符号及数值是否合理；二模型的更改是否提高了拟合优度；三模型中各个解释变量是否显著；四残差分布情况以上比较模型的一、二、三步在步骤一中已有阐述，现分析步骤一中5个不同模型的残差分布情况。分别在模型1～模型5的各方程窗口中点击View/Actual, Fitted, Residual/ Actual, Fitted, Residual Table（图3-8），可以得到各个模型相应的残差分布表（图3-9至图3-13）。可以看出，模型4的残差在前段时期内连续取负值且不断增大，在接下来的一段时期又连续取正值，说明模型设定形式不当，估计过程出现了较大的偏差。而且，模型4的表达式也说明了模型的经济意义不合理，不能用于描述我国国有工业企业的生产情况，应舍弃此模型。模型1的各期残差中大多数都落在 的虚线框内，且残差分别不存在明显的规律性。但是，由步骤一中的分析可知，模型1中除了解释变量K之外，其余变量均为通过变量显著性检验，因此，该模型也应舍弃。模型2、模型3、模型5都具有合理的经济意义，都通过了 检验和F检验，拟合优度非常接近，理论上讲都可以描述资本、劳动的投入与产出的关系。但从图3-13看出，模型5的近期误差较大，因此也可以舍弃该模型。最后将模型2与模型3比较发现，模型3的近期预测误差略小，拟合优度比模型2略有提高，因此可以选择模型2为我国国有工业企业生产函数。 图3-8 回归方程的残差分析 图3-9 模型1的残差分布图3-10 模型2的残差分布图3-11 模型3的残差分布图3-12 模型4的残差分布图3-13 模型5的残差分布

取10%的置信度还有一个变量不显著的话你要考虑下是不是存在多重共线性，异方差，自相关等的问题了。F检验是针对整个回归方程的吧，自变量的t检验也是要写出来的。就比如变量间存在多重共线性时，F值会比较大，检验也会通过，r平方也够大，但是某一个或几个自变量的系数就不会通过检验，不显著的。

百度文库。或学校图书馆网站都是有的

目的研制《糖尿病联络护士专科护理管理行为量表》，作为联络护士在糖尿病专科护理管理方面工作的评价工具；调查联络护士糖尿病相关管理工作现状，探讨其相关影响因素。方法1.应用德尔菲（Delphi）法，向12位在糖尿病专科护理和管理工作领域的专家发送2次问卷并对咨询结果进行分析，筛选优化量表条目，最终形成《糖尿病联络护士专科护理管理行为量表》；2.使用《糖尿病联络护士专科护理管理行为量表》调查粤东、粤西、粤北和珠三角区域25家医院的426例糖尿病联络护士，评价联络护士在糖尿病相关专科护理管理工作现状；3.采用单因素方差分析、两样本t检验、多元线性回归分析方法对《糖尿病联络护士专科护理管理行为量表》的结果数据统计分析。结果1.德尔菲（Delphi）法两轮咨询量表维度重要性评分的协调系数分别为(P <)和(P <)，两轮咨询量表条目的变异系数～和～，专家权威系数均值为；行为量表的重测信度系数-Spearman相关系数为，总的克朗巴赫ɑ系数为。提取公因子有4个，累计方差贡献率为；2.糖尿病联络护士专科护理管理行为平均每一条目得分为（±分）；所有维度中维度2“血糖监测及胰岛素注射管理”中平均每一条目得分最高（±分）；所有条目中得分均值最高（±分）的是条目5“检查病区血糖试纸有效期”，得分均值最低（±分）的是条目4“统计病区患者发生的低血糖事件”；3.经单因素方差分析,每周是否有固定的糖尿病护理管理时间、对糖尿病护理管理内容的了解情况、不同科室、职称、成为糖尿病联络护士的工作年限等因素在行为量表总分或某一维度得分上有差别。结论1.本研究应用Delphi法研制的《糖尿病联络护士专科护理管理行为量表》信效度良好，条目较为可靠；2.糖尿病联络护士专科护理管理工作不容乐观，整体未达到经常管理水平，相对来说，技术性的工作如血糖监测及胰岛素注射管理等相对比较到位，但仍对病区患者发生的低血糖事件缺乏重视；3.每周是否有固定的糖尿病护理管理时间、对糖尿病护理管理内容的了解情况、不同科室、职称、作为糖尿病联络护士的工作年限、月收入水平、管理体制对糖尿病联络护士护理管理工作产生不同程度的影响。[1] 丁飚,张妮娜,阚凯,卢逢娣,曾辉. 662名临床护士血糖监测相关知识认知情况分析及培训对策[J]. 护理学报. 2012(24)[2] 李意华,卢红华,钟静静. 时间管理法在护理管理工作中的应用[J]. 赣南医学院学报. 2012(05)[3] 谢晓清,姚美华. 病区胰岛素注射的质量管理[J]. 护理与康复. 2012(04)[4] 闫洋洋,穆欣,郭闯. 血液透析护理环节质量评价指标的研究[J]. 中国护理管理. 2012(04)[5] 王春枝,斯琴. 德尔菲法中的数据统计处理方法及其应用研究[J]. 内蒙古财经学院学报(综合版). 2011(04)[6] 周佩如,黄洁微. 广东省糖尿病专科护士工作现状调查及分析[J]. 护士进修杂志. 2011(13)[7] 李华,黄惠根,黄蝶卿,余桂芳,洪涛. 专科护士专职工作岗位的设立及实践[J]. 中华护理杂志. 2011(04)[8] 郑刚. 急性心肌梗死患者血糖水平与死亡率:临床研究证据说明了什么[J]. 中国心血管杂志. 2010(05)[9] 徐文聪,李华云,吴梅. 糖尿病联络小组的建立与实践[J]. 中医药临床杂志. 2010(10)[10] 彭卫群,侯晶,吴清香,郭晶晶. 糖尿病联络护士小组工作模式与实践[J]. 齐齐哈尔医学院学报. 2010(07)