# Example using sim_CI():
source("http://www2.jura.uni-hamburg.de/instkrim/kriminologie/Mitarbeiter/Enzmann/Software/sim_CI.r")

windows(record=T)
sim_CI()
sim_CI(samp=100)
sim_CI(sd=5,samp=100)
sim_CI(mean=15,sd=5,samp=100,seed=4321)
sim_CI(mean=15,sd=5,t.dist=F,samp=100,seed=4321)

# ------------------------------------------------------------------------------
# Note: The following demonstration takes time because there are 40 simulations
#       with 100,000 samples (10 to 72 cases) each:

# a) t-distribution, two-sided, sigma=1:

1-sim_CI(n=10,samp=100000,plot=F)$CR
1-sim_CI(n=18,samp=100000,plot=F)$CR
1-sim_CI(n=36,samp=100000,plot=F)$CR
1-sim_CI(n=72,samp=100000,plot=F)$CR

# b) t-distribution, one-sided, sigma=1:

1-sim_CI(n=10,one.sided=T,samp=100000,plot=F)$CR
1-sim_CI(n=18,one.sided=T,samp=100000,plot=F)$CR
1-sim_CI(n=36,one.sided=T,samp=100000,plot=F)$CR
1-sim_CI(n=72,one.sided=T,samp=100000,plot=F)$CR

# c) normal-distribution, two-sided, sigma=1:

1-sim_CI(n=10,t.dist=F,samp=100000,plot=F)$CR
1-sim_CI(n=18,t.dist=F,samp=100000,plot=F)$CR
1-sim_CI(n=36,t.dist=F,samp=100000,plot=F)$CR
1-sim_CI(n=72,t.dist=F,samp=100000,plot=F)$CR

# d) normal-distribution, one-sided, sigma=1:

1-sim_CI(n=10,t.dist=F,one.sided=T,samp=100000,plot=F)$CR
1-sim_CI(n=18,t.dist=F,one.sided=T,samp=100000,plot=F)$CR
1-sim_CI(n=36,t.dist=F,one.sided=T,samp=100000,plot=F)$CR
1-sim_CI(n=72,t.dist=F,one.sided=T,samp=100000,plot=F)$CR

# e) t-distribution, two-sided, sigma=5:

1-sim_CI(n=10,sd=5,samp=100000,plot=F)$CR
1-sim_CI(n=18,sd=5,samp=100000,plot=F)$CR
1-sim_CI(n=36,sd=5,samp=100000,plot=F)$CR
1-sim_CI(n=72,sd=5,samp=100000,plot=F)$CR

# f) t-distribution, one-sided, sigma=5:

1-sim_CI(n=10,sd=5,one.sided=T,samp=100000,plot=F)$CR
1-sim_CI(n=18,sd=5,one.sided=T,samp=100000,plot=F)$CR
1-sim_CI(n=36,sd=5,one.sided=T,samp=100000,plot=F)$CR
1-sim_CI(n=72,sd=5,one.sided=T,samp=100000,plot=F)$CR

# g) normal-distribution, two-sided, sigma=5:

1-sim_CI(n=10,sd=5,t.dist=F,samp=100000,plot=F)$CR
1-sim_CI(n=18,sd=5,t.dist=F,samp=100000,plot=F)$CR
1-sim_CI(n=36,sd=5,t.dist=F,samp=100000,plot=F)$CR
1-sim_CI(n=72,sd=5,t.dist=F,samp=100000,plot=F)$CR

# h) normal-distribution, one-sided, sigma=5:

1-sim_CI(n=10,sd=5,t.dist=F,one.sided=T,samp=100000,plot=F)$CR
1-sim_CI(n=18,sd=5,t.dist=F,one.sided=T,samp=100000,plot=F)$CR
1-sim_CI(n=36,sd=5,t.dist=F,one.sided=T,samp=100000,plot=F)$CR
1-sim_CI(n=72,sd=5,t.dist=F,one.sided=T,samp=100000,plot=F)$CR