2.3（编程，验证）模型应用p337~338
在题2.1的模型中增加x2与x3，x4的交互项后，多元线性回归模型为
y =α0 + α1 x1 + α2 x2 + α3 x3 + α4 x4 + α5 x2 x3 + α6 x2 x4 + ε
(1)（编程）同题2.2，但要去除题2.1中给出数据中编号为33的数据（异常点）
★(1) 给出程序和运行结果（数据和图形）（比较[337]表5、图5、图6）：
程序： .2 软件开发人员的薪金——基本模型 %模型：y=α0+α1*x1+α2*x2+α3*x3+α4*x4+ε clear;clc;format compact;format short g; M=dlmread('p333.txt'); %读取ASCII码文件 M(33,:)=[]; y=M(:,2); x1=M(:,3); x2=M(:,4); x3=M(:,5)==1; x4=M(:,5)==2; %教育程度 [a,aint,r,rint,stats]=regress(y,[ones(size(M,1),1) x1 x2 x3 x4 x2.*x3 x2.*x4]); 16
fprintf('%2s%4s%9s\\n','参数','估计值','置信区间');%1个汉字算1个字符 for i=1:length(a) fprintf('α %7.0f [%5.0f, %5.0f]\\n',i-1,[a(i,:),aint(i,:)]); end fprintf('\\nR2=%.3f F=%.0f p<%.4e s2=%.3e\\n',stats); subplot(121); plot(x1,r,'+'); %模型的残差ε与资历x1的关系 subplot(122); x234=(1+M(:,4)).*(M(:,5)==1)+(3+M(:,4)).*(M(:,5)==2)... + (5+M(:,4)).*(M(:,5)==3);%见p336表3 plot(x234,r,'+'); %模型的残差ε与管理-教育x2-x3,x4组合x234的关系 数据结果（比较p337表5）：
图形结果（比较p337图5、图6）： 17
(2)（验证） 6种管理—教育组合人员的“基础”薪金（资历为0）。
组合 管 理 教 育 “基础”薪金 资历 0 1 ？ 1 0 1 1 ？ 2 0 0 2 ？ 3 0 1 2 ？ 4 0 0 3 ？ 5 0 1 3 ？ 6 0 利用(1)所求得的模型求解，程序如下（与教材p337表6比较）： clear; clc; format compact; format short g; M=dlmread('p333.txt'); %读取ASCII码文件 M(33,:)=[]; %修改 y=M(:,2); x1=M(:,3); x2=M(:,4); x3=M(:,5)==1; x4=M(:,5)==2; %教育程度 [a,aint,r,rint,stats]=regress(y,[ones(size(M,1),1) x1 x2 x3 x4 x2.*x3 x2.*x4]); x1=zeros(6,1); %资历为0 x2=[0 1 0 1 0 1]'; %管理 x3=[1 1 0 0 0 0]'; %教育 x4=[0 0 1 1 0 0]'; 18
y=[ones(6,1),x1,x2,x3,x4,x2.*x3,x2.*x4]*a; fprintf('%2s%4s%4s%8s\\n', '组合', '管理', '教育', '“基础”薪金');%1个汉字算1个字符 for k=1:6 fprintf('=mm.0f\\n',k,x2(k),fix((k+1)/2),y(k)); end ☆(2) 运行程序并给出结果（比较[337]表6）：
3. 酶促反应p338~345
嘌呤霉素实验中的反应速度与底物浓度数据 底物浓度 0.02 0.06 0.11 0.22 0.56 1.10 (ppm) 反应 处理 76 47 97 107 123 139 159 152 191 201 207 200 速度 未处理 67 51 84 86 98 115 131 124 144 158 160 / 3.1（编程）线性化模型p338~340
(1) 根据给出的数据，试编写一个程序绘制出如下图所示的散点图（比较：p339图1、图2）
19
图1 y对x（经处理）的散点图220200180160140120100806040220200180160140120100806040 图2 y对x（未经处理）的散点图00.20.40.60.811.21.400.20.40.60.811.21.4
★(1) 给出程序和运行结果（比较[339]图1、图2）：
程序： clear; clc; format compact; format short g; M=dlmread('p338.txt'); %读取ASCII码文件 y=M(2,1:2:11);y2=M(2,2:2:12);x1=M(1,1:6);x2=M(1,1:5); y4=M(3,1:2:11) ;y3=M(3,2:2:10); subplot(1,2,1) plot(x1,y, 'bo'); hold on; plot(x1,y2, 'bo'); title('y对x经处理的散点图'); subplot(1,2,2) plot(x1,y4, '+'); hold on; plot(x2,y3, '+'); title('y对x未经处理的散点图'); 图形（比较p339图1、图2）： y对x经处理的散点图220200140180160140100120100806060404080120160y对x未经处理的散点图00.20.40.60.811.21.400.20.40.60.811.21.4 20