basedir=File.dirname(__FILE__)
require basedir+"/../lib/statsample"
require 'statsample/srs'
require 'statsample/multiset'
require 'gnuplot'
require 'statsample/graph/svggraph.rb'
require 'gsl'
tests=1000
sample_size=100

a=[]
r = GSL::Rng.alloc(GSL::Rng::TAUS, 1)

1.upto(tests).each {|i|
    a.push(r.ugaussian())
}



pop=a.to_vector(:scale)

pop.svggraph_histogram(10,basedir+"/images/mean_pop.svg")





s=pop.standard_deviation_population
puts "Parameters:"
puts "Mean:"+pop.mean.to_s
puts "Skew:"+pop.skew.to_s
puts "Kurtosis:"+pop.kurtosis.to_s

puts "SD:"+s.to_s
puts "SE with replacement:"+Statsample::SRS.standard_error_ksd_wr(s, sample_size, pop.size).to_s
puts "SE without replacement:"+Statsample::SRS.standard_error_ksd_wor(s, sample_size,pop.size).to_s

sd_with=[]
sd_without=[]
monte_with=Statsample::Resample.repeat_and_save(tests) {
    sample= pop.sample_with_replacement(sample_size)
    sd_with.push(Statsample::SRS.standard_error_esd_wr(sample.sds,sample_size,pop.size))
    sample.mean
}




monte_without=Statsample::Resample.repeat_and_save(tests) {
    sample= pop.sample_without_replacement(sample_size)
    sd_without.push(Statsample::SRS.standard_error_esd_wor(sample.sds,sample_size,pop.size))
    sample.mean
}





v_sd_with=sd_with.to_vector(:scale)
v_sd_without=sd_without.to_vector(:scale)

v_with=monte_with.to_vector(:scale)
v_without=monte_without.to_vector(:scale)


File.open(basedir+"/images/mean_ndp.svg","w") {|fp|
    fp.write(v_with.svggraph_normalprobability_plot().burn)
}


puts "=============="
puts "Sample distribution - with Replacement"
puts "Mean:"+v_with.mean.to_s
puts "Sd:"+v_with.sds.to_s
puts "Sd (estimated):"+v_sd_with.mean.to_s
puts "Sample distribution - without Replacement"
puts "Mean:"+v_without.mean.to_s
puts "Sd:"+v_without.sds.to_s
puts "Sd (estimated):"+v_sd_without.mean.to_s

v_without.svggraph_histogram(10,basedir+"/images/mean_wo_hist.svg")