Dear Biostars, Hi
I have about 8000 transcripts that has hit with Transcription Factors (TFs). Among them 69 are up-regulated in females (3 biological replication) and 19 up-regulated in males (with 3 biological replication). (using DESeq2, FDR 0.05)
Now I intend to draw a heat map based on my TMM-normalized expression values (3 for males and 3 for females for each transcripts) but I want to just show the DEG in it (69+19) not all 8000 transcripts; something same as Figure8b in here.
1) How can I do it? is it just by selecting those 88 transcripts TMM value and used them as input in R heatmap script?
2) As you can see in the Figure8, there is just 2 columns. but as I have 3 BR for each conditions, I have 6 columns. How I can convert them into 2 columns ?
Can I calculate the mean of each 3 TMM of each transcripts and use the mean number as input for each transcripts ?
~ Thank you in advance