Skip to contents

Proxy or delayed-outcome CMR assignment

Source: R/cmr_proxy.R
cmr_proxy.Rd

Estimate a proxy-outcome rectangle, widen it using the bridge radius zeta, and return the CMR treatment share for the primary-outcome design.

Usage

cmr_proxy(
  proxy_y,
  d,
  zeta,
  alpha = 0.05,
  method = c("auto", "bounded", "bernoulli", "maurer_pontil", "mp", "bernoulli_exact",
    "martinez_taboada_ramdas", "mtr"),
  beta = NULL,
  correction = c("bonferroni", "sidak_arms"),
  normalize = FALSE,
  lower = NULL,
  upper = NULL,
  na.rm = TRUE,
  tol = 1e-11
)

cmr_delayed_outcome(
  proxy_y,
  d,
  zeta,
  alpha = 0.05,
  method = c("auto", "bounded", "bernoulli", "maurer_pontil", "mp", "bernoulli_exact",
    "martinez_taboada_ramdas", "mtr"),
  beta = NULL,
  correction = c("bonferroni", "sidak_arms"),
  normalize = FALSE,
  lower = NULL,
  upper = NULL,
  na.rm = TRUE,
  tol = 1e-11
)

Arguments

proxy_y

Pilot proxy or delayed primary outcomes.

d

Pilot treatment indicator; treatment is 1 and control is 0.

zeta

Nonnegative standard-deviation bridge radius. Provide a scalar shared across arms or a treatment/control pair.

alpha

Target joint error level.

method

Confidence-set method. "auto" chooses exact Bernoulli bounds for 0/1 outcomes and bounded Maurer–Pontil bounds otherwise.

beta

Optional endpoint error allocation. If NULL, error is split according to correction.

correction

Endpoint error correction, either "bonferroni" or "sidak_arms".

normalize

If TRUE, normalize bounded proxy outcomes to [0, 1] before computing variances. For guarantee-bearing bounded CMR on a non-unit scale, provide known lower and upper bounds.

lower, upper

Optional lower and upper outcome bounds used when normalize = TRUE. If either is omitted, the pilot minimum and/or maximum is used with a warning; that data-dependent normalization is exploratory and does not carry the finite-sample bounded-outcome CMR guarantee.

na.rm

If TRUE, drop rows with missing proxy_y or d.

tol

Numerical tolerance for exact Bernoulli bound inversion.

Value

A list of class cmr_proxy and cmr_two_arm with treatment share pi, CMR certificate U_CMR, widened confidence rectangle, pilot summaries, zeta, bridge diagnostics, endpoint error allocation, and method metadata.

See also

Other CMR rules: binary_rectangle_corners(), cmr_multiarm(), cmr_multiarm_from_rectangle(), cmr_multiple_outcomes(), cmr_stratified(), cmr_stratified_from_rectangle(), cmr_two_arm(), cmr_unbounded(), cmr_unbounded_from_rectangle(), print.cmr_two_arm()

Examples

set.seed(12)
d <- rep(c(1, 0), each = 30)
proxy_y <- c(rbeta(30, 2, 6), rbeta(30, 4, 4))
cmr_proxy(proxy_y, d, zeta = 0.05)
#> <cmr_proxy>
#>   pi: 0.5
#>   U_CMR: 0.25
#>   method: proxy_bounded
#>   n: 60