mets (version 0.1-12)

casewise.test: Estimates the casewise concordance based on Concordance and marginal estimate using timereg and performs test for independence

Description

.. content for description (no empty lines) ..

Usage

casewise.test(conc, marg, test = "no-test")

Arguments

conc

Concordance

marg

Marginal estimate

test

Type of test for independence assumption. "conc" makes test on concordance scale and "case" means a test on the casewise concordance

Details

Uses cluster based conservative standard errors for marginal

Examples

data(prt);

prt <- prt[which(prt$id %in% sample(unique(prt$id),10000)),]
### marginal cumulative incidence of prostate cancer
times <- seq(60,100,by=10)
outm <- comp.risk(Surv(time,status==0)~+1,data=prt,prt$status,causeS=2,times=times)

cifmz <- predict(outm,X=1,uniform=0,resample.iid=1)
cifdz <- predict(outm,X=1,uniform=0,resample.iid=1)

### concordance for MZ and DZ twins
cc <- bicomprisk(Hist(time,status)~strata(zyg)+id(id),data=prt,cause=c(2,2))
cdz <- cc$model$"DZ"
cmz <- cc$model$"MZ"

cdz <- casewise.test(cdz,cifmz,test="case")
cmz <- casewise.test(cmz,cifdz,test="conc")

plot(cmz,ylim=c(0,0.7),xlim=c(60,100))
par(new=TRUE)
plot(cdz,ylim=c(0,0.7),xlim=c(60,100))