[英]R: how to apply a function to some columns into groups of data?
I appreciate your help to apply a function to a data frame grouping by some columns. 感谢您为将功能应用于按某些列分组的数据框所做的帮助。 I suppose I have to use some dplyr
function or lapply
or do.call
but I was not able to do that. 我想我必须使用某些dplyr
函数或lapply
或do.call
但我无法做到这一点。
I have the following data frame: 我有以下数据框:
dfFull <- data.frame(Cen = c("Cen01", "Cen01", "Cen01", "Cen01",
"Cen01", "Cen01", "Cen01", "Cen01",
"Cen02", "Cen02", "Cen02", "Cen02",
"Cen02", "Cen02", "Cen02", "Cen02"),
Model = c("Mod01", "Mod01", "Mod01", "Mod01",
"Mod02", "Mod02", "Mod02", "Mod02",
"Mod01", "Mod01", "Mod01", "Mod01",
"Mod02", "Mod02", "Mod02", "Mod02"),
Indiv = c(1,2,3,4,1,2,3,4,1,2,3,4,1,2,3,4),
PF = c(1,1,2,2,1,1,2,2,1,1,2,2,1,1,2,2),
Obj1 = c(0.0,0.02,0.01,0.03,0.01,
0.0,0.02,0.0,0.15,0.03,
0.02,0.08,0.1,0.06,0.02,0.09),
Obj2 = c(0.8,0.62,0.85,0.7,0.92,
0.26,0.85,0.93,0.03,0.84,
0.94,0.84,0.05,0.63,0.83,0.92))
I have to call a function (from emoa package): 我必须调用一个函数(从emoa包中):
dominated_hypervolume(matrix_points, refp)
that calculates the hypervolume
for matrix_point
using a pre-defined refp
. dominated_hypervolume(matrix_points, refp)
其计算所述hypervolume
为matrix_point
使用预先定义的refp
。 refp
is a vector ( RP <- c(1.0,1.0)
) used for every calculation. refp
是用于每次计算的向量( RP <- c(1.0,1.0)
)。 The problem relies on matrix_points: 问题取决于matrix_points:
matrix_points
is a matrix that is transposed compared to my data frame. matrix_points
是与我的数据帧相比转置的矩阵。 hypervolume
calculated using Obj1 and Obj2 of all Indiv grouped by Cen , Mod and PF columns. 我需要使用按Cen , Mod和PF列分组的所有Indiv的Obj1和Obj2计算的hypervolume
体积 。 Using small data, I know dominated_hypervolume
will do the work since I would be able to provide the proper data. 使用小数据,我知道dominated_hypervolume
可以完成工作,因为我将能够提供适当的数据。
I know it is wrong, but I was trying to do something like: 我知道这是错误的,但是我正在尝试执行以下操作:
dfFull <- dfFull %>%
group_by(Cen, Model, PF) %>%
do.call(HV =dominated_hypervolume(data.matrix(t(dfFull[,5:6]), RP)))
What I expect at the end is the below. 我期望的最终结果如下。 HV
value is just a example, not the calculated. HV
值仅是示例,而不是计算得出的值。 It is not a problem to repeat the HV
value for the lines of the individuals used in its calculation. 对计算中使用的个人行重复HV
值不是问题。
Cen Model PF Indiv Obj1 Obj2 HV
Cen01 Mod01 1 1 0.0 0.8 0.77
Cen01 Mod01 1 2 0.02 0.62 0.77
Cen01 Mod01 2 3 0.01 0.85 0.74
Cen01 Mod01 2 4 0.03 0.70 0.74
Cen01 Mod02 1 1 0.01 0.92 0.81
Cen01 Mod02 1 2 0.0 0.26 0.81
Cen01 Mod02 2 3 0.02 0.85 0.69
Cen01 Mod02 2 4 0.0 0.93 0.69
Cen02 Mod01 1 1 0.15 0.03 0.88
Cen02 Mod01 1 2 0.03 0.84 0.88
Cen02 Mod01 2 3 0.02 0.94 0.86
Cen02 Mod01 2 4 0.08 0.84 0.86
Cen02 Mod02 1 1 0.1 0.05 0.76
Cen02 Mod02 1 2 0.06 0.63 0.76
Cen02 Mod02 2 3 0.02 0.83 0.64
Cen02 Mod02 2 4 0.09 0.92 0.64
Thanks for your help. 谢谢你的帮助。
library(tidyverse)
library(emoa)
RP <- c(1.0,1.0)
dfFull %>%
nest(-Cen, -Model, -PF) %>%
mutate(HV = map_dbl(data, ~dominated_hypervolume(t(data.frame(.x$Obj1, .x$Obj2)), RP))) %>%
unnest()
# Cen Model PF HV Indiv Obj1 Obj2
# 1 Cen01 Mod01 1 0.3764 1 0.00 0.80
# 2 Cen01 Mod01 1 0.3764 2 0.02 0.62
# 3 Cen01 Mod01 2 0.2940 3 0.01 0.85
# 4 Cen01 Mod01 2 0.2940 4 0.03 0.70
# 5 Cen01 Mod02 1 0.7400 1 0.01 0.92
# 6 Cen01 Mod02 1 0.7400 2 0.00 0.26
# 7 Cen01 Mod02 2 0.1484 3 0.02 0.85
# 8 Cen01 Mod02 2 0.1484 4 0.00 0.93
# 9 Cen02 Mod01 1 0.8437 1 0.15 0.03
# 10 Cen02 Mod01 1 0.8437 2 0.03 0.84
# 11 Cen02 Mod01 2 0.1508 3 0.02 0.94
# 12 Cen02 Mod01 2 0.1508 4 0.08 0.84
# 13 Cen02 Mod02 1 0.8698 1 0.10 0.05
# 14 Cen02 Mod02 1 0.8698 2 0.06 0.63
# 15 Cen02 Mod02 2 0.1666 3 0.02 0.83
# 16 Cen02 Mod02 2 0.1666 4 0.09 0.92
Consider by
(object-oriented wrapper of tapply
for grouping) to pass grouped subets to run needed defined method, then casted into data.frame()
: 考虑by
(的面向对象的包装tapply
用于分组)来传递分组subets运行所需定义的方法,然后浇铸成data.frame()
grpcols <- c("Cen", "Model", "PF")
df_list <- by(dfFull, dfFull[grpcols], function(sub)
data.frame(Cen = sub$Cen[[1]],
Model = sub$Model[[1]],
PF = sub$PF[[1]],
HV = dominated_hypervolume(t(sub[,5:6]), RP))
)
And then process the list of data frames: 然后处理数据帧列表:
# BASE PROCESSING
final_df <- do.call(rbind, df_list)
final_df <- with(final_df, final_df[order(Cen, Model, PF),])
row.names(final_df) <- NULL
final_df <- merge(dfFull, final_df[c(grpcols, "HV")], by=c("Cen", "Model", "PF"))
final_df
# Cen Model PF Indiv Obj1 Obj2 HV
# 1 Cen01 Mod01 1 1 0.00 0.80 0.3764
# 2 Cen01 Mod01 1 2 0.02 0.62 0.3764
# 3 Cen01 Mod01 2 3 0.01 0.85 0.2940
# 4 Cen01 Mod01 2 4 0.03 0.70 0.2940
# 5 Cen01 Mod02 1 1 0.01 0.92 0.7400
# 6 Cen01 Mod02 1 2 0.00 0.26 0.7400
# 7 Cen01 Mod02 2 3 0.02 0.85 0.1484
# 8 Cen01 Mod02 2 4 0.00 0.93 0.1484
# 9 Cen02 Mod01 1 1 0.15 0.03 0.8437
# 10 Cen02 Mod01 1 2 0.03 0.84 0.8437
# 11 Cen02 Mod01 2 3 0.02 0.94 0.1508
# 12 Cen02 Mod01 2 4 0.08 0.84 0.1508
# 13 Cen02 Mod02 1 1 0.10 0.05 0.8698
# 14 Cen02 Mod02 1 2 0.06 0.63 0.8698
# 15 Cen02 Mod02 2 3 0.02 0.83 0.1666
# 16 Cen02 Mod02 2 4 0.09 0.92 0.1666
Or with dplyr handling still using the base::by()
output where do.call
, order
, merge
becomes bind_rows
, arrange
, inner_join
respectively: 或与dplyr处理仍然采用base::by()
其中输出do.call
, order
, merge
成为bind_rows
, arrange
, inner_join
分别为:
# DPLYR PROCESSING
attributes(df_list) <- NULL
final_df <- bind_rows(df_list) %>%
arrange(Cen, Model, PF) %>%
inner_join(dfFull, by=c("Cen", "Model", "PF"))
final_df
# Cen Model PF HV Indiv Obj1 Obj2
# 1 Cen01 Mod01 1 0.3764 1 0.00 0.80
# 2 Cen01 Mod01 1 0.3764 2 0.02 0.62
# 3 Cen01 Mod01 2 0.2940 3 0.01 0.85
# 4 Cen01 Mod01 2 0.2940 4 0.03 0.70
# 5 Cen01 Mod02 1 0.7400 1 0.01 0.92
# 6 Cen01 Mod02 1 0.7400 2 0.00 0.26
# 7 Cen01 Mod02 2 0.1484 3 0.02 0.85
# 8 Cen01 Mod02 2 0.1484 4 0.00 0.93
# 9 Cen02 Mod01 1 0.8437 1 0.15 0.03
# 10 Cen02 Mod01 1 0.8437 2 0.03 0.84
# 11 Cen02 Mod01 2 0.1508 3 0.02 0.94
# 12 Cen02 Mod01 2 0.1508 4 0.08 0.84
# 13 Cen02 Mod02 1 0.8698 1 0.10 0.05
# 14 Cen02 Mod02 1 0.8698 2 0.06 0.63
# 15 Cen02 Mod02 2 0.1666 3 0.02 0.83
# 16 Cen02 Mod02 2 0.1666 4 0.09 0.92
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.