StepThree {SMAHP} | R Documentation |
Post-hoc model
Description
An auxiliary function which generates p-value based on outcome and mediation model.
Usage
StepThree(X, M, C, time, status, X_sel_Y_s1, M_X_sel_s2)
Arguments
X |
An n by p matrix of exposures. |
M |
An n by p matrix of mediators. |
C |
An n by p matrix of covariates. If there are no covariates, set C = NULL. |
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_sel_s2 |
Outputs from StepTwo: A data table with selected M, X pairs and related effect size. |
Value
A list with the following components:
p_beta_m |
p-values generated from outcome model |
p_alpha_x |
p-values generated from mediation modell |
outcome_model |
coefficient estimates from outcome model |
med_results |
coefficient estimates from mediation model |
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)
res_step3 <- StepThree(X = example_dat$X, M = example_dat$M, C = example_dat$C,
time = surv_dat$time, status = surv_dat$status, X_sel_Y_s1 = res_step1$X_sel_Y_s1,
M_X_sel_s2 = M_X_sel_s2)
[Package SMAHP version 0.0.5 Index]