SAS中文论坛

标题: 怎样用SAS做简单的矩阵运算？ [打印本页]

作者: shiyiming 时间: 2004-6-5 11:03
标题: 怎样用SAS做简单的矩阵运算？
譬如：在sas中怎样表示一个矩阵。有这样两个矩阵S1,S2;怎样计算
S3=((N1-1)S1+(N2-1)S2)/N1+N2-2的值；
顺便能否介绍一下，在哪里能找到有关这方面的资料，
哪里能找到SAS/IML 编程语言的介绍。
谢谢！

作者: shiyiming 时间: 2004-6-5 19:04
标题: Re: 弱弱的问：怎样用sas做简单的矩阵运算？
[quote="wangyz":6ee6e]譬如：在sas中怎样表示一个矩阵。有这样两个矩阵S1,S2;怎样计算
S3=((N1-1)S1+(N2-1)S2)/N1+N2-2的值；
顺便能否介绍一下，在哪里能找到有关这方面的资料，
哪里能找到SAS/IML 编程语言的介绍。
谢谢！[/quote:6ee6e]

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

作者: shiyiming 时间: 2004-6-12 14:46
标题: 这是我作业的一个程序，可以参考下！
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;

欢迎光临 SAS中文论坛 (http://mysas.net/forum/)