How to make a survival plot for a gene between High expression and low expression samples?
0
1
Entering edit mode
5.6 years ago
Biologist ▴ 290

Hi,

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?

RNA-Seq R survival plot tcga • 3.2k views
ADD COMMENT
1
Entering edit mode

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 REPLY
0
Entering edit mode

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 REPLY
1
Entering edit mode

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 cbioportal.org

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

ADD REPLY
1
Entering edit mode

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 REPLY

Login before adding your answer.

Traffic: 1564 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6