CalcNpmleRSP {ICcalib}

Calculating the probabilities of positive binary exposure status at a given time point using a nonparametric risk-set calibration models

Description

For a given time point, calculate the probability of positive exposure value for multiple observations (participants). The function first fits the nonparametric risk-set calibration models at each main event time point and then calculates the probabilities of positive binary exposure status.

Usage

CalcNpmleRSP(w, w.res, point, obs.tm)

Arguments

w

A matrix of time points when measurements on the binary covariate were obtained.

w.res

A matrix of measurement results of the binary covariate. Each measurement corresponds to the time points in w

point

The time point at which the probabilities are estimated

obs.tm

Vector of observed main event time or censoring time

Details

This function calculates the NPMLE at each main event time point and then provides the estimated probabilities for positive exposure status at time point.

Value

A vector of estimated probabilities of positive exposure status at time point.

See Also

ic_np

Examples

# Simulate data set
sim.data <- ICcalib:::SimCoxIntervalCensSingle(n.sample = 200, lambda = 0.1, 
                                               alpha = 0.25, beta0 = log(0.5), 
                                               mu = 0.2, n.points = 2, 
                                               weib.shape = 1, weib.scale = 2)
# Calculate the conditional probabilities of binary covariate=1 at time one
# Unlike CalcNpmle, CalcNpmleRSP includes the calibration model fitting
probs <- CalcNpmleRSP(w = sim.data$w, w.res = sim.data$w.res, point = 1, 
                      obs.tm = sim.data$obs.tm)
summary(probs)