请问如何求得一个数在样本分布中的分位数

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我这里有一张抽样得到的表，
现在我指定了一系列的数，想要得到这些数在原样本中排名的分位数
我感觉类似于求样本分位数的反函数，不知道在SAS里该如何实现，谢谢高手解答

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：分位数 如何实现 反函数 不知道 如何 样本

You may think too much about that. According to empirical distribution function, the sample quantiles should be straightforward. note that the quantiles are not 1 to 1 at all. For example, if you want to get the percentile (in 10th) of weight =102.5 from data set sashelp.class... jingju

1. data _null_;
   
2.         a =102.5;
   
3.         set sashelp.class end =Eof;
   
4.         n ++(weight <=a);
   
5.         if Eof then do;
   
6.                 q =round(n/_n_*100,10);
   
7.                 put a=q=;
   
8.                 end;
   
9.         run;

复制代码

a=102.5 q=60
NOTE: There were 19 observations read from the data set SASHELP.CLASS.

如果抽样表中有多个数, 可以将抽样表与rank后的原表合并.

1. data big;
   
2.    do ID = 1 to 120;
   
3.        a = round(ranuni(123)*100, 0.01);
   
4.        b = round(ranuni(123)*100, 0.01);
   
5.        output;
   
6.    end;
   
7. run;
   
8. 
   
9. proc rank data=big group=100 out=big_q;
   
10.    var a b;
    
11.    ranks quantile_a quantile_b;
    
12. run;
    
13.    
    
14. data sample;
    
15.     do _n_ = 1 to 10;
    
16.       p = int(ranuni(12345)*120);
    
17.       set big point=p;
    
18.       output;
    
19.     end;
    
20.     stop;
    
21. run;
    
22. proc sql;
    
23. create table sample_q as
    
24.     select b.*
    
25.     from sample a, big_q b
    
26.     where a.id = b.id;
    
27. quit;

复制代码

还可以从rank后的原表直接抽样.

jingju11 发表于 2014-2-20 11:09
You may think too much about that. According to empirical distribution function, the sample quantile ...

多谢回复！

后来我在R里找到了ecdf函数，可以指定一张表生成经验概率密度函数，然后把这个函数赋给一个类似句柄的东西，然后就可以任意指定观测求分位数了

SAS好像没有类似的函数，是不是要自己写了？谢谢！

高人，学习了。