請教壹個關於r語言的boot strap的編程，如果已知壹堆數據服從正態分布，但是分布參數未知？

mu = mean(x)

variance = sd(x)^2*(n-1)/n

創建數據集

mu_boot = c()

建立循環

for (i in 1:1000){

y = rnorm(n,mu,sqrt(variance))

mu_boot[i] = mean(y)

variance_boot[i] = sd(y)^2*(n-1)/n

}

求置信區間

quantile(mu_boot,c(0.025,0.975))

quantile(variance_boot,c(0.025,0.975))