Hi all, I see someone had earlier posted similar question to the list. But, my problem is that? the counties in the input data don't match those on the map. Codes below:
############
library(spdep)
library(maps)
data<-read.csv("C:/Users/mitra.devkota/Desktop/MT1.csv",header=TRUE)?
names(data)
data$merge <- paste(tolower(data$State),tolower(data$NAME), sep='.')
newmap <- map("county",c("montana"), fill=T, plot=F)
?commaReplace <- function(x) {
? ? return(gsub(",",".", x))
}
newmap$names <- sapply(newmap$names, commaReplace, simplify=TRUE)
names(newmap$names) <- NULL
length(newmap$names)?
# Merge the data with the geography
index <- match(newmap$names, data$merge)
# printing index gives
index
?[1] ?1 ?2 ?3 ?4 ?5 ?6 ?7 ?8 ?9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
[26] 26 27 29 28 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50
[51] 51 52 53 54 55 56 NA
#The NA is the problem now.
Any hand is appreciated.
Mike