《应用多元分析》(第三版)各章附录中SAS程序说明等(DOC).docVIP

《应用多元解析》(第三版)各章附录中SAS程序的说明等(DOC) 《应用多元解析》(第三版)各章附录中SAS程序的说明等(DOC) PAGE/NUMPAGES 《应用多元解析》(第三版)各章附录中SAS程序的说明等(DOC) 《应用多元解析》（第三版）各章附录中 SAS程序的说明等 （王学民编） 附录1-1SAS的应用 例1-1.1的SAS程序：prociml; x={12345,24789,;g=inv(x); e=eigval(x); d=eigvec(x); h=det(x); t=trace(x); printgedht; 程序说明： “prociml”是一个矩阵运算的过程步；“x={12345,2478948 1 2 3 4 5 2 4 7 8 9 153020 ”是输入矩阵3 7 10 15 20 ，并赋值给变量x；inv(x) 是x 4 8 15 30 20 5 9 20 20 40 的逆矩阵函数，eigval(x) 是x的特色值函数，eigvec(x) 是x的特色向量函数，det(x) 是x 1 的行列式函数，trace(x)是x的迹函数，这些函数分别赋值给我们取的变量g,e,d,h,t；“print gedht”是打印语句，指定将gedht的值输出。 附录2-1SAS的应用 例2.3.3和例2.3.6的SAS程序： prociml; a={2-14,01-1,13-2}; b={5,-2,7}; c={412,19-3,2-325}; d=block(2,3,5); e=a*b; v=a*c*t(a); r=inv(d)*c*inv(d); printevr; 程序说明： “prociml”是一个矩阵运算的过程步；“a={2-14,01-1,13-2}”是输入矩阵 2 1 4 5 0 1 1 ，并赋值给变量a；“b={5,-2,7} ”是输入向量 2 ，并赋值给变量b；“c={4 1 3 2 7 4 1 2 12,19-3,2-325} ”是输入矩阵1 9 3，并赋值给变量 c；“d=block(2,3,5) ”是 2 3 25 输入对角阵 diag(2,3,5)，并赋值给变量 d；“e=a*b”是将a与b相乘，并赋值给变量 e； “v=a*c*t(a) ”是将a，c，a’三个矩阵相乘，并赋值给变量 v，其中t(a)是a的转置函数； “r=inv(d)*c*inv(d) ”是将d-1，c，d-1相乘，并赋值给变量 r，其中inv(d)是d的逆矩阵 函数；”printevr ”是打印语句，指定将 evr的值输出。 2 附录3-1SAS的应用 例3-1.1的SAS程序： proccorrdata=sasuser.examp3a1cov; varx1-x7; run; proccorrdata=sasuser.examp3a1nosimplecov; varx5x6x7; withx3x4; partialx1x2; run; 程序说明： Proc步是以proc开头的一组或几组语句，它以另一个proc步、data步或run语句结 束。SAS程序中的所有语句均以“；”结束。 “proccorr”调用corr过程，进行相关的相关解析计算。“data=sasuser.examp3a1” 规定过程解析的是sasuser逻辑库中的examp3a1数据集；选项“cov”要求输出样本协差阵； “varx1-x7”指定要解析的变量是x1至x7，缺省var语句，将对所有数值型变量进行分 析计算。上述程序中有两个proc步，第一个proc步中的语句“varx1-x7”和“run”可以省略。 第二个proc步是作偏相关解析。选项“nosimple”要求不输出简单统计量；选项“cov” 要求输出样本偏协差阵；“partialx1x2”指定在作偏相关计算时给定的变量是x1和x2； “varx5x6x7”和“withx3x4”表示计算x5、x6、x7与x3、x4之间的偏相关。该proc 步中，语句“varx5x6x7”可以省略。 附录4-1SAS的应用 例4.6.1的SAS程序： 3 procanovadata=sasuser.examp461; classg; modelx1-x4=g; manovah=g; run; 程序说明： “procanova”是一个（均衡数据的）方差解析过程；“data=sasuser.examp461”规 定过程解析的是sasuser逻辑库中的examp461数据集；class语句规定了分组变量是g；model语句规定了模型的因变量是x1，x2，x3，x4，自变量是g；manova语句要求计算别的的多元统计量。