数学实验12回归分析Word格式文档下载.docx
《数学实验12回归分析Word格式文档下载.docx》由会员分享，可在线阅读，更多相关《数学实验12回归分析Word格式文档下载.docx（18页珍藏版）》请在冰豆网上搜索。
40
42
57
54
52
50
X2
X3
X4
62
45
58
48
X5
178
185
156
166
174
164
146
（1）若x1~x5中只许选择1个变量，最好的模型是什么？
（2）若x1~x5中只许选择2个变量，最好的模型是什么？
（3）若不限制变量的个数，最好的模型是什么？
（4）对最终模型观察残查，有无异常点，若有，剔除后如何？
1．模型建立
本题不同小问需要建立不同模型，由于专业知识所限，并且提供的数据较少，难以做出精确符合现实情况的模型，因此这里采用最简单的线性回归法进行拟和,模型基本形式如下：
事实上，
中的项（高次项和交互项）对于本题目来讲意义不大，因为所给定的5个自变量和因变量之间关系比较模糊，几个变量彼此之间的联系也很难说清，因此用自变量的一次线性拟和就足以适应本题的要求。
但作为练习，还是将每种回归方法都使用到了，可以用于参考。
具体采用的各个模型将在下面单独说明，这里不再重复。
2．程序设计
由于本题需要建立多组模型，并且要在不断的调试中发现最合理的，很多命令都要在这个过程中不断使用，这里仅仅给出使用的最基本的命令。
数据
clear
A=[…];
%数据矩阵，略
n=24;
y=A（2,:
）;
%提取各个数据
x1=A（3,:
x2=A（4,:
x3=A（5,:
x4=A（6,:
x5=A（7,:
绘制散点图（大致判断影响情况）
fori=1:
5
subplot（2,3,i）,plot（A（i+2,:
）,y,'
+'
）,grid
pause
end
单参数回归（第一问）
X=[ones（n,1）,x4'
];
%这里检验的是自变量x4，实际操作时要分别检验x1~x5
[b,bint,r,rint,s]=regress（y'
X）;
%回归分析程序（
=）
b,bint,s,%输出回归系数估计值、置信区间、以及统计量
rcoplot（r,rint）%残差图
Polytool（x3'
y'
2）%检验一元多项式回归的结果，输出交互式画面
双参数回归（第二问）：
用逐步回归法找出最合理的两个变量
X5=[x1'
x2'
x3'
x4'
x5'
stepwise（X5,y'
）；
%利用输出的交互式画面，可以选出最佳的两个变量
XX=[x3'
x1'
%当得到了最佳的两个变量后（这里假设是x3\x1）
rstool（XX,y'
'
linear'
）%检验二元情况下的交互项和高次项
全部参数回归（第三问）：
%仍然用逐步回归法找出最合理的组合方式
）
第五问要求对残差进行分析，并且剔除异常点，可以在该问得到最终模型后，采用regress得到的残差值和置信区间并根据其绘制残差图，然后再进行剔除操作重新检验。
3．运行结果及分析
散点图
散点图
从左上到右下的顺序为x1~x5.可以由点的分布大致看出，除了x3自变量呈现比较明显的负相关趋势以外，对于其他的各个自变量都难以直接观测出其对于因变量的影响。
根据这种结果，可以假设自变量x3（1500m跑后心速）最直接的与锻炼耗氧能力相关，下面通过对各个自变量的单参数回归进行检验。
单参数回归
被检对象
β0
β1
β1置信区间
R^2
F
p
s^2
由单参数回归的结果可以证明X3（1500m跑后心速）可以最好的反映出y（锻炼耗氧能力）的情况。
由β1置信区间可以看出，x1、x2包含0在内，即y可能与该参数无关，所以不选择，并且两者的p值已经明显的大于
=，则不考虑x1、x2。
比较x3~x5后发现，x3的
-决定系数明显的大于x4、x5的，决定系数反映的是在因变量的总变化中自变量引起的那部分的比例，
大说明x3自变量对因变量起的决定作用最大。
并且x3的p和s^2值也都比较小，所以最终确定x3可以最好的反映出y的情况。
用Polytool检验含x3高次（2次）项的情况，参量Export如下表：
β2
回归系数估计值
置信区间下限
置信区间上限
可以同之前的仅含一次项的结果进行比较，发现各个参量的置信区间都很宽，且β2的置信区间过0。
可以认为二次项的引入是不重要的。
因此采用如下单参数模型描述y是最准确地：
其中：
双参数回归：
用stepwise作逐步回归，部分过程和最终结果如下图：
只取x3自变量的回归结果
取x3和x5自变量的回归结果
取x3和x1自变量的回归结果
根据题目要求，最终得到取双参量时的最佳结果（RMSE参量最小）是取x3（1500m跑后心速）和x1（年龄）自变量。
但事实上，实际的逐步回归过程在此时并没有结束，最终的最优结果是只取x3参量。
这说明取x3、x1参量同只取x3相比优势并不明显。
通过rstool命令检验二元情况下的交互项和高次项情况,下图是linear情况下固定单参数进行预测的结果：
项对应的系数
常数项
X3^2
X1^2
X3*X1
RMSE
Linear
Purequadratic
Interaction
Quadratic
可以看到高次项和相关项的系数都非常小，说明其对于y的影响不大。
根据rmse的结果进行比较，仍然选择linear回归方式，即只用二元自变量的一次项。
全参数回归
根据以上的分析可以验证模型建立时的猜想，本题中5个字变量和y的关系都不是很直接的，除x3外其他变量的影响很小，所以在最终完整模型中，不再考虑高次项和交互项的影响，一方面简化模型，一方面大大节省的筛选的时间。
所以采用stepwise命令，仅对五元变量x1~x5的一次项进行回归分析，结果如下：
Export参数结果：
β
置信区间
Coeff.
se
t-stat
p-val
00
00
R^2=
F=
RMSE=
P=*e-7
最终取以下三个参数得到最佳回归结果：
x3（1500m跑后心速）、x1（年龄）以及x5（跑步后心速）。
但仍需要进行一般回归分析（regress）确定常数项并观察残差，结果如下:
残差图（全部点）
可以看到10和15号数据异常，剔除，再次观察残差，结果如下
残差图（剔除10，15号点）
可以看到4号数据变为异常，再次剔除，结果仍然有异常点。
过程从略，最终经过4次剔除，去掉5个点（4,10,15,17,23）以后，得到没有异常点的模型:
去掉所有异常点之后的结果
最终得到的结果整体上优越于剔除异常点之前的结果（不再粘贴结果）。
但是事实上，由于数据点经过剔除不断的结果，模型最终的形式和实际统计到的24组数据的整体情况偏离越来越大，也就是说：
剔除异常点虽然能够一应程度上降低其对于整体情况的干扰作用，而剔除的过程也放大了其他原本正常数据点的异常性，所以异常点可能会不断产生，但是剔除的数量增加即采样数据的减少也会削弱模型反省整体性能的能力。
是一对矛盾，在数据点较少的时候尤其明显。
比较科学的做法是：
只进行1次或少次剔除，保证整体性，又去掉了最主要的异常点。
这里的最终结果采用剔除最初两个异常点（10，15号）后的结果，在此也附上完整数据（剔除之前）的结果，作为第3问的答案：
完整数据（第三题结果）：
回归参数
取值
β3
β5
R^2
其中：
一次剔除（最终结果）：
1500m跑后心速、年龄以及跑步后心速三个参数最能够反映锻炼耗氧量这个重要的身体状态指标。
三种心跳速度越快，说明耗氧量越大；
速度越慢，即时间越长，说明耗氧量越小。
习题11：
一个医药公司的新药研究部门为了掌握一种新型止痛剂的疗效，设计了一个药物试验，给24名患有同种痛病的病人使用这种新止痛剂的以下4个剂量中的某一：
2、5、7、10,（g）,并记录每个病人病痛明显减轻的时间（min）.为了了解新药的疗效与病人性别和血压有什么关系，试验过程中研究人员把柄热男性别集血压的低中高3档平均分配来进行测试。
通过比较每个病人血压的历史数据,从低到高分成3组，分别记作,,。
实验结束后，公司的记录结果见表格（略）.
请为公司建立一个模型，根据病人用药的计量、性别和血压组别，预测出服药后病痛明显减轻的时间。
病人序号
病痛减轻时间
用药剂量
性别
血压组别
35
43
55
47
10
27
26
1.模型建立
本题共提供了三种不同的自变量，设为性别-x1，血压-x2和用药计量-x3。
建立回归模型时，应该充分考虑各个自变量对于因变量的多种影响方式，这里可能涉及到高次项和交互项。
因此应该在实验过程中将所有可能的因素列于其中，分别检验，最终得到最佳的回归模型。
基本公式如下:
2.程序设计
clc,clear
x1=[000111000111000111000111];
x2=[];
x3=[222222555555777777101010101010];
y=[354355474357262728292229
升级会员 