Given a matrix of limma voom normalized genes, I want to do Survival Analysis of a particular gene X.
The design matrix applied on voom was achieved regarding recurrence or no recurrence. Yet I have three disease groups A, B and C. When I filtered patients only with disease C and did DGE over this patients.
This gene X appeared overexpressed in patients that showed cancer recurrence in group C.
Now I need to convert the matrix into a matrix underexpressed / overexpressed / normal so that I can do survival analysis to show indeed that for patients with higher expression of this gene the cancer returns.
How should I convert this matrix into this nominal shape?
Possibilities:
I understand that limma voom results in negative, neutral (near to 0) and positive, so negative means underexpressed, near to zero normal and positive mean overexpressed for all genes?
Inside the same gene X do histogram, analyze bi modality and set threshold negative/ neutral /positive by regression?
Get the log2 foldchange and log(adj.p values) of each gene and define threshold in the same way vulcano plot is designed?