interesting ! How do I read this kind of figure ? what is the x-axis ? what is 'TMB' ? ...
• 0 views
•
link
Hi there,
I am trying to find a program or code that will let me use VCF files generated from RNA-seq data (non-model organism) to create a mutational landscape plot (or rainfall plot). Any program recommendations?
Thanks.
Check out the pymaf.MafFrame.plot_oncoplot method I wrote:
from fuc import common, pymaf
common.load_dataset('tcga-laml')
f = '~/fuc-data/tcga-laml/tcga_laml.maf.gz'
mf = pymaf.MafFrame.from_file(f)
mf.plot_oncoplot()
You can convert VCF to MAF with the pymaf.MafFrame.from_vcf method:
from fuc import pymaf
mf = pymaf.MafFrame.from_vcf('annotated.vcf')
interesting ! How do I read this kind of figure ? what is the x-axis ? what is 'TMB' ? ...
Just wanted to let you know that you can now use the fuc package to draw a rainfall plot directly from VCF with the pyvcf.VcfFrame.plot_rainfall method:
import matplotlib.pyplot as plt
import seaborn as sns
from fuc import common, pyvcf
common.load_dataset('brca')
vcf_file = '~/fuc-data/brca/brca.vcf'
vf = pyvcf.VcfFrame.from_file(vcf_file)
vf.plot_rainfall('TCGA-A8-A08B',
figsize=(14, 7),
palette=sns.color_palette('Set2')[:6])
plt.tight_layout()
Log in to answer this question.