MATLAB数理统计教案解析.pptVIP

方差分析 * MATLAB和R软件 * 例8 用4种工艺生产灯泡，从各种工艺制成的灯泡中各抽出了若干个测量其寿 命，结果如下表，试推断这几种工艺制成的灯泡寿命是否有显著差异。 工艺 序号 A1 A2 A3 A4 1 1620 1580 1460 1500 2 1670 1600 1540 1550 3 1700 1640 1620 1610 4 1750 1720 1680 5 1800 解 编写程序如下： x=[1620 1580 1460 1500 1670 1600 1540 1550 1700 1640 1620 1610 1750 1720 1680 1800]; x=[x(1:4),x(16),x(5:8),x(9:11),x(12:15)]; g=[ones(1,5),2*ones(1,4),3*ones(1,3),4*ones(1,4)]; p=anova1(x,g) 求得 0.01p=0.03310.05，所以几种工艺制成的灯泡寿命有显著差异。 方差分析 * MATLAB和R软件 * 双因素方差分析:anova2 调用格式: (1)p=anova2(X) (2)p=anova2(X,reps) (3)p=anova2(X,reps,displayopt) (4)[p,table]=anova1(...) (5)[p,table,stats]=anova1(...) X=[5.5 4.5 3.5 5.5 4.5 4.0 6.0 4.0 3.0 6.5 5.0 4.0 7.0 5.5 5.0 7.0 5.0 4.5]; p=anova2(X,3) p = 0.0000 0.0001 0.7462 方差分析 * MATLAB和R软件 * 方差分析 * MATLAB和R软件 * 例9 一火箭使用了4种燃料，3种推进器作射程试验，每种燃料与每种推进器的 组合各发射火箭2次，得到结果如下： B1 B2 B3 A1 58.2,52.6 56.2,41.2 65.3,60.8 A2 49.1,42.8 54.1,50.5 51.6,48.4 A3 60.1,58.3 70.9,73.2 39.2,40.7 A4 75.8,71.5 58.2,51.0 48.7,41.4 试在水平0.05下，检验不同燃料（因素A）、不同推进器（因素B ）下的射程是 否有显著差异？交互作用是否显著？ 方差分析 * MATLAB和R软件 * 解 编写程序如下： clc,clear x0=[58.2,52.6 56.2,41.2 65.3,60.8 49.1,42.8 54.1,50.5 51.6,48.4 60.1,58.3 70.9,73.2 39.2,40.7 75.8,71.5 58.2,51.0 48.7,41.4]; x1=x0(:,1:2:5);x2=x0(:,2:2:6); for i=1:4 x(2*i-1,:)=x1(i,:); x(2*i,:)=x2(i,:); end p=anova2(x,2) 求得p=0.0035 0.0260 0.001，表明各试验均值相等的概率都为小概率，故 可拒绝均值相等假设。即认为不同燃料（因素A）、不同推进器（因素B）下的射程有显著差异，交互作用也是显著的。 回归分析 * MATLAB和R软件 * 线性回归:regress 调用格式: (1)B=regress(Y,X) (2)[B,BINT]=regress(Y,X) (3)[B,BINT,R]=regress(Y,X) (4)[B,BINT,R,RINT]= regress(Y,X) (5)[B,BINT,R,RINT,stats]= regress(Y,X) (6)[...]=regress(Y,X,alpha) 回归分析 * MATLAB和R软件 * X=[1 1 1 2 1 3 1 4 1 5 1 6 1 7 1 8 1 9 1 10]; Y=[11.1165 12.0627 13.0075 14.0352 14.9303 16.1696 17.0059 18.1797 19.0264 20.0872]; [b,bint]=regress(Y,X,0.05) b = 10.0456 1.0030 bint = 9.9165 10.1747 0.9822 1.0238 回归分析 * MATLAB和R软件 * 例10 合金的强度y 与其中的碳含量x有比较密切的关系，今从生产中收集了一批 数据如下表： x 0.10 0.11 0.12 0.13 0.14 0.15 0.16 0.17 0.18 y 42.0 41.5 45.0 45.