1、華東師范大學(xué)期中(q zhn)期末試卷(A)2012-2013學(xué)年(xunin)第二學(xué)期(xuq)課程名稱(chēng): 多元統(tǒng)計(jì)分析 學(xué)生姓名: _ 學(xué) 號(hào): _專(zhuān) 業(yè): _ 年級(jí)/班級(jí): _課程性質(zhì)：專(zhuān)業(yè)必修一二三四五六七八總分閱卷人簽名一、選擇題(每題2分，共20分).1. 程序proc cluster data=sta653 out=out653 maxc=5 list; var x1-x6; id group;run;是對(duì)數(shù)據(jù)集sta653進(jìn)行的一種統(tǒng)計(jì)過(guò)程分析，這種統(tǒng)計(jì)過(guò)程的名字是 。因子分析聚類(lèi)分析方差分析對(duì)應(yīng)分析 2. 程序proc candisc data=out653 out=can6

2、53; var x1-x6; class cluster;run; 是對(duì)數(shù)據(jù)集out653做的一種統(tǒng)計(jì)分析過(guò)程。試問(wèn)這種過(guò)程的名字是 典型判別分析因子分析典型相關(guān)分析對(duì)應(yīng)分析3. 程序proc corresp data=data902 out=out902 observed rp short; var x1-x6; id group;run;是對(duì)數(shù)據(jù)(shj)集data902做的一種統(tǒng)計(jì)分析過(guò)程(guchng)。試問(wèn)這種過(guò)程的名字是 A.聚類(lèi)分析B.因子分析C.方差分析D.對(duì)應(yīng)(duyng)分析4. 程序proc cancorr data=d1031 all vprefix=PHYS wpre

3、fix=EXER vname=生理指標(biāo) wname=訓(xùn)練指標(biāo); var weight waist pulse; with chins situps jumps;是對(duì)數(shù)據(jù)集d1031做的一種統(tǒng)計(jì)分析過(guò)程。試問(wèn)這種過(guò)程的名字是 A典型判別分析B因子分析C. 典型相關(guān)分析D.對(duì)應(yīng)分析 5. 程序proc princomp data=data711 prefix=z out=o711; var x1-x8;run;是對(duì)數(shù)據(jù)集data711做的一種統(tǒng)計(jì)分析過(guò)程。試問(wèn)這種過(guò)程的名字是 逐步回歸分析主成分分析因子分析對(duì)應(yīng)分析6. 程序 proc factor data=data808 score n=2 o

4、ut=o88; var x1-x5;run;是對(duì)數(shù)據(jù)集data808做的一類(lèi)統(tǒng)計(jì)分析過(guò)程。試問(wèn)這種過(guò)程的名字是 判別分析主成分分析因子分析對(duì)應(yīng)分析7. 程序 proc varclus data=d672 maxc=8 summary outtree=o672;run;是對(duì)數(shù)據(jù)(shj)變量的相關(guān)矩陣數(shù)據(jù)(shj)集d672做的一類(lèi)(y li)統(tǒng)計(jì)分析過(guò)程。試問(wèn)這種過(guò)程的名字是 變量聚類(lèi)分析觀測(cè)值聚類(lèi)分析列聯(lián)表分析對(duì)應(yīng)分析8. 程序 proc discrim data= d511 simple wcov pcov wsscp psscp distance list; class group; v

5、ar x1-x4;run;是對(duì)數(shù)據(jù)集d511做的一類(lèi)統(tǒng)計(jì)分析過(guò)程。試問(wèn)這種過(guò)程的名字是 變量聚類(lèi)分析觀測(cè)值聚類(lèi)分析列聯(lián)表分析判別分析9. 程序 proc discrim data=ds511 testdata=d511test list testlist; class group; var x1-x4;run;是對(duì)數(shù)據(jù)集d511做的一類(lèi)統(tǒng)計(jì)分析過(guò)程。試問(wèn)下面說(shuō)法正確的是 這是一個(gè)判別分析過(guò)程，其中數(shù)據(jù)集ds511是建立判別函數(shù)的輸入數(shù)據(jù)集，數(shù)據(jù)集d511test是輸出數(shù)據(jù)集。這是一個(gè)判別分析過(guò)程，其中數(shù)據(jù)集ds511是建立判別函數(shù)的輸入數(shù)據(jù)集，數(shù)據(jù)集d511test也是幫助建立判別函數(shù)的輸入數(shù)

6、據(jù)集。這是一個(gè)判別分析過(guò)程，其中數(shù)據(jù)集ds511是建立判別函數(shù)的輸入數(shù)據(jù)集，數(shù)據(jù)集d511test是檢驗(yàn)判別函數(shù)判別正確性的輸出數(shù)據(jù)集。這是一個(gè)判別分析過(guò)程，其中數(shù)據(jù)集ds511是建立判別函數(shù)的輸入數(shù)據(jù)集，數(shù)據(jù)集d511test是檢驗(yàn)判別函數(shù)判別正確性的輸入數(shù)據(jù)集。 10. 程序 proc reg data= d431; model y1 y2=x1-x5/selection=stepwise sle=0.05 sls=0.05;run; proc reg data=d431; model y1 y2=x3-x5;run;是對(duì)數(shù)據(jù)(shj)集d431做的一類(lèi)(y li)統(tǒng)計(jì)分析過(guò)程。試問(wèn)下面說(shuō)

