Displaying 1 result from an estimated 1 matches for "msacmsa".
Did you mean:
msacm32
2006 Oct 07
3
merge and polylist
...;NA"))
nc=read.shape("astae.shp", dbf.data=TRUE, verbose=TRUE)
mappolys=Map2poly(nc)
submap <- subset(mappolys, nc$att.data$NAME!="Honolulu, HI")
nc$att.data=subset(nc$att.data, nc$att.data$NAME!="Honolulu, HI")
nc$att.data[,1]=as.numeric(paste(nc$att.data$MSACMSA))
#attributes(nc$att.data)
nc$att.data=merge(nc$att.data, dat, by.x="AREA", by.y="cod", all.x=TRUE,
sort=FALSE)
#attributes(nc$att.data)
tmp=file("tmp")
write.polylistShape(submap, nc$att.data, "tmp")
#___________________________________________________...