MATLAB在统计中的应用.pptVIP

MATLAB在数理统计中的应用 汪晓银 华中农业大学 x=[21,29,24,22,25,30,27,26,20,25,25,23,29,31,24,26,20,21,24,22,28,25,21,26]; g=[ones(1,8),2*ones(1,10),3*ones(1,6)]; [p,d,f]=anova1(x,g) finv(0.95,2,21) 列出方差分析表 X=[19, 23, 21,13; 21,24,27,20; 20,18,19,15; 22,25,27,22]; [p,d,f]=anova1(X) 双因素试验的方差分析（不交互） 例2.1《产量分析》某5名工人在4台机床上分别各工作了一天，其日产量(单位：件)服从正态分布，试作方差分析，估计各个总体的未知参数μi.和μ.j 及μ，如有必要，试求出两两机床机床日产量的总体均值差的双侧0.95置信区间。观测值如下： 双因素试验的方差分析（不交互） x=[53 56 45 52 49 47 50 47 47 53 57 63 54 57 58 45 52 42 41 48]; [p,t,st]=anova2(x) 注意 row指行，column指的是列 双因素试验的方差分析（不交互） 双因素试验的方差分析（交互） 例3.1《火箭试验》设火箭的射程(单位：km)在其它条件基本相同时与燃料的种类及推进器的型号有关，现有4种燃料及3种推进器搭配试验的数据如下表。试作方差分析 双因素试验的方差分析（交互） 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,t,st]=anova2(x,2) 回归分析 回归分析 x1=[1.5 1.8 2.4 3.0 3.5 3.9 4.4 4.8 5.0]; y=[4.8 5.7 7.0 8.3 10.9 12.4 13.1 13.6 15.3]; x=[ones(9,1),x1] [b,bint,r,rint,stats]=regress(y,x) 回归分析 x1=[1.5 1.8 2.4 3.0 3.5 3.9 4.4 4.8 5.0]; y=[4.8 5.7 7.0 8.3 10.9 12.4 13.1 13.6 15.3]; x=[ones(9,1),x1];[b,bint,r,rint,stats]=regress(y,x); y1=0.2569*ones(9,1)+2.9303*x1;y0=0.2569+2.9303*2 sse=sum((y-y1).^2);llx=sum(x1.*x1)-9*mean(x1)^2 delta=tinv(0.975,7)*sqrt((sse/7)*(1+1/9+(2-mean(x1))^2/llx)) a=y0-delta,b=y0+delta;[a,b] 例2.8《药品疗效》为了解某种药品对某种疾病的疗效是否与患者的年龄有关，查到了91名老年患者、100名中青年患者、109名儿童患者的疗效记录，整理后列表如下： 疗效 老年 中青年 儿童 行求和 ────────────────── 显著 32 38 58 128 一般 45 44 28 117 较差 14 18 23 55 ────────────────── 列求和 91 100 109 300 试检验年龄与疗效相互独立。 x=[32,38,58;45,44,28;14,18,23]; x1=x(1,:);x2=x(2,:);x3=x(3,:); y1=x(:,1);y2=x(:,2);y3=x(:,3); h1=sum(x1);h2=sum(x2);h3=sum(x3); le1=sum(y1);le2=sum(y2);le3=sum(y3); n=h1+h2+h3; ph1=h1/