ICSchedule: Generates interval-censored Failure Time data with schedule...

Description Usage Arguments Details Value Examples

Description

Generates interval-censored Failure Time data with schedule visits. The failure time follows the Weibull regression model

Usage

1
2
ICSchedule(n, shape, scale, visits, p, X = as.matrix(0),
  beta = as.matrix(0))

Arguments

n

number of observations to be generated

shape

the shape parameter of Weibull distribution

scale

the scale parameter of Weibull distribution

visits

the vector with the time of visits

p

the probability of a observation fail each visit

X

the matrix with the values of covariates. The number of rows of X must be equal n

beta

the vector with the value of the regression coefficients. The lenght of beta must be equal the number of covariaties

Details

First the function generates the failure time using [rweibull()]. The shape is the same and the second parameter is lambda. Lambda is equal to:

lambda = scale*exp(X*beta)

If there is no value to X and beta lambda = scale. Then is generated from the function [rbinom()] for each observation to know what visits really happened. Then is created the vector with the result. If the failure time is greater than the maximum of the visited time the, the first time of that observation is the maximum of the time visited, the second time is [Inf] and the failure time is right-censored. If the failure time less than the minimum of the visited times, the first time of that observation is [-Inf], the second time is the minimum of the visited time and the failured time is left-censored. Orthewise the failure time is between two visited times, than the first time and second time are these visited times and the failure is interval-censored.

Value

The function return a data.frame on what:

L

the lower limit of the interval generated

R

the upper limit of the interval generated

event

the indicator of interval censored data in the same way that is used in [Surv()] function from survival package.

X

the value of the covariates

Examples

1
ICSchedule(n = 50, shape = 2, scale = 2, visits = c(1,2,3,4), p = 0.1)

felipeamorimbr/ICdata documentation built on May 30, 2019, 1:59 p.m.