Description Usage Arguments Value References See Also Examples
Simulation of continuous longitudinal outcome and competing risks data Currently, only the simulation in Elashoff et al(2008) is implemented.
1 2 3 4 5 6 7 8 9 10 11 12 13 | SimDataC(
k_val,
p1_val,
p1a_val,
p2_val,
g_val,
truebeta,
truegamma,
randeffect,
yfn,
cfn,
mfn
)
|
k_val |
The number of subjects in study. |
p1_val |
The dimension of fixed effects in longitudinal measurements. |
p1a_val |
The dimension of random effects in longitudinal measurements. |
p2_val |
The dimension of fixed effects in competing risks failure time data. |
g_val |
The number of type of failure in competing risks data. |
truebeta |
True values for beta, the longitudinal coefficients. |
truegamma |
True values for gamma, the survival coefficients. |
randeffect |
True values for random effects in longitudinal and competing risks parts,namely in the order of σ,σ_b,ν_2,σ_u. |
yfn |
Filename of genereated Y matrix for longitudinal measurements in long format. |
cfn |
Filename of genereated C matrix for competing risks failure time data. |
mfn |
Filename of genereated M vector to indicate the number of longitudinal measurements per subject. |
Files with names yfn, cfn and mfn.
censoring_rate | Censoring rate of the survival data. |
rate1 | Censoring rate of competing risk 1. |
rate2 | Censoring rate of competing risk 2. |
yfn | Filename of genereated Y matrix for longitudinal measurements. |
cfn | Filename of genereated C matrix for competing risks failure time data. |
mfn | Filename of genereated M vector to indicate the number of longitudinal measurements per subject. |
Elashoff, Robert M., Gang Li, and Ning Li. "A joint model for longitudinal measurements and survival data in the presence of multiple failure types." Biometrics 64.3 (2008): 762-771.
1 2 3 4 5 6 7 8 9 10 11 12 13 | # A toy example testint data generations
require(JMcmprsk)
set.seed(123)
yfn=tempfile(pattern = "", fileext = ".txt")
cfn=tempfile(pattern = "", fileext = ".txt")
mfn=tempfile(pattern = "", fileext = ".txt")
k_val=30;p1_val=4;p1a_val=1; p2_val=2;g_val=2;
truebeta=c(10,-1,1.5,0.6);truegamma=c(0.8,-1,0.5,-1); randeffect=c(5,0.5,0.5,0.5);
#writing files
SimDataC(k_val, p1_val, p1a_val, p2_val, g_val,truebeta,
truegamma, randeffect, yfn, cfn, mfn)
jmc_0(p=4,yfn,cfn,mfn,point=6)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.