How to use prctile in grpstats?

Question

Lei 2013-5-15

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/75902-how-to-use-prctile-in-grpstats

Hello, I am using the grpstats function to get some statistic based on groups; however, I have trouble to use prctile for the whichstats. As the prctile requires a both input X and p. The p I want is [ 5 95]. But for the input x, I don't know how to handle with this as I have predefined X in the grpstats. Here is my grpstats input example:

[A,B,C,D,E,F] = grpstats(sample(:,2),sample(:,7),{'gname','numel','mean','var','skewness',@iqr});

Sample(:,2) is my input data, and sample(:,7) is the group. I am wondering how can I use the prctile in the grpstats properly. Thanks.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Tom Lane 2013-5-16

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/75902-how-to-use-prctile-in-grpstats#answer_85696

在 MATLAB Online 中打开

The trick is to use an anonymous function that has the desired percentiles built in. For example, the following two calls to grpstats yield the same results:

x = randn(1000,2);
g = randi(5,1000,1);
grpstats(x,g, @median)
grpstats(x,g, @(x) prctile(x,50))

The following yields a 3-D array where the middle page, a(:,:,2), has the median and the other pages have the 5% and 90% points:

a = grpstats(x,g, @(x) prctile(x,[5 50 90]))

3 个评论
显示 1更早的评论隐藏 1更早的评论

Tom Lane 2013-5-16

That's called an anonymous function. It's a way to define a function on the fly, at the command line, without entering code into a file. So "@(x)prctile(x,50)" is a function of x that is computed by invoking the prctile function with x as the first input and 50 as the second.

Lei 2013-5-16