# Package ID: knb-lter-vcr.26.20 Cataloging System:https://pasta.lternet.edu.
# Data set title: Extratropical Storms (1885-1996) Count by  Month (USA).
# Data set creator: Dr. Bruce Hayden -  
# Metadata Provider:    - Virginia Coast Reserve Long-Term Ecological Research Project 
# Contact:    - Information Manager LTER Network Office  - tech-support@lternet.edu
# Contact: Dr. Bruce Hayden -    - bph@virginia.edu
# Contact:    - Information manager - Virginia Coast Reserve Long-Term Ecological Research Project   - jporter@lternet.edu
# Metadata Link: https://portal.lternet.edu/nis/metadataviewer?packageid=knb-lter-vcr.26.20
# Stylesheet v2.11 for metadata conversion into program: John H. Porter, Univ. Virginia, jporter@virginia.edu 

inUrl1  <- "https://pasta.lternet.edu/package/data/eml/knb-lter-vcr/26/20/0f736eb14b6c0d55349a7f62b2dc4388" 
infile1 <- tempfile()
try(download.file(inUrl1,infile1,method="curl"))
if (is.na(file.size(infile1))) download.file(inUrl1,infile1,method="auto")

                   
 dt1 <-read.csv(infile1,header=F 
          ,skip=22
            ,sep=","  
                ,quot='"' 
        , col.names=c(
                    "SEQ",     
                    "YEAR",     
                    "MONTH",     
                    "LONG",     
                    "LAT55",     
                    "LAT525",     
                    "LAT50",     
                    "LAT475",     
                    "LAT45",     
                    "LAT425",     
                    "LAT40",     
                    "LAT375",     
                    "LAT35",     
                    "LAT325",     
                    "LAT30",     
                    "LAT275",     
                    "LAT25"    ), check.names=TRUE)
               
unlink(infile1)
		    
# Fix any interval or ratio columns mistakenly read in as nominal and nominal columns read as numeric or dates read as strings
                
if (class(dt1$SEQ)=="factor") dt1$SEQ <-as.numeric(levels(dt1$SEQ))[as.integer(dt1$SEQ) ]               
if (class(dt1$SEQ)=="character") dt1$SEQ <-as.numeric(dt1$SEQ)
if (class(dt1$YEAR)=="factor") dt1$YEAR <-as.numeric(levels(dt1$YEAR))[as.integer(dt1$YEAR) ]               
if (class(dt1$YEAR)=="character") dt1$YEAR <-as.numeric(dt1$YEAR)
if (class(dt1$MONTH)=="factor") dt1$MONTH <-as.numeric(levels(dt1$MONTH))[as.integer(dt1$MONTH) ]               
if (class(dt1$MONTH)=="character") dt1$MONTH <-as.numeric(dt1$MONTH)
if (class(dt1$LONG)=="factor") dt1$LONG <-as.numeric(levels(dt1$LONG))[as.integer(dt1$LONG) ]               
if (class(dt1$LONG)=="character") dt1$LONG <-as.numeric(dt1$LONG)
if (class(dt1$LAT55)=="factor") dt1$LAT55 <-as.numeric(levels(dt1$LAT55))[as.integer(dt1$LAT55) ]               
if (class(dt1$LAT55)=="character") dt1$LAT55 <-as.numeric(dt1$LAT55)
if (class(dt1$LAT525)=="factor") dt1$LAT525 <-as.numeric(levels(dt1$LAT525))[as.integer(dt1$LAT525) ]               
if (class(dt1$LAT525)=="character") dt1$LAT525 <-as.numeric(dt1$LAT525)
if (class(dt1$LAT50)=="factor") dt1$LAT50 <-as.numeric(levels(dt1$LAT50))[as.integer(dt1$LAT50) ]               
if (class(dt1$LAT50)=="character") dt1$LAT50 <-as.numeric(dt1$LAT50)
if (class(dt1$LAT475)=="factor") dt1$LAT475 <-as.numeric(levels(dt1$LAT475))[as.integer(dt1$LAT475) ]               
if (class(dt1$LAT475)=="character") dt1$LAT475 <-as.numeric(dt1$LAT475)
if (class(dt1$LAT45)=="factor") dt1$LAT45 <-as.numeric(levels(dt1$LAT45))[as.integer(dt1$LAT45) ]               
if (class(dt1$LAT45)=="character") dt1$LAT45 <-as.numeric(dt1$LAT45)
if (class(dt1$LAT425)=="factor") dt1$LAT425 <-as.numeric(levels(dt1$LAT425))[as.integer(dt1$LAT425) ]               
if (class(dt1$LAT425)=="character") dt1$LAT425 <-as.numeric(dt1$LAT425)
if (class(dt1$LAT40)=="factor") dt1$LAT40 <-as.numeric(levels(dt1$LAT40))[as.integer(dt1$LAT40) ]               
if (class(dt1$LAT40)=="character") dt1$LAT40 <-as.numeric(dt1$LAT40)
if (class(dt1$LAT375)=="factor") dt1$LAT375 <-as.numeric(levels(dt1$LAT375))[as.integer(dt1$LAT375) ]               
if (class(dt1$LAT375)=="character") dt1$LAT375 <-as.numeric(dt1$LAT375)
if (class(dt1$LAT35)=="factor") dt1$LAT35 <-as.numeric(levels(dt1$LAT35))[as.integer(dt1$LAT35) ]               
if (class(dt1$LAT35)=="character") dt1$LAT35 <-as.numeric(dt1$LAT35)
if (class(dt1$LAT325)=="factor") dt1$LAT325 <-as.numeric(levels(dt1$LAT325))[as.integer(dt1$LAT325) ]               
if (class(dt1$LAT325)=="character") dt1$LAT325 <-as.numeric(dt1$LAT325)
if (class(dt1$LAT30)=="factor") dt1$LAT30 <-as.numeric(levels(dt1$LAT30))[as.integer(dt1$LAT30) ]               
if (class(dt1$LAT30)=="character") dt1$LAT30 <-as.numeric(dt1$LAT30)
if (class(dt1$LAT275)=="factor") dt1$LAT275 <-as.numeric(levels(dt1$LAT275))[as.integer(dt1$LAT275) ]               
if (class(dt1$LAT275)=="character") dt1$LAT275 <-as.numeric(dt1$LAT275)
if (class(dt1$LAT25)=="factor") dt1$LAT25 <-as.numeric(levels(dt1$LAT25))[as.integer(dt1$LAT25) ]               
if (class(dt1$LAT25)=="character") dt1$LAT25 <-as.numeric(dt1$LAT25)
                
# Convert Missing Values to NA for non-dates
                


# Here is the structure of the input data frame:
str(dt1)                            
attach(dt1)                            
# The analyses below are basic descriptions of the variables. After testing, they should be replaced.                 

summary(SEQ)
summary(YEAR)
summary(MONTH)
summary(LONG)
summary(LAT55)
summary(LAT525)
summary(LAT50)
summary(LAT475)
summary(LAT45)
summary(LAT425)
summary(LAT40)
summary(LAT375)
summary(LAT35)
summary(LAT325)
summary(LAT30)
summary(LAT275)
summary(LAT25) 
                # Get more details on character variables
                
detach(dt1)