图的置信区间

【问题标题】：Confidence interval for graph图的置信区间
【发布时间】：2020-02-04 18:30:41
【问题描述】：

我正在尝试使用我之前编写和使用但现在让我头疼的代码绘制图表。我认为这可能是因为包 Rmisc 和 dplyr 之间的交互，如果我在 dplyr 之后调用包 Rmisc，至少我会收到通知。有办法解决吗？

  library(Rmisc)
  library(dplyr)
  Desc <- Acc %>%
  group_by(Epoch, Probability, Session) %>%
  summarise(mean = mean(Acc), CI = CI(Acc, ci=0.95)[1]- CI(Acc, ci=0.95)[2])

预期：

Epoch Probability Session mean    CI 
 1         1          1     0.89    0.002
 2         1          2     0.87    0.001

...

相反，我得到了这个：

 mean    CI
 0.96     0.002

【问题讨论】：

标签： r group-by confidence-interval

【解决方案1】：

我不确定 Rmisc 的这个问题是否是因为 dplyr 或 Rmisc 的安装在没有你的 'sessionInfo()' 的情况下很难判断，但是下面的技巧应该适合你。请注意，CI 函数与 Rmisc 中的函数完全相同。只需在控制台中执行 Rmisc::CI 即可自行查看。如果您觉得这有帮助，请投票。

library(dplyr)
same_CI <- function (x, ci = 0.95) 
{
  a <- mean(x)
  s <- sd(x)
  n <- length(x)
  error <- qt(ci + (1 - ci)/2, df = n - 1) * s/sqrt(n)
  return(c(upper = a + error, mean = a, lower = a - error))
}
Desc <- Acc %>%
      group_by(Epoch, Probability, Session) %>%
      summarise(mean = mean(Acc), 
                up = same_CI(Acc)[1],
                low = same_CI(Acc)[3],
                CI = up - low)

更新

请注意，您是从置信区间上限中减去平均值，这可能是错误所在。还有一件事要确保您在 Acc 数据集中获得了正确的级别。

【讨论】：