8

私はRを初めて使用し、2010年から2020年までのアラバマ郡の予測人口増加をプロットするプロジェクトをいじっていました。これが私のコードです:

dat <- read.table("C:/Users/rasmus/Documents/countyinfo.txt", sep="\t", header=TRUE)

library(ggplot2)
library(maps)
library(ggmap)

mdat <- map_data('county')
str(mdat)
al1 = get_map(location = c(lon = -86.304474, lat = 32.362563), 
              zoom = 7, maptype = 'roadmap')
al1MAP = ggmap(al1) + 
         geom_point(data=dat,inherit.aes = FALSE, 
                    aes(y=Lat, x=Lon, map_id=County, size=Growth), col="red") + 
         borders("state", colour="red", alpha=0.8) + 
         borders("county", colour="blue", alpha=0.5)
al1MAP

さて、私は2つの質問があります。

1) 州境は奇妙なことをしているようだ. 郡のオーバーレイがオフになっているスクリーンショットを次に示します。

ここに画像の説明を入力

2) これはアラバマ州だけの話なので、州境以外の国境をオフにしたいのですが、その方法がわかりません。xlimとを試してみると思いylimますが、それらをアラバマ国境ポリゴンに制限する方法がわかりません。

4

1 に答える 1

12

borders()いくつかの状態の座標の関数で接続されているようです。

この問題を解決するには、 を使用しstateて境界線を別のデータ フレームとして保存し、を使用して境界線をマップにmap_data()追加します。さらに、ポイントが 1 つのリージョンのみに従って接続されるようにします。stategeom_path()geom_path()group=region

アラバマ郡だけの境界線を表示するには、region="alabama"function に引数を追加する必要がありますborders()。

al1 = get_map(location = c(lon = -86.304474, lat = 32.362563), 
      zoom = 6, maptype = 'roadmap')
mdat <- map_data('state')
ggmap(al1) + 
  geom_path(data=mdat,aes(x=long,y=lat,group=region),colour="red",alpha=0.8)+
  borders("county", colour="blue", alpha=0.5,region="alabama")

ここに画像の説明を入力

于 2013-03-28T08:41:05.447 に答える