data b;                                                                                                                                 
set a;                                                                                                                                 
x=substr(left(date),5,2);                                                                                                               
if x in (12,01,02) then y=1;                                                                                                            
if x in (03,04,05) then y=2;                                                                                                            
if x in (06,07,08) then y=3;                                                                                                            
if x in (09,10,11) then y=4;                                                                                                            
year=substr(left(date),1,4);                                                                                                            
if x =12 then year=year+1;                                                                                                              
jidu=trim(left(stkcd))||trim(left(year))||trim(left(y));                                                                                
run;                                                                                                                                    
proc transpose data=b out=c prefix=V;                                                                                                   
    var msmvttl;                                                                                                                        
    by jidu;                                                                                                                           
run;                                                                                                                                    
proc transpose data=b out=c1 prefix=p;                                                                                                  
    var pb;                                                                                                                             
    by jidu;                                                                                                                           
run;                                                                                                                                    
data d;                                                                                                                                 
merge c c1;                                                                                                                             
by jidu;                                                                                                                                
summsm = sum(of v1-v3);                                                                                                                 
sumpb = sum(of p1-p3);                                                                                                                  
lagmsm=lag(summsm);                                                                                                                     
lagsumpb=lag(sumpb);                                                                                                                    
keep jidu summsm sumpb lagmsm lagsumpb;                                                                                                
run;                                                                                                                                    
                                                                                                                                       
                                                                                                                                       
                                                                                                                                       
                                                                                                                                       
data e;                                                                                                                                 
set d;                                                                                                                                 
z =dif(substr(jidu,1,1));                                                                                                               
if z^=0  then lagmsm=. ;                                                                                                               
if z^=0  then  lagsumpb=. ;                                                                                                            
run;                                                                                                                                    
                                                                                                                                       
                                                                                                                                       
data fin;                                                                                                                              
merge b e;                                                                                                                              
by jidu;                                                                                                                                
keep stkcd date msmvttl pb lagmsm lagsumpb;                                                                                             
run;                                                                                                                                    
data fa;                                                                                                                                
set fin;                                                                                                                                
gengxin =  (msmvttl/lagmsm)*lagsumpb;                                                                                                   
run;     

data fin;                                                                                                                              
merge b e;                                                                                                                              
by jidu;                                                                                                                                
keep stkcd date msmvttl pb lagmsm lagsumpb;                                                                                             
run;                                                                                                                                    
data fa;                                                                                                                                
set fin;                                                                                                                                
gengxin =  (msmvttl/lagmsm)*lagsumpb;                                                                                                   
run;     

大神，感谢你的编程，我程序的数据有20多万个观测，而且参数也比这个多得多，不知道你能不能把你的想法回复我一下，你编得程序有的地方我不是很懂，万分感谢

初学SAS方法比较笨。。。。。
data b
把月份按季度归类 把stkcd变成stkcd年季度的格式
data d
算出每个月的值，在用lag返回上个月的值
data e
由于是返回上一个值，需要把下一个stkcd的第一条记录设为空

最后就把上个月的值合并到这个月后面就算出来了

第一jidu=trim(left(stkcd))||trim(left(year))||trim(left(y)); ？
第二if z^=0  then lagmsm=. ; 中 z^=0 是什么意思？这一步是什么意思？
第三z =dif(substr(jidu,1,1));？
这三步是什么意思啊，有点儿疑惑