我最近发现,在用“ [”子集对象(即数据帧)之后,结果对象可能是同一行代码中“ [”的子集(我应该早就意识到了!)。这是一个例子:
# Create a data frame
df1 <- as.data.frame(matrix(1:9, nrow = 3))
# Take a look at the data frame
df1
V1 V2 V3
1 1 4 7
2 2 5 8
3 3 6 9
# If I want the value which is on the 3rd row and 2nd column
df1[3,2]
[1] 6
# But I could also
df1[,2][3]
[1] 6
关于第二种选择的几句话。df[,2]
返回原子向量,然后是的子集df[,2][3]
。
以下数据框将有助于说明我的问题。这是一个简单的数据框,其中包含26名学生的姓名,他们各自的系以及一个数值。添加种子编号以提高可重复性。
set.seed(123)
df2 <- data.frame(name = letters, dept = sample(c("econ", "stat", "math"), 26, replace = TRUE), value = runif(26, 0, 100))
head(df2)
name dept value
1 a econ 54.40660
2 b math 59.41420
3 c stat 28.91597
4 d math 14.71136
5 e math 96.30242
6 f econ 90.22990
我想知道谁在econ
部门中价值最低。我尝试的第一件事是:
df2[df2$dept == "econ" & df2$value == min(df2$value),]
[1] name dept value
<0 rows> (or 0-length row.names)
我花了一段时间才了解自己在做错什么,但我最终意识到问题是我的代码假定总体价值最低的人也来自econ
部门,事实并非如此(这就是答案)这R
给了我)。实际上,总值最低的人来自stat
部门。
i <- which(df$value == min(df$value))
df[i,]
name dept value
9 i stat 2.461368
当然,我可以通过以下方式轻松找到问题的答案:
df_econ <- df2[df2$dept == "econ",]
df_econ
name dept value
1 a econ 54.40660
6 f econ 90.22990
15 o econ 14.28000
17 q econ 41.37243
18 r econ 36.88455
19 s econ 15.24447
df_econ[df_econ$value == min(df_econ$value),]
name dept value
15 o econ 14.28
但是我想知道是否可以通过与[
操作员一起使用“嵌套”子集来获得相同的结果。我的意思是用这样的代码:
df2[df2$dept == "econ",][... ,]
我目前尚不知道如何引用该value
列,因为第一个子设置操作的结果数据帧df2[df2$dept == "econ",]
是不同于的数据帧df2
。我也知道该value
列是第三列,但是我不知道如何使用列索引而不是它们的名称来设置子集条件。
感谢您的帮助。
以下是一些选项:
library(dplyr)
# also in @bramtayl's answer:
df2 %>% filter(dept == "econ") %>% filter(value==min(value))
# or
df2 %>% filter(dept == "econ") %>% slice(which.min(value))
# or...
library(data.table)
setDT(df2)[dept == "econ"][value==min(value)]
# or
setDT(df2)[dept == "econ"][which.min(value)]
这些软件包提供了方便的链接方式,这些方式在笨拙的R中不可用,例如笨拙的
subset(subset(df2, dept=="econ"), value == min(value))
可能还有其他软件包,但是这两个软件包最近被广泛使用。
评论。如果您只是浏览数据,建议您在以下dept
级别进行汇总:
# dplyr:
df2 %>% group_by(dept) %>% slice(which.min(value))
# data.table:
df2[, .SD[which.min(value)], by=dept]
dept name value
1: econ o 14.280002
2: math t 13.880606
3: stat i 2.461368
本文收集自互联网,转载请注明来源。
如有侵权,请联系 [email protected] 删除。
我来说两句