1

我正在使用 R 包 rnoaa（连同它需要的其他包）来收集历史天气数据。我编写了这个嵌套循环来收集所有数据集，但是当我运行它时我不断收到错误。它似乎运行了第二次罚款

循环：

require('triebeard')
require('bindr')
require('colorspace')
require('mime')
require('curl')
require('openssl')
require('R6')
require('urltools')
require('httpcode')
require('stringr')
require('assertthat')
require('bindrcpp')
require('glue')
require('magrittr')
require('pkgconfig')
require('rlang')
require('Rcpp')
require('BH')
require('plogr')
require('purrr')
require('stringi')
require('tidyselect')
require('digest')
require('gtable')
require('plyr')
require('reshape2')
require('lazyeval')
require('RColorBrewer')
require('dichromat')
require('munsell')
require('labeling')
require('viridisLite')
require('data.table')
require('rjson')
require('httr')
require('crul')
require('lubridate')
require('dplyr')
require('tidyr')
require('ggplot2')
require('scales')
require('XML')
require('xml2')
require('jsonlite')
require('rappdirs')
require('gridExtra')
require('tibble')
require('isdparser')
require('geonames')
require('hoardr')
require('rnoaa')
install.package('ncdf4')

install.packages("devtools")
library(devtools)
install_github("rnoaa", "ropensci")
library(rnoaa)


list <- buoys(dataset='wlevel')
lid <- data.frame(list$id)
foo <- for(range in 1990:2017){
for(bid in lid){
bid_range <- buoy(dataset = 'wlevel', buoyid = bid, year = range)
bid.year.data <- data.frame(bid.year$data)
write.csv(bid.year.data, file='cwind/bid_range.csv')
} 
}

响应：

Using c1990.nc
Using 
Error: length(url) == 1 is not TRUE

它保存了第一个数据集，但它没有在文件名中应用 for，它只是将其命名为 bid_range.csv。

4

1 回答 1

1

这个错误信息表明在 1990 年没有给定站 id 的任何数据。因为您使用的是 for 循环，一旦出现错误，它就会停止。

这里我介绍使用tidyverse下载 NOAA 浮标数据。以下许多功能来自purrr包，它是tidyverse.

# Load packages
library(tidyverse)
library(rnoaa)

第 1 步：创建一个包含 id 和 year 的所有组合的“Grid”

expand函数 from可以创建不同值的tidyr组合。

data_list <- buoys(dataset = 'wlevel')

data_list2 <- data_list %>%
  select(id) %>%
  expand(id, year = 1990:2017)

第 2 步：创建一个在没有数据时不会中断的“安全”版本。也使这个函数适用于 map2 函数

因为我们将使用它的和参数来循环使用函数的map2所有组合id和year使用。我们将参数序列修改为 create 。我们还使用该函数来创建. 现在，当它遇到错误时，它将存储错误消息并移动到下一个而不是中断。map2.x.ybuoy_modifysafelysafebuoy_modify

# Modify the buoy function
buoy_modify <- function(buoyid, year, dataset, ...){
  buoy(dataset, buoyid = buoyid, year = year, ...)
}

# Creare a safe version of buoy_modify
buoy_safe <- safely(buoy_modify)

第 3 步：应用 buoy_safe 函数

wlevel_data <- map2(data_list2$id, data_list2$year, buoy_safe, dataset = "wlevel")

# Assign name for the element in the list based on id and year
names(wlevel_data) <- paste(data_list2$id, data_list2$year, sep = "_")

完成此步骤后，所有数据都下载到wlevel_data. 中的每个元素wlevel_data都有两个部分。$result下载成功则显示数据，否则显示NULL. $error显示NULL下载是否成功，否则显示错误信息。

第 4 步：访问数据

transpose可以将列表“由内而外”。所以现在wlevel_data2有两个元素：结果和错误。我们可以存储这两个并访问数据。

# Turn the list "inside out"
wlevel_data2 <- transpose(wlevel_data)

# Get the error message
wlevel_error <- wlevel_data2$error

# Get he result
wlevel_result <-  wlevel_data2$result

# Remove NULL element in wlevel_result
wlevel_result2 <- wlevel_result[!map_lgl(wlevel_result, is.null)]
于 2017-09-15T07:38:58.697 回答