保留相邻的最大（小）值

有如图的数据（SAS数据格式见附件）:


现在high列中1和-1不是交叉出现的数据点，并根据p_sm的值重新定义high。

具体来讲，我要在high列两个或更多相邻为1的，取取p_sm最大的一个，将high重新定义为1，其他的定义为. (missing); 类似，我要在high列两个或更多相邻为-1的，取出p_sm最小的一个，将high重新定义为-1，其他的定义为 . (missing)

比如，在附件的数据里（以从第一个到第19个观测值为例），要将第2,13,15,19个观测值改为 “.”， 其他不变。

请问在SAS如何实现。

1. * method 1;
   
2. data want;
   
3.         retain min max;
   
4.         do until (last.high);
   
5.                 set tmp1.have;
   
6.                 by high notsorted;
   
7.                 if first.high then do;
   
8.                         min=p_sm;
   
9.                         max=p_sm;
   
10.                 end;
    
11.                 else do;
    
12.                         min=min(min,p_sm);
    
13.                         max=max(max,p_sm);
    
14.                 end;
    
15.         end;
    
16.         do until (last.high);
    
17.                 set tmp1.have;
    
18.                 by high notsorted;
    
19.                 if high=1 and p_sm^=max then high=.;
    
20.                 if high=-1 and p_sm^=min then high=.;
    
21.                 output;
    
22.         end;
    
23. run;

复制代码

2# elek.me 一直没彻底明白notsorted的使用，看了这个代码，彻底明白了。

1. * method 2;
   
2. data tmp1;
   
3.         set tmp1.have;
   
4.         obs=_n_;
   
5. run;
   
6. proc means data=tmp1 noprint idmin;
   
7.         by high notsorted;
   
8.         var p_sm;
   
9.         id  p_sm obs;
   
10.         output out=tmp2(drop=_freq_ _type_) min(p_sm)=min max(p_sm)=max;
    
11. run;
    
12. proc means data=tmp1 noprint;
    
13.         by high notsorted;
    
14.         var p_sm;
    
15.         id  p_sm obs;
    
16.         output out=tmp3(drop=_freq_ _type_) min(p_sm)=min max(p_sm)=max;
    
17. run;
    
18. data want;
    
19.         merge tmp1 tmp2 tmp3;
    
20.         by obs;
    
21.         if high=1 and max^=p_sm then high=.;
    
22.         if high=-1 and min^=p_sm then high=.;
    
23.         drop min max obs;
    
24. run;

复制代码

baoaibaobao 发表于 2011-1-24 15:58
2# elek.me 一直没彻底明白notsorted的使用，看了这个代码，彻底明白了。

呵呵。 其实挺简单。
by 和 class 在sas里的关系，有些微妙。

elek.me 发表于 2011-1-24 16:11

1. * method 2;
   
2. data tmp1;
   
3.         set tmp1.have;
   
4.         obs=_n_;
   
5. run;
   
6. proc means data=tmp1 noprint idmin;
   
7.         by high notsorted;
   
8.         var p_sm;
   
9.         id  p_sm obs;
   
10.         output out=tmp2(drop=_freq_ _type_) min(p_sm)=min max(p_sm)=max;
    
11. run;
    
12. proc means data=tmp1 noprint;
    
13.         by high notsorted;
    
14.         var p_sm;
    
15.         id  p_sm obs;
    
16.         output out=tmp3(drop=_freq_ _type_) min(p_sm)=min max(p_sm)=max;
    
17. run;
    
18. data want;
    
19.         merge tmp1 tmp2 tmp3;
    
20.         by obs;
    
21.         if high=1 and max^=p_sm then high=.;
    
22.         if high=-1 and min^=p_sm then high=.;
    
23.         drop min max obs;
    
24. run;

复制代码

楼主的数据里 date 这个变量，看上去是唯一的而且是排序的，如果确实是的话，就没必要用obs这个变量了，程序可以更简练。

baoaibaobao 发表于 2011-1-24 15:58
2# elek.me 一直没彻底明白notsorted的使用，看了这个代码，彻底明白了。

确实，又学了一招

支持elek.me
多多学习

我尝试着写了一个很笨拙的程序来解决这个问题，由于太笨拙了不怒好意思拿出手，过完年再给ELEK师兄看看，呵呵

1. data have2;
   
2. set have;
   
3. by high notsorted;
   
4. if first.high then group+1;
   
5. run;
   
6. data have3(drop=max min cnt type);
   
7.   max=-999;min=999;
   
8.   do _n_=1 by 1 until(last.group);
   
9.      set have2;
   
10.       by group notsorted;
    
11.       max=max(p_sm,max);
    
12.           min=min(p_sm,min);
    
13.   end;
    
14.   cnt=_n_;
    
15.   type=high;
    
16.   do _n_=1 to _n_;
    
17.    set have2;
    
18.       if cnt=1 then output;
    
19.          else if type=1 then do;
    
20.                  high=ifn(p_sm=max,1,.);output;
    
21.                   end;       
    
22.              else do;high=ifn(p_sm=min,1,.);output;end;
    
23.   end;
    
24. run;
    
25. 
    
26.   
    
27. 
    
28.    

复制代码