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)