楼主: yayacuiliu
4361 14

[原创博文] 连续数据的个数计算 [推广有奖]

11
cathy3212 发表于 2011-8-18 15:23:40
  1. data aa;
  2.         input Reporter Partner Product year;
  3. cards;
  4. 1 003      30303      1995
  5. 1 004      99284      1995
  6. 1 003      30303      1996
  7. 1 003      30303      1997
  8. 1 003      30303      2000
  9. ;
  10. run;
  11. proc sort data=aa out=aa1;by reporter partner product year;run;/*这里把你的所有要区别的组别放在前面,最后加year*/
  12. data aa2;
  13.         set aa1;
  14.         retain p 0;
  15.         if year-lag(year)>1 then p+1;
  16.         if reporter>lag(reporter) then p+1;
  17.         if partner>lag(partner) then p+1;
  18.         if product>lag(product) then p+1;       
  19. run;
  20. proc sql;
  21.         create table aa3 as select distinct
  22.         reporter
  23.         ,partner
  24.         ,product
  25.         ,count(p) as cnt
  26.         from aa2
  27.         group by p;
  28. quit;
复制代码
已有 1 人评分学术水平 热心指数 收起 理由
yayacuiliu + 1 + 1 很热心,回答很有帮助。

总评分: 学术水平 + 1  热心指数 + 1   查看全部评分

12
yayacuiliu 发表于 2011-8-18 22:40:56
cathy3212 发表于 2011-8-18 15:23
对我很有帮助,谢谢你!!

13
jasonscut 在职认证  发表于 2011-8-18 23:16:45
这样好像可以:

1)  proc sort data=data1 out=data2;
     by Reporter Partner Product  year;     run;
2) proc transpose data=data2 out=data3 prefix=year_;
    by Reporter Partner Product;
    var year;
   run;
3) data data4; set data3;
   array year(13) year1_year13;  (you may change to the maximum num. of year_1-year_max)
   do i=1 to 13;
   if year(i+1)=year(i)+1 and year(i+2)^=year(i+1)+1 then var1=1;
  else  if year(i+1)=year(i)+1 and year(i+2)=year(i+1)+1 and year(i+3)^=year(i+2)+1 then var1=2;
....


第三步里面的if condition 要看你具体的怎么去count而作修改,  可能未必是我上面的情况.


但是这样去做应该是可以行得通的.


Good luck!






北美统计金融博士

14
cathy3212 发表于 2011-8-19 08:43:39
yayacuiliu 发表于 2011-8-18 22:40
对我很有帮助,谢谢你!!
不客气,对你有帮助就好。

15
jingju11 发表于 2011-8-23 21:25:42
  1. proc sort out =have1; by reporter partner product year; run;
  2. data have2;
  3.         set have1;
  4.         by reporter partner product;
  5.         if dif(year) ^=1 then n ++1;
  6.         if first.product then n =1;
  7. run;
  8. data have3; drop year n;
  9.         do until(last.n);
  10.                 set have2;
  11.                 by reporter partner product n;
  12.                 duration =sum(duration, 1);
  13.         end;        
  14. run;
复制代码
不考虑缺失值。如果有,应预先删除。
京剧

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群
GMT+8, 2026-1-1 10:34