求助，怎么用 SAS 对 极端值 处理？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

在用SAS处理数据时需要对极端值进行处理，对 最大和最小的5%数据 进行重新赋值，分别 等于95%和5%的那个值，请问

这个程序要怎么写啊？用排序吗，还是有其他更简单的方法？

        谢谢各位了！！！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏3 回帖

关键词：极端值 SAS处理数据 处理数据 程序

用proc means或者proc univariate求出95%和5%的值，然后进行处理

试了下，好像只能算 p1 和 p99的，换成p0.5和 p99.5 就跑不出来啊？

ccgzl 发表于 2010-6-27 12:32
试了下，好像只能算 p1 和 p99的，换成p0.5和 p99.5 就跑不出来啊？

用proc univariate

data Crackman;
do i=1 to 100;
x=10.67+3.178*rannor(4);
output;
end;
run;
proc univariate data=Crackman ;
var x;
output out=Result pctlpts=0.5,99.5  pctlpre=per;
run;

多谢了，又学到了很多，呵呵

这个macro很好用

来源：http://ihome.cuhk.edu.hk/~b121456/tools.html

1. /**********************************************************
   
2. Purpose: trim or winsorize SAS dataset to remove the impact from extreme values
   
3. 
   
4. Input
   
5. dsetin  : dataset to winsorize/trim
   
6. byvar   : define subset to winsorize/trim,e.g. 'date'. type 'none' for the whole dataset
   
7. type    : 'delete' or 'winsor' (delete will trim, winsor will winsorize
   
8. vars    : subsetting variables to winsorize/trim on; type 'none for no byvar
   
9. pctl    : the percenagte of left and right tails to trim/winsorize
   
10. 
    
11. Output
    
12. dsetout : dataset to output with winsorized/trimmed values
    
13. ************************************************************/
    
14. 
    
15. %macro winsor(dsetin=, dsetout=, byvar=none, vars=, type=winsor, pctl=5 95);
    
16. 
    
17. %if &dsetout = %then %let dsetout = &dsetin;
    
18.    
    
19. %let varL=;
    
20. %let varH=;
    
21. %let xn=1;
    
22. 
    
23. %do %until ( %scan(&vars,&xn)= );
    
24.     %let token = %scan(&vars,&xn);
    
25.     %let varL = &varL &token.L;
    
26.     %let varH = &varH &token.H;
    
27.     %let xn=%EVAL(&xn + 1);
    
28. %end;
    
29. 
    
30. %let xn=%eval(&xn-1);
    
31. 
    
32. data xtemp;
    
33.     set &dsetin;
    
34.     run;
    
35. 
    
36. %if &byvar = none %then %do;
    
37. 
    
38.     data xtemp;
    
39.         set xtemp;
    
40.         xbyvar = 1;
    
41.         run;
    
42. 
    
43.     %let byvar = xbyvar;
    
44. 
    
45. %end;
    
46. 
    
47. proc sort data = xtemp;
    
48.     by &byvar;
    
49.     run;
    
50. 
    
51. proc univariate data = xtemp noprint;
    
52.     by &byvar;
    
53.     var &vars;
    
54.     output out = xtemp_pctl PCTLPTS = &pctl PCTLPRE = &vars PCTLNAME = L H;
    
55.     run;
    
56. 
    
57. data &dsetout;
    
58.     merge xtemp xtemp_pctl;
    
59.     by &byvar;
    
60.     array trimvars{&xn} &vars;
    
61.     array trimvarl{&xn} &varL;
    
62.     array trimvarh{&xn} &varH;
    
63. 
    
64.     do xi = 1 to dim(trimvars);
    
65. 
    
66.         %if &type = winsor %then %do;
    
67.             if not missing(trimvars{xi}) then do;
    
68.               if (trimvars{xi} < trimvarl{xi}) then trimvars{xi} = trimvarl{xi};
    
69.               if (trimvars{xi} > trimvarh{xi}) then trimvars{xi} = trimvarh{xi};
    
70.             end;
    
71.         %end;
    
72. 
    
73.         %else %do;
    
74.             if not missing(trimvars{xi}) then do;
    
75.               if (trimvars{xi} < trimvarl{xi}) then delete;
    
76.               if (trimvars{xi} > trimvarh{xi}) then delete;
    
77.             end;
    
78.         %end;
    
79. 
    
80.     end;
    
81.     drop &varL &varH xbyvar xi;
    
82.     run;
    
83. 
    
84. %mend winsor;

复制代码

BraveMadMan 发表于 2010-6-28 02:45
这个macro很好用

来源：http://ihome.cuhk.edu.hk/~b121456/tools.html

好人啊！正需要呢！！！

其实也可以求出5%的样本数是a，然后用winsorized=a

BraveMadMan 发表于 2010-6-28 02:45
这个macro很好用

来源：http://ihome.cuhk.edu.hk/~b121456/tools.html

是直接复制粘贴吗？为什么用过了没有效果？