3

R「マップ」パッケージを使用して、郡ごとに米国のヒートマップ/コロプレス/テーマプロットを作成しようとしています。残念ながら、いくつかの郡がデータベースから欠落しているようです。

library(maps)
data(county.fips)
which(county.fips[,1] == 35006)
which(county.fips[,2] == 'new mexico,cibola')

エントリが見つかりませんでした

> library(maps)
> data(county.fips)
> which(county.fips[,1] == 35006)
integer(0)
> which(county.fips[,2] == 'new mexico,cibola')
integer(0)

この fips コードが存在することは知っています。(http://en.wikipedia.org/wiki/Cibola_County,_New_Mexico http://en.wikipedia.org/wiki/List_of_counties_in_New_Mexico )

完全なプロットを作成する方法を知っている人はいますか? これで見つけたものはすべてここにリストされていました：http://grokbase.com/t/r/r-sig-geo/0964515wmd/missing-county-in-maps

ご協力いただきありがとうございます！！

4

2 に答える 2

1

ああ、問題はマップではなく、county.fips にあります。Cibola County はそのデータセットに表示されないため、データ {variable-to-plot, FIPS code} をリストに一致させようとすると、Cibola は常に欠落していました。

シボラ郡を表示するための回避策があります... county.fips データセットを含まない郡の明示的な行を書きます。

countynames <- match(county.fips[,"fips"], mydatavector$fips)
cibola_value <- mydatavector[which(mydatavector$fips == 35006)]
m <- map("county")
map("county", countynames, col=rainbow(100)[mydatavector$var_to_plot], fill=TRUE, add=TRUE)
map("county", "new mexico,cibola", col=rainbow(100)[cibola_value], fill=TRUE, add=TRUE)

注: map 関数では、米国内のすべての郡をプロットするのに問題があることもわかりましたが、呼び出しを 2 つのセットに分割するとうまくいきました。

于 2013-01-23T16:52:25.543 に答える
0

シボラ郡は「county」データベースにありません。地図を描くには十分なはずです。例を参照してください。

map.text('county', 'new mexico')
于 2012-12-18T22:09:58.750 に答える