แบบฝึกหัดจำลองสถานการณ์ขนาดเล็กเพื่อแสดงให้เห็นว่าคำตอบของ @soakley ทำงานได้หรือไม่:
# Set the number of trials, M
M=10^6
# Set the true mean for each trial
mu=rep(0,M)
# Set the true standard deviation for each trial
sd=rep(1,M)
# Set counter to zero
count=0
for(i in 1:M){
# Control the random number generation so that the experiment is replicable
set.seed(i)
# Generate one draw of a normal random variable with a given mean and standard deviation
x=rnorm(n=1,mean=mu[i],sd=sd[i])
# Estimate the lower confidence bound for the population mean
lower=x-9.68*abs(x)
# Estimate the upper confidence bound for the population mean
upper=x+9.68*abs(x)
# If the true mean is within the confidence interval, count it in
if( (lower<mu[i]) && (mu[i]<upper) ) count=count+1
}
# Obtain the percentage of cases when the true mean is within the confidence interval
count_pct=count/M
# Print the result
print(count_pct)
[1] 1
ออกจากหนึ่งล้านการทดลองสุ่มช่วงความเชื่อมั่นรวมถึงความจริงหมายถึงหนึ่งล้านครั้ง, ที่อยู่, เสมอ ว่าไม่ควรเกิดขึ้นในช่วงความเชื่อมั่นกรณีเป็น95%ช่วงความเชื่อมั่น
ดังนั้นสูตรดูเหมือนจะไม่ทำงาน ... หรือฉันทำผิดรหัส?
แก้ไข:ผลลัพธ์เชิงประจักษ์ที่เหมือนกันเมื่อใช้ ;
อย่างไรก็ตามเป็นสำหรับ - ดังนั้นจึงค่อนข้างใกล้เคียงกับช่วงความมั่นใจ 95%(μ,σ)=(1000,1)
0.950097≈0.95(μ,σ)=(1000,1000)