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)