* Swishex.do
* 11/6/97

clear
pause on

capture log close
set more off
set matsize 200

log using swishex.log, replace
di %dD_M_Y date("${S_DATE}", "dmy") "    $S_TIME"
use dta/master01, replace
gen double Guesprob = 0
quietly summarize lnpce, detail
replace Guesprob = 1 if lnpce > _result(10)
#delimit ;
eq main : haz docdist nursdist phardist father mother gender
lnpce crop_siz urban rural ;
eq regime : Guesprob lnpce members maxfe1 maxfe2 highfe
  urban rural KZN ECape NorP reticul elect ; 
switchr main regime [pweight=rsweight] , cl(clustnum) noisyn(200) sige;

#delimit cr
quietly summarize Guesprob, detail
local av = _result(3)
di "Average of probability vector is: `av'"
di "Log-likelihood is $S_E_ll"