program define z_rvrfy
*! version 1.1.1 02Feb96

   version 4.0
   local n1 : rownames(`1')
   local n2 : rownames(`2')
   local c = ("`n1'" == "`n2'")
   local n2 = .5*log((1+`1'[1,2])/(1-`1'[1,2])) + /*
         */    invnorm((1+`3'[1,3])/2)/sqrt(`3'[1,1]-3)
   local n1 = .5*log((1+`2'[1,2])/(1-`2'[1,2]))
   local c = `c' * (float(`n2') == float(`n1'))
   if `c' { exit }
   di in red "output from z_r not found"
   error 499
end