StepTwo {SMAHP}R Documentation

Select M-X pairs via SIS

Description

An auxiliary function which selects M-X paris using SIS.

Usage

StepTwo(
  X,
  M,
  time,
  status,
  X_sel_Y_s1,
  M_X_s1,
  M_sel_Y_s1,
  SIS_thres = "n/log(n)"
)

Arguments

X

An n by p matrix of exposures.

M

An n by p matrix of mediators.

time

A vector of survival time of samples.

status

A vector of status indicator: 0=alive, 1=dead.

X_sel_Y_s1

Outputs from StepOne: A vector of selected X for Y.

M_X_s1

Outputs from StepOne: A list of selected X for M.

M_sel_Y_s1

Outputs from StepOne: A vector of selected M for Y.

SIS_thres

SIS thresholds. Default is "n/log(n)". Other options include "n-1", "n/2log(n)", "2n/log(n)", "3n/log(n)".

Value

A data table with selected M, X pairs and related effect size.

Examples


data(example_dat)
surv_dat <- example_dat$surv_dat
res_step1 <- StepOne(X = example_dat$X, M = example_dat$M, time = surv_dat$time,
status = surv_dat$status, model_option = "MCP")
M_X_sel_s2 <- StepTwo(X = example_dat$X, M = example_dat$M, time = surv_dat$time,
status = surv_dat$status, X_sel_Y_s1 = res_step1$X_sel_Y_s1, M_X_s1 = res_step1$M_X_s1,
M_sel_Y_s1 = res_step1$M_sel_Y_s1)


[Package SMAHP version 0.0.5 Index]