stat_density2d：删除包含非有限值的行

Question

嗨〜我看到了一些关于此的讨论，但似乎没有多少答案。所以我在城市地图上绘制密度/热图。

df = data.frame(lon=rnorm(1000, mean=-87.62, sd=0.01), 
                lat=rnorm(1000, mean=41.88, sd=0.01))


map = get_googlemap('chicago', zoom=14, color='bw')

ggmap(map, extent='device') +
  stat_density2d(
    aes(x=lon, y=lat, fill=..level.., alpha =..level..),
    size = 2, bins = 4, 
    data = df,
    geom = "polygon"
  ) +
  scale_fill_gradient(low='green', high='red') +
  scale_alpha(guide=FALSE)

它绘制得很好，但显示：

Warning message:
Removed 96 rows containing non-finite values (stat_density2d).

使用真实的城市犯罪数据，它实际上报告删除了大多数行。所以它让我担心，情节可能并不反映真实的分布。

这是什么意思？我应该担心情节不能反映数据吗？谢谢！

Answer 1

我认为这意味着你的一些数据超出了地图的范围。

运行上面的代码后，这个数字与地图警告中排除的行数相匹配。

bb <- attr(map, "bb")

sum(df$lat < bb$ll.lat | 
    df$lat > bb$ur.lat | 
    df$lon < bb$ll.lon | 
    df$lon > bb$ur.lon)

Answer 2

这也可能是由于您的绘图中使用的变量有一些缺失值。检查is.na（dataframe $ variable）会有所帮助。