当前课程知识点:系统工程原理 > 第四章 系统预测 > 第五节 回归分析预测法应用 > 回归分析预测法应用
接下来第五节
我们通过两个例子来具体地理解一下
回归分析预测的应用
首先看一元线性回归分析预测的例子
这个表格里面给出来了
阿尔卑斯山和苏格兰的17个地方
沸点以及大气压的数据
其中第二列
是17个地点的沸点数据
第三列是
以英寸汞柱为单位的
大气压的数据
把这一列数据做一个转换以后
得到最后一列的以
100乘以log英寸汞柱为单位的气压的数据
那么令
因变量是大气压
自变量是沸点x
我们把它先描在散点图上一看
大体上呈现出一个
线性相关的关系
y是随着x的递增而递增的
我们采用
一元线性回归模型
去进行预测
按照最小二乘法
可以求解到a和b
两个回归系数
分别是-42.138
和0.895
那么
预测方程就是y= -42.138+0.895x
然后进行统计检验
根据
所给出来的样本点值可以计算
y和x的相关系数rxy
计算得到rxy=0.997
而给出显著性水平是0.01的时候
α取0.01 然后查相关系数表 得到临界值rc
在样本点数是17的情况下
等于0.606
很显然这个检验是通过了的
于是这个一元线性回归方程是可以用于预测的
第二个例子
我们来看一下多元线性回归分析预测的例子
这个问题的背景是
有一个部队维修厂
它的固定资产x1 维修工人数x2
和维修的总产值y
它的统计数据在这个表格里面给出来了
总共10个点
现在请建立以x1、x2为自变量的
总产值的回归预测模型
现在我们看到的这个表
除了刚才给出的信息以外
还增加了很多列
这都是在计算回归系数的时候需要用到的一些
我们先把它算好放在这里
根据上一节我们所介绍的
回归系数的最小二乘法进行估计
= (XTX)-1 XTY
可以最终得到
a0 a1 a2
分别是最后这个结果向量中的有关的元素
那么这样的话回归模型就已经得到了
就是y=-106.7218+0.498921x1+1.34047x2
我们现在
只有一个因变量和两个自变量
因此还是可以画图的
把它画在一个三维坐标系上来观察一下
是不是大体上呈现一个线性相关关系呢?
现在我们看到的这个图
是不太看得清楚的
把这个坐标旋转一下
可以看到y与x1、x2
确实大体呈现出一个线性相关关系
但是
光这样观察散点图
还是不足以
进行检验
所以接下来要进行的是
上一节我们介绍的一系列的检验
第一个检验是标准离差检验
按照标准离差的计算公式可以算到S等于4.3776
标准离差S除以y的
样本平均
等于0.0160
低于阈值10%
因此标准离差检验是通过的
第二个检验是相关系数检验
根据样本点值计算出来相关系数r
等于0.9897
取α=0.05
自变量有两个
n就是样本点数是10个
去查相关系数表
等于0.666
这是阈值
计算出来的相关系数r是超过这个阈值的
因此相关系数检验通过
第三个检验是F检验
将样本点数据代入F统计量的公式
算到F=167.17
然后取α=0.05
自由度分别为
2和7
2就是自变量的个数
7就是n-k-1
就是
样本点个数减去自变量个数再减1
那么查到临界值等于4.74
计算出来的F值是超过这个临界值的
因此 F检验通过
第四个检验是t检验
这里面要进行两次t检验
是说
相对于两个自变量而言 分别要进行t检验
我们首先看t1
它的计算式里面a1
其实就等于
x1所对应的那个回归系数
的估计值0.498921
S是标准离差
4.3776
c11呢?
现在我们回过头来看矩阵C
我们目前看到的
就是
回归系数矩阵它的估计式以及它的计算过程
C矩阵就是(XTX)-1
那么c11呢
就是在这个矩阵上的第2行第2列元素
为什么是第2行第2列元素
因为只有这个元素
才是
对应于第一个自变量的
c22呢
是C矩阵上的第3行第3列元素
因为这个元素是对应于第二个自变量的
t1和t2两个量值都可以计算出来
接下来我们取α=0.05
去查t表得到了阈值是2.365
发现
t检验不能通过
如果取α=0.1的话
阈值等于1.895
检验通过了
第五个检验是DW检验
根据这个公式可以计算出来
DW值是1.9636
去查DW表
得到du等于1.32
dl等于0.879
然后去
对查我们的DW检验的表格
发现检验是通过的
接下来可以给结论了
预测的区间估计
如果取
x10等于350 x20等于190
进行一个系统预测的话
把x10等于350 x20等于190 分别代入
预测方程
可以得到y0=322.5899
这就是一个点估计
又已知S=4.3776
在90%的置信度下的
预测上下限分别为y0 + 1.645S
和y0 - 1.645S
结果分别为
预测上限为329.7911
预测下限为315.3887
最后给出的结论是
预测模型为y = -106.7218+0.498921x1+1.34047x2
其中
还需要给出一些其他的信息
S标准离差等于4.3776
复相关系数r等于0.9897
F量值等于167.17
DW
等于1.9636
当x10等于350 x20等于190时
90%的置信区间为
315.3887到329.7911
这样的话
一个多元线性回归分析预测就做完了
-课程概述
--课程概述
-为何学习系统工程
--为何学习系统工程
--人类基因工程
-第一节 系统的基本概念
--系统的基本概念
-第一节测试 系统的基本概念
-第二节 系统的分类
--系统的分类
-第二节测试 系统的分类
-第三节 系统思想的发展
--系统思想的发展
-第三节测试 系统思想的发展
-第四节 系统工程的概念
--系统工程的概念
-第四节测试 系统工程的概念
-第五节 系统工程的特点
--系统工程的特点
--扩展阅读
-第五节测试 系统工程的特点
-第六节 系统工程的发展与应用
--扩展资料
-第六节测试 系统工程的发展与应用
-第一节 还原论与整体论
--还原论与整体论
-第一节测试 还原论与整体论
-第二节 霍尔三维结构
--霍尔三维结构
-第二节测试 霍尔三维结构
-第三节 霍尔三维结构的应用
-第四节 切特兰德的调查学习法
-第四节测试 切特兰德的调查学习法
-第五节 调查学习模式的应用
-第六节 综合集成工程
--综合集成工程
-第六节测试 综合集成工程
-第一节 系统分析概述
--系统分析概述
-第一节测试 系统分析概述
-第二节 系统分析应用
--系统分析应用
-第三节 系统目标分析的目标手段法
--目标手段法
-第三节测试 系统目标分析的目标手段法
-第四节 系统环境分析的SWOT法
--SWOT法
-第四节测试 系统环境分析的SWOT法
-第五节 系统功能分析的IDEF0方法
--IDEF0法
-第五节测试 系统功能分析的IDEF0方法
-第六节 系统结构分析的ISM法
--ISM方法
-第六节测试 系统结构分析的ISW法
-第七节 ISM方法的应用
--ISM应用
--ISM-扩展阅读
-第一节 系统预测概述
--系统预测概述
-第一节测试 系统预测概述
-第二节 时间序列分析预测
--时间序列分析预测
-第二节测试 时间序列分析预测
-第三节 时间序列分析预测应用
--时间序列预测应用
-第四节 回归分析预测法
--回归分析预测法
-第四节测试 回归分析预测法
-第五节 回归分析预测法应用
-第六节 马尔可夫预测法
--马尔可夫预测法
-第六节测试 马尔可夫预测法
-第七节 马尔可夫预测法应用
--马尔可夫预测应用
-第一节 系统评价基本概念
-第一节测试 系统评价基本概念
-第二节 系统评价的一般步骤
-第三节 指标体系的建立
--指标体系的建立
-第三节测试 指标体系的建立
-第四节 指标权重确定方法
--指标权重确定方法
-第四节测试 指标权重确定方法
-第五节 指标数量化方法
--指标数量化方法
-第五节测试 指标数量化方法
-第六节 指标综合方法
--指标综合方法
-第六节测试 指标综合方法
-第七节 层次分析法及应用
--层次分析法
-第七节测试 层次分析法及应用
-第八节 主成分分析及应用
--主成分分析及应用
-第八节测试 主成分分析及应用
-第九节 数据包络分析法
--数据包络分析法
-第九节测试 数据包络分析法
-第一节 系统决策概述
--系统决策概述
-第一节测试 系统决策概述
-第二节 完全不确定型决策方法
--完全不确定型决策
-第二节测试 完全不确定型决策方法
-第三节 风险型决策方法
--风险型决策方法
--决策树应用案例
-第三节测试 风险型决策方法
-第四节 贝叶斯决策方法
--贝叶斯决策方法
-第四节测试 贝叶斯决策方法
-引言
--引言
-第七章 网络计划技术--引言
-第一节 网络图的组成与绘制
-第一节测试 网络图的组成与绘制
-第二节 网络图的参数计算
--网络图的参数计算
-第二节测试 网络图的参数计算
-第三节 网络图关键路线的确定
-第三节测试 网络图关键路线的确定
-第四节 任务按期完成概率分析与计算
-第四节测试 任务按期完成概率分析与计算
-期末考核