IPCWdf: Inverse probability of censoring weighting estimator for the...

In gsoutinho/survrec: Estimation of the multivariate distributions for censored gap times

Description

Provides estimates for the bivariate distribution function based on the Inverse Probability of Censoring Weighting estimator, IPCW.

Usage

IPCWdf(object, x, y, covariate, cov.value, bw, window = "gaussian")

Arguments

object	An object of class multidf.
x	The first time for obtaining estimates for the bivariate distribution function
y	The second time for obtaining estimates for the bivariate distribution function.
covariate	Name of the quantitative covariate.
cov.value	The value of the quantitative covariate.
bw	A single numeric value to compute a kernel density bandwidth. Use "dpik" for the KernSmooth package based selector or "np" for the 'npudensbw' function of the np package.
window	A character string specifying the desired kernel. See details below for possible options. Defaults to "gaussian" where the gaussian density kernel will be used.

Value

Vector with the IPWC estimates for the bivariate distribution function.

Author(s)

Gustavo Soutinho and Luis Meira-Machado.

References

de Una-Alvarez, J. and Meira-Machado, L. (2008). A simple estimator of the bivariate distribution function for censored gap times, Statistics and Probability Letters 78, 2440-2445.

See Also

KMWdf, LDMdf, LINdf and WCHdf.

Examples

b3state2 <- multidf(time1=bladder4state$y1, event1=bladder4state$d1, 
                    time=bladder4state$y1+bladder4state$y2,status=bladder4state$d2, 
                    size=bladder4state$size)
b3size <- multidf(time1=bladder3$t1, event1=bladder3$d1, 
                   time=bladder3$t2, status=bladder3$d2, size=bladder4state$size)

head(b3state2[[1]])
library(KernSmooth)
IPCWdf(object=b3state2, x=13, y=15, covariate="size", cov.value=3, window = "gaussian")
IPCWdf(object=b3state2, x=13, y=15, covariate="size", bw=2, cov.value=3, window = "gaussian")

IPCWdf(object=b3size, x=13, y=15, covariate="size", cov.value=3, window = "gaussian")
IPCWdf(object=b3size, x=13, y=15, covariate="size", bw=2, cov.value=3, window = "gaussian")

gsoutinho/survrec documentation built on Dec. 20, 2021, 1:46 p.m.