多元回归分析第一次上机

1、第一次上机题目1. 对某地区农村的6名2周岁男婴的身高、胸围、上半臂围进行测量，得样本数 据如表1所示。根据以往资料，该地区城市 2周岁男婴的这三个指标的均值%=（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 7860.616.52 7658.112.53 9263.214.54 8159.014.05 816

2、0.815.56 8459.514.0proc 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 In; 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); print dd d2 t2 f; p

3、O=1-probf(f,p,n-p); print pO;fa=finv(0.95,p,n-p); beta=probf(fa,p,n-p,t2); print fa beta;run;解:Al0.0372fi 0.12&37S 0.0773292 -0.126378 0.516S014 -0.323064 0.0773292 -0.3230S4 0.3076572DDD2T2F14,014822 70.07411 420.44466 84.088332PO0.0021553FA BETA9.27S6282 5.5885E-9由图可知检验的 P值小于0.05，故H0不相容，即农村男婴 三个指标

4、的均值与城市男婴的三个指标均值不同。2. 根据经验，今天和昨天的湿度差 X1及今天的压温差（气压与温度之差）X2是预报明天下雨或不下雨的两个重要因素。现有一批已收集的数据资料，如表2所示。今测得x1=8.1，x2=2.0，试问预报明天下雨还是预报明天不下雨？用距 离判别来得到你所需要的结论。雨天非雨天X1 （湿度差）X2 （压温差）X1 （湿度差）X2 （压温差）-0.2-6.910.4-2.00.414.65.07.30.012.7--15.4-1.710.9-12.5-2.5-2.613.11.51.

5、32.6-2.810.0表2湿度差与压温差数据data ds511;input x1-x2 group $;cards;-1.93.2A-6.910.4A5.22A52.5A7.30 A6.812.7A0.9-15.4A-12.5-2.5 A1.51.3A3.86.8A0.20.2B-0.17.5B0.414.6B2.78.3B2.10.8 B-4.64.3 B-1.710.9 B-2.613.1 B2.612.8B-2.810Bdata d511test;input x1-x2 group $;cards;8.1 2.0 .proc discrim data=ds511

6、testdata=d511test list testlist distance;class group;var x1-x2;run;Squared DIstan匚呂 ta groupFromgroupABAB01.D1074L010740X?F Stat ist icsT NDF-Z, DDF-17 for Squered Distance to groupFrom groupABA02.3864?B2JS6470Prob Mahdi knob i 3 D1 stanceforSqua red D i st ance to groupFromgroupABLOOOO0.1220(JJ2201

7、.0000L i near Duscr immantFunct i on foirgroupVari ab1eAB-0.81562-0.039580.19590Constant-0.06472xl0.033430.0469SPosterior Probabi Iity of Membership in &roupFrom groupCI assi fledi nto group1AA0.533S0.46612Ae0.2139DJ8G13AAD.eaea8.30314AA0.G7780.32225AA0.78310.2169SA0*0.34440.G55S7AA1).86733.1)4278AA

8、0.5&240.447B3AAD.BGOB0.39921UAA0.60940.488111BAD.B7B1CJ239EB0.40780.582213BBD.ldflB0.801214EB0.42350.5715!5eAD.eeee0,313215BB0.443S0.556117BB0,28970.790318BB0.19950.8005BB0,27580.724220EB0.28040.719G$ Miscl ass i f i ed observationPosterior rooabi i it of wieift&ersrup in sroupFromClassifiedObs groupinto groupAB1A$0J3370.2S03* N i sclossi f ied observat(on