I'm learning the dataRetrieval package. Following the example in Section
1.1.2 of the vignette (whatNWISdata) I prepared this script:
------
library("dataRetrieval")
siteNumbers <-
c("14207920","14208000","14208200","14208300","14208500","14208600","14208700","14208850","14209000","14209100","14209250","14209500","14209600","14209670","14209700","14209710","14209750","14209775","14209790","14209900","14210000","14210005","14210020","14210025","14210030","14210100","14210150","14210152","14210160","14210200","14210255","14210400","14210480","14210490","14210500","14210530","14210535","14210600","14210650","14210676","14210750","14210760","14210765","14210800","14210830","14210850","14210900","14211000","14211004","14211005","14211006","14211008","14211010","14211023","14211494")
dataAvailable <- whatNWISdata(siteNumbers, service="all",
parameterCD="all", statCD="all")
-----
The vignette says that for service, parameterCD, and statCD the default is
"all", but the package wants that explicitly. So that's what I
did. Yet, my
syntax is still off:
> source("R-scripts/get-site-data-list.R")
Error: All components of query must be named
Please show me what I'm missing.
Rich
You failed to name the first parameter, siteNumbers? Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Sat, Sep 5, 2020 at 4:44 PM Rich Shepard <rshepard at appl-ecosys.com> wrote:> I'm learning the dataRetrieval package. Following the example in Section > 1.1.2 of the vignette (whatNWISdata) I prepared this script: > ------ > library("dataRetrieval") > > siteNumbers <- > c("14207920","14208000","14208200","14208300","14208500","14208600","14208700","14208850","14209000","14209100","14209250","14209500","14209600","14209670","14209700","14209710","14209750","14209775","14209790","14209900","14210000","14210005","14210020","14210025","14210030","14210100","14210150","14210152","14210160","14210200","14210255","14210400","14210480","14210490","14210500","14210530","14210535","14210600","14210650","14210676","14210750","14210760","14210765","14210800","14210830","14210850","14210900","14211000","14211004","14211005","14211006","14211008","14211010","14211023","14211494") > > dataAvailable <- whatNWISdata(siteNumbers, service="all", > parameterCD="all", statCD="all") > ----- > > The vignette says that for service, parameterCD, and statCD the default is > "all", but the package wants that explicitly. So that's what I did. Yet, my > syntax is still off: > > > source("R-scripts/get-site-data-list.R") > Error: All components of query must be named > > Please show me what I'm missing. > > Rich > > ______________________________________________ > R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. >[[alternative HTML version deleted]]
This worked: > dataAvailable <- whatNWISdata(siteNumber=siteNumbers) > str(dataAvailable) 'data.frame':??? 2565 obs. of? 24 variables: ?$ agency_cd???????? : chr? "USGS" "USGS" "USGS" "USGS" ... ?$ site_no?????????? : chr? "14207920" "14208000" "14208000" "14208000" ... ?$ station_nm??????? : chr? "POOP CREEK NEAR BIG BOTTOM, OR" "CLACKAMAS RIVER AT BIG BOTTOM, OREG." "CLACKAMAS RIVER AT BIG BOTTOM, OREG." "CLACKAMAS RIVER AT BIG BOTTOM, OREG." ... ?$ site_tp_cd??????? : chr? "ST" "ST" "ST" "ST" ... ?$ dec_lat_va??????? : num? 45 45 45 45 44.9 ... ?$ dec_long_va?????? : num? -122 -122 -122 -122 -122 ... ?$ coord_acy_cd????? : chr? "U" "U" "U" "U" ... ?$ dec_coord_datum_cd: chr? "NAD83" "NAD83" "NAD83" "NAD83" ... ?$ alt_va??????????? : chr? " 2760.00" " 2040.00" " 2040.00" " 2040.00" ... ?$ alt_acy_va??????? : chr? " 20" " 20" " 20" " 20" ... ?$ alt_datum_cd????? : chr? "NGVD29" "NGVD29" "NGVD29" "NGVD29" ... ?$ huc_cd??????????? : chr? "17090011" "17090011" "17090011" "17090011" ... ?$ data_type_cd????? : chr? "pk" "dv" "pk" "sv" ... ?$ parm_cd?????????? : chr? NA "00060" NA NA ... ?$ stat_cd?????????? : chr? NA "00003" NA NA ... ?$ ts_id???????????? : num? 0 114288 0 0 0 ... ?$ loc_web_ds??????? : chr? NA NA NA NA ... ?$ medium_grp_cd???? : chr? "wat" "wat" "wat" "wat" ... ?$ parm_grp_cd?????? : chr? NA NA NA NA ... ?$ srs_id??????????? : num? 0 1645423 0 0 0 ... ?$ access_cd???????? : num? 0 0 0 0 0 0 0 0 0 0 ... ?$ begin_date??????? : Date, format: "1966-05-05" "1920-04-01" "1921-01-03" ... ?$ end_date????????? : Date, format: "1983-01-06" "1970-09-29" "1970-01-23" ... ?$ count_nu????????? : num? 17 18444 50 26 10 ... ?- attr(*, "comment")= chr? "#" "#" "# US Geological Survey" "# retrieved: 2020-09-05 20:18:46 -04:00\t(caas01)" ... ?- attr(*, "queryTime")= POSIXct, format: "2020-09-05 17:18:46" ?- attr(*, "url")= chr "https://waterservices.usgs.gov/nwis/site/?seriesCatalogOutput=true&sites=14207920,14208000,14208200,14208300,14"| __truncated__ ?- attr(*, "header")=List of 12 ? ..$ date?????????????????????? : chr "Sun, 06 Sep 2020 00:18:45 GMT" ? ..$ server???????????????????? : chr "Apache-Coyote/1.1" ? ..$ strict-transport-security? : chr "max-age=31536000" ? ..$ vary?????????????????????? : chr "Accept-Encoding" ? ..$ content-encoding?????????? : chr "gzip" ? ..$ content-type?????????????? : chr "text/plain;charset=UTF-8" ? ..$ cache-control????????????? : chr "max-age=900" ? ..$ expires??????????????????? : chr "Sun, 06 Sep 2020 00:33:46 GMT" ? ..$ x-ua-compatible??????????? : chr "IE=edge,chrome=1" ? ..$ access-control-allow-origin: chr "*" ? ..$ x-frame-options??????????? : chr "deny" ? ..$ transfer-encoding????????? : chr "chunked" ? ..- attr(*, "class")= chr? "insensitive" "list" -- David. On 9/5/20 4:44 PM, Rich Shepard wrote:> I'm learning the dataRetrieval package. Following the example in Section > 1.1.2 of the vignette (whatNWISdata) I prepared this script: > ------ > library("dataRetrieval") > > siteNumbers <- > c("14207920","14208000","14208200","14208300","14208500","14208600","14208700","14208850","14209000","14209100","14209250","14209500","14209600","14209670","14209700","14209710","14209750","14209775","14209790","14209900","14210000","14210005","14210020","14210025","14210030","14210100","14210150","14210152","14210160","14210200","14210255","14210400","14210480","14210490","14210500","14210530","14210535","14210600","14210650","14210676","14210750","14210760","14210765","14210800","14210830","14210850","14210900","14211000","14211004","14211005","14211006","14211008","14211010","14211023","14211494") > > dataAvailable <- whatNWISdata(siteNumbers, service="all", > parameterCD="all", statCD="all") > ----- > > The vignette says that for service, parameterCD, and statCD the > default is > "all", but the package wants that explicitly. So that's what I did. > Yet, my > syntax is still off: > >> source("R-scripts/get-site-data-list.R") > Error: All components of query must be named > > Please show me what I'm missing. > > Rich > > ______________________________________________ > R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code.
Name all your arguments (the vignette gets this wrong), including siteNumber=siteNumbers. Also, the vignette uses 'statCd', not 'statCD'. dataAvailable <- whatNWISdata(siteNumber=siteNumbers, service="all", statCd="all") gives some results. -Bill On Sat, Sep 5, 2020 at 4:44 PM Rich Shepard <rshepard at appl-ecosys.com> wrote:> I'm learning the dataRetrieval package. Following the example in Section > 1.1.2 of the vignette (whatNWISdata) I prepared this script: > ------ > library("dataRetrieval") > > siteNumbers <- > c("14207920","14208000","14208200","14208300","14208500","14208600","14208700","14208850","14209000","14209100","14209250","14209500","14209600","14209670","14209700","14209710","14209750","14209775","14209790","14209900","14210000","14210005","14210020","14210025","14210030","14210100","14210150","14210152","14210160","14210200","14210255","14210400","14210480","14210490","14210500","14210530","14210535","14210600","14210650","14210676","14210750","14210760","14210765","14210800","14210830","14210850","14210900","14211000","14211004","14211005","14211006","14211008","14211010","14211023","14211494") > > dataAvailable <- whatNWISdata(siteNumbers, service="all", > parameterCD="all", statCD="all") > ----- > > The vignette says that for service, parameterCD, and statCD the default is > "all", but the package wants that explicitly. So that's what I did. Yet, my > syntax is still off: > > > source("R-scripts/get-site-data-list.R") > Error: All components of query must be named > > Please show me what I'm missing. > > Rich > > ______________________________________________ > R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. >[[alternative HTML version deleted]]
On Sat, 5 Sep 2020, Bert Gunter wrote:> You failed to name the first parameter, siteNumbers?Bert,>> siteNumbers <- >> c("14207920","14208000","14208200","14208300","14208500","14208600","14208700","14208850","14209000","14209100","14209250","14209500","14209600","14209670","14209700","14209710","14209750","14209775","14209790","14209900","14210000","14210005","14210020","14210025","14210030","14210100","14210150","14210152","14210160","14210200","14210255","14210400","14210480","14210490","14210500","14210530","14210535","14210600","14210650","14210676","14210750","14210760","14210765","14210800","14210830","14210850","14210900","14211000","14211004","14211005","14211006","14211008","14211010","14211023","14211494")This is where I defined that vector; the MUA split it into two lines. The vector's name is the same as in the vignette example>> dataAvailable <- whatNWISdata(siteNumbers, service="all", parameterCD="all", statCD="all")Have I not specified the list of site numbers properly here? Thanks, Rich
On Sat, 5 Sep 2020, David Winsemius wrote:> This worked: >> dataAvailable <- whatNWISdata(siteNumber=siteNumbers)David/Bert, Mea culpa! The vignette example used a single siteNumber which is also the variable name and I missed that last point. Thank you very much. Regards, Rich
On Sat, 5 Sep 2020, Bill Dunlap wrote:> Name all your arguments (the vignette gets this wrong), including > siteNumber=siteNumbers. > Also, the vignette uses 'statCd', not 'statCD'.Bill, Thanks very much. I missed these.> dataAvailable <- whatNWISdata(siteNumber=siteNumbers, service="all", > statCd="all") > gives some results.Much appreciated, Rich