江湖救急：如何替换SAS数据集中的观测值

苹果叶 发表于 2014-5-2 18:11
如果可以多设几个值，就这样吧。如果纯合子只能1,2 就比较麻烦了。

恩，谢谢哈。纯合子只能是0、2，还有什么好的办法吗

进行式 发表于 2014-5-2 18:37
恩，谢谢哈。纯合子只能是0、2，还有什么好的办法吗

就是只能用0、1、2、3个数字替换

进行式 发表于 2014-5-2 18:40
就是只能用0、1、2、3个数字替换

1. proc import datafile = 'C:\Users\XXX\Desktop\gene.txt'
   
2.   dbms = dlm
   
3.   out  = gene
   
4.   replace;
   
5.   delimiter = '09'x;
   
6.   getnames  = no;
   
7.   datarow   = 2;
   
8. run;
   
9. 
   
10. proc transpose data = gene out = trans_gene;
    
11.   var var2-var2181;
    
12. run;
    
13. 
    
14. proc format;
    
15.   value $C2N
    
16.     'AA'  = '1'
    
17.     'CC'  = '2'
    
18.     'GG'  = '3'
    
19.     '--'  = '--'
    
20.     other = '0'
    
21.   ;
    
22. run;
    
23. 
    
24. data trans_gene1;
    
25.   set trans_gene;
    
26.   array col{20};
    
27.   do i = 1 to 20;
    
28.     col{i} = put(col{i},C2N.);
    
29.   end;
    
30.   drop i;
    
31. run;
    
32. 
    
33. %macro tranwrd;
    
34.   %do i = 1 %to 20;
    
35. 
    
36.     proc sql noprint;
    
37.       select distinct col&i. into :values separated by '' from trans_gene1
    
38.           order by col&i.;
    
39.     quit;
    
40. 
    
41.     %if %sysfunc(compress("&values",'-')) ne '012' %then %do;
    
42.           data trans_gene1;
    
43.             set trans_gene1;
    
44.              if compress("&values",'-') = '013' then col&i. = tranwrd(col&i.,'3','2');
    
45.         else if compress("&values",'-') = '023' then col&i. = tranwrd(col&i.,'3','1');
    
46.           run;
    
47.     %end;
    
48. 
    
49.   %end;
    
50. %mend;
    
51. 
    
52. %tranwrd;
    
53. 
    
54. proc transpose data = trans_gene1 out = final;
    
55.   id _name_;
    
56.   var col1-col20;
    
57. run;
    
58. 
    
59. data final;
    
60.   merge gene final;
    
61.   drop _name_;
    
62.   rename var1 = snpname;
    
63. run;

复制代码

我的方法比较笨，先转置，把3改成0 或者 2， 再转置回来。  

苹果叶 发表于 2014-5-2 20:00
我的方法比较笨，先转置，把3改成0 或者 2， 再转置回来。

恩，十分感谢您的帮助。可是突然发现这种思路始终还是分不开纯合子，因为只能是0、1、2三个数字代替。就比如AG型里纯合为AA=0，GG=2，AC型里AA=0，CC=2，而CG型里CC=0，GG=2，这种方法没法区别每一行的CC和GG。所以我在想是不是可以通过字符串来修改？无论在哪一行里，第一个识别出的两个相同字符为，不同为0，第二个识别出的相同字符为2.但是具体的程序不知道怎么写？再次感谢您的帮助！

进行式 发表于 2014-5-2 21:14
恩，十分感谢您的帮助。可是突然发现这种思路始终还是分不开纯合子，因为只能是0、1、2三个数字代替。就比 ...

按你的想法，其实我的程序也可以分出来了，AA=1, CC=2 这两种是永远不变的，变的其实只是GG

比如某一行有 AA，GG，CG， 那么AA=1，CG=0，则GG=2
另一行有 CC, GG, AG, 那么 CC =2 ，AG=0，则GG=1

进行式 发表于 2014-5-2 21:14
恩，十分感谢您的帮助。可是突然发现这种思路始终还是分不开纯合子，因为只能是0、1、2三个数字代替。就比 ...

按你的想法，其实我的程序也可以分出来了，AA=1, CC=2 这两种是永远不变的，变的其实只是GG

比如某一行有 AA，GG，CG， 那么AA=1，CG=0，则GG=2
另一行有 CC, GG, AG, 那么 CC =2 ，AG=0，则GG=1

苹果叶 发表于 2014-5-2 22:38
按你的想法，其实我的程序也可以分出来了，AA=1, CC=2 这两种是永远不变的，变的其实只是GG

比如某一行 ...

恩，多谢了！我再试试

进行式 发表于 2014-5-3 09:38
恩，多谢了！我再试试

还有个小问题想向您请教，我想保留一个数据集中含某一特定观测值的变量，改怎么写程序呢？谢谢哈