We lost variableInfo on readNWISqw?
pCode <- c("00665","80154") #Phos and Sed
sitesVA <- readNWISdata(stateCd="VA", parameterCd=pCode,
service="site", seriesCatalogOutput=TRUE)
sitesVA.1 <- filter(sitesVA, parm_cd %in% pCode) %>%
filter(data_type_cd == "qw") %>%
filter(count_nu > 500) %>%
filter(duplicated(site_no)) %>%
arrange(desc(count_nu))
comboData <- readNWISqw(sitesVA.1$site_no[1],
pCode, reshape = TRUE) %>%
filter(!is.na(result_va_00665) & !is.na(result_va_80154))
variableInfo <- attr(comboData, "variableInfo")
variableInfo
NULL