r-如何找到并绘制不同的可信区间

Question

如果我使用"发现"的数据，我如何找到80%、90%和95%的可信区间，并将它们绘制在置信曲线中？

通过使用summary(discoveries)，我的平均值a=3.1，样本量n=310，标准偏差s=2.5

我能做下面的操作找到90%吗？如果是这样，我如何绘制我从几个%得到的数据？(80%、90%、95%等(

a <- 3.1
s <- 2.5
n <- 310
error <- qnorm(0.95)*s/sqrt(n)
left <- a-error
right <- a+error
left
right

Answer 1

您可以创建这样的函数：

mu_interval<-function(data, prob){ 
lower<-mean(data)+qnorm((1-prob)/2,length(data)-1)*sqrt(var(data)/length(data)) 
upper<-mean(data)+qnorm((1+prob)/2,length(data)-1)*sqrt(var(data)/length(data)) 
interval<-c(lower,upper)
interval
}