Skip to content

dataRetrieval query error

7 messages · Bert Gunter, David Winsemius, Bill Dunlap +1 more

Rich Shepard
# 
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:
Error: All components of query must be named

Please show me what I'm missing.

Rich
Bert Gunter
# 
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:

  
  


  


      

    

    

      
      

        


  

        

      David Winsemius
    

    

      
      #
    

  


  

      
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:

        

            

      
      
      
    

        

      
      
    

  
  


  


      

    

    

      
      

        


  

        

      Bill Dunlap
    

    

      
      #
    

  


  

      
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:

            

      
      
      
    

        

      
      
    

        

  
  


  


      

    

    

      
      

        


  

        

      Rich Shepard
    

    

      
      #
    

  


  

      
On Sat, 5 Sep 2020, Bert Gunter wrote:

        

            

      
      
      
    

        
Bert,

            

      
      
      
    

        
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

            

      
      
      
    

        
Have I not specified the list of site numbers properly here?

Thanks,

Rich

  
  


  


      

    

    

      
      

        


  

        

      Rich Shepard
    

    

      
      #
    

  


  

      
On Sat, 5 Sep 2020, David Winsemius wrote:

        

            

      
      
      
    

        
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

  
  


  


      

    

    

      
      

        


  

        

      Rich Shepard
    

    

      
      #
    

  


  

      
On Sat, 5 Sep 2020, Bill Dunlap wrote:

        

            

      
      
      
    

        
Bill,

Thanks very much. I missed these.

            

      
      
      
    

        
Much appreciated,

Rich