求计算状态转移

如图所示，我已经按照每年每月，按照monret指标分好五组。现在想计算出的是一个组间股票转移的情况。具体的步骤为，在t月份，按照当期的monret分成五组0,1，2,3,4;在t+1月，计算还有多少比率的股票仍然在t时期所在组；第t+2月，计算还有多少比率的股票仍然在t时期所在组；一直计算到t+4期。然后按照每个月分别在未来1,2,3,4月的情况，求出这种转移率的平均

最后希望能呈现大概如下的样子
thanks

1. proc import out=input
   
2.         datafile="C:\tmp\zhuanyi.xls"
   
3.         dbms=excel replace;
   
4. run;
   
5. 
   
6. proc sql;
   
7.         create table out as
   
8.         select distinct a.year,a.month,
   
9.         sum((a.year=b.year)*(a.month=b.month)) as tot,
   
10.         sum(((b.year-a.year)*12+b.month-a.month=1)*(a.rflag=b.rflag)) as cnt_1,
    
11.         sum(((b.year-a.year)*12+b.month-a.month=2)*(a.rflag=b.rflag)) as cnt_2,
    
12.         sum(((b.year-a.year)*12+b.month-a.month=3)*(a.rflag=b.rflag)) as cnt_3,
    
13.         sum(((b.year-a.year)*12+b.month-a.month=4)*(a.rflag=b.rflag)) as cnt_4
    
14.         from input a
    
15.         left join input b on a.stkcd=b.stkcd
    
16.         group by a.year,a.month
    
17.         order by a.year,a.month;
    
18. quit;
    
19. 
    
20. proc sql;
    
21.         select max(year),max(month) into :max_year,:max_mth from input;
    
22. quit;
    
23. 
    
24. data out1;
    
25.         retain time tot rate_1-rate_4;
    
26.         keep time tot rate_1-rate_4;
    
27.         set out;
    
28.         time=compress(year||put(month,Z2.));
    
29.         array cnt(*) cnt_1-cnt_4;
    
30.         array rate(*) rate_1-rate_4;
    
31.         do i=1 to dim(cnt);
    
32.                 if (&max_year-year)*12+&max_mth-month<=i-1 then rate(i)=.;else
    
33.                 rate(i)=cnt(i)/tot;
    
34.         end;
    
35.         format rate_1-rate_4 percent12.2;
    
36. run;

复制代码

rflag就是根据每年每月对monret分好组的

freerunning_sky 发表于 2014-8-22 14:02

不好意思，可能我没有表达清楚，我再想想怎么表达

1. /* 和楼上的方法一样，加了两个变量，看起来好懂些。　
   
2. sum（condition）函数（）内部的条件满足时值为 1，不满足时为0. */
   
3. 
   
4. 
   
5. PROC IMPORT OUT= WORK.zhuanyi_all
   
6.    DATAFILE= "C:\temp\zhuanyi.xls"
   
7.                DBMS = EXCEL REPLACE;
   
8.               SHEET = "Sheet1$";
   
9. RUN;
   
10. 
    
11. data test_all;
    
12.    set zhuanyi_all;
    
13.    ym = mdy(month, 1, year);   
    
14.    id = cats(stkcd, '_', put(rflag, 1.));
    
15.    format ym yymmdd10.;
    
16. run;
    
17. 
    
18. proc sql;
    
19.     create table t1  as
    
20.     select  a.ym, sum( intck('month', a.ym, b.ym)=0 ) as t0,
    
21.                   sum( intck('month', a.ym, b.ym)=1 ) as t1,
    
22.                   sum( intck('month', a.ym, b.ym)=2 ) as t2,
    
23.                   sum( intck('month', a.ym, b.ym)=3 ) as t3,
    
24.                   sum( intck('month', a.ym, b.ym)=4 ) as t4
    
25.     from test_all a, test_all b
    
26.     where a.id=b.id
    
27.     group by a.ym;
    
28. quit;

复制代码