如何只保留想要的观测

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

求教大神，假如有这样一组数据：
ID    SIG
1      正常
1      正常
1      异常有临床意义
1      异常有临床意义
1      异常无临床意义
1      正常
1      正常
2      正常
2      正常
2     异常有临床意义
2     正常
2     正常
2     正常

2     异常有临床意义
2     正常

对于每个ID，只想保留SIG第一次出现异常有临床意义到到不再出现异常有临床意义的最后一次正常，比如对于ID为1的数据来说，要保留3~6行，对于ID为2的数据来说要保留3~8行，如何编程

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

没想到什么好办法，写着写着就变的复杂了。

data have;
input ID    SIG $20.;
cards;
1      正常
1      正常
1      异常有临床意义
1      异常有临床意义
1      异常无临床意义
1      正常
1      正常
2      正常
2      正常
2     异常有临床意义
2     正常
2     正常
2     正常
2     异常有临床意义
2     正常
;run;
data want_1;
        set have;
        by id;
        if first.id then do;
                rows=0;
                flg=0;
        end;
        rows+1;
        if sig='异常有临床意义' then flg+1;
        if flg>0;

run;

proc sql ;
        create table work.want_2 as select t1.id,min(t1.rows) as normal_min_rows
                from work.want_1 t1
                        ,(select id,max(rows) as max_rows
                                from work.want_1
                                where sig='异常有临床意义'
                                group by id) t2
        where t1.id=t2.id
        and t1.rows>t2.max_rows
        and t1.sig='正常'
        group by t1.id;
quit;

data want;
        merge work.want_1 work.want_2;
        by id;
        if rows<=normal_min_rows;
        keep id sig;
run;

要方便的进行倒序查找，建议使用多重set或hash技巧：
data have;
input ID    SIG $20.;
cards;
1      正常
1      正常
1      异常有临床意义
1      异常有临床意义
1      异常无临床意义
1      正常
1      正常
2      正常
2      正常
2     异常有临床意义
2     正常
2     正常
2     正常
2     异常有临床意义
2     正常
3      正常
3     异常有临床意义
3     异常有临床意义
3      正常
3     异常无临床意义
;run;

data want(drop=tmp:);
  set have;
  by id notsorted;

  retain csfn;    *此前存在cs记录的标识;
  if first.id then csfn = .;
  if sig = '异常有临床意义' then csfn = 1;

  do i = rec to 1 by -1;    *此前cs记录已经转变为正常的标识;
    set have(rename=(id=tmpid sig=tmpsig)) nobs=rec point=i;
    if id=tmpid and tmpsig = '异常有临床意义' then leave;
    else if id=tmpid and tmpsig = '正常' and _N_>i then csfn = .;
  end;
run;

1. data test1;
   
2.     set test;
   
3.         by id;
   
4.         n+1;
   
5.         if first.id then n=1;
   
6. run;
   
7. 
   
8. proc sql;
   
9.     * 找出第一次和最后一次异常有临床意义的位置;  
   
10.     create table  test2 as select id,  max(n) as last, min(n) as first from test1 where sig="异常有临床意义" group by id;
    
11.         create table test3 as select * from test1 a full join test2 b on a.id=b.id where a.sig="正常" and a.n>b.last;
    
12.         * 找出最后一次异常有临床意义后的第一个正常位置;
    
13.         create table test4 as select distinct id, first,  min(n) as normal_first from test3 group by id;
    
14. quit;
    
15. 
    
16. data total;
    
17.     merge test1(in=a) test4;
    
18.         by id;
    
19.         if a;
    
20.     if first<=n<=normal_first;
    
21. run;

复制代码