Hi All, So I have been doing some survival analysis using the extremely helpful tutorial Survival analysis of TCGA patients integrating gene expression (RNASeq) data, but being new to R I was wondering if anyone could help me with some changes that I'd like to make. Instead of just choosing one gene to check the p-value and plot, I instead want to write a code that would divide the genes into up and down-regulated, and create a loop or something equivalent that does survival analysis on and extracts the p-value for every gene and then lists the top 20 up-regulated and top 20 down-regulated genes after it is done. I've tried a lot of code but I'm afraid that I'm really not good at this and any help to achieve my goal would be hugely appreciated! Thanks so much