我正在处理一个包含足球数据的嵌套数据框,并为每个数据框拟合了一个模型。使用goalmodel::predict_result函数,我希望基于模型预测每场比赛的结果。predict_result函数需要三个参数:模型、主队的向量和客队的向量。我如何在嵌套数据框中引用列名?这是我的嵌套数据框的一个例子:
library(goalmodel)library(tidyverse)nested_df# A tibble: 4 x 3# Groups: League [4] League data model <chr> <list<df[,133]>> <list> 1 F1 [380 x 133] <goalmodl>2 E0 [380 x 133] <goalmodl>3 SP1 [380 x 133] <goalmodl>4 D1 [308 x 133] <goalmodl>
如果我想预测这个数据框的第一个元素F1的结果,我会写:
predict_result(nested_df$model[[1]], team1=nested_df$data[[1]]$HomeTeam, team2=nested_df$data[[1]]$AwayTeam, return_df = TRUE)
这会返回一个包含所需结果的数据框。我尝试使用purrr::map来迭代上述函数,代码如下:
map(nested_df,~predict_result(.x$model, team1=.x$data[[.]]$HomeTeam, team2=.x$data[[.]]$AwayTeam, return_df = TRUE))
但它不起作用,显示的错误是:
Error in .x$data : $ operator is invalid for atomic vectors
我非常感谢任何帮助和建议,提前谢谢。—–为了可重现性,这里是一个例子:
df <- tibble(League = c("F1","E0","SP1","D1"), HomeTeam = c("TeamA","TeamB","TeamC","TeamD"), AwayTeam = c("TeamE","TeamF","TeamG","TeamH"), FTHG = c(0,1,2,0), FTAG = c(0,1,0,2))nested_df <- df %>% group_by(League)%>% nest()%>% mutate(model = map(data,~goalmodel::goalmodel(goals1 = .x$FTHG, goals2 = .x$FTAG, team1 = .x$HomeTeam, team2 = .x$AwayTeam, ,rs=TRUE)))nested_df# A tibble: 4 x 3# Groups: League [4] League data model <chr> <list<df[,4]>> <list> 1 F1 [1 x 4] <goalmodl>2 E0 [1 x 4] <goalmodl>3 SP1 [1 x 4] <goalmodl>4 D1 [1 x 4] <goalmodl>
我希望使用goalmodel::predict_result函数来迭代nested_df中的所有四个数据及其各自的goalmodl对象。提前谢谢。
回答:
您应该使用map2
函数。
map2(.x = nested_df$model, .y = nested_df$data, .f = ~ predict_result(model = .x, team1 = .y$HomeTeam, team2 = .y$AwayTeam, return_df = TRUE))
您会得到一个数据框列表。
或者也可以这样做:
nested_df <- nested_df %>% mutate(pred = map2(.x = model, .y = data, .f = ~ predict_result(model = .x, team1 = .y$HomeTeam, team2 = .y$AwayTeam, return_df = TRUE)))
您会在nested_df
中得到一个列表列。使用unnest
您可以将其转换为数据框。