0

私はいくつかの経度位置データを持っており、いくつかの異なる調査対象者のそれぞれについて、時間の経過に伴うその変化を示したいと考えています。また、その上に限界タイルのような密度プロットを作成したいのですが、この密度がマップ上のどこにあるかを示す必要があるため、その上に地理マップを重ねる必要があります。

私のデータは次のようになります。

SO <- data.frame(date = rep(seq(as.Date("2000/1/1"), by = "day", length.out = 365), 3),
             julian = rep(seq(1,365,1),3),          
             ind = c(rep(1,365), rep(2,365), rep(3,365)),
             longitude = c(rnorm(365, mean = 90, sd =5), rnorm(365, mean = 85, sd =2), rnorm(365, mean = 92, sd =3)))

これまでのところ、平滑化された線を当てはめた経度をプロットし、座標を反転して x 軸に経度を残すことができました。私のコードは現在次のようになっています。

ggplot(SO, aes(x = julian, y = longitude)) +
          geom_point(aes(color=factor(ind)), size = 0.1) +
          stat_smooth(aes(group= factor(ind)), se = FALSE)+
          stat_smooth(aes(color=factor(ind))) +
          coord_flip() + scale_x_reverse()

ただし、タイル密度をプロットしてマップをオーバーレイすることに行き詰まりました。結果は次のようになります。

望ましい結果

マップを密度プロットに重ねる方法を考え出すことができれば、それはすでに非常に役立ちます。どうもありがとうございました。

4

0 に答える 0