7、法正確的是 這是一類(lèi)回歸分析過(guò)程，前一個(gè)回歸分析是逐步回歸分析，后一個(gè)是有兩個(gè)(lin )因變量和三個(gè)自變量的一般回歸分析。這是一類(lèi)回歸分析過(guò)程，后一個(gè)回歸分析是逐步回歸分析，前一個(gè)是有兩個(gè)因變量和三個(gè)自變量的一般回歸分析。這是一類(lèi)判別分析過(guò)程，前一個(gè)判別分析是逐步判別分析，后一個(gè)是有兩個(gè)因變量和三個(gè)自變量的一般判別分析。這是一類(lèi)判別分析過(guò)程，前一個(gè)判別分析是逐步判別分析，后一個(gè)是有兩個(gè)因變量和三個(gè)自變量的一般判別分析。 二、填空題(每空2分，共20分) 1 設(shè)為服從某一總體的p維隨機(jī)向量，期望為EX,方差為而為服從另一總體的q維隨機(jī)向量，期望為EY,方差為 X和Y的協(xié)方差為。試寫(xiě)出X和Y的兩

8、種相關(guān)矩陣 2. 設(shè)為來(lái)自p元正態(tài)總體的隨機(jī)樣本，又設(shè)為來(lái)自q元正態(tài)總體的隨機(jī)樣本， 記，又記,其中是階單位矩陣，是元素全是1的維向量，用寫(xiě)出方差和協(xié)方差的矩估計(jì)形式 的矩估計(jì)形式是 的矩估計(jì)形式是 3設(shè)為來(lái)自(li z)p元正態(tài)總體(zngt)的隨機(jī)樣本，記試寫(xiě)出因子分析(fnx)和主成分分析中的主因子和主成分公式？ 因子分析中的主因子公式 ；主成分分析中主成分公式 ；4設(shè)為來(lái)自p元正態(tài)總體的隨機(jī)樣本，又設(shè)是階投影矩陣，其秩滿足，記，試寫(xiě)出矩陣和的分布函數(shù)？矩陣的分布函數(shù) ；矩陣的分布函數(shù) ；5 試寫(xiě)出關(guān)聯(lián)分析和相干分析的數(shù)學(xué)定義？關(guān)聯(lián)分析的數(shù)學(xué)定義 ；相干分析的數(shù)學(xué)定義 ；三、改錯(cuò)題(每題

9、2分，共20分) 閱讀如下程序，找出其中的10個(gè)錯(cuò)誤 %let n=5; %let m=10;data ds511; input x1-x4 group $; cards;13.85 2.79 7.80 49.60 A22.31 4.67 12.31 47.80 A28.82 4.63 16.18 62.15 A15.29 3.54 7.50 43.20 A28.79 4.90 16.12 58.10 A 2.18 1.06 1.22 20.60 B 3.85 0.80 4.06 47.10 B11.40 0.00 3.50 0.00 B 3.66 2.42 2.14 15.10 B12.10

10、 0.00 5.68 0.00 B;data d511test; input x1-x4 group $; cards; 8.85 3.38 5.17 26.10 .28.60 2.40 1.20 127.0 .20.70 6.70 7.60 30.20 . 7.90 2.40 4.30 33.20 . 3.19 3.20 1.43 9.90 .12.40 5.10 4.43 24.60 .16.80 3.40 2.31 31.30 .15.00 2.70 5.02 64.00 .;data a;set ds511-d511test;u1=.;u2=.;if group=A then do;u

11、1=100; u2=0;end;if group=B then do; u1=0;u2=100; end;if group=. then do;u2=.; u1=.;end;data aaa;set a ;x5=(x1*x2)*(1/2);x6=(x1*x3)*(1/2);x7=(x1*x4)*(1/2);x8=(x2*x3)*(1/2);x9=(x2*x4)*(1/2);x10=(x3*x4)*(1/2);x11=(x1*x3*x4)*(1/3);x12=(x1*x2*x3*x4)*(1/4);array mx&m x1-x&m;do i=1 to &m;mxi=round(mxi,0.01

12、);end;test=_n_;proc factor data=aaa n=&n out=bbb outstat=bbbccc ;class x1-x&m;data bbb000;set bbb;array mfactor&n factor1-factor&n;array mf&n f1-f&n; test=_n_;ttt=0;do i=1 to &n;mfi=sign(mfactori);ttt=ttt+(mfi+1)/2)*2*(&n-i);end; drop factor1-factor&m;%macro fx;%do iii=1 %to &n;proc means out=bbb000

