悬赏求代码 - 经管之家

悬赏求代码

hamsik11

4920

收藏 2013-08-07

悬赏 200 个论坛币已解决

如下面数据所示，原始数据格式内容如下所示，每一列为一人名，每一行数据表示两个人之间建立联系，

Tom	Jack
Jack	Mike
Tom	Hellen
Dora	Hellen
Ross	Jack
Will	Paul
Paul	Howard
Paul	Andy
Andy	Kim
Andy	Tom

要求把数据转成一个矩阵，上面变量名为人名，第一列也为人名，如果两个人之间建立直接联系，则标识为1，没有则为0，最终数字矩阵为一对称矩阵。因不懂IML，往高手帮忙，谢谢。

names	Andy	Dora	Hellen	Howard	Jack	Kim	Mike	Paul	Ross	Tom	Will
Andy	1	0	0	0	0	1	0	1	0	1	0
Dora	0	1	0	0	0	0	0	0	0	0	0
Hellen	0	0	1	0	0	0	0	0	0	1	0
Howard	0	0	0	1	0	0	0	1	0	0	0
Jack	0	0	0	0	1	0	1	0	1	1	0
Kim	1	0	0	0	0	1	0	0	0	0	0
Mike	0	0	0	0	1	0	1	0	0	0	0
Paul	1	0	0	1	0	0	0	1	0	0	0
Ross	0	0	0	0	1	0	0	0	1	0	0
Tom	1	0	1	0	1	0	0	0	0	1	0
Will	0	0	0	0	0	0	0	0	0	0	1

最佳答案

ziyenano 查看完整内容

扫码加我拉你入群

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

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

全部回复

ziyenano

2013-8-7 13:06:11

data have;
input name1 $ name2 $;
datalines;
Tom       Jack
Jack       Mike
Tom       Hellen
Dora       Hellen
Ross       Jack
Will       Paul
Paul       Howard
Paul       Andy
Andy       Kim
Andy       Tom
;
run;

proc sql;
create table name_id as
select name,monotonic() as id from
(select name1 as name from have
union
select name2 as name from have);

create table name_id1 as
select t.*,t1.id as id1
from
have t
inner join
name_id t1
on t.name1=t1.name;

create table name_id2 as
select t.*,t1.id as id2
from
name_id1 t
inner join
name_id t1
on t.name2=t1.name;
quit;

proc iml;
use name_id;
read all var {name} into name;
close name_id;
use name_id2;
read all var{id1,id2} into id;
close name_id2;
res=j(nrow(name),nrow(name),0);
res[id[,1],id[,2]]=1;
create want from res[colnames=name rownames=name];
quit;

扫码加我拉你入群

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

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

playmore

2013-8-7 14:46:06

data have;
input name1 $ name2 $;
datalines;
Tom Jack
Jack Mike
Tom Hellen
Dora Hellen
Ross Jack
Will Paul
Paul Howard
Paul Andy
Andy Kim
Andy Tom
;
run;

data have;
set have;
id=1;
run;

proc sort data=have;
by name1;
run;

proc transpose data=have out=want;
by name1;
id name2;
var id;
run;

data want(drop=i);
set want;
array numvar _NUMERIC_;
do i=1 to dim(numvar);
if numvar(i)=. then numvar(i)=0;
end;
run;

扫码加我拉你入群

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

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

akakak123

2013-8-7 15:50:01

data a;
input name $ @@;
cards;
Andy Dora Hellen Howard Jack Kim Mike Paul Ross Tom Will
;
run;
data b;
set a;
length Andy $1 Dora $1 Hellen $1 Howard $1 Jack $1 Kim $1 Mike $1 Paul $1 Ross $1 Tom $1 Will $1 ;
run;


