a <- c(4,6,3,9,7,5,5,8,6,5,7,4,9,7,3,5,7,3,7,8,6,4,9,3,8,6,2,7)
b <-c(1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3)
pf <- factor(b)
tapply(a,pf,quantile)

新手学习中，如有错误请各位指正

谢谢，请问有没有办法算不是标准的分位数，比如说20%分位数？

tapply(a,pf,function(x) {return(quantile(x,probs=c(0,0.2,0.4,0.6,0.8,1)))})
probs 任意自己设定

自己编写一下关于quantile的函数就可以实现了，比如我要实现20%分位数统计：编写一个qt20的函数，然后应用到dataorz的代码中即可。代码如下：

聪明！