数据计算程序求助（百分位数）

n（观测号）	a	v（根据n和a计算）
1	0.1
2	0.1
3	0.1
4	0.1
5	0.2
6	0.2	1-5样本区间a的95%百分位数
7	0.2	2-6样本区间a的95%百分位数
8	0.2	3-7样本区间a的。。
9	0.2	4-8的
10	0.7	5-9的
11	0.7	。。
12	0.7	。。
13	0.7	。。
14	0.7	。。
15	0.7	。。
16	0.7	12-16样本区间a的95%的百分位数

根据n和a计算v。。。。


谢谢！

1. data b;
   
2.         retain x1-x6;
   
3.         set a;
   
4.         array xx[6] x1-x6;
   
5.         do i=6 to 2 by -1;
   
6.                 xx[i]=xx[i-1];
   
7.         end;
   
8.         xx[1]=a;
   
9.         if _N_ >= 6 then do;
   
10.                 v=pctl(95,of x2-x6);
    
11.                 sum=sum(of x2-x6);
    
12.                 max=max(of x2-x6);
    
13.                 min=min(of x2-x6);
    
14.         end;
    
15.         drop x1-x6 i;       
    
16. run;

复制代码

求教：聚合函数对于临时数组怎么写啊？
一时写不出来没用临时数组......
我只会用宏把临时数组显式化来写在聚合函数里面......

data have;
   input n a;
   cards;
1 0.1
2 0.1
3 0.1
4 0.1
5 0.2
6 0.2
7 0.2
8 0.2
9 0.2
10 0.7
11 0.7
12 0.7
13 0.7
14 0.7
15 0.7
16 0.7
;

%macro lags(var,n);
   lag(&var)
   %do i=2 %to &n;
       ,lag&i(&var)
   %end;
%mend;

data wanted;
    set have;
    v=pctl(95,%lags(a,5));       
    if _n_<=5 then call missing(v);
run;


data test;
    do n=1 to 10000;
           a=ranuni(0);
           output;
        end;
run;

data testl;
     set test;
    v=pctl(95,%lags(a,100));       
    if _n_<=100 then call missing(v);
run;

2# soporaeternus

1. data b;
   
2.   set a;
   
3.   array t{5} _temporary_;
   
4. if _n_ >5 then v =pctl(95, of t[*] );
   
5. t[5 -mod(_n_, 5)] =a;
   
6. run;

复制代码

只要是前五个数，其在数组里的位置是无所谓的。如果要根据记录的位置求白分位，n值反而多余。
为了300币。。。当然交流第一。京剧

16 0.7 12-16样本区间a的95%的百分位数

本文来自: 人大经济论坛 SAS专版 版，详细出处参考：http://www.pinggu.org/bbs/viewthread.php?tid=964816&page=1&from^^uid=1032463

按照规律，16 0.7 12-16样本区间a的95%的百分位数，
12-16应该为11-15才对。

京剧前辈，
的确，我做移动窗口方法太繁琐，在求和，求百分位数等方法上，不需要刻意保留之前数据的顺序，看了您的方法，简单高效，受益匪浅。多谢！
现在还有个问题，在我运行if _n_ >5 then v =pctl(95, of t[*] );这句的时候，还是报错
ERROR: ARRAYNAME[*] 说明要求基于数组的变量。
我用非临时数组的话，就没有这个问题，就如我之前问的那样。我不知道对于临时数组应该怎样用聚合函数
还是您使用的是9.2，这是9.2的新特性.......
望赐教，十分感谢

6# soporaeternus
不敢当。HOPELWELL，你知道其中的缘故吗？我用的是9.22
使用pobel的test数据（较小的数据集），结果如下

580  data test;
581      do n=1 to 100;
582             a=ranuni(1);
583             output;
584          end;
585  run;
NOTE: The data set WORK.TEST has 100 observations and 2 variables.
NOTE: DATA statement used (Total process time):
      real time           0.01 seconds
      cpu time            0.01 seconds

586
587
data b;
589     set test;
590     array t{5} _temporary_;
591     if _n_ >5 then v =pctl(95, of t

);
592     t[5 -mod(_n_, 5)] =a;
593  run;
NOTE: There were 100 observations read from the data set WORK.TEST.
NOTE: The data set WORK.B has 100 observations and 3 variables.
NOTE: DATA statement used (Total process time):
      real time           0.03 seconds
      cpu time            0.03 seconds

不知道为什么总是贴不好。

怎么求百分位，就是一个变量的所有值的百分位怎么求？