sas多元统计例题程序.ppt

《sas多元统计例题程序.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、procprincompdata=duo.innovationprefix=zout=o;run;procplotdata=o;plotz2*z1$region='*'/href=0vref=0;run;procsortdata=o;bydescendingz1;run;procprintdata=o;varz1z2;run;quit;主成分分析在SAS中用princomp过程：plot过程已经很熟悉了。href=0表示在横坐标z1=0处画一条垂线，vref=0表示在纵坐标z2=0处画一条垂线。$region

2、=‘*’表示每个点在图上用*表示，并且在*后显示该样本点的region变量的值。主成分分析有一个princomp过程就足够了。prefix=z表示，在输出数据集中（o中），主成分变量是z1、z2、…sort和print过程也是很熟悉的过程了。procstandarddata=vregex1out=svmean=0std=1;varx1-x3y;run;procprincompdata=svprefix=zout=opcr;varx1x2x3;run;procprintdata=opcr;varz1z2y;ru

3、n;procregdata=opcr;modely=z1z2;run;quit;procregdata=vregex1outest=out;modely=x1-x3/pcomit=1,2;run;quit;procprintdata=out;run;首先对数据标准化然后对标准化后的数据进行主成分分析。这也就相当于是对相关系数进行的主成分分析。作回归分析，自变量取第1、2主成分，因变量为y。直接作主成分回归分析。pcomit=1,2表示分别作两个回归，分别是剔除1个主成分，和剔除2个主成分，所做的主成分回归。主

4、成分回归的结果:由于刚才我们分析了y与z1、z2的回归，z1、z2的回归系数非0都是99%显著的。而z1、z2又有99%的累积贡献率。因此我们可以认为我们对y与z1、z2、z3之间的回归，取y对z1、z2的回归最佳。y对z1、z2的回归，还原为x1、x2、x3的系数后就是上图中红线圈起的第2行。这是看_PCOMIT_=1知道的，这表示，这一行是主成分回归，而且删除了最不重要的一个主成分后，回归的结果。而y与z1回归，再还原为x1、x2、x3的系数，就是上图中最后一行。因这一行_PCOMIT_=2。procpr

5、incompdata=文件名n=？out=文件名1outstat=文件名2covprefix=？;run;SAS程序data=文件名prefix=n=说明指定主成分的个数指定用于分析的数据文件指定主成分的代号out=outstat=文件名自己取，保存原始数据和主成分得分文件名自己取，保存主成分分析过程中的统计量特别注意，分号表示一个语句的结束，不能遗漏。cov指定基于协差阵分析，若省略，基于相关系数矩阵procfactordata=文件名r=vn=？out=文件名1outstat=文件名2;run;SAS程序

6、data=文件名r=vn=说明指定公共因子的个数指定用于分析的数据文件指定旋转的方法，是rotate=varimax的缩写，也可以用rotate=quartimax(r=q)out=outstat=文件名自己取，保存原始数据和因子得分文件名自己取，保存因子分析过程中的统计量特别注意，分号表示一个语句的结束，不能遗漏。6procfactordata=tmp1.innovationr=vn=4out=out;run;dataa1;setout;f=(5.8617*factor1+2.6989*factor2+2.

7、0995*factor3+1.5842*factor4)/14;run;dataa2;seta1;keepregionffactor1factor2factor3factor4;run;procsortdata=a2;bydescendingf;run;procprintdata=a2;run;Procclustermethod=选项data=文件名outtree=文件名1standard;varvariable-list;idvariable;run;Proctreedata=文件名1horizontalg

8、raphics;idvariable;run;Method＝选项single最短距离法complete最长距离法median中间距离法centroid重心法average类平均法ward离差平方和法（Ward法）聚类分析SAS程序Procfastclusmaxcluster=data=文件名seep=文件名mean=文件名out=文件名;varvariable-list;run;Maxcluster