如何使用 ggplot2 中使用的 mean_CI_boot 计算自举置信区间？

Question

我有一个 2 x 2 阶乘数据集，我使用mean_cl_boot函数绘制了置信区间。我想使用适当的函数在 R 中计算它。我怎样才能做到这一点？

我的数据集示例如下：

df <- data.frame(
      fertilizer = c("N","N","N","N","N","N","N","N","N","N","N","N","P","P","P","P","P","P","P","P","P","P","P","P","N","N","N","N","N","N","N","N","N","N","N","N","P","P","P","P","P","P","P","P","P","P","P","P"), 
      level = c("low","low","high","high","low","low","high","high","low","low","high","high","low","low","high","high","low","low","high","high","low","low","high","high","low","low","high","high","low","low","high","high","low","low","high","high","low","low","high","high","low","low","high","high","low","low","high","low"), 
      repro = c(0,90,2,4,0,80,1,90,2,33,56,0,99,100,66,80,1,0,2,33,0,0,1,2,90,5,2,2,5,8,0,1,90,2,4,66,0,0,0,0,1,2,90,5,2,5,8,55)
    )    

我知道有从图表中提取CI点的方式，但我并不想这样做。我想使用计算这个的函数。