Question: Integrating methylation in survival analysis
gravatar for Pin.Bioinf
5 months ago by
Pin.Bioinf230 wrote:

Hello, I am trying to do a survival analysis with TCGA data following this tutorial: Survival analysis of TCGA patients integrating gene expression (RNASeq) data , but I want to integrate methylation data instead of rnaseq data. What I am doing is selecting the gene of interest and averaging its CpG beta values per column: then, I create a vector for that gene where, for samples where avg beta values are >0.3 , the value is a 1, and a 0 if its <0.3. Then, my survival function is as follows:

s <- survfit(Surv(as.numeric(as.character(all_clin$new_death))[ind_clin],all_clin$death_event[ind_clin])~event_r[ind_clin])
s1 <- tryCatch(survdiff(Surv(as.numeric(as.character(all_clin$new_death))[ind_clin],all_clin$death_event[ind_clin])~event_r[ind_clin]), error = function(e) return(NA))

where event_r is a vector of 1s and 0s for the corresponding clinical samples in that gene of interest.

Do you think averaging the beta-values for the whole gene is a good approach to integrate methylation with survival? Any other suggestions? Thank you

rna-seq survival methylation R • 200 views
ADD COMMENTlink written 5 months ago by Pin.Bioinf230
gravatar for Kevin Blighe
5 months ago by
Kevin Blighe37k
Republic of Ireland
Kevin Blighe37k wrote:

I think that you should treat the methylation probes as separate entities.

ADD COMMENTlink written 5 months ago by Kevin Blighe37k

Okay, I will try that approach. Thank you!

ADD REPLYlink written 5 months ago by Pin.Bioinf230

But would averaging the beta values be recommended, can't we consider the continuous values for each probe?

ADD REPLYlink written 8 days ago by Researcher10

I replied to you in the other thread: A: Survival analysis with gene expression

ADD REPLYlink written 8 days ago by Kevin Blighe37k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2055 users visited in the last hour