ICA.BinBin.Grid.Full {Surrogate}R Documentation

Assess surrogacy in the causal-inference single-trial setting in the binary-binary case when monotonicity for SS and TT is assumed using the full grid-based approach

Description

The function ICA.BinBin.Grid.Full quantifies surrogacy in the single-trial causal-inference framework (individual causal association and causal concordance) when both the surrogate and the true endpoints are binary outcomes. This method provides an alternative for ICA.BinBin and ICA.BinBin.Grid.Sample. It uses an alternative strategy to identify plausible values for π\pi. See Details below.

Usage

ICA.BinBin.Grid.Full(pi1_1_, pi1_0_, pi_1_1, pi_1_0, pi0_1_, pi_0_1, 
Monotonicity=c("General"), pi_1001=seq(0, 1, by=.02), 
pi_1110=seq(0, 1, by=.02), pi_1101=seq(0, 1, by=.02),
pi_1011=seq(0, 1, by=.02), pi_1111=seq(0, 1, by=.02), 
pi_0110=seq(0, 1, by=.02), pi_0011=seq(0, 1, by=.02), 
pi_0111=seq(0, 1, by=.02), pi_1100=seq(0, 1, by=.02), 
Seed=sample(1:100000, size=1))

Arguments

pi1_1_

A scalar that contains P(T=1,S=1Z=0)P(T=1,S=1|Z=0), i.e., the proability that S=T=1S=T=1 when under treatment Z=0Z=0.

pi1_0_

A scalar that contains P(T=1,S=0Z=0)P(T=1,S=0|Z=0).

pi_1_1

A scalar that contains P(T=1,S=1Z=1)P(T=1,S=1|Z=1).

pi_1_0

A scalar that contains P(T=1,S=0Z=1)P(T=1,S=0|Z=1).

pi0_1_

A scalar that contains P(T=0,S=1Z=0)P(T=0,S=1|Z=0).

pi_0_1

A scalar that contains P(T=0,S=1Z=1)P(T=0,S=1|Z=1).

Monotonicity

Specifies which assumptions regarding monotonicity should be made: Monotonicity=c("General"), Monotonicity=c("No"), Monotonicity=c("True.Endp"), Monotonicity=c("Surr.Endp"), or Monotonicity=c("Surr.True.Endp"). When a general analysis is requested (using Monotonicity=c("General") in the function call), all settings are considered (no monotonicity, monotonicity for SS alone, for TT alone, and for both for SS and TT. Default Monotonicity=c("General").

pi_1001

A vector that specifies the grid of values that should be considered for πpi1001\pi_{pi_1001}. Default pi_1001=seq(0, 1, by=.02).

pi_1110

A vector that specifies the grid of values that should be considered for πpi1110\pi_{pi_1110}. Default pi_1110=seq(0, 1, by=.02).

pi_1101

A vector that specifies the grid of values that should be considered for πpi1101\pi_{pi_1101}. Default pi_1101=seq(0, 1, by=.02).

pi_1011

A vector that specifies the grid of values that should be considered for πpi1011\pi_{pi_1011}. Default pi_1011=seq(0, 1, by=.02).

pi_1111

A vector that specifies the grid of values that should be considered for πpi1111\pi_{pi_1111}. Default pi_1111=seq(0, 1, by=.02).

pi_0110

A vector that specifies the grid of values that should be considered for πpi0110\pi_{pi_0110}. Default pi_0110=seq(0, 1, by=.02).

pi_0011

A vector that specifies the grid of values that should be considered for πpi0011\pi_{pi_0011}. Default pi_0011=seq(0, 1, by=.02).

pi_0111

A vector that specifies the grid of values that should be considered for πpi0111\pi_{pi_0111}. Default pi_0111=seq(0, 1, by=.02).

pi_1100

A vector that specifies the grid of values that should be considered for πpi1100\pi_{pi_1100}. Default pi_1100=seq(0, 1, by=.02).

Seed

The seed to be used to generate πr\pi_r. Default Seed=sample(1:100000, size=1).

Details

In the continuous normal setting, surroagacy can be assessed by studying the association between the individual causal effects on SS and TT (see ICA.ContCont). In that setting, the Pearson correlation is the obvious measure of association.

When SS and TT are binary endpoints, multiple alternatives exist. Alonso et al. (2014) proposed the individual causal association (ICA; RH2R_{H}^{2}), which captures the association between the individual causal effects of the treatment on SS (ΔS\Delta_S) and TT (ΔT\Delta_T) using information-theoretic principles.

The function ICA.BinBin.Grid.Full computes RH2R_{H}^{2} using a grid-based approach where all possible combinations of the specified grids for the parameters that are allowed that are allowed to vary freely are considered. When it is not assumed that monotonicity holds for both SS and TT, the computationally less demanding algorithm ICA.BinBin.Grid.Sample may be preferred.

Value

An object of class ICA.BinBin with components,

Pi.Vectors

An object of class data.frame that contains the valid π\pi vectors.

R2_H

The vector of the RH2R_H^2 values.

Theta_T

The vector of odds ratios for TT.

Theta_S

The vector of odds ratios for SS.

H_Delta_T

The vector of the entropies of ΔT\Delta_T.

Author(s)

Wim Van der Elst, Paul Meyvisch, Ariel Alonso & Geert Molenberghs

References

Alonso, A., Van der Elst, W., & Molenberghs, G. (2014). Validation of surrogate endpoints: the binary-binary setting from a causal inference perspective.

Buyse, M., Burzykowski, T., Aloso, A., & Molenberghs, G. (2014). Direct estimation of joint counterfactual probabilities, with application to surrogate marker validation.

See Also

ICA.ContCont, MICA.ContCont, ICA.BinBin, ICA.BinBin.Grid.Sample

Examples

## Not run:  # time consuming code part
# Compute R2_H given the marginals, 
# assuming monotonicity for S and T and grids
# pi_0111=seq(0, 1, by=.001) and 
# pi_1100=seq(0, 1, by=.001)
ICA <- ICA.BinBin.Grid.Full(pi1_1_=0.2619048, pi1_0_=0.2857143, pi_1_1=0.6372549, 
pi_1_0=0.07843137, pi0_1_=0.1349206, pi_0_1=0.127451,  
pi_0111=seq(0, 1, by=.01), pi_1100=seq(0, 1, by=.01), Seed=1)

# obtain plot of R2_H
plot(ICA, R2_H=TRUE)

## End(Not run)

[Package Surrogate version 3.3.0 Index]