R/data_functions.R

In sathish-deevi/Roads-Package: Find the Modes and Assess the Modality of Complex and Mixture Distributions, Especially with Big Datasets

#####################################
#		Load in Data		#
#####################################

Load.data<-function(...){
	
	temp <- tempfile()
	download.file("http://www.fhwa.dot.gov/planning/processes/tools/nhpn/2015/nhpnv14-05shp.zip"
		,temp, mode="wb")

	unzip(temp, "NHPNLine.shp")
	unzip(temp, "NHPNLine.shx")
	unzip(temp, "NHPNLine.dbf")

	#unzip(temp, "NHPNLine.prj")
	#unzip(temp, "NHPNLine.shp")

	unlink(temp)

	roads=readShapeLines("NHPNLine.shp",
		delete_null_obj=TRUE)

	if(!dim(rds)==NULL){
		print("Data loaded successfully.")
	}
	else{
		stop("Data not loaded.")
		print("Suggest downloading data manually and using Read.data function.")
	}

	return(roads)
}


#####################################
#		Read in Data		#
#####################################

Read.data<-function(...){
	roads=readShapeLines("NHPNLine.shp",
		delete_null_obj=TRUE)

	if(!dim(rds)==NULL){
		print("Data loaded successfully.")
	return(roads)
}


#####################################
#		  Select States			 #
#####################################

states.rds<-function(states,...){
    
	if(!is.list(states)){
		stop("Argument states is not in list format.")


    if(!isTRUE(
        all(
         suppressWarnings(
           unlist(lapply(states,function (X){ !is.na(as.numeric(noquote(X)))} ))
            ))
            )) 
    stop("Please use numbers to identify your states.
            Consider using the ... function.")



    state.data<-vector("list",length(states))
    
    if(fortify=="TRUE"){
        if(!require("maptools",character.only = TRUE)){
            {install.packages("maptools",dep=TRUE)
             if(!require(x,character.only = TRUE)) 
                stop("Please install maptools.")
            }
        }
        state.data<-lapply(states,function(x){
            tmp.d<-roads[roads@data$STFIPS=="x",]
            f.tmp.d<-fortify(tmp.d)
            return(f.tmp.d)
            }
        )
    }
        
    if(fortify=="FALSE"){
        state.data<-lapply(states,function(x){
            roads[roads@data$STFIPS=="x",]
            }
        )
    }


    if(tmp.d%in%ls){rm(tmp.d)}
    if(f.tmp.d%in%ls){rm(f.tmp.d)}
    names(state.data)<-states
    return(state.data)    
}



#####################################
#		  Select Counties			 #
#####################################


county.rds<-function(counties,...){

	if(!is.list(counties)){
		stop("Argument states is not in list format.")

length(counties)
vector("list",length)
#WArds<-Wards[Wards@data$CTFIPS=="33",]


#Wards[Wards@data$CTFIPS=="33",]
#WArds2<-fortify(WArds)


if(!is.list(counties)){
		stop("Argument counties is not in list format.")

    if(!isTRUE(
        all(
         suppressWarnings(
           unlist(lapply(counties,function (X){ !is.na(as.numeric(noquote(X)))} ))
            ))
            )) 
    stop("Please use numbers to identify your counties.
            Consider using the ... function.")

        
    county.data<-vector("list",length(counties))
    
    if(fortify=="TRUE"){
        if(!require("maptools",character.only = TRUE)){
            {install.packages("maptools",dep=TRUE)
             if(!require(x,character.only = TRUE)) 
                stop("Please install maptools.")
            }
        }
        county.data<-lapply(counties,function(x){
            tmp.d<-roads[roads@data$CTFIPS=="x",]
            f.tmp.d<-fortify(tmp.d)
            return(f.tmp.d)
            }
        )
    }
        
    if(fortify=="FALSE"){
        county.data<-lapply(counties,function(x){
            roads[roads@data$CTFIPS=="x",]
            }
        )
    }


    if(tmp.d%in%ls){rm(tmp.d)}
    if(f.tmp.d%in%ls){rm(f.tmp.d)}
    names(county.data)<-counties
    return(county.data)    
}