Question: plotTracks - empty plot and custom genomes
0
gravatar for agata88
6 days ago by
agata88790
Poland
agata88790 wrote:

Hi all!

I am having trouble with Gvis library and plotTracks function. Here is my code:

ideogram <- IdeogramTrack(chromosome = "chr4",genome="hg38")
axis <- GenomeAxisTrack()
plotTracks(list(ideogram, axis),from=189233986,to=189234447)

Unfortunately, it is not plotting anything. Do you have any idea why?

My second question include custom genome. If new genome is not available through UCSC (or is modified), how can I use it in this example? When I use my custom hg38_clean.fa file I have an error:

ideogram <- IdeogramTrack(chromosome = "chr4",genome="../hg38_clean.fa")
axis <- GenomeAxisTrack()
plotTracks(list(ideogram, axis),from=189233986,to=189234447)

    'hg38_clean.fa' is not a valid UCSC genome.

Many thanks for any suggestions,

Best, Agata

gviz R • 65 views
ADD COMMENTlink modified 6 days ago by benformatics1.1k • written 6 days ago by agata88790
0
gravatar for benformatics
6 days ago by
benformatics1.1k
ETH Zurich
benformatics1.1k wrote:

When I run your first block of code I am able to produce the following plot as expected (R version 3.6.0; Gviz version 1.28.0): gviz plot

An earlier R install gave me an error: There doesn't seem to be any cytoband data available for genome 'hg38' at UCSC or the service is temporarily down. Trying to fetch the chromosome length data. So perhaps you need to update your packages.

In regards to your second question see the following posts:

In general what you want to do should be possible... but acquiring the Ideogram information is going to be relatively difficult or you will have to generate it yourself.

ADD COMMENTlink written 6 days ago by benformatics1.1k

I am also running R version 3.6.0 and Gviz 1.28.0, so I don't know why this is not working for me... For second question, thank you for tips! I will dig a little more. Best, Agata

ADD REPLYlink written 6 days ago by agata88790
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1105 users visited in the last hour