r语言 - 图形的置信区间

Question

我正在尝试使用以前编写和使用的代码绘制图形，但现在让我头疼。我认为这可能是因为软件包 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

Answer 1

我不确定 Rmisc 的这个问题是否是因为没有您的"sessionInfo(("很难分辨 dplyr 或 Rmisc 的安装已损坏，但以下技巧应该适合您。请注意，函数 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 数据集中获得正确的水平。