如何在SAS/IML中提取某些特定的分位数

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

请问，现在有一批数据，如何在SAS/IML中提取某些特定的分位数，请高手指点。谢谢！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：分位数 高手指点 如何

在proc iml;这个statement前添加ods trace on;
run;这个statement后添加ods trace off;
从新run; 然后看log, 找出你要分位数.
最后用ods select ***分位数名称***;

PROC IML;
   scores={90,62,66,68,70,72,73,74,78,78,78,79,80,81,82,
           82,82,84,84,85,85,85,85,87,88,89,89,89,89,61};
   PRINT scores;
   quartiles=QUARTILE(scores);
   PRINT quartiles;
/** Start Percentile **/
START percentile(x,pct);
   CALL SORT(x,1);
   count=NROW(x);
   pct=pct/100*(count+1);
   fpct=FLOOR(pct);
   cpct=CEIL(pct);
   IF fpct=cpct THEN percentile=x[pct];
   ELSE percentile=x[fpct]+(pct-fpct)* (x[cpct]-x[fpct]);
/* SAS PCTLDF */
RETURN(percentile);
FINISH;
percentile=PERCENTILE(scores,25);
/** End Percentile **/


CALL SORT(scores,1); /* sorts the scores in ascending order */
count=NROW(scores);  
quantile=scores[CEIL(count*.25)];
PRINT quartiles, "25th Percentile:"
      percentile, "25th Quantile:" quantile;
QUIT;

非常感谢，但我以前记得有非常简单的程序来实现，现在找不到了，能否有更简单的呢？

good answer....jingju