如何找到和绘制不同的可信区间？

问题描述

如果我使用来自“发现”的数据，我如何找到例如 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

解决方法

你可以创建一个这样的函数：

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
     }