***Jan Ditzen - November 2013
***	File creates 5 funnels, one for each SR
***	Attention, file is overwritten! 
/*
Changelog:
	17.02.2015 - added PEESE 
				 added seed to produce identical results
*/

clear all
set seed 345623

cd "C:\Users\Martin Paldam\Documents\AAA"
capture mkdir "`c(pwd)'\datasets" , public
*******************************************************************************************************
** Settings
*******************************************************************************************************
****Parameters:
local ci = 0.05 	// set confidence level
local r = 1			// number of funnels
local N = 500		// number of regressions in each funnel
local obs = 20		// number of observations of the regressions
local mint = 2 		// minimum acceptable t value (needed to determine point a)

****Variables:
local x_mean = 0 	// mean of x
local e_mean = 0	// mean of error term
local x_sd = 2		// standard deviation of x
local e_sd = 10		// standard deviation of error term

local beta = 1		// set true parameter for b1

local J = 50		// J Set
*******************************************************************************************************
*******************************************************************************************************

**Include programs
qui include programs_meta.do

	
local j = `J'
local m = `N'*`J'
set obs 1

foreach sr in SR0 SR1 SR2 SR3 SR4{	
	local fat_pet "`sr'_fat_pet_bf = r(`sr'_fat_pet_bf) `sr'_fat_pet_cbf = r(`sr'_fat_pet_cbf) `sr'_fat_pet_bm = r(`sr'_fat_pet_bm) `sr'_fat_pet_cbm = r(`sr'_fat_pet_cbm)"
	local peese "`sr'_peese_bf = r(`sr'_peese_bf) `sr'_peese_cbf = r(`sr'_peese_cbf) `sr'_peese_bp = r(`sr'_peese_bp) `sr'_peese_cbp = r(`sr'_peese_cbp)"
	local `sr' "`sr'_bbar = r(`sr'_bbar) `sr'_b_sd = r(`sr'_b_sd)  `sr'_t_mean = r(`sr'_t_mean) `fat_pet' `peese'  "
	gen `sr' = "`sr'"
}

simulate `SR0' `SR1' `SR2' `SR3' `SR4' ,  reps(`r') : build_funnel `beta' , x_mean(`x_mean') x_sd(`x_sd') e_mean(`e_mean')  e_sd(`e_sd') /*
			*/ n(`N') j(`j') ci(`ci') obs(`obs') mint(`mint')  save