SIS_combo(X, M, C, time, status, p_adjust_option ="BH", thres =0.05)
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.
time: A vector of survival time of samples.
status: A vector of status indicator: 0=alive, 1=dead.
p_adjust_option: The method for multiple correction. Option include q-value, holm, hochberg, hommel, bonferroni, BH, BY, and fdr. Default is BH.
thres: Threshold for determining significance.
Returns
A list which includes the final p-value matrix (p_final_matrix), adjusted p-value matrix (p_adjusted_matrix) and mediation-exposure matrix (p_med_matrix).