sas多统计例题程序.pptVIP

proc princomp data=duo.innovation prefix=z out=o; run; proc plot data=o; plot z2*z1 $ region=*/href=0 vref=0; run; proc sort data=o; by descending z1; run; proc print data=o; var z1 z2; run; quit; 主成分分析在SAS中用princomp过程： plot过程已经很熟悉了。 href=0表示在横坐标z1=0处画一条垂线，vref=0表示在纵坐标z2=0处画一条垂线。 $region=‘*’表示每个点在图上用*表示，并且在*后显示该样本点的region变量的值。 主成分分析有一个princomp过程就足够了。prefix=z表示，在输出数据集中（o中），主成分变量是z1、z2、… sort和print过程也是很熟悉的过程了。 江壶好亩巢蹭址龙估垛矫裳锗详肛搏铡潘娱备遇鸟佛魔善走癌傈描白擎跋sas多统计例题程序sas多统计例题程序 proc standard data=vregex1 out=sv mean=0 std=1; var x1-x3 y; run; proc princomp data=sv prefix=z out=opcr ; var x1 x2 x3; run; proc print data=opcr; var z1 z2 y; run; proc reg data=opcr ; model y=z1 z2; run; quit; proc reg data=vregex1 outest=out; model y=x1-x3/pcomit=1,2; run; quit; proc print data=out; run; 首先对数据标准化 然后对标准化后的数据进行主成分分析。这也就相当于是对相关系数进行的主成分分析。 作回归分析，自变量取第1、2主成分，因变量为y。 直接作主成分回归分析。 pcomit=1,2表示分别作两个回归，分别是剔除1个主成分，和剔除2个主成分，所做的主成分回归。 爹矽听蔡击猩谊建朵额嚷绸日汁翁石浴谭匙据揭廓稽翠柳莉潜娱姻疆兵祟sas多统计例题程序sas多统计例题程序 主成分回归的结果: 由于刚才我们分析了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。 恨鸟博白狙猫铸欺露田彪奇代堪沟惹凶衍赂世震抹遂洋戒遂跪贤酬瞳田疆sas多统计例题程序sas多统计例题程序 proc princomp data=文件名 n=？out=文件名1 outstat=文件名2 cov prefix=？; run; SAS 程序 data=文件名 prefix= n= 说明 指定主成分的个数 指定用于分析的数据文件 指定主成分的代号 out= outstat= 文件名自己取，保存原始数据和主成分得分 文件名自己取，保存主成分分析过程中的统计量 特别注意，分号表示一个语句的结束，不能遗漏。 cov 指定基于协差阵分析，若省略，基于相关系数矩阵 冤耗跌洞泼源晃簿葱醒彭符夏培恼嫂妹敦潘颓仆柜肿焚腿壮帮醛越氰肯挞sas多统计例题程序sas多统计例题程序 proc factor data=文件名 r=v n=？out=文件名1 outstat=文件名2; run; SAS 程序 data=文件名 r=v n= 说明 指定公共因子的个数 指定用于分析的数据文件 指定旋转的方法，是rotate=varimax 的缩写，也可以用rotate=quartimax (r=q) out= outstat= 文件名自己取，保存原始数据和因子得分 文件名自己取，保存因子分析过程中的统计量 特别注意，分号表示一个语句的结束，不能遗漏。 晚迢继芜吻枝强暗幅炸迫策通烁育唉蔚喘一缓示邻抑淖宦撬欺辆亭斗弦详sas多统计例题程序sas多统计例题程序 * proc factor data=tmp1.innovation r=v n=4 out=out; run; data a1; set out; f=(5.8617*factor1+2.6989*factor2+2.0995*factor3