rm(list = ls()) 								#THIS LINE CLEARS THE WORKSPACE
setwd("")									#SET TO THE FOLDER THAT HAS YOUR SUPPLEMENTARY FILES IN IT
library(data.table)								#INSTALL IF YOU DON'T HAVE IT ALREADY
BBS <- rbindlist(lapply(list.files("",full.names=TRUE),read.csv))		#NEEDS THE FULL PATH TO THE BBS FOLDER
BBS$countrynum <- NULL
BBS$SpeciesTotal <- NULL
BBS$StopTotal <- NULL
BBS$rteno <- paste(BBS$statenum,formatC(BBS$Route, width=3, flag="0"),sep="")

BBS <- melt(BBS, id.vars = c("Year","statenum","Route","RPID","Aou","rteno"),variable.name = "interval",value.name="count")

routeinfo <- read.csv("weather.csv")
routeinfo <- routeinfo[c("Year","ObsN","rteno","RunType")]
routeinfo$firstyear <- c(1L,routeinfo$ObsN[-1]!=routeinfo$ObsN[-nrow(routeinfo)])
BBS <- merge(BBS,routeinfo,by=c("rteno","Year"))
route_data <- read.csv("routes.csv")
route_data <- route_data[c("rteno",...)]
BBS <- merge(BBS,route_data,by="rteno")
BBS <- BBS[BBS$RunType > 0,]
BBS <- BBS[BBS$RPID < 102,]