如何绘制占用概率?

How to plot occupation probability?

我目前遇到以下问题。我收集了蝙蝠箱的数据并检查了箱子里是否有蝙蝠(0=不存在,1=存在)。现在我想在 y 轴上绘制占用概率,在 x 轴上绘制蝙蝠箱周围 500 米半径范围内的树木覆盖率。

我找不到任何提示,也不知道如何操作。

这是数据的一部分(所以你看:随着树木覆盖率的增加,蝙蝠箱被蝙蝠使用,这就是我想要绘制的(我想也许 geom_smooth 会很好)

Tree Cover (%)   Occupation

11,25806         0 

9,84961          0    

9,77807          1

14,89228         0

14,38696         0

14,40911         0

14,24987         0

18,43799         1

20,28749         0

20,13795         1

19,87961         0

22,14027         1

46,17592         1

46,82814         1 

28,36949         0

28,40777         1

我感谢每一个想法!

这似乎是 glm 的经典问题。您可以直接在 ggplot:

中绘制 glm 结果
ggplot(d, aes(Tree_Cover, Occupation)) +
  geom_smooth(method = 'glm', method.args = list(family = 'binomial')) +
  geom_point(shape = '|', size = 8, alpha = 0.3)

y 轴表示蝙蝠箱被占用的概率,具体取决于树木覆盖率。垂直条对应于原始观察结果。