[理學(xué)]第2章1均向量的統(tǒng)計推斷-資料下載頁
2024-10-19 00:51本頁面
【正文】 ? proc corr csscp out=ssa noprint。 ? var ? %do i=1 %to amp。m。 ? amp。amp。varamp。i ? %end。 ? 。 ? by amp。groupvar。 ? proc iml。 ? use amp。dataname。 ? read all var {amp。groupvar} into gv。 ? min=gv[,]。 ? num=min。 ? nu=0。 ? ss=j(amp。m,amp。m,0)。 ? %do i=1 %to amp。groupnum。 ? ssamp。i=j(amp。m,1,0)。 ? %end。 ? use ssa。 ? %do i=1 %to amp。groupnum。 ? %do k=1 %to amp。m。 ? read all var {amp。amp。varamp。k} where ( _TYPE_=CSSCP amp。 amp。groupvar=num) into ssamp。iamp。k。 ? ssamp。i=ssamp。i||ssamp。iamp。k。 ? %end。 ? ssamp。i=ssamp。i[,2:%eval(amp。m+1)]。 ? read all var {amp。var1} into namp。i where ( _TYPE_=N amp。 amp。groupvar=num)。 ? nu=nu+namp。i。 ? vamp。i=ssamp。i/(namp。i1)。 ? ss=ss+ssamp。i。 ? lvamp。i=log(det(vamp。i))。 ? num=num+1。 ? %end。 ? vv=ss/(nuamp。groupnum)。 ? lvv=log(det(vv))。 ? mm=(nuamp。groupnum)*lvv。 ? %do i=1 %to amp。groupnum。 ? mm=mm(namp。i1)*lvamp。i。 ? %end。 ? d=(2*amp。m*amp。m+3*amp。m1)/(6*(amp。m+1)*(amp。groupnum1))。 ? e=0。 ? %do i=1 %to amp。groupnum。 ? e=e+1/(namp。i1)。 ? %end。 ? d1=d*(e1/(nuamp。groupnum))。 ? df=amp。m*(amp。m+1)*(amp。groupnum1)/2。 ? chi=mm*(1d1)。 ? pp=1probchi(chi,df)。 ? print d1。 ? print chi。 ? print df。 ? print pp。 ? quit。 ? %mend。 ? %testcov(a2,x1/x2/x3,g,3)
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號-1