SAS矩阵

·现在想把一个数据库中变量的值都作为行做一个对应矩阵(观测中有该元素，对应值为1，无则为0)，原数据如下：

ID	y1	y2	y3
1	A	V	E
2	G		X
3	A	C	R

变换后数据结果如下：

ID	A	V	E	G	X	C	R
1	1	1	1	0	0	0	0
2	0	0	0	1	1	0	0
3	1	0	0	0	0	1	1

请问如何操作?谢谢！

here is my solution.

data a;
     input ID y1 $2. y2 $2. y3 $2.;
     cards;
1 A V E
2 G   X
3 A C R
run;

%macro trans_data;
     proc sql noprint;
           select y into: value1-  from
           (select distinct y1 as y from a where y1 is not null
                union select distinct y2 as y from a where y2 is not null
                union select distinct y3 as y from a where y3 is not null)
           ;
     quit;
     data a_new;
           set a;
           array yvar y1-y3;
           %do i=1 %to &sqlobs.;
                &&value&i=0;
           %end;
           do over yvar;
                select(yvar);
                     %do i=1 %to &sqlobs.;
                           when("&&value&i")  &&value&i=1;
                     %end;
                     otherwise;
                end;
           end;
           drop y1-y3;
     run;
%mend;
%trans_data;

It works, but is not efficient for a big data.
data a;
input ID y1 $2. y2 $2. y3 $2. ;
cards;
1 A V E
2 G   X
3 A C R
run;
data b (keep=ID A V E G X C R);
set a;
array aa{*} A V E G X C R;
array bb{*} y1-y3;
do i=1 to dim(aa);
        aa=0;
        do j=1 to dim(bb);
                if bb[j]=vname(aa) then do;
                        aa=1;
                        leave ;
                end;
        end;
end;
run;

teqel 发表于 2015-8-13 03:21
It works, but is not efficient for a big data.
data a;
input ID y1 $2. y2 $2. y3 $2. ;

Your code is not only inefficient, but also inflexible. It's not good idea to hard-code the data value (A, V, E, ..., R).

chiant 发表于 2015-8-13 05:17
here is my solution.

data a;

朋友你的代码确实灵活，但是用distinct这种需要排序的statement外加连用3个select from，要是遇到大一些的数据，也是挺“呵呵”的……

sniperhgy 发表于 2015-8-13 13:25
朋友你的代码确实灵活，但是用distinct这种需要排序的statement外加连用3个select from，要是遇到大一些的 ...

ok, 那您有本事给个大数据的solution呗

teqel 发表于 2015-8-13 03:21
It works, but is not efficient for a big data.
data a;
input ID y1 $2. y2 $2. y3 $2. ;

谢谢您！

chiant 发表于 2015-8-12 09:47
here is my solution.

data a;

好厉害，看不懂！谢谢

sniperhgy 发表于 2015-8-13 13:25
朋友你的代码确实灵活，但是用distinct这种需要排序的statement外加连用3个select from，要是遇到大一些的 ...

确实是大数据~

chiant 发表于 2015-8-13 16:31
ok, 那您有本事给个大数据的solution呗

说话请不要带节奏哦^_^

而且不光是处理大数据上，还有就是3个变量就需要union三次，要是10个变量就来十次？

我发下我的解决方案，欢迎斧正：

1. data a;
   
2.      input ID y1 $1. y2 $1. y3 $1.;
   
3.      cards;
   
4. 1 AVE
   
5. 2 G X
   
6. 3 ACR
   
7. ;
   
8. run;
   
9. 
   
10. %macro CreateCountTable;
    
11.   proc contents
    
12.     noprint
    
13.     data = a
    
14.     out = var_of_a(keep = NAME where = (NAME ne "ID"));
    
15.   run;
    
16. 
    
17.   data _null_;
    
18.     set var_of_a end = eof;
    
19.         call symput(strip("VAR") || strip(_N_), strip(NAME));
    
20.         if eof then
    
21.       call symput(strip("COUNT_OF_VAR"), strip(_N_));
    
22.   run;
    
23. 
    
24.   data _NULL_;
    
25.     if 0 then set A nobs = n;
    
26.     call symput('obs_of_a',trim(left(put(n, 8.))));
    
27.     stop;
    
28.   run;
    
29. 
    
30.   %do i = 1 %to &obs_of_a.;
    
31.     data _NULL_;
    
32.       set a(firstobs = &i. obs = &i.);
    
33.           %do j = 1 %to &COUNT_OF_VAR.;
    
34.         call symput(strip("VALUE") || strip(&j.), strip(&&VAR&j.));
    
35.       %end;
    
36.     run;
    
37.     %put _user_;
    
38.     data part_&i.;
    
39.       %do j = 1 %to &COUNT_OF_VAR.;
    
40.         %if &&VALUE&j. ne %str() %then
    
41.           &&VALUE&j. = 1;;
    
42.       %end;
    
43.     run;
    
44.   %end;
    
45. 
    
46.   data wanted;
    
47.     set part_1 - part_&obs_of_a.;
    
48. 
    
49.     array myArr _numeric_;
    
50.         do over myArr;
    
51.       myArr = coalesce(myArr, 0);
    
52.     end;
    
53.   run;
    
54. 
    
55.   proc datasets;
    
56.     delete part_1 - part_&obs_of_a.;
    
57.   run;
    
58.   quit;
    
59. %mend;
    
60. 
    
61. %CreateCountTable;

复制代码