プロットに二項信頼区間を追加する必要があります。
私の手順は次のとおりです。
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)
私はすべてを正しく行ったと思いますが、出力プロットは正しくないようです:
何か提案はありますか？
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)
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")
出力は次のようになります。
