nrsp.GTDL: Normally-transformed randomized survival probability...

In carrascojalmar/GTDL: The Generalized Time-Dependent Logistic Family

Normally-transformed randomized survival probability residuals for the GTDL model

Description

Normally-transformed randomized survival probability residuals for the GTDL model

Usage

nrsp.GTDL(t, formula, pHat, censur)

Arguments

t	non-negative random variable representing the failure time and leave the snapshot failure rate, or danger.
formula	The structure matrix of covariates of dimension n x p.
pHat	Estimate of the parameters from the GTDL model.
censur	Censoring status 0=censored, a=fail.

Value

Normally-transformed randomized survival probability residuals

References

• Li, L., Wu, T., e Cindy, F. (2021). Model diagnostics for censored regression via randomized survival probabilities. Statistics in Medicine, 40, 1482–1497.

• de Oliveira, L. E. F., dos Santos L. S., da Silva, P. H. F., Fabio, L. C., Carrasco, J. M. F.(2022). Análise de resíduos para o modelo logístico generalizado dependente do tempo (GTDL). Submitted.

Examples

### Example 1

require(survival)
data(lung)
lung <- lung[-14,]
lung$sex <- ifelse(lung$sex==2, 1, 0)
lung$ph.ecog[lung$ph.ecog==3]<-2
t1 <- lung$time
formula1 <- ~lung$sex+factor(lung$ph.ecog)+lung$age
censur1 <- ifelse(lung$status==1,0,1)
start1 <- c(0.03,0.05,-1,0.7,2,-0.1)
fit.model1 <- mle2.GTDL(t = t1,start = start1,
           formula = formula1,
           censur = censur1)
r1 <- nrsp.GTDL(t = t1,formula = formula1 ,pHat = fit.model1$Coefficients[,1],
             censur = censur1)
r1

### Example 2

data(tumor)
t2 <- tumor$time
formula2 <- ~tumor$group
censur2 <- tumor$censured
start2 <- c(1,-0.05,1.7)
fit.model2 <- mle2.GTDL(t = t2,start = start2,
                       formula = formula2,
                       censur = censur2)
r2 <- nrsp.GTDL(t = t2,formula = formula2, pHat = fit.model2$Coefficients[,1],
            censur = censur2)
r2

carrascojalmar/GTDL documentation built on May 18, 2022, 1:12 p.m.