SimDataC: Data simulation of continuous outcomes and competing risks

In JMcmprsk: Joint Models for Longitudinal and Competing Risks Data

Description

Simulation of continuous longitudinal outcome and competing risks data Currently, only the simulation in Elashoff et al(2008) is implemented.

Usage

SimDataC(
  k_val,
  p1_val,
  p1a_val,
  p2_val,
  g_val,
  truebeta,
  truegamma,
  randeffect,
  yfn,
  cfn,
  mfn
)

Arguments

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.

Value

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.

References

• 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.

See Also

SimDataO

Examples

# 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)

JMcmprsk documentation built on March 22, 2021, 9:07 a.m.