Title: | Model Free Approach to Quantifying Surrogacy |
---|---|
Description: | Identifies an optimal transformation of a surrogate marker such that the proportion of treatment effect explained can be inferred based on the transformation of the surrogate and nonparametrically estimates two model-free quantities of this proportion. Details are described in Wang et al (2020) <doi:10.1093/biomet/asz065>. |
Authors: | Xuan Wang [aut], Layla Parast [aut, cre], Ming Yang [aut], Lu Tian [aut], Tianxi Cai [aut] |
Maintainer: | Layla Parast <[email protected]> |
License: | GPL-3 |
Version: | 1.0 |
Built: | 2024-11-16 02:47:32 UTC |
Source: | https://github.com/cran/OptimalSurrogate |
Simulated data with continuous surrogate marker
Simulated data with discrete surrogate marker
PTE estimation with continuous surrogate marker
pte_cont(sob, yob, aob, var = TRUE, conf.int = TRUE, rep = 500)
pte_cont(sob, yob, aob, var = TRUE, conf.int = TRUE, rep = 500)
sob |
CONTINUOUS surrogate marker |
yob |
outcome of interest (continuous or binary) |
aob |
treatment assignment (1: treatment; 0: control) |
var |
whether variance should be calculated (TRUE/FALSE) |
conf.int |
whether 95% confidence intervals should be calculated (TRUE/FALSE) |
rep |
number of resampling replications (default is 500) |
Estimates |
Estimates of the treatment effect on the priamry outcome, delta, the treatment effect on the transformation of the surrogate, delta.gs, two versions of the proportion of treatment effect explained by the surrogate, pte1 and pte2; if var = TRUE, standard error estimates are also provided (se); if conf.int = TRUE, 95% confidence intervals are also provided |
Transformed.S |
the transformed surrogate, g(s), for each value of the surrogate, s; if var = TRUE, standard error estimates are also provided (se); if conf.int = TRUE, 95% confidence intervals are also provided |
data(marker_cont) out <- pte_cont( sob = marker_cont$sob, yob = marker_cont$yob, aob = marker_cont$aob, var = FALSE, conf.int = FALSE) out x <- as.numeric(names(out$Transformed.S)) plot(x, out$Transformed.S, ylim = range(out$Transformed.S), type = "l", las = 1, xlab = "Surrogate Marker", ylab = "Optimal Transformation") out <- pte_cont( sob = marker_cont$sob, yob = marker_cont$yob, aob = marker_cont$aob, var = TRUE, conf.int = TRUE, rep = 1000) out$Estimates x <- as.numeric(rownames(out$Transformed.S)) plot(x, out$Transformed.S[, "est"], ylim = range(out$Transformed.S[, -2]), type = "l", las = 1, xlab = "Surrogate Marker", ylab = "Optimal Transformation") lines(x, out$Transformed.S[, "lower"], lty = 2) lines(x, out$Transformed.S[, "upper"], lty = 2)
data(marker_cont) out <- pte_cont( sob = marker_cont$sob, yob = marker_cont$yob, aob = marker_cont$aob, var = FALSE, conf.int = FALSE) out x <- as.numeric(names(out$Transformed.S)) plot(x, out$Transformed.S, ylim = range(out$Transformed.S), type = "l", las = 1, xlab = "Surrogate Marker", ylab = "Optimal Transformation") out <- pte_cont( sob = marker_cont$sob, yob = marker_cont$yob, aob = marker_cont$aob, var = TRUE, conf.int = TRUE, rep = 1000) out$Estimates x <- as.numeric(rownames(out$Transformed.S)) plot(x, out$Transformed.S[, "est"], ylim = range(out$Transformed.S[, -2]), type = "l", las = 1, xlab = "Surrogate Marker", ylab = "Optimal Transformation") lines(x, out$Transformed.S[, "lower"], lty = 2) lines(x, out$Transformed.S[, "upper"], lty = 2)
PTE estimation with discrete surrogate marker
pte_disc(sob, yob, aob, var = TRUE, conf.int = TRUE, rep = 500)
pte_disc(sob, yob, aob, var = TRUE, conf.int = TRUE, rep = 500)
sob |
DISCRETE surrogate marker |
yob |
outcome of interest (continuous or binary) |
aob |
treatment assignment (1: treatment; 0: control) |
var |
whether variance should be calculated (TRUE/FALSE) |
conf.int |
whether 95% confidence intervals should be calculated (TRUE/FALSE) |
rep |
number of resampling replications (default is 500) |
Estimates |
Estimates of the treatment effect on the priamry outcome, delta, the treatment effect on the transformation of the surrogate, delta.gs, two versions of the proportion of treatment effect explained by the surrogate, pte1 and pte2; if var = TRUE, standard error estimates are also provided (se); if conf.int = TRUE, 95% confidence intervals are also provided |
Transformed.S |
the transformed surrogate, g(s), for each value of the surrogate, s; if var = TRUE, standard error estimates are also provided (se); if conf.int = TRUE, 95% confidence intervals are also provided |
data(marker_disc) out <- pte_disc( sob = marker_disc$sob, yob = marker_disc$yob, aob = marker_disc$aob, var = FALSE, conf.int = FALSE) out out <- pte_disc( sob = marker_disc$sob, yob = marker_disc$yob, aob = marker_disc$aob, var = TRUE, conf.int = TRUE, rep = 1000) out
data(marker_disc) out <- pte_disc( sob = marker_disc$sob, yob = marker_disc$yob, aob = marker_disc$aob, var = FALSE, conf.int = FALSE) out out <- pte_disc( sob = marker_disc$sob, yob = marker_disc$yob, aob = marker_disc$aob, var = TRUE, conf.int = TRUE, rep = 1000) out