原创：如何用SAS的SVD

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

最近做项目刚刚开发的测试代码。
%let n=5;
%let cutoff=0.7;
proc iml;
reset deflib=work;
a={22  10   2   3   7,
      14   7  10   0   8,
      -1  13  -1 -11   3,
      -3  -2  13  -2   4,
       9   8   1  -2   4,
       9   1  -7   5  -1,
       2  -6   6   5   1,
       4   5   0  -2   2};
create tb from a;
append from a ;
quit;
PROC STDIZE data=tb out=tb_std;run;
proc princomp data=tb_std out=ww;run;
proc iml;
reset deflib=work;
use tb_std;
read all into tb_std;
print tb_std;
call svd(u,q,v,tb_std);
reset fuzz;
print  u q v ;
xx=tb_std*v;
print xx;/*the result xx matrix is consistent with data set ww**/
xxx=u*diag(q);/*the result xx matrix is consistent with data set ww**/
print xxx;
call svd(tu,tq,tv,t(tb_std));
reset fuzz;
print  tu tq tv ; /*tu with rank is equal to v with rank; so are tq and tv(u)*/
sq=sum(q);
print sq;
retain_eig={0};
idx={0};
do i=1 to &n.;
eig=q;
retain_eig=retain_eig+eig;
rate_eig=retain_eig/sq;
if rate_eig >= &cutoff. then do;idx=i; new_eig=q[1:idx];stop;end;
end;
print retain_eig;print idx;print new_eig;
v1=v[,1:idx];print v1;
u1=u[,1:idx];print u1;
q1=diag(q[1:idx]);print q1;
create u_matrix  from u1 ;
append from u1;
create v_matrix  from v1 ;
append from v1;
/***valid eigvalue and single value***/
call eigen(eig_value,eig_vect,t(tb_std)*tb_std);
print eig_value;
single=sqrt(abs(eig_value));
print single;
quit;

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：sas的 SVD 如何用 consistent princomp

这么好的博文没人看？

nkwilling,这应该sas自带的svd，与你何敢，这不算你开发吧。

保存一下。
最近做项目刚刚开发的测试代码。
%let n=5;
%let cutoff=0.7;
proc iml;
reset deflib=work;
a={22  10   2   3   7,
      14   7  10   0   8,
      -1  13  -1 -11   3,
      -3  -2  13  -2   4,
       9   8   1  -2   4,
       9   1  -7   5  -1,
       2  -6   6   5   1,
       4   5   0  -2   2};
create tb from a;
append from a ;
quit;
PROC STDIZE data=tb out=tb_std;run;
proc princomp data=tb_std out=ww;run;
proc iml;
reset deflib=work;
use tb_std;
read all into tb_std;
print tb_std;
call svd(u,q,v,tb_std);
reset fuzz;
print  u q v ;
xx=tb_std*v;
print xx;/*the result xx matrix is consistent with data set ww**/
xxx=u*diag(q);/*the result xx matrix is consistent with data set ww**/
print xxx;
call svd(tu,tq,tv,t(tb_std));
reset fuzz;
print  tu tq tv ; /*tu with rank is equal to v with rank; so are tq and tv(u)*/
sq=sum(q);
print sq;
retain_eig={0};
idx={0};
do i=1 to &n.;
eig=q;
retain_eig=retain_eig+eig;
rate_eig=retain_eig/sq;
if rate_eig >= &cutoff. then do;idx=i; new_eig=q[1:idx];stop;end;
end;
print retain_eig;print idx;print new_eig;
v1=v[,1:idx];print v1;
u1=u[,1:idx];print u1;
q1=diag(q[1:idx]);print q1;
create u_matrix  from u1 ;
append from u1;
create v_matrix  from v1 ;
append from v1;
/***valid eigvalue and single value***/
call eigen(eig_value,eig_vect,t(tb_std)*tb_std);
print eig_value;
single=sqrt(abs(eig_value));
print single;
quit;

爱萌 发表于 2013-12-26 08:59
nkwilling,这应该sas自带的svd，与你何敢，这不算你开发吧。

版主为什么编辑别人的帖子呢？

程序运行后出现warning和error！

dxystata 发表于 2013-12-27 17:41
版主为什么编辑别人的帖子呢？

你自己没有看到吗，有error,而且SVD是sas自身携带的，和当时标题不符合，这是不道德的一种表现，
我只是做了我该做的