* Simulation of contingency table chisquares

cap log close
cap prog drop sim
prog def sim
* 1=# of obs, 2=# of reps, 3=square table dimension.
local obs `1'
local reps `2'
local rc `3'
tempvar row col
drop _all
set obs `obs'
gen int `row'=0
gen int `col'=0
quietly {
	local i 1
	while `i'<=`reps' {
		replace `row'=1+int(`rc'*uniform())
		replace `col'=1+int(`rc'*uniform())
		tab `row' `col', chi
		noi di `obs',`rc',_result(4)
		local i=`i'+1
	}
}
end

set more 1
log using sim.sto, replace
set seed 655
sim  5 500 3
sim 10 500 3
sim 20 500 3
sim 45 500 3
log close