怎样用SAS做简单的矩阵运算？

shiyiming

譬如：在sas中怎样表示一个矩阵。有这样两个矩阵S1,S2;怎样计算
S3=((N1-1)S1+(N2-1)S2)/N1+N2-2的值；
顺便能否介绍一下，在哪里能找到有关这方面的资料，
哪里能找到SAS/IML 编程语言的介绍。
谢谢！

shiyiming

[quote="wangyz":6ee6e]譬如：在sas中怎样表示一个矩阵。有这样两个矩阵S1,S2;怎样计算
S3=((N1-1)S1+(N2-1)S2)/N1+N2-2的值；
顺便能否介绍一下，在哪里能找到有关这方面的资料，
哪里能找到SAS/IML 编程语言的介绍。
谢谢！[/quote:6ee6e]

<!-- m --><a class="postlink" href="http://v8doc.sas.com/sashtml/">http://v8doc.sas.com/sashtml/</a><!-- m -->

shiyiming

proc iml;
     x={0.17152   1.0199   0.16937  0.72754   0.41796 -0.5262   -1.9801,
       -1.6353    1.1552    0.8272  -0.02896   0.07701 -0.97931  0.58415,
       -0.64938  -1.45253 -0.4396   0.58909   0.00851  0.25128   1.69263,
        0.00179   0.22683   0.6019   1.077     0.0518    0.13931  -2.09863,
        0.04473   2.17419  -0.55693 -0.43061  0.01369 -0.61954 -0.62553,
        0.32969   1.90226   0.26921 -0.33414 -0.41633 -0.43786 -1.31283,
      -0.59835  -0.33368 -0.39323 -0.39323 -0.13517 -0.39323   2.24688};

     A={0.84178   0.48551   0.06594,
       -0.10124   0.17863   0.97860,
       -0.49755  -0.79680   0.04653,
        0.99136   0.10427  -0.02384,
        0.07030   0.95477   0.21091,
        0.38673   0.89230   0.13798,
      -0.88658  -0.29619   0.29450};

      R={1.00000   0.06793  -0.78649   0.87017   0.56464   0.74251  -0.84905,
         0.06793   1.00000  -0.04224  -0.10542   0.37211   0.25343   0.32772,
        -0.78649  -0.04224   1.00000  -0.57959  -0.73505  -0.93801   0.75294,
         0.87017  -0.10542  -0.57959   1.00000   0.15357   0.48447  -0.91902,
         0.56464   0.37211  -0.73505   0.15357   1.00000   0.87190  -0.25461,
         0.74251   0.25343  -0.93801   0.48447   0.87190   1.00000  -0.58797,
        -0.84905   0.32772   0.75294  -0.91902  -0.25461  -0.58797   1.00000};
     F=A`*inv(R)*x;
     print F;
run;