? Set the panel
?
SETPANEL ; Group = id ; Pds = ti $
?
? Define the base case list
?
NAMELIST ; xp  = your list of variables $ (do not include HEALTHY) $
NAMELIST ; xfe = your 3 variables $  (for example, age, income, hhkids)
?
? Compare fixed effects logit estimators
?
LOGIT     ; Lhs = doctor ; Rhs = xfe ; Panel ; Table = logit_c $
LOGIT     ; Lhs = doctor ; Rhs = xfe ; Panel ; FEM ; Table = logit_u $
MAKETABLE ; logit_c,logit_u $
?
? Compare quadrature and simulation based RE estimators
?
PROBIT ; Lhs = ... ; Rhs = one,... ; Panel ; Random $ (B&M)
PROBIT ; Lhs = ... ; Rhs = one,... ; Panel
       ; RPM ; Fcn = one(n) ; Halton ; Pts = 50 $
?
? Delta method for a nonlinear function
?
NAMELIST ; xp = ... $ your specification $
PROBIT   ; lhs = doctor ; rhs=xp $
WALD     ; if [ female=1 ]  ;  parameters = b 
         ; covariance = varb ; labels = kreg_b
         ; fn1 = n01(b1'xp)/phi(b1'xp) $

? Boostrapping a nonlinear function
?
PROCEDURE $
PROBIT    ;  quietly ; lhs=doctor ; rhs = xp $
CREATE    ;  imr = n01(b'xp)/phi(b'xp) $
CALC      ;  meanimr = female'imr/sum(female) $
ENDPROC   $
EXECUTE   ; n=100 ; bootstrap = meanimr  ; histogram $

? Bootstrapping a test statistic
?
PROBIT    ; lhs = doctor ; rhs = xp $ (Obtains the 'true' coefficients.)
CREATE    ; ysim = (b'xp + rnn(0,1)) > 0 $ (Simulates the homoscedastic data)
?
? LM test
?
PROBIT    ; Lhs = ysim ; Rhs = xp $ (this computes the restricted estimates)
PROBIT    ; Lhs = ysim ; Rhs = xp ; Het ; Hfn = female ; start = b,0 ; Maxit = 0 $
PROCEDURE $
PROBIT    ;  quietly ; lhs=ysim ; rhs = xp $
PROBIT    ;  quietly ; lhs=ysim ; rhs=x
          ;  het ; hfn=female;start=b,0 ; maxit=0 $
ENDPROC   $
EXECUTE   ; n=100 ; bootstrap = lmstat $
HISTOGRAM ; rhs = bootstrp $

? Bootstrapping a Parameter Vector
?
PROCEDURE $
LOGIT     ; quiet ; Lhs = healthy ; Rhs = xp ; Prob = p $
CREATE    ; scale = p*(1-p) $
CALC      ; avgscale = xbr(scale) $
MATRIX    ; ape = avgscale * b $
ENDPROC   $
EXECUTE   ; n = 50 ; bootstrap = ape $
LOGIT     ; quiet ; Lhs = healthy ; Rhs = xp ; Marginal  $