如何在R数据框中找到组的百分等级？
百分位一词是指低于或高于百分位值的百分比。例如，如果我们有一个位于第50个百分位数的值，则可以说有50％的值低于或高于该值。此处的值50称为百分等级。要找到R数据帧中组的百分数等级，我们可以使用dplyr包的mutate函数。
示例
考虑以下数据帧-
Group<-sample(LETTERS[1:4],20,replace=TRUE) Response<-rpois(20,5) df1<-data.frame(Group,Response)df1
输出结果
Group Response 1 D 5 2 B 7 3 D 5 4 C 4 5 D 5 6 C 5 7 A 10 8 D 3 9 B 2 10 D 0 11 B 4 12 D 5 13 A 3 14 A 6 15 D 2 16 A 7 17 A 6 18 C 2 19 A 9 20 C 3
示例
加载dplyr软件包-
library(dplyr)
查找组的响应百分等级-
示例
df1%>%group_by(Group)%>%mutate(Percentile_Rank=rank(Response)/length(Response)) # A tibble: 20 x 3 # Groups: Group [4]
输出结果
Group Response Percentile_Rank <chr> <int> <dbl> 1 D 5 0.786 2 B 7 1 3 D 5 0.786 4 C 4 0.75 5 D 5 0.786 6 C 5 1 7 A 10 1 8 D 3 0.429 9 B 2 0.333 10 D 0 0.143 11 B 4 0.667 12 D 5 0.786 13 A 3 0.167 14 A 6 0.417 15 D 2 0.286 16 A 7 0.667 17 A 6 0.417 18 C 2 0.25 19 A 9 0.833 20 C 3 0.5
示例
Class<-sample(c("I","II","III"),20,replace=TRUE)
Y<-rnorm(20,25,3.27)
df2<-data.frame(Class,Y)df2输出结果
Class Y 1 III 32.88152 2 III 23.35048 3 III 19.78199 4 III 26.05137 5 I 26.16563 6 III 20.30466 7 I 22.93382 8 II 30.03620 9 I 16.89365 10 I 27.33329 11 I 27.46550 12 III 27.59028 13 II 27.40766 14 III 23.29442 15 II 28.69237 16 II 31.25723 17 II 22.58002 18 III 22.48583 19 I 26.08357 20 III 24.51681
查找类别的响应百分等级-
示例
df2%>%group_by(Class)%>%mutate(Percentile_Rank=rank(Y)/length(Y)) # A tibble: 20 x 3 # Groups: Class [3]
输出结果
Class Y Percentile_Rank <chr> <dbl> <dbl> 1 III 32.9 1 2 III 23.4 0.556 3 III 19.8 0.111 4 III 26.1 0.778 5 I 26.2 0.667 6 III 20.3 0.222 7 I 22.9 0.333 8 II 30.0 0.8 9 I 16.9 0.167 10 I 27.3 0.833 11 I 27.5 1 12 III 27.6 0.889 13 II 27.4 0.4 14 III 23.3 0.444 15 II 28.7 0.6 16 II 31.3 1 17 II 22.6 0.2 18 III 22.5 0.333 19 I 26.1 0.5 20 III 24.5 0.667