sas取均值问题

有一笔资料，希望做到：将第五个数据和第十个数据取一个均值，之后第十五笔和第二十笔取均值……以此类推，最终输出均值。

data.xlsx (67.4 KB)

2018-1-7 17:39:23 上传

data

静态的 凑合

1. data a;
   
2. x=1;output;
   
3. x=2;output;
   
4. x=3;output;
   
5. x=4;output;
   
6. x=5;output;
   
7. x=3;output;
   
8. x=5;output;
   
9. x=7;output;
   
10. x=2;output;
    
11. x=1;output;
    
12. x=3;output;
    
13. x=8;output;
    
14. x=7;output;
    
15. x=4;output;
    
16. x=3;output;
    
17. x=5;output;
    
18. x=4;output;
    
19. x=3;output;
    
20. x=1;output;
    
21. x=5;output;
    
22. x=8;output;
    
23. x=9;output;
    
24. x=22;output;
    
25. x=33;output;
    
26. x=32;output;
    
27. x=43;output;
    
28. x=56;output;
    
29. x=44;output;
    
30. x=6;output;
    
31. x=9;output;
    
32. x=13;output;
    
33. x=15;output;
    
34. x=17;output;
    
35. x=15;output;
    
36. x=15;output;
    
37. x=14;output;
    
38. x=13;output;
    
39. x=12;output;
    
40. x=12;output;
    
41. x=10;output;
    
42. run;
    
43. 
    
44. 
    
45. data b1(keep=i x) b2(keep=i x rename=(x=x_));
    
46. set a;
    
47. do i=1 to 99;
    
48.   if _n_=10*i then output b1;
    
49.   if _n_=5*(2*i-1) then output b2;
    
50. end;
    
51. run;
    
52. 
    
53. data want;
    
54. merge b1 b2;
    
55. by i;
    
56. if n(x,x_)=2 then av_x=(x+x_)/2;
    
57. run;
    
58. 
    
59. proc sql noprint;
    
60. create table c as select avg(av_x) as avg from want where av_x ne .;
    
61. quit;

复制代码

lovexialulu 发表于 2018-1-7 20:11
静态的 凑合

您好，能不能麻烦稍微解释一下呢

滞点 发表于 2018-1-8 10:03
您好，能不能麻烦稍微解释一下呢

b1 放的第10 20 30 笔值；
b2 放的第5 15 25 笔值；
merge b1 b2;横向得到第5 10；15 20；25 30；单个block的值

/*分组，按5个观察分组*/
data        tmp;
        do i=1 to 181;
                ee=int(i/5);
                output;
        end;
run;
/*按分组统计均值*/
proc summary data=tmp nway;
        var i;
        class ee;
        output out=tmp1(drop=_: ee) mean=;
run;

流水不朽 发表于 2018-1-9 10:18
/*分组，按5个观察分组*/
data        tmp;
        do i=1 to 181;

您好，可否麻烦帮忙解释一下呢，有点不太明白呀

滞点 发表于 2018-1-14 20:14
您好，可否麻烦帮忙解释一下呢，有点不太明白呀

第一步是：模拟测试数据，使用 int 来求整，这样，只会按照分组间隔产生 新变量（5是你要求的间隔）
第二步是：按照新变量进行 分组求均值

看能否明白？
然后实际应用到你的数据集的时候，用_n_替代 i即可。

你可以step1 留下obs为5的倍数的观测，保存为data1；
step2 直接merge data1 和data1(firstobs=2 rename=(value=value1)) 不需要by  留奇数观测 ，取平均就好了

_n_    value   value1

1  value5 value10
2  value10 value15
3  value15 value20
4 value20 value25
5 value25 value30
...



然后我也不知道你需要什么

DATA A;
        SET DATA_1;  /*DATA_1为原始数据*/
        N=_N_;
RUN;
/*取能被5整除的行*/
DATA B;
        SET A;
        IF MOD(N,5)=0;
RUN;
DATA C;
        SET B;
        OBS=_N_;
RUN;
DATA D;
        SET C;
        FUZHU=LAG(OBS);
        IF MOD(OBS,2)=0 THEN NEW_OBS=FUZHU;
        ELSE NEW_OBS=OBS;
RUN;
/*求平均*/
PROC SQL;
        CREATE TABLE FINAL AS
                SELECT
                        T1.NEW_OBS,
                        T1.TIME,
                        MEAN(PM) AS GROUP_MEAN
                        FROM D AS T1
                        GROUP BY NEW_OBS;
QUIT;