Best way to draw heatmap using Gfold deferentially expressed no replicate data
0
0
Entering edit mode
3.4 years ago
lookamrita • 0

Hi, I have two sample (Control and Test) RNA seq data (Ref: Human) after analyzing the files using Star and GATK pipeline , I Have used Gfold to check the differential expression analysis. Now I have Gfold output file of differential expression analysis. I want to draw heatmap using it. There are 475 genes which are having more than 2 log2fc, I want that genes to use in heatmap from the Gfold output file. What is the best way to draw heatmap using it? New in RNA-seq. Need help. Thanks in advance

RNA-Seq next-gen R • 914 views
0
Entering edit mode

Please see as a starting point:

Heatmap based with FPKM values

There are plenty of posts here on Biostars. Please use the search function.

0
Entering edit mode

Even though I have log2fc value, shall I use FPKM only. Please see the example values and columns in my file. This is how it looks like:

Gene    GFOLD(2)    E-FDR   log2fdc 1stRPKM 2ndRPKM
7SK 3.63629 1   2.7294  0.0064312   0.0753654
A2ML1   8.08695 1   7.09134 0.00184317  0.47231
AADAC   2.57739 1   2.49758 0.108541    0.610544
AATK    2.68001 1   2.63272 0.0563924   0.34059
ABCA6   4.80114 1   4.7294  0.0358213   0.941189

0
Entering edit mode

You can use whatever you think is best. I linked this one for example code that you can adapt to feed in your FCs.
{"extraction_info": {"found_math": true, "script_math_tex": 0, "script_math_asciimath": 0, "math_annotations": 0, "math_alttext": 0, "mathml": 0, "mathjax_tag": 0, "mathjax_inline_tex": 0, "mathjax_display_tex": 0, "mathjax_asciimath": 1, "img_math": 0, "codecogs_latex": 0, "wp_latex": 0, "mimetex.cgi": 0, "/images/math/codecogs": 0, "mathtex.cgi": 0, "katex": 0, "math-container": 0, "wp-katex-eq": 0, "align": 0, "equation": 0, "x-ck12": 0, "texerror": 0, "math_score": 0.18887880444526672, "perplexity": 6155.071765790617}, "config": {"markdown_headings": false, "markdown_code": true, "boilerplate_config": {"ratio_threshold": 0.18, "absolute_threshold": 10, "end_threshold": 15, "enable": true}, "remove_buttons": true, "remove_image_figures": true, "remove_link_clusters": true, "table_config": {"min_rows": 2, "min_cols": 3, "format": "plain"}, "remove_chinese": true, "remove_edit_buttons": true, "extract_latex": true}, "warc_path": "s3://commoncrawl/crawl-data/CC-MAIN-2022-40/segments/1664030337490.6/warc/CC-MAIN-20221004085909-20221004115909-00020.warc.gz"}