simGLcox: Simulation of two-stage recurrent events data based on...

In kkholst/mets: Analysis of Multivariate Event Times

Simulation of two-stage recurrent events data based on Cox/Cox or Cox/Ghosh-Lin structure

Description

Simulation of two-stage recurrent events data based on Cox/Cox or Cox/Ghosh-Lin structure. type=3 will generate Cox/Cox twostage mode, type=2 will generate Ghosh-Lin/Cox model. If the variance is var.z=0, then generates data without any dependence or frailty. If model="twostage" then default is to generate data from Ghosh-Lin/Cox model, and if type=3 then will generate data with marginal Cox models (Cox/Cox). Simulation based on linear aproximation of hazard for two-stage models based on grid on time-scale. Must be sufficientyly fine.

Usage

simGLcox(
  n,
  base1,
  drcumhaz,
  var.z = 0,
  r1 = NULL,
  rd = NULL,
  rc = NULL,
  fz = NULL,
  fdz = NULL,
  model = c("twostage", "frailty", "shared"),
  type = NULL,
  share = 1,
  cens = NULL,
  nmin = 100,
  nmax = 1000
)

Arguments

n	number of id's
base1	baseline for cox/ghosh-lin models
drcumhaz	baseline for terminal event
var.z	variance of gamma frailty
r1	relative risk term for baseline
rd	relative risk term for terminal event
rc	relative risk term for censorings
fz	possible transformation (function) of frailty term
fdz	possible transformation (function) of frailty term for death
model	twostage, frailty, shared (partly shared two-stage model)
type	type of simulation, default is decided based on model
share	to fit patly shared random effects model
cens	censoring rate for exponential censoring
nmin	default 100, at least nmin or number of rows of the two-baselines max(nmin,nrow(base1),nrow(drcumhaz)) points in time-grid from 0 to maximum time for base1
nmax	default 1000, at most nmax points in time-grid

Details

Must specify baselines of recurrent events and terminal event and possible covariate effects.

References

Scheike (2024), Twostage recurrent events models, under review.

kkholst/mets documentation built on June 14, 2025, 9:19 a.m.