1

デバッグできないように見えるデータを使用しようとすると、いくつかのエラーが発生します。

Rスクリプトはこちらhttps://dl.dropbox.com/u/28231177/This%20Should%20Work.R
データはこちらhttps://dl.dropbox.com/u/28231177/my_data.csv

実行したときにエラーを含む最後の数行を次に示します。

pds <- fortify(sf_map)
# Using OBJECTID to define regions.
pds$OBJECTID <- as.integer(pds$OBJECTID)
# Error in `$<-.data.frame`(`*tmp*`, "OBJECTID", value = integer(0)) : 
# replacement has 0 rows, data has 16249


### Make the map

p1 <- ggplot(my_data, aes(map_id = zip))
p1 <- p1 + geom_map(aes(fill=vol, map_id = zip), map = pds)
p1 <- p1 + expand_limits(x = pds$lon, y = pds$lat) + coord_equal()
p1 + xlab("Basic Map with Default Elements")
# Error in unit(x, default.units) : 'x' and 'units' must have length > 0
4

1 に答える 1

2

設定すると、コードが機能します

pds <- fortify(sf_map, region = "ID")

そして線を取り除く

pds$OBJECTID <- as.integer(pds$OBJECTID)

また、@mdsumner が述べたようにlongnotを使用する必要があります。lon

これは少し異なる解決策です（geom_polygonnotを使用geom_map）

library(maptools)
library(ggplot2)

## 'fortify' needs 'library(gpclib)' locally available
gpclibPermit()

## Import the shapefile
sf_map <- readShapeSpatial("sfzipcodes", ID = "ID")

## Import the data
my_data <- read.csv("my_data.csv")
my_data <- unique(my_data)

## Merge the data
sf_df <- fortify(sf_map, region='ID')
sf_df <- merge(sf_df, my_data, by.x="id", by.y="zip", all=FALSE)
sf_df <- sf_df[order(sf_df$group, sf_df$order), ]

## Make the map
p2 <- ggplot(sf_df, aes(x=long, y=lat, group=group, fill=vol))
p2 <- p2 + geom_polygon() + coord_equal()
p2 <- p2 + xlab("Basic Map with Default Elements")
p2
于 2012-08-05T09:32:42.433 に答える