全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 SAS专版
1439 3
2010-01-15
有关sas软件中如何将矩阵转化为sas数据集的问题
proc iml;
A={2.684 2.876 2.793 0.536 2.010,
  2.694 2.924 2.823 0.563 2.487,
  2.737 2.981 2.846 0.633 2.709,
  2.698 2.965 2.850 0.608 2.789,
  2.601 2.878 2.775 0.549 2.855,
  2.647 2.893 2.840 0.571 2.913,
  2.677 2.942 2.836 0.595 2.965,
  2.679 2.967 2.882 0.605 3.052,
  2.714 3.029 2.903 0.649 3.090,
  2.711 3.027 2.915 0.636 3.124};
B={1.763,
   2.127,
   2.248,
   2.267,
   2.270,
   2.274,
   2.318,
   2.405,
   2.441,
   2.462};
X=STANDARD(A);
Y=STANDARD(B);
Z=t(X)*X;
val=eigval(Z);
vet=eigvec(Z);
print X Y Z val vet;
F1=X*vet[,1];
F2=X*vet[,2];
F3=X*vet[,3];
Q=(val[1,]+val[2,]+val[3,])/(10*5);
print F1 F2 F3 Q;
quit;run;
这是我在sas软件上编的一个程序,现在我想做B对F1 F2 F3的线性回归。这需要先把四个矩阵转化为一个数据集的形式。请教各位高人,帮帮忙,查了很多书都没有找到啊,谢谢了
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2010-1-15 16:18:59
不懂~~~~~~~~~
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2010-1-16 11:23:53
zhangxl121 发表于 2010-1-15 16:17
有关sas软件中如何将矩阵转化为sas数据集的问题
proc iml;
A={2.684 2.876 2.793 0.536 2.010,
  2.694 2.924 2.823 0.563 2.487,
  2.737 2.981 2.846 0.633 2.709,
  2.698 2.965 2.850 0.608 2.789,
  2.601 2.878 2.775 0.549 2.855,
  2.647 2.893 2.840 0.571 2.913,
  2.677 2.942 2.836 0.595 2.965,
  2.679 2.967 2.882 0.605 3.052,
  2.714 3.029 2.903 0.649 3.090,
  2.711 3.027 2.915 0.636 3.124};
B={1.763,
   2.127,
   2.248,
   2.267,
   2.270,
   2.274,
   2.318,
   2.405,
   2.441,
   2.462};
X=STANDARD(A);
Y=STANDARD(B);
Z=t(X)*X;
val=eigval(Z);
vet=eigvec(Z);
print X Y Z val vet;
F1=X*vet[,1];
F2=X*vet[,2];
F3=X*vet[,3];
Q=(val[1,]+val[2,]+val[3,])/(10*5);
print F1 F2 F3 Q;
quit;run;
这是我在sas软件上编的一个程序,现在我想做B对F1 F2 F3的线性回归。这需要先把四个矩阵转化为一个数据集的形式。请教各位高人,帮帮忙,查了很多书都没有找到啊,谢谢了
You may use

create DSN_A from a;
  append from a;

See example attached below.
*****************************************;

proc iml;
A={2.684 2.876 2.793 0.536 2.010,
  2.694 2.924 2.823 0.563 2.487,
  2.737 2.981 2.846 0.633 2.709,
  2.698 2.965 2.850 0.608 2.789,
  2.601 2.878 2.775 0.549 2.855,
  2.647 2.893 2.840 0.571 2.913,
  2.677 2.942 2.836 0.595 2.965,
  2.679 2.967 2.882 0.605 3.052,
  2.714 3.029 2.903 0.649 3.090,
  2.711 3.027 2.915 0.636 3.124};
  create DSN_A from a;
  append from a;
  quit;

  proc print data=DSN_A;
  run;
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2010-1-16 22:39:11
类似是做pls的吧
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群