如何根据观测的奇偶数拆分数据集

data aaa1;
input x$ y$;
cards;
a b
c d
e f
g h
l j
;
run;

data aaa2;
input x$ y$;
cards;
a b
c d
e f
g h
;
run;

根据数据集中观测的个数拆分数据集
如果为观测为奇数如aaa1，则数据集aaa1拆分为
x1 y1 x2 y2
a b g h
c d l j
e f

如果为观测为偶数如aaa2，则数据集aaa2拆分为
x1 y1 x2 y2
a b e f
c d g h

1. data aa1;
   
2. input x $ y $;
   
3. cards;
   
4. a b
   
5. c d
   
6. e f
   
7. g h
   
8. i j
   
9. ;
   
10. run;
    
11. data aa2;
    
12. input x $ y $;
    
13. cards;
    
14. a b
    
15. c d
    
16. e f
    
17. g h
    
18. ;
    
19. run;
    
20. %macro m(data);
    
21. data m;
    
22. set &data nobs=nobs;
    
23. n=nobs;
    
24. n2=_n_;
    
25. run;
    
26. proc sql;
    
27. create table m2 as
    
28. select x as x1,y as y1
    
29. from m
    
30. where n2>(n+1)/2;
    
31. quit;
    
32. data m2;
    
33. set m2;
    
34. n2=_n_;
    
35. run;
    
36. data data_&data;
    
37. merge m m2;
    
38. by n2;
    
39. if n2>(n+1)/2 then delete;
    
40. keep x y x1 y1;
    
41. run;
    
42. proc print data=data_&data;
    
43. run;
    
44. %mend;
    
45. %m(aa1);
    
46. %m(aa2);

复制代码

就喜欢马 发表于 2015-2-1 12:44

data m3 改成data _&data 这样的比较好 要不覆盖前面结果了

sushe1527 发表于 2015-2-1 20:52
data m3 改成data _&data 这样的比较好 要不覆盖前面结果了

是哦，

1. data aaa1;
   
2. set aaa1;
   
3. n=_n_;
   
4. run;
   
5. proc sql;
   
6. select count(*) into:num from aaa1;
   
7. quit;
   
8. data aaa(drop=n);
   
9. merge aaa1(where=(n<=int((&num.+1)/2)) rename=(x=x1 y=y1))  aaa1(where=(n>int((&num.+1)/2)) rename=(x=x2 y=y2) );
   
10. run;

复制代码

data _null_;
        if 0 then set aaa1 nobs=nobs;
         cut1=ceil(nobs/2); cut2=cut1+1;
        call execute(cat('data want;
        merge aaa1(rename=(x=x1 y=y1) obs=',cut1,')
        aaa1(rename=(x=x2 y=y2) firstobs=',cut2,');
        run;'));
        stop;
run;