imsig_survival {imsig}

Survival analysis based on relative abundance of immune infiltration estimated by ImSig

Description

Patients are split into two groups based on their immune cell abundance (median aundance value) and a regular survival analyis is carried out.

Usage

imsig_survival(exp, cli, time = "time", status = "status", r = 0.6)

Arguments

exp

Dataframe of transcriptomic data (natural scale) containing genes as rows and samples as columns. Note: Gene names should be set as row names and duplicates are not allowed. Missing values are not allowed within the expression matrix. Check example- head(example_data): example_data.

cli

Clinical metadata containting the event data (dead or alive) and time to event data. Samples names should be in rownames and same as that in the expression file. Check head() of example_cli for an example clinical data.

time

Column name of time-to-event parameter.

status

Column name of event (dead or alive) parameter.

r

Use a value between 0 and 1. Default is 0.6. This is a user defined correlation cut-off to perform feature selection (feature_select). Feature selection aids to enrich the prediction of relative abundance of immune cells by filtering off poorly correlated ImSig genes. To get an idea of what cut-off to use check the results of (gene_stat) and choose a cut-off that displays high median correlation and maintains a high proportion of genes after feature selection.

Value

Hazard Ratio

See Also

feature_select, example_data, example_cli

Examples

survival = imsig_survival (exp = example_data, cli = example_cli)
head(survival)