连续出现次数

做法是 ， 判断r是否大于等于15   大于等于15取1   生成 r15   然后con是同一id连续出现大于等于15的次数。

id	year	qtr	r	r15	con
1	1997	1	7
1	1997	2	8
1	1997	3	9
1	1997	4	10
1	1998	1	11
1	1998	2	12
1	1998	3	13
1	1998	4	14
1	1999	1	15	1	4
1	1999	2	16	1	4
1	1999	3	17	1	4
1	1999	4	16	1	4
2	1996	3	15	1	1
2	1996	4	14
2	1997	1	13
2	1997	2	12
2	1997	3	14
2	1997	4	16	1	4
2	1998	1	18	1	4
2	1998	2	20	1	4
2	1998	3	22	1	4
2	1998	4	7
2	1999	1	17	1	2
2	1999	2	16	1	2

一个data步：
proc sort data=a;
by id year qtr;
run;

data b;
  do tmp=1 to nobs;
           set a nobs=nobs point=tmp;
           con=.;
           if r<15 then output;
           else do;
              con=1;
              _start=tmp;
                  _end=tmp;
                  _endflag=0;
              do while(_endflag=0 and _end<nobs);
                     _end=_end+1;
                         _id=id;
                         set a point=_end;
                         if r>=15 and _id=id and _end<=nobs then con=con+1;
                         else _endflag=1;
                   end;
                   do i=_start to ifn(_end=nobs and r>=15,_end,_end-1);
                            set a point=i;
                            output;
                        end;
                        if _end=nobs and r>=15 then stop;
                        tmp=_end-1;
                end;
   end;
   stop;
   drop _:;
run;
      

data a;
input id year qtr r;
cards;
1 1997 1 7
1 1997 2 8
1 1997 3 9
1 1997 4 10
1 1998 1 11
1 1998 2 12
1 1998 3 13
1 1998 4 14
1 1999 1 15
1 1999 2 16
1 1999 3 17
1 1999 4 16
2 1996 3 15
2 1996 4 14
2 1997 1 13
2 1997 2 12
2 1997 3 14
2 1997 4 16
2 1998 1 18
2 1998 2 20
2 1998 3 22
2 1998 4 7
2 1999 1 17
2 1999 2 16
;
proc sort data=a;
by id year qtr;
run;
data b;
set a;
if r ge 15 then r15=1;
run;
proc sql;
create table mice as
select id,year,sum(r15) as con
from b
group by id,year;
create table abc as
select *
from b,mice
where b.id=mice.id and b.year=mice.year;
quit;
data result;
set abc;
if r15=. then con=.;
run;

proc print data=result;
run;

连续出现>15，还是说只需要知道大于15的个数？

littlepig818 发表于 2015-6-15 09:31
data a;
input id year qtr r;
cards;

这个应该是4    而不是1333

一般方法：
data c;
   set a;
   by id year qtr;
   if first.id then con=.;
   if r<15 then con=.;
   else con+1;
run;
   
proc sort data=c;
   by id descending year descending qtr;
run;

data d;
   set c;
   by id descending year descending qtr;
   retain _con;
   if missing(_con) ne missing(con) then _con=con;
   else con=_con;
   drop _con;
run;

proc sort data=d;
   by id year qtr;
run;