SimDataO: Data simulation of ordinal outcomes and competing risks

Description Usage Arguments Value References See Also Examples

View source: R/RcppExports.R

Description

Simulation of ordinal longitudinal outcome and competing risks data Currently, only the simulation in Li et al(2010) is implemented.

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
SimDataO(
  k_val,
  p1_val,
  p1a_val,
  p2_val,
  g_val,
  truebeta,
  truetheta,
  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.

truetheta

True values for theta, subset of the non-proportional odds 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

See Also

SimDataC

Examples

 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=50;p1_val=3;p1a_val=1; p2_val=2;g_val=2;
truebeta=c(-1,1.5,0.8);truetheta=c(-0.5,1);truegamma=c(0.8,-1,0.5,-1); randeffect=c(1,0.5,0.5);
#writing files
SimDataO(k_val, p1_val, p1a_val, p2_val, g_val,
      truebeta, truetheta, truegamma, randeffect, yfn,  cfn,  mfn)

jmo_0(p=3,s=1, yfn,cfn,mfn,point=10,do.trace = TRUE)

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