Question: How to make a survival plot for a gene between High expression and low expression samples?
0
gravatar for Biologist
8 months ago by
Biologist150
Biologist150 wrote:

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?

plot survival rna-seq tcga R • 406 views
ADD COMMENTlink written 8 months ago by Biologist150
1

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 8 months ago by roy.granit790

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 8 months ago by Biologist150
1

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 REPLYlink modified 8 months ago • written 8 months ago by roy.granit790
1

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 8 months ago by Biologist150
Please log in to add an answer.

Help
Access

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