求向量quantile in sas iml？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

<P>急着请教：如何在sas iml 里计算一个向量的%25 quantile? funtion PCTL 可以用，但必须是用逗号隔开的几各数，不能求一个向量的。 请朋友们帮忙！谢谢！</P>

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：quantile quant Tile Ant PCT 朋友 如何

在等待答案中。学习学习

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;