求助如何按条件统计缺省值

如图，求每个stkcd的所有缺省值跟连续的缺省值。
比如对于002067，总共的dcost缺省为2，07与12缺失。而连续缺省值，07年取1,12年取1，对于002066，则08,09,10的连续缺省，值是3

小鳄鱼a 发表于 2014-8-31 08:48
不好意思，我试了试,totalmiss全部为0

不好意思，
第10行后加入一行
if missing(dcost) then totalMiss + 1;

data a;
        set 你的数据；
        where Dcost=.;
run;

proc sort data=a out=b;
        by stkcd year;
run;

data c;
        set b;
        if first.stkcd then zong_quesheng=0;
        zong_quesheng+1;
        if last.stkcd;
run;

data d;
        set b;
        lag_year=lag(year);
        if first.stkcd then lianxu_quesheng=1;
        if year=lag_year+1 then lianxu_quesheng+1;
        if last.stkcd;       
run;

data final;
        merge c d;
        by stkcd year;
run;


-----------------------
随手写来，如果报错请告知

quekehanmu 发表于 2014-8-30 09:51
data a;
        set 你的数据；
        where Dcost=.;

first.stkcd没有定义   传一份数据吧

小鳄鱼a 发表于 2014-8-30 10:22
first.stkcd没有定义   传一份数据吧

你把“你的数据” 改为你导入的sas数据名就可以了

小鳄鱼a 发表于 2014-8-30 10:22
first.stkcd没有定义   传一份数据吧

data a;
        set 你的数据；
        where Dcost=.;
run;

proc sort data=a out=b;
        by stkcd year;
run;

data c;
        set b;
        by stkcd year;
        if first.stkcd then zong_quesheng=0;
        zong_quesheng+1;
        if last.stkcd;
run;

data d;
        set b;
        by stkcd year;
        lag_year=lag(year);
        if first.stkcd then lianxu_quesheng=1;
        if year=lag_year+1 then lianxu_quesheng+1;
        if last.stkcd;      
run;

data final;
        merge c d;
        by stkcd year;
run;

1. /* try this */
   
2. data want;
   
3.     array dvalue{9} _temporary_;
   
4.      i = 1;
   
5.     do until (last.stkcd);
   
6.        set temp;
   
7.        by stkcd;
   
8.        dvalue[i] = dcost;  
   
9.        i+1;
   
10.     end;
    
11. 
    
12.         flag = 0;
    
13.     do i = 1 to 9;                     
    
14.        if missing(dvalue[i]) then do;
    
15.             if flag = 0 then do;  flag = 1; nmiss = 1;  s1 = i;  s2 = i;
    
16.                             end;   
    
17.        else if flag = 1 then  do; nmiss + 1; s2+1;
    
18.                              end;
    
19.        end;
    
20. 
    
21.        if (^missing(dvalue[i]) and flag=1) or (i=9 and flag=1)then do;
    
22.             do j = s1 to s2;
    
23.                dvalue[j] =nmiss;      
    
24.             end;
    
25.               flag=0;
    
26.         end;  
    
27.     end;
    
28. 
    
29.    i = 1;
    
30.    do until (last.stkcd);
    
31.       set temp;
    
32.       by stkcd;
    
33.       dcost= dvalue[i];
    
34.       output;
    
35.    i + 1;
    
36.    end;
    
37.    drop i j flag s1 s2 nmiss;
    
38. run;

复制代码

yongyitian 发表于 2014-8-30 11:03

不好意思，是产生两个新的变量来分别计算dcost的总的缺失值跟连续的缺失值，而不是用连续的缺失值将dcost覆盖掉

小鳄鱼a 发表于 2014-8-30 14:36
不好意思，是产生两个新的变量来分别计算dcost的总的缺失值跟连续的缺失值，而不是用连续的缺失值将dcost ...

1. data want;
   
2.     array dvalue{9} _temporary_;
   
3.     array cmiss{9}  _temporary_ ;
   
4.     do i = 1 to 9; cmiss[i]=0; end;
   
5.                i = 1;
   
6.        totalMiss = 0;
   
7.     do until (last.stkcd);
   
8.        set temp;
   
9.        by stkcd;
   
10.        dvalue[i] = dcost;  
    
11.        i+1;
    
12.     end;
    
13. 
    
14.         flag = 0;
    
15.     do i = 1 to 9;                     
    
16.        if missing(dvalue[i]) then do;
    
17.             if flag = 0 then do;  flag = 1; nmiss = 1;  s1 = i;  s2 = i;
    
18.                             end;   
    
19.        else if flag = 1 then  do; nmiss + 1; s2+1;
    
20.                              end;
    
21.        end;
    
22. 
    
23.        if (^missing(dvalue[i]) and flag=1) or (i=9 and flag=1)then do;
    
24.             do j = s1 to s2;
    
25.                cmiss[j] = nmiss;      
    
26.             end;
    
27.               flag=0;
    
28.         end;  
    
29.     end;
    
30. 
    
31.    i = 1;
    
32.    do until (last.stkcd);
    
33.       set temp;
    
34.       by stkcd;
    
35.         continue_miss = cmiss[i];
    
36.       output;
    
37.    i + 1;
    
38.    end;
    
39.    drop i j flag s1 s2 nmiss;
    
40. run;

复制代码

yongyitian 发表于 2014-8-30 23:08

不好意思，我试了试,totalmiss全部为0