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

Question

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

通过使用summary(discoveries)i 意味着 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

标签： rstatisticsconfidence-intervalcredible-interval

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
     }