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)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 \(\rho_{\Delta}\).
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 \(\rho_{\Delta}\).
A scalar that specifies the variance of the true endpoint in the control treatment condition that should be considered in the computation of \(\rho_{\Delta}\). Default 1.
A scalar that specifies the variance of the true endpoint in the experimental treatment condition that should be considered in the computation of \(\rho_{\Delta}\). Default 1.
A scalar that specifies the variance of the surrogate endpoint in the control treatment condition that should be considered in the computation of \(\rho_{\Delta}\). Default 1.
A scalar that specifies the variance of the surrogate endpoint in the experimental treatment condition that should be considered in the computation of \(\rho_{\Delta}\). Default 1.
A scalar or vector that contains the correlation(s) between the counterfactuals T0 and T1 that should be considered in the computation of \(\rho_{\Delta}\). Default seq(-1, 1, by=.001).
A scalar or vector that contains the correlation(s) between the counterfactuals T0 and S1 that should be considered in the computation of \(\rho_{\Delta}\). Default seq(-1, 1, by=.001).
A scalar or vector that contains the correlation(s) between the counterfactuals T1 and S0 that should be considered in the computation of \(\rho_{\Delta}\). Default seq(-1, 1, by=.001).
A scalar or vector that contains the correlation(s) between the counterfactuals S0 and S1 that should be considered in the computation of \(\rho_{\Delta}\). Default seq(-1, 1, by=.001).
The number of runs that should be conducted. Default 50000.
An object of class ICA.ContCont with components,
An object of class numeric that contains the total number of matrices that can be formed as based on the user-specified correlations in the function call.
A data.frame that contains the positive definite matrices that can be formed based on the user-specified correlations. These matrices are used to compute the vector of the \(\rho_{\Delta}\) values.
A scalar or vector that contains the individual causal association (ICA; \(\rho_{\Delta}\)) value(s).
A data.frame that contains the ICA (\(\rho_{\Delta}\)), \(\sigma_{\Delta_{T}}\), and \(\delta\).
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.
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
# NOT RUN {
# 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)
# }
Run the code above in your browser using DataCamp Workspace