# R_example_consistency.txt
# Population is exponential with mean 2
# Generates a random samples with increasing samplesize and plots 
# a 95% confidence interval for each sample. 
# nci is number of confidence intervals generated
# this code has sample size start at 10 and increase by 10 for each CI

nci=200

x=c(1:nci)
ybar=c(1:nci)
s=c(1:nci)
size=c(1:nci)

samplesize=10;
for (i in x) {
  y=rexp(samplesize,0.5)
  ybar[i]=mean(y)
  s[i]=sd(y)
  size[i]=samplesize
  samplesize=samplesize+10
}

yl=ybar-1.96*s/sqrt(size)
yu=ybar+1.96*s/sqrt(size)
plot(x,ybar,ylim=c(1,3))
segments(x,yl,x,yu)
abline(h=2.0)