【发布时间】:2019-07-16 13:08:05
【问题描述】:
我想提取具有最小值的列和行的名称,按物种名称分组。例如,如果 DB 是下面的示例,则结果将是物种 Ag 的 V1(列名)和 ow(行名),物种 Au 的结果是 V1 和 pre,物种 Bri 的结果是 V1 和 post。
species <- rep(c("Ag", "Au", "Bri"), each=4)
period<- rep(c("ow", "pre", "fp", "post"), l=12)
V1 <- c(1,3,5,10,18,3,15,6,3,3,46,1)
V2 <- c(2,4,5,2,8,13,5,98,54,34,56,11)
df<-data.frame(species, period, V1, V2)
我首先尝试了列名,但没有成功:
AIC.min<- df %>%
group_by(species) %>%
mutate(min.model = names(df[2:3])[apply(df[2:3],1,function(x) which(x==min(x)))])
我还尝试使用过滤器和选择,以及聚合而不是应用以下类似问题,但没有成功。没有一个对我有用,在这里我收到一条错误消息“名称错误(models.AIC.all[2:8])[apply(models.AIC.all[2:8], 1, function(x) which(x == : 无效的下标类型'list'"
任何人都可以建议将是令人难以置信的棒!
【问题讨论】:
-
您是否意识到您正在创建一个
matrix,并且所有数字都转换为字符串? (如果没有,请使用data.frame代替cbind。) -
你需要
lapply(split(df[3:4], df$species), function(x) which(x == min(as.matrix(x)), arr.ind = TRUE)) -
请加AIC.min:
object 'AIC.min' not found. -
如果您只需要列名,请尝试使用
max.colstack(lapply(split(df[3:4], df$species), function(x) names(x)[max.col(-x, "first")]))