SAS如何实现

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

数据大概是这样的：


rsid  SNPid  AlleleA  AlleleB  
1      1         A          T          0
1      1         A          T          1
1      1         A          T          0
1      1         A          T          1
1      1         A          T          0
1      1         A          T          0
...
...
...
2      2         C          G        1
2      2         C          G        0
2      2         C          G        0
2      2         C          G        0
2      2         C          G        1
2      2         C          G        0


想新生成变量 m1 and m2。rsid=1时 ，m1重复 1 to 3， m2 每隔观测加1。rsid=2时，m1还是重复1 to 3，m2 从10开始，然后也是每隔三个观测加1；

rsid  SNPid  AlleleA  AlleleB    m1  m2
1      1         A          T          0         1     1
1      1         A          T          1         2     1
1      1         A          T          0         3     1
1      1         A          T          1         1     2
1      1         A          T          0         2     2
1      1         A          T          0         3     2
...
...
...
2      2         C          G        1          1     10
2      2         C          G        0           2     10
2      2         C          G        0           3     10
2      2         C          G        0           1     11
2      2         C          G        1           2     11
2      2         C          G        0           3     11








请大家帮忙想想如何编写code，谢谢

扫码加我 拉你入群

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

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

关键词：如何实现 code 生成变量 RSI sid 如何

1. 每个rsid的数量是否都是3的倍数？
2. m2 从10 开始的原因，是这样定义， 还是你给的数据中省略号的部分会是3,4,5...， 到了rsid=2的时候正好是10？

wwang111 发表于 2016-7-28 11:35
1. 每个rsid的数量是否都是3的倍数？
2. m2 从10 开始的原因，是这样定义， 还是你给的数据中省略号的部分 ...

您好，谢谢您的回复！
问题1答复：每个rsid的数量是3的倍数
问题2答复：m2 从10 开始的原因，是这样定义，无论当rsid=1时，m2=10、11或者12……

我还是不太明白你说的，“无论当rsid=1时，m2=10、11或者12……”，这句话应该怎么理解？

wwang111 发表于 2016-7-28 12:35
我还是不太明白你说的，“无论当rsid=1时，m2=10、11或者12……”，这句话应该怎么理解？

rsid=1 与rsid=2时的m2是相互独立的

Try this:

data wanted;
m2+1;
  do m1=1 to 3;
   set test;
   if rsid ne 1 and rsid ne lag(rsid) then m2=10;
   output;
end;
run;

wwang111 发表于 2016-7-28 12:42
Try this:

data wanted;

非常感谢！