data c;
set b;
if name='Tom' then Jack=1;
else Jack=0;
if name='Jack' then Mike=1;
else Mike=0;
if name='Tom' then Hellen=1;
else Hellen=0;
if name='Dora' then Hellen=1;
if name='Ross' then Jack=1;
if name='Will' then  Paul=1;
else Paul=0;
if name='Paul' then Howard=1;
else Howard=0;
if name='Paul' then Andy=1;
else Andy=0;
if name='Andy' then Kim=1;
else Kim=0;
if name='Andy' then Tom=1;
else Tom=0;
if name='Tom' then Jack=1;
if name='Mike' then Jack=1;
if name='Hellen' then Tom=1;
if name='Hellen' then Dora=1;
else Dora=0;
if name='Jack' then Ross=1;
else Ross=0;
if name='Paul' then  Will=1;
else Will=0;
if name='Howard' then Paul=1;
if name='Andy' then Paul=1;
if name='Kim' then Andy=1;
if name='Tom' then Andy=1;
if name='Andy' then andy=1;
if name='Dora' then Dora=1;
if name='Hellen' then hellen=1;
if name='Howard' then howard=1;
if name='Jack' then jack=1;
if name='Kim' then kim=1;
if name='Mike' then Mike=1;
if name='Paul' then paul=1;
if name='Ross' then ross=1;
if name='Tom' then tom=1;
if name='Will' then will=1;
run;

扫码加我拉你入群

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

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

hamsik11

2013-8-7 16:48:04

playmore 发表于 2013-8-7 14:46
data have;
input name1 $ name2 $;
datalines;

这个不行你看Tom和andy在有关联，那tom-andy，andy-tom对应的都应该是1

扫码加我拉你入群

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

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

hamsik11

2013-8-7 16:48:33

playmore 发表于 2013-8-7 14:46
data have;
input name1 $ name2 $;
datalines;

这个不行你看Tom和andy在有关联，那tom-andy，andy-tom对应的都应该是1

扫码加我拉你入群

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

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

点击查看更多内容…

hamsik11

2013-8-7 16:49:23

akakak123 发表于 2013-8-7 15:50
data a; ...

这样不行，没扩展性，倒腾半天，不如我在excel里面填数字再导入sas算了

扫码加我拉你入群

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

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

playmore

2013-8-7 17:08:15

hamsik11 发表于 2013-8-7 16:48
这个不行你看Tom和andy在有关联，那tom-andy，andy-tom对应的都应该是1

嗯，刚看了下
我的方法适合单向的关联
你要是做多项的话
把我得到的结果，用IML对矩阵沿对角线做个镜像吧，一边为1则别一边也为1就好
否则改下初始表格也行，如下所示：

data have(drop=name0);
set have;
output;
name0=name1;
name1=name2;
name2=name0;
output;
run;

扫码加我拉你入群

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

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

hamsik11

2013-8-7 17:29:23

playmore 发表于 2013-8-7 17:08
嗯，刚看了下
我的方法适合单向的关联
你要是做多项的话

捣腾出来了,以前没用过Proc IML，....

data have ;
input name1 :$8. name2 :$8. ;
dummy = 1 ;
datalines;
Tom Jack
Jack Mike
Tom Hellen
Dora Hellen
Ross Jack
Will Paul
Paul Howard
Paul Andy
Andy Kim
Andy Tom
;
;;;;

/* create list of unique names that are in either column */
/* assumes names in name1 and name2 may be different lists */
proc sql ;
create table unique_names as
select distinct name
from
( select distinct name1 as name
from have )
union (select distinct name2 as name
from have )
order by name
;
quit;

/* merge with original data to create complete list */
proc sort data=have ;
by name1 name2;
run;

data have2 ;
merge have (in=a) unique_names (in=b rename=(name=name1));
by name1; /* left side */
if a or b;
run;

proc sort data=have2 ;
by name2 name1;

data havenew ;
merge have2 (in=a) unique_names (in=b rename=(name=name2));
by name2; /*right side */
if a or b;
run;

/* resort */
proc sort data=havenew;
by name1 name2 ;
run;

/* read column names into macro var in sorted order */
/* used to sort column output from proc transpose */
proc sql noprint;
select name into :namelist separated by ' '
from unique_names
order by name
;
quit;

/* get number of names - used for array processing later on */
proc sql noprint;
select count ( distinct name ) into :namecnt
from unique_names
;
quit;

