#计算包含s1gma2=4的区间数 sum(UCL 4) #计算经验的覆盖率（置信水平） mean(UCL 4) Code replicate函数要重复执行的代码放在{}中，参数expr可以调用一个函数： FCode calCI <-function(n,alpha){ x <-rnorm(n,mean =0,sd 2) return((n-1)*var(x)/qchisq(alpha,df n-1)) ] UCL<-replicate(1000,expr=calCI(n=20,alpha=.05)) mean(UCL>4) Coda Previous Next First Last Back Forward 14#Oéù¹sigma^2=4´mÍ sum(UCL > 4) #Oé²CX«(ò&Y²) mean(UCL > 4) ↓Code replicateºÍá­Eâ1ìËò3{}•, ÎÍexprå±N^òáºÍ: ↑Code calCI <- function(n,alpha){ x <- rnorm(n, mean = 0, sd = 2) return((n-1) * var(x) / qchisq(alpha, df = n-1)) } UCL<-replicate(1000,expr=calCI(n=20,alpha=.05)) mean(UCL>4) ↓Code Previous Next First Last Back Forward 14