Log2 Scale Deeptools

Greetings,

I hope everyone is doing well.

So I used computeMatrix from deeptools to calculate the coverage of certain samples. The bigwig signal is in BPM (generated by bamCoverage). I want to transform the values in the matrix to a log2 scale to get the enrichment profile or heatmap in a log2 scale rather than the original BPM.

I’ve tried importing the matrix generated by compute matrix into Rstudio, converted the values into log2 scale, and then rezipped the file and tried using plotHeatmap with no success.

 Traceback (most recent call last):
  File "/home/bol7asan/anaconda3/envs/pyngs/bin/plotHeatmap", line 12, in <module>
    main(args)
  File "/home/bol7asan/anaconda3/envs/pyngs/lib/python3.7/site-packages/deeptools/plotHeatmap.py", line 815, in main
    hm.read_matrix_file(matrix_file)
  File "/home/bol7asan/anaconda3/envs/pyngs/lib/python3.7/site-packages/deeptools/heatmapper.py", line 779, in read_matrix_file
    chrom, start, end, name, score, strand = region[0:6]
ValueError: not enough values to unpack (expected 6, got 1)

I’ve also attached a picture of the matrix after modification, I just changed the values after the strand column.

enter image description here

I would appreciate it if someone can help with this.

Read more here: Source link