/* transpose rows to columns */
proc transpose data=havenew out=want ;
id name2;
by name1 ;
run;

/* reorder proc transpose output */
/* assign 1 when name1 var matches column name */

data want;
attrib name1 &namelist label='';
set want (drop=_name_);
where name1 ne ' ' ;
run;

data havesorted (drop=i);
set want;
/* array nm {&namecnt} &namelist ;*/
array nm

&namelist;
do i = 1 to &namecnt ;
/* deal with missing values (set to 1 or 0 ) */
if nm{i} = . then do;
/* set to 1 if the person's name is the same as the column name */
if vname(nm)=name1 then nm{i}=1;
else nm{i} = 0;
end;
end;
run;

proc iml;
use havesorted;
read all var _NUM_ into m;
n_row=nrow(m);
n_col=ncol(m);
do i=1 to n_col;
do j=1 to n_row;
if m[j,i]=1 then m[i,j]=1;
end;
end;
print m;
create MyData from m;
append from m;
close MyData;
quit;

扫码加我拉你入群

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

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

10楼

hamsik11

2013-8-7 19:14:45

ziyenano 发表于 2013-8-7 18:45
data have;
input name1 $ name2 $;
datalines;
这个很大神，最后proc iml有点小问题出不来结果能再看看改一下吗？谢谢拉

扫码加我拉你入群

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

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

11楼

hamsik11

2013-8-7 19:15:06

ziyenano 发表于 2013-8-7 18:45
data have;
input name1 $ name2 $;
datalines;
这个很大神，最后proc iml有点小问题出不来结果能再看看改一下吗？谢谢拉

扫码加我拉你入群

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

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

12楼

ziyenano

2013-8-7 20:43:11

hamsik11 发表于 2013-8-7 19:14
这个很大神，最后proc iml有点小问题出不来结果能再看看改一下吗？谢谢拉
最后加一句：
create want from res[colname=name rowname=name];
append from res[rowname=name];

扫码加我拉你入群

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

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

13楼

hamsik11

2013-8-7 21:58:06

ziyenano 发表于 2013-8-7 20:43
最后加一句：
create want from res[colname=name rowname=name];
append from res[rowname=name];
高手谢谢，不过还有问题啊，譬如“Andy”列，按照数据，Andy只跟kim，tom，paull三个直接联结，加上对角线，应该只有4个1，可是生成的“Andy”列却有7个1.......

扫码加我拉你入群

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

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

14楼

ziyenano

2013-8-7 22:10:42

hamsik11 发表于 2013-8-7 21:58
高手谢谢，不过还有问题啊，譬如“Andy”列，按照数据，Andy只跟kim，tom，paull三个直接联结，加上对角 ...
下标表示的有点问题，这样应该就没问题了：
proc iml;
use name_id;
read all var {name} into name;
close name_id;
use name_id2;
read all var{id1,id2} into id;
close name_id2;
res=j(nrow(name),nrow(name),0);
do i=1 to nrow(name);
res[id[i,1],id[i,2]]=1;
end;
create want from res[rowname=name colname=name];
append from res[rowname=name];
quit;

扫码加我拉你入群

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

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

15楼

jingju11

2013-8-7 23:34:18

如果没有IML, 数组也有类似的功能.值得注意的是数组和矩阵只是形似..京剧
复制代码

扫码加我拉你入群

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

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

16楼

yongyitian

2013-8-8 08:05:05

复制代码

扫码加我拉你入群

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

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

17楼

hamsik11

2013-8-8 10:42:32

ziyenano 发表于 2013-8-7 22:10
下标表示的有点问题，这样应该就没问题了：
proc iml;
use name_id;
比较悲剧的是我处理一些大一点的数据，譬如10万个用户是，在res=j(nrow(name),nrow(name),0);这一步出现错误提示
ERROR: (execution) Unable to allocate sufficient memory. At least 2147483647 more bytes required.

请问这个有办法解决吗？

扫码加我拉你入群

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

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

18楼

hamsik11

2013-8-8 10:48:31

