Hi,
I understand that bioconductor packages such as edgeR and DESeq require readcounts for further downstream analysis.I have FPKM values from cufflinks.Can I just use the values to produce a heatmap?
Also,I want to get the top expressed genes from normals vs Tumors.Any suggestions.
Best,
Ron
•
written
by
Ron
2 answers
If you have bam files, run cuffdiff. The cuffdiff output will have differentiallly expressed isoforms, genes, FPKM values for genes and isoforms etc.
The cuffdiff output can be read in to cummeRbund package which has many visualisations including heatmaps on FPKM values.
If you want to use edgeR or DESeq, use htseq-count to get the raw counts per gene/exon and feed that data into edger/DESeq.
Old question but this may be of help for those dealing with the same issue:
For clustering, you need to transform your data into a z-score. There is a Bioconductor package, zFPKM do this on FPKM reads.
Please open your own question rather than piggybacking on another if you have one. This helps us keep the site organized and allows other to find your question more easily.
@rhasanvandj,
Do a differential expression analysis, Set a cut-off to retrieve DE genes, Use a package for heatmap plotting. Please check this tuturial from Kevin Blighe.
Log in to answer this question.
-
MA plot from Differential Expression between One Sample vs another
written by Ron
-
CAGE seq & RNA-seq
written by Ron
-
Remove batch effect for differential expression
written by Ron
-
RNA-Seq analysis for individual genes
written by sumithra.das •
-
Differential expression analysis of transcriptome data
written by vijay
-
GSEA Preranked analysis on RNA-SEQ
written by Ron
-
Rank Genes By FPKM for only One Sample
written by Ron
-
FPKM for Normal vs Tumor
written by Ron
-
FPKM values / Data analysis
written by Biogeek
-
Unsupervised Clustering/Heatmap Of Paired Sample Rna-Seq Expression
written by bstephenwhite •
