您为使用 magrittr 管道和地图函数的第一步选择了一个具有挑战性的问题!我会尽力给你一个有用的答案,但我也建议你在练习时找到一些更容易使用的数据。了解管道 %>% 的好地方是 Hadley Wickham 书中的 "Pipes" 章节。关于iteration 的章节也很好地介绍了map_* 函数。一旦你对概念有了更坚定的理解,你就可以回到更复杂的问题。我认为 Hadley 对这些工具的解释比我以往任何时候都好,所以我不会在这里详细介绍它们,而是重点解释为什么您的代码不起作用,以及为什么我的代码起作用。
对代码的分析
Map 函数允许使用一些有用的快捷方式,您已经发现了其中一个 - 即,如果您将向量或列表作为函数参数传递,它们会自动转换为提取器函数。所以,你在正确的轨道上!
要记住的是,map 函数返回一个与输入向量具有相同长度和名称的向量。您的输入向量是 jsonData,它有 5 个名称为 [1] "copyright" "allPlays" "currentPlay" "scoringPlays" "playsByInning" 的元素。当您运行jsonData %>% map("playEvents") %>% map("hitData") 时,正在提取数据,但 R 仍返回一个包含五个元素且名称与原始向量相同的向量。如果你看一下下面的代码,你会发现你的代码确实是剥离了最上层,但长度保持不变,这不是很有帮助:
> unlist(map(jsonData, class))
copyright allPlays currentPlay scoringPlays playsByInning
"character" "data.frame" "list" "integer" "data.frame"
> unlist(map(jsonData %>% map("playEvents"), class))
copyright allPlays currentPlay scoringPlays playsByInning
"NULL" "list" "data.frame" "NULL" "NULL"
> unlist(map(jsonData %>% map("playEvents") %>% map("hitData"), class))
copyright allPlays currentPlay scoringPlays playsByInning
"NULL" "NULL" "data.frame" "NULL" "NULL"
最终的输出,以及您试图与上面对bind_rows 的调用相结合的内容是:
> jsonData %>% map("playEvents") %>% map("hitData")
$copyright
NULL
$allPlays
NULL
$currentPlay
launchSpeed launchAngle totalDistance trajectory hardness location coordinates.coordX coordinates.coordY
1 NA NA NA <NA> <NA> <NA> NA NA
2 81.3 61.92 187.5 popup medium 6 75.78 167.97
$scoringPlays
NULL
$playsByInning
NULL
显然这不是你想要的。经过一番修补,我想出了以下解决方案。
我自己的策略
图书馆:
library(jsonlite)
library(purrr)
library(dplyr)
library(readr)
library(stringr)
library(magrittr)
我使用稍微不同的方法来下载和解析 JSON,因为我需要查看结构。我会把它包括在内,以防你觉得它有用:
url <- paste0("http://statsapi-prod-alt-968618993.us-east-1.elb.amazonaws",
".com/api/v1/game/565711/playByPlay")
url %>% read_file() %>% prettify() %>% write_file("bball.json")
jsonData <- fromJSON("bball.json")
我首先提取并清理 hitData 数据帧。我知道它们都可以在playEvents 中找到,所以我可以使用$ 语法跳过几个步骤。第一次调用map 从列表playEvents 的每个元素中提取hitData。 hitData 数据框是嵌套的(它们包含其他数据框),因此第二次使用 jsonlite::flatten 调用 map 会使它们变平。函数 safely 确保 R 在遇到数据帧以外的内容时不会抛出错误(只有 46 个元素包含 hitData)。许多hitData 数据帧包含充满NAs 的行,因此对map 的第三次调用使用匿名函数(再次在safely 中)来摆脱这些。第四次调用map 然后从每个元素的result 变量中提取数据帧,该变量由safely 创建(以及我们不需要的error 变量):
hitdata_list <- jsonData$allPlays$playEvents %>%
map("hitData") %>%
map(safely(jsonlite::flatten)) %>%
map(safely(~.$result[complete.cases(.$result),])) %>%
map("result")
现在我有一个 hitData 数据框列表。正如我上面提到的,80 个条目中只有 46 个包含 hitData,所以我需要一种从 atBatIndex 获取相应值的方法。当hitdata_list 中的元素包含数据框时,我可以通过使用TRUE 生成逻辑向量来做到这一点,否则FALSE。我使用map_lgl 来返回一个逻辑向量而不是一个列表:
lgl_index <- map_lgl(hitdata_list, ~ !is.null(.))
atbatindex_vec <- jsonData$allPlays$atBatIndex[lgl_index]
然后我使用stringr 函数从URL 中获取game_pk。我不确定它是否适用于每个 URL,但在这种情况下它可以正常工作:
game_pk_vec <- str_match(url, "/(\\d+)/")[2] %>%
as.integer()
最后,我将atBatIndex 和game_pk 组合到一个tibble 中,然后使用bind_cols 将该tibble 与hitData 数据组合起来。 hitData 数据帧仍在列表中,所以我需要先将它们与 bind_rows 结合起来。 set_colnames 函数来自 magrittr 包,并按照它所说的做。我需要设置列名,因为在展平 hitData 数据框时创建了一些复合名称:
hitdata_df <- tibble(game_pk = game_pk_vec, atBatIndex = atbatindex_vec) %>%
bind_cols(bind_rows(hitdata_list)) %>%
set_colnames(str_extract(names(.), "\\w+$"))
我唯一没有做的是提取pitchNumber。调用 jsonData$allPlays$playEvents %>% map("pitchNumber") 返回序列 1 到 n 的列表,其中每个向量的长度 > 1。我假设您只需要每个序列中的最终数字,但我不确定所以我会省自己的力气。你可以做我用atBatIndex 做的事情来获取相关元素,然后提取你需要的东西。这是最终的数据框:
# A tibble: 46 x 10
game_pk atBatIndex launchSpeed launchAngle totalDistance trajectory hardness location coordX coordY
<chr> <int> <dbl> <dbl> <dbl> <chr> <chr> <chr> <dbl> <dbl>
1 565711 4 76.6 2.74 188. ground_ball medium 9 178. 145.
2 565711 5 101. 15.4 328. line_drive hard 8 145. 62.2
3 565711 6 103. 29.4 382. line_drive medium 9 237. 79.4
4 565711 8 109. 15.6 319. line_drive hard 9 181. 102.
5 565711 9 75.8 47.8 239. fly_ball medium 7 99.8 103.
6 565711 10 91.6 44.1 311. fly_ball medium 8 140. 69.3
7 565711 12 79.1 23.4 246. line_drive medium 7 52.3 126.
8 565711 13 67.3 -21.3 124. ground_ball medium 6 108. 156.
9 565711 14 89.9 -21.6 7.41 ground_ball medium 6 108. 152.
10 565711 15 110. 27.7 420. fly_ball medium 9 250. 69.0
# … with 36 more rows