mets (version 1.3.11)

casewise: Estimate Casewise Concordance from prodlim Objects

Description

Estimates the casewise concordance based on concordance and marginal estimates derived from prodlim objects. Unlike test_casewise, this function does not perform hypothesis testing but focuses on estimation and plotting.

Usage

casewise(conc, marg, cause.marg)

Value

An object of class "casewise" containing:

casewise

Matrix with time, casewise concordance, and standard errors.

marg

Matrix with time, marginal CIF, and standard errors.

concordance

Matrix with time, concordance, and standard errors.

timer

Time points used.

P1, se.P1

Extracted concordance values and SEs.

Arguments

conc

Concordance object from prodlim (output of bicomprisk with prodlim=TRUE).

marg

Marginal cumulative incidence object from prodlim (output of prodlim).

cause.marg

Specifies which cause should be used for the marginal CIF based on the Event object.

Author

Thomas Scheike

See Also

test_casewise, bicomprisk

Examples

 ## Reduce Ex.Timings
library(prodlim)
data(prt);
prt <- force_same_cens(prt,cause="status")

### marginal cumulative incidence of prostate cancer
outm <- prodlim(Hist(time,status)~+1,data=prt)

times <- 60:100
cifmz <- predict(outm,cause=2,time=times,newdata=data.frame(zyg="MZ"))
cifdz <- predict(outm,cause=2,time=times,newdata=data.frame(zyg="DZ"))

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

cdz <- casewise(cdz,outm,cause.marg=2) 
cmz <- casewise(cmz,outm,cause.marg=2)

plot(cmz,ci=NULL,ylim=c(0,0.5),xlim=c(60,100),legend=TRUE,col=c(3,2,1))
par(new=TRUE)
plot(cdz,ci=NULL,ylim=c(0,0.5),xlim=c(60,100),legend=TRUE)
summary(cdz)
summary(cmz)