请教高手一个IML的问题

shiyiming

我现在用IML的I函数产生了一个单位矩阵：
proc iml;
i3=i(3);
print i3;
quit;
结果如下：

1 0 0
0 1 0
0 0 1
；
我现在想把这个结果导成SAS数据集，如下格式：

x y z
1 0 0
0 1 0
0 0 1

怎么实现？多谢高手指点！

shiyiming

[code:j31eg3ux]proc iml;
        i3=i(3);
        print i3;
        create result from i3;
        append from i3;
quit;
[/code:j31eg3ux]

by the way，可以直接存储矩阵！

shiyiming

多谢指点，非常感谢！能否接着问一下：
如果数据集result已经有了，如下：
x y z
1 0 0
0 1 0
0 0 1
;
我现在想把这个SAS数据集转换成对应矩阵的三个列向量X1;Y1;Z1.
即：X1=(1,0,0);Y1=(0,1,0);Z1=(0,0,1).请问怎么实现？

shiyiming

proc iml;
use result;
read all into x;
x1=x[1:3,1];
x2=x[1:3,2];
x3=x[1:3,3];
print x1 x2 x3;
quit;

shiyiming

非常感谢，算是解决了SAS数据集和矩阵之间的数据转换。

shiyiming

非常佩服！受教了！什么是IML？

shiyiming

徐大人还不知道SAS有这个模块啊。专门用来做矩阵分析的，很好很强大。不过很多数据处理的功能其实base sas也可以搞定。只是如果你有时候必须要和矩阵打交道的时候，那只好用iml或者类似的如matlab之类的工具了。