如何找到和plot不同的可信區間？

Question

如果我使用“發現”中的數據，我如何在置信曲線中找到例如 80%、90% 和 95% 的可信區間和 plot 它們？

通過使用summary(discoveries) ，我的意思是 a = 3.1，樣本量 n = 310 和標准偏差 s =2.5

我可以執行以下操作來找到 90% 嗎？ 如果是這樣，我怎么能 plot 我從幾個％的數據得到？ （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

您可以像這樣創建一個 function：

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
     }