Биномиальные доверительные интервалы с графикой

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) 

(Это объединяет в себе две полученных заявления для построения точек и засовов ошибок одновременно.)

В качестве альтернативы можно построить пропорции и их КУ:

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") 

Выход будет: