プロットの二項式信頼区間をプロットする

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) 

を試してみてください（これと同時に、ポイントとエラーバーをプロットするために、あなたの2プロットの文を兼ね備えています。）

代わりにあなたはプロポーションとそのCIをプロットすることができます：

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

Answer 2

ggplot2を使用するオプションを考慮しましたか？
geom_smoothは、線形モデル（「lm」）からの予測の95％信頼水準間隔を与えます。

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

出力は次のようになります。