mmmultPost: Bayesian mixed effects model with a DP prior on by-subject...

Description Usage Arguments Value Note Author(s) See Also

Description

An internal function to dpgrowmult

Usage

1
2
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 Γ prior on the DP concentration parameter. The rate parameter is set of 1.

ratebeta

The rate parameter for the Γ prior on the DP concentration parameter. Default value is 1.

strength.mm

The shape and rate parameters for the Γ prior on the CAR precision parameter, τ_{γ}

Value

res A list object containing MCMC runs for all model parameters.

Note

Intended as an internal function for dpgrowmult

Author(s)

Terrance Savitsky tds151@gmail.com

See Also

dpgrowmm


growcurves documentation built on May 2, 2019, 7:03 a.m.