《多元回归分析第一次上机》由会员分享,可在线阅读,更多相关《多元回归分析第一次上机(5页珍藏版)》请在金锄头文库上搜索。
1、第一次上机题目第一次上机题目1.对某地区农村的 6 名 2 周岁男婴的身高、胸围、上半臂围进行测量,得样本数据如表 1 所示。根据以往资料,该地区城市 2 周岁男婴的这三个指标的均值=,现欲用在多元正态性假定下检验该地区农村男婴是否与城市0/(90, 58, 16)男婴由相同的均值。编号身高(cm)胸围(cm)上半臂围(cm)17860.616.527658.112.539263.214.548159.014.058160.815.568459.514.0data d321;input num x1-x3;cards;1 78 60.6 16.52 76 58.1 12.53 92 63.2
2、14.54 81 59.0 14.05 81 60.8 15.56 84 59.5 14.0;proc iml;n=6; p=3;m0=90 58 16;use d321;xa=x1 x2 x3;read all var xa into x; print x;ln=6 1 ; print ln;x0=(ln*x)/n; print x0;xm=x0-m0; print xm;mm=i(6)-j(6,6,1)/n;a=x*mm*x; print a;ai=inv(a); print ai;dd=xm*ai*xm; d2=(n-1)*dd;t2=n*d2;f=(n-p)*t2/(n-1)*p);p
3、rint dd d2 t2 f;p0=1-probf(f,p,n-p);print p0;fa=finv(0.95,p,n-p);beta=probf(fa,p,n-p,t2);print fa beta;run;解:由图可知检验的 P 值小于 0.05,故 H0不相容,即农村男婴三个指标的均值与城市男婴的三个指标均值不同。2.根据经验,今天和昨天的湿度差 X1 及今天的压温差(气压与温度之差)X2是预报明天下雨或不下雨的两个重要因素。现有一批已收集的数据资料,如表2 所示。今测得 x1=8.1,x2=2.0,试问预报明天下雨还是预报明天不下雨?用距离判别来得到你所需要的结论。雨天非雨天X1(
4、湿度差)X2(压温差)X1(湿度差)X2(压温差)-1.9-6.95.25.07.36.80.9-12.51.53.83.210.42.02.50.012.7-15.4-2.51.36.80.2-0.10.42.72.1-4.6-1.7-2.62.6-2.80.27.514.68.30.84.310.913.112.810.0表 2 湿度差与压温差数据data ds511;input x1-x2 group $;cards;-1.9 3.2 A-6.9 10.4 A5.2 2 A5 2.5 A7.3 0 A6.8 12.7 A0.9 -15.4 A-12.5 -2.5 A1.5 1.3 A3.
5、8 6.8 A0.2 0.2 B-0.1 7.5 B0.4 14.6 B2.7 8.3 B2.1 0.8 B-4.6 4.3 B-1.7 10.9 B-2.6 13.1 B2.6 12.8 B-2.8 10 B;data d511test;input x1-x2 group $;cards;8.1 2.0 .;proc discrim data=ds511 testdata=d511test list testlist distance;class group;var x1-x2;run;线性判别函数为Y1=-0.06472+0.03343X1+0.04699X2Y2=-0.81562-0.03959X1+0.19590X2回判结果给出对来自 A 的 2 号和 6 号,B 的 1 号和 5 号样品判错了,对待判样品判断结果为:样品属于 A,即明天是雨天。