sim.ci.stdmean2(.05, 20, 20, 1.5, 3, 4, .75, 5000)
# Should return (within sampling error):
# Coverage Lower Error Upper Error Ave CI Width Ave Est
# Unweighted Standardizer 0.9058 0.0610 0.0332 1.342560 0.7838679
# Group 1 Standardizer 0.9450 0.0322 0.0228 1.827583 0.7862640
Run the code above in your browser using DataLab