jingju11 发表于 2013-8-7 23:34
如果没有IML, 数组也有类似的功能.值得注意的是数组和矩阵只是形似..京剧
多谢多谢只是这个也存在一个问题，在处理大量用户数据的时候，在catx函数那一步出现，字符截断的问题，在 array mat[&nobs, &nobs] _temporary_(%eval(&nobs**2)*0) ;也出现内存不足的问题，我已经把设置MAXMEMQUERY=MAX了.....

请问jingju老师这个有办法解决吗？

扫码加我拉你入群

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

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

19楼

hamsik11

2013-8-8 10:49:15

jingju11 发表于 2013-8-7 23:34
如果没有IML, 数组也有类似的功能.值得注意的是数组和矩阵只是形似..京剧
多谢多谢只是这个也存在一个问题，在处理大量用户数据的时候，在catx函数那一步出现，字符截断的问题，在 array mat[&nobs, &nobs] _temporary_(%eval(&nobs**2)*0) ;也出现内存不足的问题，我已经把设置MAXMEMQUERY=MAX了.....

请问jingju老师这个有办法解决吗？

扫码加我拉你入群

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

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

20楼

南方de梦

2013-8-8 10:51:50

data have ;
      input name1 :$8. name2 :$8. ;
      dummy = 1 ;
datalines;
Tom Jack
Jack Mike
Tom Hellen
Dora Hellen
Ross Jack
Will Paul
Paul Howard
Paul Andy
Andy Kim
Andy Tom
;
;;;;

/* create list of unique names that are in either column */
/* assumes names in name1 and name2 may be different lists */
proc sql ;
      create table unique_names as
      select distinct name
      from
      ( select distinct name1 as name
      from have )
      union (select distinct name2 as name
      from have )
      order by name
;
quit;

/* merge with original data to create complete list */
proc sort data=have ;
      by name1 name2;
      run;

data have2 ;
      merge have (in=a) unique_names (in=b rename=(name=name1));
      by name1; /* left side */
      if a or b;
run;

proc sort data=have2 ;
      by name2 name1;

data havenew ;
      merge have2 (in=a) unique_names (in=b rename=(name=name2));
      by name2; /*right side */
      if a or b;
run;

/* resort */
proc sort data=havenew;
      by name1 name2 ;
run;

/* read column names into macro var in sorted order */
/* used to sort column output from proc transpose */
proc sql noprint;
      select name into :namelist separated by ' '
      from unique_names
      order by name
;
quit;

/* get number of names - used for array processing later on */
proc sql noprint;
      select count ( distinct name ) into :namecnt
      from unique_names
;
quit;

/* transpose rows to columns */
proc transpose data=havenew out=want ;
      id name2;
      by name1 ;
run;

/* reorder proc transpose output */
/* assign 1 when name1 var matches column name */

data want;
attrib name1 &namelist label='';
set want (drop=_name_);
where name1 ne ' ' ;
run;

data havesorted (drop=i);
      set want;
      /* array nm {&namecnt} &namelist ;*/
      array nm
&namelist;
      do i = 1 to &namecnt ;
            /* deal with missing values (set to 1 or 0 ) */
            if nm{i} = . then do;
                     /* set to 1 if the person's name is the same as the column name */
                     if vname(nm)=name1 then nm{i}=1;
                     else nm{i} = 0;
            end;
      end;
run;

proc iml;
use havesorted;
read all var _NUM_ into m;
n_row=nrow(m);
n_col=ncol(m);
do i=1 to n_col;
do j=1 to n_row;
   if m[j,i]=1 then m[i,j]=1;
end;
end;
print m;
create MyData from m;
append from m;
close MyData;
quit;

扫码加我拉你入群

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

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

21楼

jingju11

2013-8-8 11:28:08

what about the solution on 16th floor? Transpose procedure might work (but not sure). In practice, 100K*100K table is useless in any way. unable to calcuate, unable to print, unable to query...
But any way, I think the matrix is additive by row of same names(?), plus to mirror the elements by diagonal from  right to left corner...
JingJu

