Description Usage Arguments Details Value See Also Examples
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.
1 | CalcNpmleRSP(w, w.res, point, obs.tm)
|
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 |
point |
The time point at which the probabilities are estimated |
obs.tm |
Vector of observed main event time or censoring time |
This function calculates the NPMLE at each main event time point and then provides the estimated probabilities for positive
exposure status at time point
.
A vector of estimated probabilities of positive exposure status at time point
.
1 2 3 4 5 6 7 8 9 10 | # 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)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.