Patients are split into two groups based on their immune cell abundance (median aundance value) and a regular survival analyis is carried out.
imsig_survival(exp, cli, time = "time", status = "status", r = 0.6)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.
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.
Column name of time-to-event parameter.
Column name of event (dead or alive) parameter.
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.
Hazard Ratio
# NOT RUN {
survival = imsig_survival (exp = example_data, cli = example_cli)
head(survival)
# }
Run the code above in your browser using DataLab