program define gr_log
version 6.0

local or = `2'
local xunits = `3'
local b1 = ln(`or')

* make summary of logistic data from equation
set obs `xunits'
generate pgty = 1 - 1/(1 + exp(score))