Question: How to make a survival plot for a gene between High expression and low expression samples?
gravatar for Biologist
2.2 years ago by
Biologist190 wrote:


I have downloaded TCGA data for lung cancer. It is counts data. I want to make survival plot for a specific gene between high expression and low expression samples. I want to make both Overall survival and Disease free survival plots. It should look something like this Kaplan-Meier survival analysis of OS (P < 0.001, log-rank) and DFS (P < 0.001, log-rank) rate in 144 patients based on the expression level of SNHG20 expression

How to divide the samples into high and low expression samples based on expression of single gene? I have counts data. Is there any cutoff for that?

plot survival rna-seq tcga R • 1.4k views
ADD COMMENTlink written 2.2 years ago by Biologist190

You could plot the data and see how it spreads and use this do define the cutoff (e.g. if you notice there is clear separation into 2 groups). Another approach might be to use z-score and take +-1-2 score as high and low..

ADD REPLYlink written 2.2 years ago by roy.granit830

Based on z-score, so +1 and +2 are high, -1 and -2 are low? And for this survival do I also need to consider normal samples? or only tumor samples?

ADD REPLYlink written 2.2 years ago by Biologist190

Yes, you can try a cutoff on 1 or 2 and see how it looks like ( for a two sided distribution the critical Z score values when using a 95% confidence level are -1.96 and +1.96). This task could be easily accomplished using

I would plot cancer and normal on different charts (not sure there is point in showing a survival chart for normal individuals).

ADD REPLYlink modified 2.2 years ago • written 2.2 years ago by roy.granit830

Ok. Thankyou. then I will convert counts data to logCPM and then to Z-score. From that samples with values less than -1.96 are low and samples with more than +1.96 are high and then make a survival plot for that.

ADD REPLYlink written 2.2 years ago by Biologist190
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: 1956 users visited in the last hour