/* Assess the bootstrap models:
1. Load all of the boostrapped models
2. Make sure the models differ
3. Assess the variation we see in the parameters (including none)
*/

local maxrep 25
local directory psid

local list : dir "../FEM_CPP_settings/`directory'/models/models_rep1" files "*est"
local n : word count `list'

local y
local z

forvalues i = 1/`n' {
	local counter = 0
	forvalues j = 2/`maxrep' {
		local filenm : word `i' of `list'
		qui tex_equal, file1("../FEM_CPP_settings/`directory'/models/models_rep1/`filenm'") file2("../FEM_CPP_settings/`directory'/models/models_rep`j'/`filenm'")
		local outcome = r(comparison)
		if `outcome' == 0 {
			* files differ
			local counter = `counter' + 1 
		}
	}
	if `counter' == 0 {
		di "`filenm' does not differ. THIS IS NOT BOOTSTRAPPED!"
		local y `y' `i'
	}
	if `counter' > 0 {
		local cntp1 = `counter' + 1
		di "`filenm' differs in `cntp1' files out of `maxrep'."
		local z `z' `i'
	}
}

di "These models ARE bootstrapped: "
foreach model of local z {
	local filnm : word `model' of `list'
	di "`filnm'"
}

di "These models ARE NOT bootstrapped: "
foreach model of local y {
	local filnm : word `model' of `list'
	di "`filnm'"
}


capture log close