Question

What is zdelta/2 when delta is 0.01?

Answer 1

Model:
y_{hi} = 1 with pr=exp(x_{hi}'beta_h)/(1+exp(x_{hi}'beta_h). beta_h is nvar x 1.
h=1,...,length(lgtdata) units or "respondents" for survey data.

beta_h= ZDelta[h,] + u_h.
Note: here ZDelta refers to Z%*%Delta, ZDelta[h,] is hth row of this product.
Delta is an nz x nvar array.

u_h ~ N(0,V_{beta}).
Priors:
delta= vec(Delta) ~ N(vec(Deltabar),V_{beta} (x) ADelta^{-1})
V_{beta} ~ IW(nu,V)

Lists contain:

lgtdatalist of lists with each cross-section unit MNL data
lgtdata[[h]]$y n_h vector of binary outcomes (0,1)
lgtdata[[h]]$X n_h by nvar design matrix for hth unit
Deltabarnz x nvar matrix of prior means (def: 0)
ADelta prior prec matrix (def: .01I)
nu d.f. parm for IW prior on norm comp Sigma (def: nvar+3)
V pds location parm for IW prior on norm comp Sigma (def: nuI)
sbeta scaling parm for RW Metropolis (def: .2)
R number of MCMC draws
keep MCMC thinning parm: keep every keepth draw (def: 1)
Value

a list containing:

Deltadraw R/keep x nz*nvar matrix of draws of Delta
betadraw nlgt x nvar x R/keep array of draws of betas
Vbetadraw R/keep x nvar*nvar matrix of draws of Vbeta
llike R/keep vector of log-like values
reject R/keep vector of reject rates over nlgt units
Note

Some experimentation with the Metropolis scaling paramter (sbeta) may be required.