如何计算随访记录中每个人每种药的使用时间？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

有这样一个类似的数据集，id代表患者编号，n代表随访次数，treat是用药方案，v_date 是随访日期。
怎么统计每个人用每种方案的分别的使用时间呢？
比如第一个人用a方案的时间应该为第3次减去第1次的时间，即（2012-5-3） -（2012-1-1）谢谢!
three是原始数据集。
data three;
input id n treat$ v_date : yymmdd10.;
format v_date yymmdd10.;
cards;
1        1        a        2012-1-1
1        2        a        2012-3-1
1        3        c        2012-5-3
1        4        c        2012-7-28
1        5        a        2012-10-1
1        6        a        2013-1-3
2        1        a        2009-3-9
2        2        a        2009-7-3
2        3        c        2009-12-1
2        4        a        2010-2-3
2        5        c        2010-5-4
2        6        c        2010-8-26
2        7        c        2010-12-3

;
run;
想得到每个人的这些观测

1	1	a	2012-1-1
1	3	c	2012-5-3
1	5	a	2012-10-1
1	6	a	2013-1-3
2	1	a	2009-3-9
2	3	c	2009-12-1
2	4	a	2010-2-3
2	5	c	2010-5-4
2	7	c	2010-12-3

或者更好的就是产生新变量，每个人只有一条记录，如：

id

n

treat

v_date

treat_1

v_date_1

treat_2

v_date_2

treat_3

v_date_3

treat_4

v_date_4

treat_5

v_date_5

treat_6

v_date_6

1

1

a

2012-1-1

a

2012-1-1

c

2012-5-3

a

2012-10-1

a

2013-1-3

2

1

a

2009-3-9

a

2009-3-9

c

2009-12-1

a

2010-2-3

c

2010-5-4

c

2010-12-3

非常感谢！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：使用时间 format treat Three Input treat 记录 如何 统计

有人能帮忙看看吗？是说的不够清楚吗？就是想算随访中的暴露的时间。
谢谢各位啊！

1. data three;
   
2. input id n treat $ v_date : yymmdd10.;
   
3. format v_date yymmdd10.;
   
4. cards;
   
5. 1        1        a        2012-1-1
   
6. 1        2        a        2012-3-1
   
7. 1        3        c        2012-5-3
   
8. 1        4        c        2012-7-28
   
9. 1        5        a        2012-10-1
   
10. 1        6        a        2013-1-3
    
11. 2        1        a        2009-3-9
    
12. 2        2        a        2009-7-3
    
13. 2        3        c        2009-12-1
    
14. 2        4        a        2010-2-3
    
15. 2        5        c        2010-5-4
    
16. 2        6        c        2010-8-26
    
17. 2        7        c        2010-12-3
    
18. ;run;
    
19. 
    
20. proc sort data=three out=three_sort;
    
21.     by id n treat;
    
22. run;
    
23. 
    
24. data three_a;
    
25.     set three_sort;
    
26.        by id;
    
27.        retain tg;
    
28.        treat_lag = lag(treat);
    
29.        if first.id then tg=0;  
    
30.        if treat ^=treat_lag then tg+1;
    
31.    drop treat_lag;
    
32. run;
    
33. 
    
34. data three_b;
    
35.     set three_a;
    
36.        by id tg;
    
37.        if first.tg or last.id;
    
38.    drop tg;
    
39. run;
    
40. proc print data=three_b; title 'three_b'; run;

复制代码

yongyitian 发表于 2013-8-20 11:15

非常感谢楼上。