Description Usage Arguments Details Value Author(s) References See Also Examples
This function simulates a dataset based on the negative binomial mixed-effect independent/AR(1) model with two treatment groups described in Zhao et al (2013). The group mean can be different at each time point, but no other covariates are allowed.
See fitParaIND
, fitParaAR1
for details of the model explanations.
1 2 3 4 5 |
gdist |
The distribution of the random effect term G[i]. If If If If If the generated values are negative, they are truncated to zero. If |
n |
The number of patients. It must be an even number. |
sn |
The number of repeated measures per patient. Generated datasets are balanced design. |
th |
If If If If |
u1 |
A vector of length sn, specifying the mean of the treatment group 1 E(Y[ij]) = |
u2 |
A vector of length sn, specifying the mean of the treatment group 2 E(Y[ij]) = |
a |
The dispersion parameter α of the negative binomial mixed-effect independent model. See description in lmeNB. |
d |
If |
othrp |
If |
The generated datasets have equal number of scans per person.
The number of patients in the two groups are the same.
If gdist=="GN"
, datasets are generated from:
othrp$p.mx
*N(mean
=othrp$u.n
,s.d
=othrp$s.n
) + (1-othrp$p.mx
)*gamma(scale=th,shape)
,
where shape
of the gamma distribution is chosen to ensure E(G[i])=1.
id |
The vector of length |
vn |
The vector of length |
gp |
The vector of length |
y |
The vector of length |
g |
The vector of length |
Gpara |
The record of the distribution and parameter specifications used to generate the dataset |
Zhao, Y. and Kondo, Y.
Detection of unusual increases in MRI lesion counts in individual multiple sclerosis patients. (2013) Zhao, Y., Li, D.K.B., Petkau, A.J., Riddehough, A., Traboulsee, A., Journal of the American Statistical Association.
The main function to fit the Negative Binomial mixed-effect model:
lmeNB
,
The functions to fit related models:
fitParaIND
,
fitParaAR1
,
fitSemiIND
,
fitSemiAR1
,
The subroutines of index.batch
to compute the conditional probability index:
jCP.ar1
,
CP1.ar1
,
MCCP.ar1
,
CP.ar1.se
,
CP.se
,
jCP
,
1 2 3 4 5 | ## Not run:
## See the examples in help files of fitParaIND, fitParaAR1, fitSemiIND, fitSemiAR1 and lmeNB
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.