Heads up:
We’re moving the GATK website, docs and forum to a new platform. Read the full story and breakdown of key changes on this blog.
If you happen to see a question you know the answer to, please do chime in and help your fellow community members. We encourage our fourm members to be more involved, jump in and help out your fellow researchers with their questions. GATK forum is a community forum and helping each other with using GATK tools and research is the cornerstone of our success as a genomics research community.We appreciate your help!

Test-drive the GATK tools and Best Practices pipelines on Terra

Check out this blog post to learn how you can get started with GATK and try out the pipelines in preconfigured workspaces (with a user-friendly interface!) without having to install anything.

PlotDenoisedCopyRatios (v4.0.1.2) error "need finite 'ylim' values"

stefan_stojanovicstefan_stojanovic Member
edited February 4 in Ask the GATK team
Hello, I'm getting the following error from PlotDenoisedCopyRatios tool, version

> 2019-01-31T09:56:23.782315172Z Stderr: Error in plot.window(...) : need finite 'ylim' values
> 2019-01-31T09:56:23.782317446Z Calls: source ... SetUpPlot -> plot -> plot.default -> localWindow -> plot.window

I assume this is to do with the inputs denoised and standardized copy ratios, which are produced by DenoiseReadCounts (v4.0.1.2 as well). Below is chunk of these files. There are "Infinity" values in LOG2_COPY_RATIOS column so I assume this is causing PlotDenoisedCopyRatios to fail, but maybe I am wrong. There is also number of NaN values in this column.

chr22 50776501 50776750 0.336554
chr22 50776751 50777000 Infinity
chr22 50777751 50778000 0.250825

chr22 50807751 50808000 NaN
chr22 50808001 50808250 NaN
chr22 50808251 50808500 NaN

I ran both tools with default parameter values, here are command lines.

gatk --java-options -Xmx13000m PlotDenoisedCopyRatios --output . --output-prefix cnv_NORMAL --denoised-copy-ratios cnv_NORMAL.denoisedCR.tsv --sequence-dictionary Homo_sapiens_assembly38.dict --standardized-copy-ratios cnv_NORMAL.standardizedCR.tsv

gatk --java-options -Xmx13000m DenoiseReadCounts --input cnv_NORMAL.counts.hdf5 --denoised-copy-ratios cnv_NORMAL.denoisedCR.tsv --standardized-copy-ratios cnv_NORMAL.standardizedCR.tsv --count-panel-of-normals PoN_4.0_WGS_for_public.pon.hdf5

My question is how do I solve this issue. Do I need to change some parameters or is this issue with input files?
Post edited by bshifaw on

Best Answers


  • Hi @bshifaw, I managed to find a workaround which is to use larger interval bins to collect read counts (bin_length=1000). I am still getting the same error with smaller bins (bin_length=250) which are used in best practice FullSomaticPipeline example (github.com/gatk-workflows/gatk4-somatic-with-preprocessing). It is possible my BAM files are not suitable for CNV calling (they are around 2GB each), but DenoiseReadCounts finished without errors so I'm not sure. Hope this helps your team.
Sign In or Register to comment.