用plotrix

Question

我需要在我的情节中添加二项式置信区间。

以下是我的步骤：

library(binom)
library(plotrix)
x <- c(1:6)
y <- c(68, 69, 70, 75, 75, 87)

CI <- binom.confint(y, 265, conf.level = 0.95, methods = "exact")
plot(x, y)
plotCI(x, y, ui = CI$upper, li = CI$lower, add = TRUE)

我认为我做的一切都正确，但我的输出图似乎不正确：

你有什么建议吗？

Answer 1

binom.confint会返回比例的置信区间，而不是总数（如果您通过打印检查了CI对象，则可能会注意到这一点） 。尝试

plotCI(x,y,ui=CI$upper*CI$n,li=CI$lower*CI$n)

（这将您的两个绘图语句组合在一起绘制点和误差线。）

或者，您可以绘制比例及其CI：

plotCI(x,y/CI$n,ui=CI$upper,li=CI$lower)

Answer 2

您是否考虑过使用ggplot2的选项？
geom_smooth为您提供线性模型预测的95％置信水平区间（“lm”）。

data<-data.frame(y=c(20.7, 18, 21.4, 15.3, 27.3, 20),x=c(1:6))
library(ggplot2)
g<-ggplot(data,aes(x,y))
g+geom_point()+geom_smooth(method="lm")

输出将是：