在 ggplot2 / R 中添加指数 geom_smooth

Question

我正在尝试使用 ggplot2 生成一些示例图形，我选择的示例之一是 生日问题，这里使用的是从 Revolution 计算演示中“借用”的代码在奥斯康。

birthday<-function(n){
    ntests<-1000
    pop<-1:365
    anydup<-function(i){
        any(duplicated(sample(pop,n,replace=TRUE)))
        }
    sum(sapply(seq(ntests), anydup))/ntests
    }

x<-data.frame(x=rep(1:100, each=5)) 
x<-ddply(x, .(x), function(df) {return(data.frame(x=df$x, prob=birthday(df$x)))})
birthdayplot<-ggplot(x, aes(x, prob))+
        geom_point()+geom_smooth()+
        theme_bw()+
        opts(title = "Probability that at least two people share a birthday in a random group")+
        labs(x="Size of Group", y="Probability")

这里我的图表是我所描述的指数图表，但 geom_smooth 不太适合数据。我尝试过黄土方法，但这并没有太大改变。谁能建议如何添加更好的平滑度？

谢谢

保罗。

Answer 1

平滑例程无法足够快地对 x 的低值的突然变化做出反应（并且它无法知道 prob 的值被限制为 0- 1 个范围）。由于变异性如此之低，一个快速的解决方案是减少在每个点进行平滑的值的跨度。查看该图中的红线：

birthdayplot + geom_smooth(span=0.1, colour="red")

Answer 2

问题在于概率遵循逻辑曲线。如果您更改生日函数以返回原始的成功和失败而不是概率，则可以拟合适当的平滑线。

birthday<-function(n){
  ntests<-1000
  pop<-1:365
  anydup<-function(i){
    any(duplicated(sample(pop,n,replace=TRUE)))
  }
  data.frame(Dups = sapply(seq(ntests), anydup) * 1, n = n)
}
x<-ddply(x, .(x),function(df) birthday(df$x))

现在，您必须添加点作为摘要，并指定逻辑回归作为平滑类型。

ggplot(x, aes(n, Dups)) +
  stat_summary(fun.y = mean, geom = "point") +
  stat_smooth(method = "glm", family = binomial)