https://communities.sas.com/t5/SAS-Procedures/The-PCTL-function-in-PROC-IML/m-p/64210#M18253 The IML manual says that the data step function PCTL will work in PROC<BR /> IML. Great, except I can make it work in a useful way.<BR /> <BR /> If x is a 10x1 vector, then<BR /> <BR /> median=pctl(50,x);<BR /> <BR /> doesn't really work the way I would hope it works. It gives me a 10x1<BR /> vector in return, the median of each value in the original matrix x.<BR /> How non-useful!<BR /> <BR /> It does work properly if I write<BR /> <BR /> median=pctl(50,x[1],x[2],x[3],x[4],x[5],x[6],x[7],x[8],x[9],x[10]);<BR /> <BR /> but this isn't a useful way to do things for general problems where x<BR /> can be any size.<BR /> <BR /> So, is there a way to have the PCTL function work in a useful way in<BR /> PROC IML? Wed, 03 Dec 2008 20:39:14 GMT Paige 2008-12-03T20:39:14Z https://communities.sas.com/t5/SAS-Procedures/The-PCTL-function-in-PROC-IML/m-p/64210#M18253 The IML manual says that the data step function PCTL will work in PROC<BR /> IML. Great, except I can make it work in a useful way.<BR /> <BR /> If x is a 10x1 vector, then<BR /> <BR /> median=pctl(50,x);<BR /> <BR /> doesn't really work the way I would hope it works. It gives me a 10x1<BR /> vector in return, the median of each value in the original matrix x.<BR /> How non-useful!<BR /> <BR /> It does work properly if I write<BR /> <BR /> median=pctl(50,x[1],x[2],x[3],x[4],x[5],x[6],x[7],x[8],x[9],x[10]);<BR /> <BR /> but this isn't a useful way to do things for general problems where x<BR /> can be any size.<BR /> <BR /> So, is there a way to have the PCTL function work in a useful way in<BR /> PROC IML? Wed, 03 Dec 2008 20:39:14 GMT https://communities.sas.com/t5/SAS-Procedures/The-PCTL-function-in-PROC-IML/m-p/64210#M18253 Paige 2008-12-03T20:39:14Z https://communities.sas.com/t5/SAS-Procedures/The-PCTL-function-in-PROC-IML/m-p/64211#M18254 You can get the median (50th pctl) and the 25th and 75th percentiles by using the MEDIAN or QUARTILE modules in IMLMLIB:<BR /> <BR /> proc iml;<BR /> a = {2 3 2 5 6 8 9 4 12 5}`;<BR /> med = median(a);<BR /> print med;<BR /> <BR /> q = quartile(a);<BR /> print q[rowname={"min" "Q1" "median" "Q3" "max"}];<BR /> quit;<BR /> <BR /> For more general percentiles, you can use the SORT call and the RANK function to rank your data. Fri, 29 May 2009 20:12:03 GMT https://communities.sas.com/t5/SAS-Procedures/The-PCTL-function-in-PROC-IML/m-p/64211#M18254 Rick_SAS 2009-05-29T20:12:03Z