13、 noprint;var u1 u2;class f&iii;output out=ccc mean=m&iii.u1 m&iii.u2; data ccc&iii; set ccc;if _n_=1 then do; call symput(mu1,m&iii.u1); call symput(mu2,m&iii.u2);delete;end; keep m&iii.u1 m&iii.u2;%end;proc iml;x=1,2;%do iii=1 %to &n;xxx=j(2*(&iii-1),1,1)x j(2*(&n-&iii),1,1);de= de | xxx; %end;yyy1

14、=j(2*(&n),&n,1);yyy2=j(2*(&n),&n,1);mu1=&mu1;mu2=&mu2; %do iii=1 %to &n;use ccc&iii;read all into ddd&iii;yyy1 ,&iii=ddd&iiide,&iii,1-mu1;yyy2 ,&iii=ddd&iiide,&iii,2-mu2;%end;zzz1=mu1+yyy1,+;zzz2=mu2+yyy2,+;mv1=zzz1:;mv2=zzz2:;zzz1=zzz1-zzz1:;zzz2=zzz2-zzz2:;uuu1=(mu1+mv1)/2;uuu2=(mu2+mv2)/2;%do iii

15、=1 %to &n;c=design(de,&iii);AAA=c*ginv(c*c)*c-j(nrow(c),nrow(c),1)/nrow(c);uuu1=uuu1+aaa*zzz1;uuu2=uuu2+aaa*zzz2;%end;mvv1=uuu1:;mvv2=uuu2:;uuu1=uuu1-uuu1:;uuu2=uuu2-uuu2:;www1=(mvv1+mv1)/2;www2=(mvv2+mv2)/2;%do iii=1 %to &n;c=design(de,&iii);AAA=c*ginv(c*c)*c-j(nrow(c),nrow(c),1)/nrow(c);www1=www1+

16、aaa*uuu1;www2=www2+aaa*uuu2;%end;mvvv1=www1:;mvvv2=www2:;www1=www1-www1:;www2=www2-www2:;zzz1=(mvvv1+mvv1)/2;zzz2=(mvvv2+mvv2)/2;%do iii=1 %to &n;c=design(de,&iii);AAA=c*ginv(c*c)*c-j(nrow(c),ncol(c),1)/nrow(c);zzz1=zzz1+aaa*www1;zzz2=zzz2+aaa*www2;%end;mvvvv1=zzz1:;mvvvv2=zzz2:;zzz1=zzz1-zzz1:;zzz2

17、=zzz2-zzz2:;www1=(mvvvv1+mvvv1)/2;www2=(mvvvv2+mvvv2)/2;%do iii=1 %to &n;c=design(de,&iii);AAA=c*ginv(c*c)*c-j(nrow(c),nrow(c),1)/nrow(c);www1=www1+aaa*zzz1;www2=www2+aaa*zzz2;%end;mvvvvv1=www1:;mvvvvv2=www2:;www1=www1-www1:;www2=www2-www2:;zzz1=(mvvvvv1+mvvvv1)/2;zzz2=(mvvvvv2+mvvvv2)/2;%do iii=1 %

18、to &n;c=design(de,&iii);AAA=c*ginv(c*c)*c-j(nrow(c),nrow(c),1)/nrow(c);zzz1=zzz1+aaa*www1;zzz2=zzz2+aaa*www2;%end;mvvvvvv1=zzz1:;mvvvvvv2=zzz2:;zzz1=zzz1-zzz1:;zzz2=zzz2-zzz2:;www1=(mvvvvvv1+mvvvvv1)/2;www2=(mvvvvvv2+mvvvvv2)/2;%do iii=1 %to &n;c=design(de,&iii);AAA=c*ginv(c*c)*c-j(nrow(c),nrow(c),1

19、)/nrow(c);ddd&iiide,&iii,1=aaa*zzz1;ddd&iiide,&iii,2=aaa*zzz2;www1=www1+aaa*zzz1;www2=www2+aaa*zzz2;%end;nnn=0;%do iii=1 %to &n; nnn=nnn+(de,&iii-1)*2*(&n-&iii); ddd=ddd|round(ddd&iii,0.1); %end;www1=round(www1,0.01);www2=round(www2,0.01);rf=(www1www2);www= www1|www2 |nnn|rf;print www1 mvvvvvv1 mvvv

20、vv1, mu1 www2 mvvvvvv2 mvvvvv2 mu2;print www;print ddd;create www from www;append from www;close www;quit;%mend;%fx;data www;set www;col=col1-col2; rf=(col1col2);proc sort; by col3;%macro test;data zzz;%do jjj=1 %to 18;data bbbxxx;set bbb000;col3=ttt; if test=&jjj then call symput(&ttt,ttt);data www

21、; set www;if col3=&ttt thendo; %put &jjj;put rf rf rf rf rf col3 col1 col2;call symput(rf, rf);call symput(gua, col3);call symput(energe1, col1);call symput(energe2, col2); end;data ccc;jjj=&jjj; rf=&rf; gua=&gua; energe1=&energe1; energe2=&energe2;data zzz set zzz ccc;%end;%mend;%test;data www;set www;energe1=col1; energe2=col2; gua=&col3;keep rf gua energe1-en