View source: R/ICA.Sample.ContCont.R
ICA.Sample.ContCont | R Documentation |
The function ICA.Sample.ContCont
quantifies surrogacy in the single-trial causal-inference framework. It provides a faster alternative for ICA.ContCont
. See Details below.
ICA.Sample.ContCont(T0S0, T1S1, T0T0=1, T1T1=1, S0S0=1, S1S1=1, T0T1=seq(-1, 1, by=.001),
T0S1=seq(-1, 1, by=.001), T1S0=seq(-1, 1, by=.001), S0S1=seq(-1, 1, by=.001), M=50000)
T0S0 |
A scalar or vector that specifies the correlation(s) between the surrogate and the true endpoint in the control treatment condition that should be considered in the computation of |
T1S1 |
A scalar or vector that specifies the correlation(s) between the surrogate and the true endpoint in the experimental treatment condition that should be considered in the computation of |
T0T0 |
A scalar that specifies the variance of the true endpoint in the control treatment condition that should be considered in the computation of |
T1T1 |
A scalar that specifies the variance of the true endpoint in the experimental treatment condition that should be considered in the computation of |
S0S0 |
A scalar that specifies the variance of the surrogate endpoint in the control treatment condition that should be considered in the computation of |
S1S1 |
A scalar that specifies the variance of the surrogate endpoint in the experimental treatment condition that should be considered in the computation of |
T0T1 |
A scalar or vector that contains the correlation(s) between the counterfactuals T0 and T1 that should be considered in the computation of |
T0S1 |
A scalar or vector that contains the correlation(s) between the counterfactuals T0 and S1 that should be considered in the computation of |
T1S0 |
A scalar or vector that contains the correlation(s) between the counterfactuals T1 and S0 that should be considered in the computation of |
S0S1 |
A scalar or vector that contains the correlation(s) between the counterfactuals S0 and S1 that should be considered in the computation of |
M |
The number of runs that should be conducted. Default |
Based on the causal-inference framework, it is assumed that each subject j has four counterfactuals (or potential outcomes), i.e., T_{0j}
, T_{1j}
, S_{0j}
, and S_{1j}
. Let T_{0j}
and T_{1j}
denote the counterfactuals for the true endpoint (T
) under the control (Z=0
) and the experimental (Z=1
) treatments of subject j, respectively. Similarly, S_{0j}
and S_{1j}
denote the corresponding counterfactuals for the surrogate endpoint (S
) under the control and experimental treatments, respectively. The individual causal effects of Z
on T
and S
for a given subject j are then defined as \Delta_{T_{j}}=T_{1j}-T_{0j}
and \Delta_{S_{j}}=S_{1j}-S_{0j}
, respectively.
In the single-trial causal-inference framework, surrogacy can be quantified as the correlation between the individual causal effects of Z
on S
and T
(for details, see Alonso et al., submitted):
\rho_{\Delta}=\rho(\Delta_{T_{j}},\:\Delta_{S_{j}})=\frac{\sqrt{\sigma_{S_{0}S_{0}}\sigma_{T_{0}T_{0}}}\rho_{S_{0}T_{0}}+\sqrt{\sigma_{S_{1}S_{1}}\sigma_{T_{1}T_{1}}}\rho_{S_{1}T_{1}}-\sqrt{\sigma_{S_{0}S_{0}}\sigma_{T_{1}T_{1}}}\rho_{S_{0}T_{1}}-\sqrt{\sigma_{S_{1}S_{1}}\sigma_{T_{0}T_{0}}}\rho_{S_{1}T_{0}}}{\sqrt{(\sigma_{T_{0}T_{0}}+\sigma_{T_{1}T_{1}}-2\sqrt{\sigma_{T_{0}T_{0}}\sigma_{T_{1}T_{1}}}\rho_{T_{0}T_{1}})(\sigma_{S_{0}S_{0}}+\sigma_{S_{1}S_{1}}-2\sqrt{\sigma_{S_{0}S_{0}}\sigma_{S_{1}S_{1}}}\rho_{S_{0}S_{1}})}},
where the correlations \rho_{S_{0}T_{1}}
, \rho_{S_{1}T_{0}}
, \rho_{T_{0}T_{1}}
, and \rho_{S_{0}S_{1}}
are not estimable. It is thus warranted to conduct a sensitivity analysis.
The function ICA.ContCont
constructs all possible matrices that can be formed based on the specified vectors for \rho_{S_{0}T_{1}}
, \rho_{S_{1}T_{0}}
, \rho_{T_{0}T_{1}}
, and \rho_{S_{0}S_{1}}
, and retains the positive definite ones for the computation of \rho_{\Delta}
.
In contrast, the function ICA.ContCont
samples random values for \rho_{S_{0}T_{1}}
, \rho_{S_{1}T_{0}}
, \rho_{T_{0}T_{1}}
, and \rho_{S_{0}S_{1}}
based on a uniform distribution with user-specified minimum and maximum values, and retains the positive definite ones for the computation of \rho_{\Delta}
.
The obtained vector of \rho_{\Delta}
values can subsequently be used to examine (i) the impact of different assumptions regarding the correlations between the counterfactuals on the results (see also plot Causal-Inference ContCont
), and (ii) the extent to which proponents of the causal-inference and meta-analytic frameworks will reach the same conclusion with respect to the appropriateness of the candidate surrogate at hand.
The function ICA.Sample.ContCont
also generates output that is useful to examine the plausibility of finding a good surrogate endpoint (see GoodSurr
in the Value section below). For details, see Alonso et al. (submitted).
Notes
A single \rho_{\Delta}
value is obtained when all correlations in the function call are scalars.
An object of class ICA.ContCont
with components,
Total.Num.Matrices |
An object of class |
Pos.Def |
A |
ICA |
A scalar or vector that contains the individual causal association (ICA; |
GoodSurr |
A |
Wim Van der Elst, Ariel Alonso, & Geert Molenberghs
Alonso, A., Van der Elst, W., Molenberghs, G., Buyse, M., & Burzykowski, T. (submitted). On the relationship between the causal-inference and meta-analytic paradigms for the validation of surrogate markers.
MICA.ContCont
, ICA.ContCont
, Single.Trial.RE.AA
,
plot Causal-Inference ContCont
# Generate the vector of ICA values when rho_T0S0=rho_T1S1=.95,
# sigma_T0T0=90, sigma_T1T1=100,sigma_ S0S0=10, sigma_S1S1=15, and
# min=-1 max=1 is considered for the correlations
# between the counterfactuals:
SurICA2 <- ICA.Sample.ContCont(T0S0=.95, T1S1=.95, T0T0=90, T1T1=100, S0S0=10,
S1S1=15, M=5000)
# Examine and plot the vector of generated ICA values:
summary(SurICA2)
plot(SurICA2)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.