信頼区間９５%が真値を含む様子の図

データの９５%信頼区間が

場合によっては真値を含まないケースを図示するので

以前、どこかを参考に書いたはずの絵が見つからないので

ggplot2で描いてみた。

もっと良い方法があるのかもしれないけど

ggplot2の練習もかねて、これで一応、良しとする。

dataset <- NULL
setmean <- 5
for (i in 1:100){
sampledata <- rnorm(100, setmean, 1)
testresu <- t.test(sampledata)
check <- NULL
if (testresu$conf.int[2] < setmean | testresu$conf.int[1] > setmean){
check <- 1
} else{
check <- 2
}
dataset <- rbind(dataset, c(i, testresu$estimate, testresu$conf.int[1], testresu$conf.int[2], check))
}
colnames(dataset) <- c("No", "ave", "min", "max", "check")
dataset <- as.data.frame(dataset)

ggplot() +
geom_line(aes(x=c(0,101), y =c(setmean, setmean)), size=2, colour="#00FF00") +
geom_point(data = subset(dataset, dataset$check == 1), aes(x=No, y=ave),
colour = "#FF0000") +
geom_errorbar(data = subset(dataset, dataset$check == 1),
aes(x=No, ymin = min, ymax = max), width =.2,
colour = "#FF0000") +
geom_point(data = subset(dataset, dataset$check == 2), aes(x=No, y=ave)) +
geom_errorbar(data = subset(dataset, dataset$check == 2),
aes(x=No, ymin = min, ymax = max), width =.2,
colour = "#000000") +
xlab("") + ylab("") +
theme(plot.title = element_text(size=20),
legend.title = element_text(size=15),
legend.text = element_text(size=15),
axis.text.x=element_blank(),
axis.text.y=element_text(size=15),
axis.title.y = element_text(size=15),
panel.border= element_rect(colour="#000000", fill=NA, size=1.5),
panel.background = element_rect(fill="#FFFFFF"),
panel.grid.major = element_blank())