2016-08-26 7 views
0

内のすべての緯度と経度が、私は今、私は予測列で「1」すべてのグループに表示一度MAP

latitude longitude Class prediction 
9.7   21.757  244732  1 
12.21  36.736  112206  0 
-15.966  126.844 133969  1 

のようなデータセットをしようとしていますし、彼らの緯度と経度を取りき、後で私が表示したいです1つの地図上のすべての点。

実際に私が書いたコードは予測列にそれぞれ '1'をとり、毎回緯度と経度をとり、毎回地図上に1点を表示します。しかし、予測が1の緯度と経度をすべて集めて、すべての点を1つのマップに表示したい。

library(ggplot2) 
library(ggmap) #install.packages("ggmap") 

#data set name testData1 
for (i in 1:100){ 

    if (testData1$prediction[i]==1) { 

    lat <- testData1$latitude[i] 
    lon <- testData1$longitude[i] 
    df <- as.data.frame(cbind(lon,lat)) 

    # getting the map 
    mapgilbert <- get_map(location = c(lon = mean(df$lon), lat = mean(df$lat)), zoom = 4, 
          maptype = "satellite", scale = 2) 

    # plotting the map with some points on it 
    ggmap(mapgilbert) + 
     geom_point(data = df, aes(x = lon, y = lat, fill = "red", alpha = 0.8), size = 5, shape = 21) + 
     guides(fill=FALSE, alpha=FALSE, size=FALSE) 

    } 

} 
+0

通常、私は5分で答えを得ます。それは1時間でした。 .wht hpndをアクティブなユーザーに!助けてください –

+0

あなたはSOの合計3つの質問(これを含む)をお願いしました。答えのスピードで正確な予測モデルを構築するのに十分なものではありません。 – hrbrmstr

+0

サンプルサイズが予測するには低すぎます! 。 。 –

答えて

2

私はあなたがものをオーバーコンプリートしていると思う。

ggmap(mapgilbert) + 
    geom_point(data = subset(df, prediction == 1), aes(x = lon, y = lat, fill = "red", alpha = 0.8), size = 5, shape = 21) + 
    guides(fill = FALSE, alpha = FALSE, size = FALSE) 
+0

私はdfを作成しました。ここで、V1 = lat、V2 = long、V3 =予測です。私はコードを変更しました。それは完全なデータのサブセットです。しかし、私は地図上のポイントを見ることができませんでした。私が見ているのは地図だけです。マップ上の座標無し+ geom_point(データ=部分集合(df、V3 == 2)、aes(x = V1、y = V2、塗りつぶし= "赤"、アルファ= 0.8)、サイズ= 9、形= 21)+ ガイド(塗りつぶし= FALSE、アルファ= FALSE、サイズ= FALSE) –

+0

表示したい座標の範囲と、表示されているマップの範囲外の可能性があります。 – beetroot

関連する問題