*! Version 1.00
* Calculation of Wald and Atkinson's R (partial correlation) after logistic 
* Joseph Hilbe, CRC/STB  5/5/92 

program define lwald
version 3.0
   if "$S_E_cmd"~="logistic" {
         error 301
   }
   local in "opt"
   local if "opt"
   local weight "fweight"
   parse "`*'"
   parse "$S_E_vl", parse(" ")
   if "$S_E_wgt"!="" {
       if "`weight'"=="" {
            local wgt "[{S_E_wgt}{S_E_exp}]"
            local subttl ", weighted statistics"
       }
       else local subttl ", unweighted statistics" 
    }   
    if "`if'`in'`all'"=="" {
         local if "$S_E_if"
         local in "$S_E_in"
    }
    qui {
    local depvar "$S_E_depv"
    local llo=(_result(6) + (-2*_result(2)))/2
    tempvar cof secof wald parcor
        gen `cof'   =.
        gen `secof' =.
        gen `wald'  =.
        gen `parcor'=.
        local i = 0
        mac shift

     noisily di _n(2)
     noisily di in gr "Wald Statistics and Atkinson's Partial Correlations (R)"
     noisily di in gr _n _col(1) "No.  Var       Wald     Prob(Chi)    Partial Cor"
     noisily di in gr _dup(48) "="
     
       while "`1'"!= "" {
          local i=`i'+1
            replace `cof'   = _b[`1'] in `i'
            replace `secof' = _se[`1'] in `i'
            replace `wald'  = (`cof'/`secof')^2 in `i'
            replace `parcor'= sqrt((`wald'-2)/(2*abs(`llo'))) 
            replace `parcor'=-`parcor' if (`cof')<0
            replace `parcor'= 0 if `wald'<2
            #delimit ;
            noisily di in gr `i' %9.4f _col(5) "`1'" %9.4f in ye _col(13) 
                             `wald'[`i'] %9.4f _col(23) chiprob(1,`wald'[`i']) 
                             %9.4f _col(37) `parcor'[`i'];
            #delimit cr
            mac shift
       }
    }
end