Avatar

timh

  • Total activity 10
  • Last activity
  • Member since
  • Following 0 users
  • Followed by 0 users
  • Votes 0
  • Subscriptions 3

Activity overview

Latest activity by timh
  • Avatar

    timh commented,

    Thanks Genevieve. All right, sounds like I will automate something based on annotation stats. The --max-attempts function is usefull but the output file, which then contains multiple models, requir...

  • Avatar

    timh created a post,

    VQSR: how is insufficient variance inferred for annotations?

    Answered

    Hi there. Some of the annotations I use in VQSR do not always vary significantly and hence I get the following error: ```A USER ERROR has occurred: Positive training model failed to converge. One o...

  • Avatar

    timh commented,

    Hi. No I didn't change the R script.

  • Avatar

    timh commented,

    update: this appears to be the case for only a number of my datasets, works fine for others. So I am guessing the problem is not VQSR but my data, will have a closer look. Thanks

  • Avatar

    timh commented,

    Here you are: [exact command works without "-an AS_MQRankSum"] java -jar ~/Programs/gatk-4.1.4.1/gatk-package-4.1.4.1-local.jar VariantRecalibrator -R ref.fa -V test.4.1.4.1.vcf.gz -AS --resource:r...

  • Avatar

    timh created a post,

    postive AS_MQRankSum and VQSR

    Hi. I have a technical question about AS_MQRankSum. As is mentioned in the info page in the tool section, positive scores for this measure are not really informative. I wonder how VQSR and ApplyVQS...

  • Avatar

    timh created a post,

    AS_MQRanksum craches R script VQSR

    Hi. I am running GATK 4.1.4.1. and when I am running allele specific VQSR and include "-an AS_MQRankSum" I get the following error: 13:02:54.929 INFO VariantRecalibrator - Executing: Rscript /home/...