Skip to contents

Performs CARE analysis with bootstrapping for robust estimation.

Usage

CARE2_boot(
  gamma.exp_sel,
  gamma.out_sel,
  se.exp_sel,
  se.out_sel,
  nx,
  ny,
  nrep = 1000,
  algorithm = "Lasso",
  random_start = 0,
  biascorrect = "no",
  etamean = 0.5,
  pthr = 5e-05
)

Arguments

gamma.exp_sel

Vector of genetic associations with the exposure

gamma.out_sel

Vector of genetic associations with the outcome

se.exp_sel

Vector of standard errors for the exposure associations

se.out_sel

Vector of standard errors for the outcome associations

nx

Sample size for exposure

ny

Sample size for outcome

nrep

Number of bootstrap repetitions

algorithm

Algorithm to use: "Lasso" or "CD" (coordinate descent)

random_start

Number of random starting points

biascorrect

Bias correction method: "no", "direct", "rerand", "rerand2", or "rerand3"

etamean

Mean of the random variable for re-randomization

pthr

P-value threshold for selecting significant SNPs

Value

A list containing CARE analysis results

Examples

if (FALSE) { # \dontrun{
results <- CARE2_boot(beta_exp, beta_out, se_exp, se_out, 10000, 10000)
} # }