扫码加我拉你入群

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

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

22楼

hamsik11

2013-8-8 11:35:20

jingju11 发表于 2013-8-8 11:28
what about the solution on 16th floor? Transpose procedure might work (but not sure). In practice, 1 ...
yongyitian高手用了select distinct name1 into : NameList，按照我的数据量的话这个肯定也超出sas允许的字符串长度了......

扫码加我拉你入群

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

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

23楼

jingju11

2013-8-8 11:38:46

很显然名称的问题不是关键.一个字符串放不下,多用几个好了.但是内存的问题是难点.京剧

扫码加我拉你入群

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

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

24楼

ziyenano

2013-8-8 11:46:21

hamsik11 发表于 2013-8-8 10:42
比较悲剧的是我处理一些大一点的数据，譬如10万个用户是，在res=j(nrow(name),nrow(name),0);这一步出现 ...
好吧，10W*10W，每个单元算一个字节，10^10/1024/1024/1024=9.31G；
这种稀疏矩阵耗内存；
楼主是需要通过矩阵形式实现什么算法么？
我不大清楚SAS运行的机制，如果运算都是放在内存中进行，估计是爱莫能助了，
即使换了服务器，你的数据量再上一个等级，服务器也撑不住。
新建一个库，libname test "E:\";保证库的空间足够大。
proc iml;
reset deflib=test;
默认数据放在test库中；
或者用楼上几位的data步，数据也存储在test库中。
试试看吧。

扫码加我拉你入群

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

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

25楼

jingju11

2013-8-8 21:05:30

说的不错.通过实验发现,当维度在10K*10K时, 数据尺寸1.2G. 也就是说, 如此推算, 100K的方阵大概需要100G以上(当然这不是内存).而且时间的消耗也相当令人心烦(通过逐步叠加的方法,10K*10K的表格需要15分钟).
京剧

扫码加我拉你入群

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

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

26楼

playmore

2013-8-9 08:22:09

hamsik11 发表于 2013-8-8 11:35
yongyitian高手用了select distinct name1 into : NameList，按照我的数据量的话这个肯定也超出sas允许的 ...
看了下你的问题
如果你的数据量大，矩阵全存起来肯定不是好办法
SAS有解稀疏矩阵的办法，尤其适用于对称矩阵
是用索引的办法存储矩阵，在资源占用和速度肯定有优化

所以归根结底，
如果你不需要做复杂的矩阵运算
你之前的表就够用了
不用再想着变成矩阵形了

扫码加我拉你入群

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

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

27楼

hamsik11

2013-8-9 09:58:47

谢谢各位了学习了不少感谢中....

另外，我怎么给各位发论坛币，哥们儿准备给大家发点小小的酬劳............

扫码加我拉你入群

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

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

栏目导航

SAS专版

文献求助专区

会计与财务管理

金融学（理论版）

学道会

经管文库（原现金交易版）

热门文章

【同程商旅】中国企业出海差旅研究报告

“十四五”能源发展成就报告

智算无界AIDC的超越和重构2025

当社科基础理论重大理论发现的时候

【10+指标】2007-2024年上市公司污染物排放 ...

【24重磅,自用整理!】2000-2024上市公司投资 ...

2025年我国医药航空冷链发展现状与趋势展望 ...

是相信人工智能？还是否定人工智能？相信就 ...

ibm-AI时代的银行业-以AI驭险，更须为AI设防 ...

CNNIC-生成式人工智能应用发展报告（2025）

推荐文章

AI狂潮席卷学术圈，不会编程也能打造专属智 ...

10月重磅来袭｜《打造Coze/Dify专属学术智能 ...

最快1年拿证，学费不足5W！热门美国人工智能 ...

关于如何利用文献的若干建议

关于学术研究和论文发表的一些建议

关于科研中如何学习基础知识的一些建议 (一 ...

一个自编的经济学建模小案例 --写给授课本科 ...

AI智能体赋能教学改革: 全国AI教育教学应用 ...

2025中国AIoT产业全景图谱报告-406页

关于文献求助的一些建议

最佳答案

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群