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 this list is equal to the number of terms under prior option set to "mmigrp"
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 equal to the number terms with prior option set to "mmcar"
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}$