gravatar for Kevin Blighe

2 hours ago by

Republic of Ireland

Hi,

You just need to construct a standard survival regression model of the form (in R):

library(survival)
Surv(Months, DFS) ~ gene1 + gene2 + gene3

NB - 'DFS' = disease-free survival

I have a tutorial here: Survival analysis with gene expression

###

For a quick example, try something like:

library(survminer)
ggsurvplot(survfit(Surv(Months, DFS) ~ gene1 + gene2 + gene3, data = survplotdata),
    data = survplotdata,
    risk.table = TRUE,
    pval = TRUE,
    break.time.by = 500,
    ggtheme = theme_minimal(),
    risk.table.y.text.col = TRUE,
    risk.table.y.text = FALSE)

Kevin



Source link