#include <oxstd.h>
#include <oxdraw.h>
#include <oxprob.h>
#import  <per.h>
decl iter_num=10000;
decl samp_num=1000;
get_boot(decl data)
{
	decl bootstrap=zeros(samp_num,1);
	for(decl i=0;i<samp_num;i++)
		bootstrap[i][0]=data[(int(ranu(1,1)*samp_num))][0];		
	return bootstrap;
}
main()
{
	ranseed("GM");
	decl t_data=rant(samp_num,1,1);
	decl median_bs=zeros(iter_num,1);
	for(decl j=0;j<iter_num;j++)
	{
		decl boot=get_boot(t_data);
		median_bs[j][0] = quantilec(boot);
		percentage(j,iter_num);
	}
	println(" The median 95% interval ");
	println("( ",quantilec(median_bs,0.025)[0][0]," , ",quantilec(median_bs,0.975)[0][0]," )");	
	DrawDensity(0,median_bs',"bootstrap",TRUE, FALSE, FALSE);
	DrawDensity(1,median_bs',"bootstrap",TRUE, FALSE, FALSE);
	DrawDensity(1,t_data',"t_dis",TRUE, FALSE, FALSE);
	DrawDensity(2,t_data',"t_dis",TRUE, FALSE, FALSE);
	ShowDrawWindow();
}