CONICS icon indicating copy to clipboard operation
CONICS copied to clipboard

How can we use tpm values provided in GEO datasets for CNV or tumor cell prediction

Open biobug16 opened this issue 4 years ago • 2 comments

Hi @diazlab and @soerenmueller I want to predict tumor cells in Suva’s 2014 IDH wt Glioma dataset, for which I only have log2tpm values. How can we convert log2tpm to log2cpm or is there a way to use log2tpm itself with CONICSmat?

Thanks

biobug16 avatar Jan 30 '20 14:01 biobug16

Hi @biobug16

here is some example code for the Patel et al dataset

a=read.table("Glioblastoma_expressed_genes.txt",sep="\t",header=T,row.names=1,check.names = F) a=((2^a)-1) a=log2(a/10+1)

soerenmueller avatar Feb 04 '20 17:02 soerenmueller

Hi @soerenmueller Thank you for your reply. Can you please explain your answer a bit. I guess after reading, you first calculated the antilog for the log2tpm by using 2^a and then deducted 1 from it, followed you dividing it by 10 and then again you did log2 of it.

May be I am missing something, but are you sure it will be equivalent to log2CPM?

Correct me if I am getting it wrong.

Thanks again!

biobug16 avatar Feb 04 '20 19:02 biobug16