Usage
mmmultPost(y, X, Z, Wcases, Mmats, Omegas, omegapluses, ngs, subjects, typet,
niter, nburn, nthin, shapealph, ratebeta, strength.mm)
Arguments
y
An N x 1 response (of subject-measure
cases)
X
Fixed effects design matrix
Z
Random effects design matrix. Assumed grouped
by subjects
Wcases
A list object containing N x 1
multiple membership (MM) weight matrices; one for each MM
term.
Mmats
A list object containing S[i] x G[i]
matrices, where S[i] are the number of effects for
MM term i and G[i] is the number of
(unique) groups for the S[i] effects. The length
of t
Omegas
A list object containing S[i] x S[i]
unnormalized adjacency matrices for MM term i,
each with entries equal to 1 where two effects
communicate and 0, otherwise. Diagonal elements are
zero. The length of this list is
omegapluses
A list object containing S x 1
vector of row sums of Omega[[i]]. The length of
omegaplus should equal Omega.
ngs
A numeric vector containing the number of
total groups for each MM term under prior option
"mmcar".
subjects
An N x 1 set of subject
identifiers.
typet
A numeric vector of length equal to the
number of MM terms, where each entry specifies the prior
formulation for that effect term (block). Prior
formulation options are 1 = "mmcar", 2 = "mmi", 3 =
"mmigrp", 4 = "mmdp".
niter
The number of MCMC iterations
nburn
The number of MCMC burn-in iterations to
discard
nthin
The step increment of MCMC samples to
return
shapealph
The shape parameter for the $\Gamma$
prior on the DP concentration parameter. The rate
parameter is set of 1.
ratebeta
The rate parameter for the $\Gamma$
prior on the DP concentration parameter. Default value is
1.
strength.mm
The shape and rate parameters for the
$\Gamma$ prior on the CAR precision parameter,
$\tau_{\gamma}$