R中的置信区间问题

Question

是否可以使用function创建一个方程，该方程将产生正态分布并同时产生所述数据的95％置信区间？我知道我可以使用rnorm(n,mean,sd)生成随机正态分布，但是如何获得输出来告诉我置信区间？

我尝试了sample_CI<- function(n,j,k){list(g<-rnorm(n, mean=j, sd=k), confint(g, level=.95))}。

感谢所有帮助。

Answer 1

可能正是您要找的东西

sample_CI <- function(n,j,k){
  error <- qnorm(0.975)*k/sqrt(n)
  left <- j-error
  right <- j+error
  paste("[",round(left,4)," ; ",round(right,4),"]")
}
sample_CI(1000,2,4)