Description Usage Arguments Value Note Author(s) References See Also Examples
Employs a Dirichlet Process (DP) prior on the set of by-subject random effect parameters
under repeated waves of measurements to allow the number of random effect parameters specified
per subject, q
, to be equal to the number of measurement waves, T
.
Random effects are grouped by subject and
all q
parameters receive the DP prior. The resulting joint marginal
distribution over the data is a DP mixture.
1 2 |
y |
A univariate continuous response, specified as an N x 1 matrix or vector,
where |
subject |
The objects on which repeated measures are conducted that serves
as the random effects grouping factor. Input as an N x 1 matrix or vector of
subject-measure cases in either integer or character formt; e.g.
|
trt |
An integer or character matrix/vector of length |
time |
A univariate vector of length |
n.random |
The desired number of subject random effect terms, |
n.fix_degree |
The desired polynomial order in time to use for
generating time-based fix effects.
The fixed effects matrix will be constructed as,
|
formula |
Nuisance fixed and random effects may be entered in
|
random.only |
A Boolean variable indicating whether the input formula contains
random (for fixed) effects in the case that only
one set are entered. If excluded and |
data |
a |
n.iter |
Total number of MCMC iterations. |
n.burn |
Number of MCMC iterations to discard. |
n.thin |
Gap between successive sampling iterations to save. |
shape.dp |
Shape parameter under a c ~ G(shape.dp, 1)
prior on the concentration parameter of the DP (prior
on the set of random effects parameters, b_1, ..., b_n ~ DP(c,G_0)
where |
rate.dp |
Rate parameter under a c ~ G(shape.dp, rate.dp) prior on the concentration parameter of the DP. |
plot.out |
A boolean variable indicating whether user wants to return plots with output results. Defaults to |
option |
Modeling option, of which there are two: 1. |
S3 dpgrow
object, for which many methods are available to return and
view results. Generic functions applied
to an object, res
of class dpgrow
, includes:
summary(res) |
returns |
print(summary(res)) |
prints contents of summary to console. |
plot(res) |
returns results plots, including the set of subject random effects values and credible intervals, a sample of by-subject growth curves, mean growth curves split by each treatment and control, as well as selected trace plots for number of clusters and for precision parameters for the likehilood and random effects. Lastly, a trace plot for the deviance statistic is also included. |
samples(res) |
contains ( |
resid(res) |
contains the model residuals. |
The intended focus for this package are data where both number of subjects and number of
repeated measures are limited. A DP prior
is placed on the by-subject random effects to borrow strength across subjects for
each estimation of each subject's growth curve. The
imposition of the DP prior also allows the resulting posterior distributions
over the subject random effects to be non-Gaussian.
The dpgrow
function is very similar to dpgrowmm
;
only the latter includes a separate set of random effects not grouped
by subject (e.g. for treatment dosages allocated to subjects) mapped
back to subject-time cases through a multiple membership design matrix.
The dpgrowmult
function generalizes dpgrowmm
by allowing
more than one multiple membership effects term.
See Savitsky and Paddock (2011) for detailed model constructions.
Terrance Savitsky tds151@gmail.com Susan Paddock paddock@rand.org
S. M. Paddock and T. D. Savitsky (2012) Bayesian Hierarchical Semiparametric Modeling of Longitudinal Post-treatment Outcomes from Open-enrollment Therapy Groups, submitted to: JRSS Series A (Statistics in Society).
T. D. Savitsky and S. M. Paddock (2011) Visual Sufficient Statistics for Repeated Measures data with growcurves for R, submitted to: Journal of Statistical Software.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 | ## Not run:
## extract simulated dataset
library(growcurves)
data(datsim)
## attach(datsim)
## run dpgrow mixed effects model; returns object of class "dpgrow"
shape.dp = 4
res = dpgrow(y = datsim$y, subject = datsim$subject,
trt = datsim$trt, time = datsim$time,
n.random = datsim$n.random,
n.fix_degree = 2, n.iter = 10000,
n.burn = 2000, n.thin = 10,
shape.dp = shape.dp, option = "dp")
plot.results = plot(res) ## ggplot2 plot objects, including growth curves
summary.results = summary(res) ## parameter credible intervals, fit statistics
samples.posterior = samples(res) ## posterior sampled values
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.