* remember to set seed
clear
*set seed 5734
set obs 400
gen nbhd=_n
gen greenspace=uniform()>.3
gen U_nbhd=uniform()>.5
gen U_nbhd2=uniform()>.2
expand 25
sort nbhd
by nbhd: gen id=_n
gen male=uniform()<.4
gen BMI=invnorm(uniform())+1*U_nbhd
gen dep=invnorm(uniform())*.5+U_nbhd2
gen U_ind=uniform()<.5
gen sbp=.3*greenspace+1*U_nbhd+.2*male*.2*BMI+.5*dep+.1*U_ind+invnorm(uniform())
reg sbp greenspace male BMI dep
xtset nbhd id
xtmixed  sbp greenspace male BMI dep || nbhd: