如何在 r 中使用 quantile() 函数
在统计学中,分位数是将分类数据集分为相等组的值。
R 中的quantile()函数可用于计算数据集的样本分位数。
该函数使用以下基本语法:
分位数(x, 概率 = seq(0, 1, 0.25), na.rm = FALSE)
金子:
- x :向量名称
- probs :概率的数值向量
- na.rm :是否删除 NA 值
以下示例展示了如何在实践中使用此功能。
示例 1:计算向量的分位数
以下代码展示了如何在 R 中计算向量的分位数:
#define vector of data data = c(1, 3, 3, 4, 5, 7, 8, 9, 12, 13, 13, 15, 18, 20, 22, 23, 24, 28) #calculate quartiles quantile(data, probs = seq(0, 1, 1/4)) 0% 25% 50% 75% 100% 1.0 5.5 12.5 19.5 28.0 #calculate quintiles quantile(data, probs = seq(0, 1, 1/5)) 0% 20% 40% 60% 80% 100% 1.0 4.4 8.8 13.4 21.2 28.0 #calculate deciles quantile(data, probs = seq(0, 1, 1/10)) 0% 10% 20% 30% 40% 50% 60% 70% 80% 90% 100% 1.0 3.0 4.4 7.1 8.8 12.5 13.4 17.7 21.2 23.3 28.0 #calculate random quantiles of interest quantile(data, probs = c(.2, .5, .9)) 20% 50% 90% 4.4 12.5 23.3
示例 2:计算数据框中列的分位数
以下代码显示如何计算数据框中特定列的分位数:
#create data frame
df <- data. frame (var1=c(1, 3, 3, 4, 5, 7, 7, 8, 12, 14, 18),
var2=c(7, 7, 8, 3, 2, 6, 8, 9, 11, 11, 16),
var3=c(3, 3, 6, 6, 8, 4, 4, 7, 10, 10, 11))
#calculate quartiles of column 'var2'
quantile(df$var2, probs = seq(0, 1, 1/4))
0% 25% 50% 75% 100%
2.0 6.5 8.0 10.0 16.0
我们还可以使用sapply()函数一次计算多列的分位数:
#calculate quartiles of every column
sapply(df, function (x) quantile(x, probs = seq(0, 1, 1/4)))
var1 var2 var3
0% 1.0 2.0 3
25% 3.5 6.5 4
50% 7.0 8.0 6
75% 10.0 10.0 9
100% 18.0 16.0 11
示例 3:按组计算分位数
以下代码演示了如何使用dplyr包中的函数通过分组变量计算分位数:
library (dplyr)
#define data frame
df <- data. frame (team=c('A', 'A', 'A', 'A', 'B', 'B', 'B', 'B', 'C', 'C', 'C') ,
points=c(1, 3, 3, 4, 5, 7, 7, 8, 12, 14, 18))
#define quantiles of interest
q = c(.25, .5, .75)
#calculate quantiles by grouping variable
df %>%
group_by(team) %>%
summarize(quant25 = quantile(points, probs = q[1]),
quant50 = quantile(points, probs = q[2]),
quant75 = quantile(points, probs = q[3]))
# A tibble: 3 x 4
team quant25 quant50 quant75
1 to 2.5 3 3.25
2 B 6.5 7 7.25
3 C 13 14 16
其他资源
以下教程展示如何使用quantile()函数计算其他常见分位数值: