[Date Prev][Date Next] [Thread Prev][Thread Next] [Date Index] [Thread Index]

Re: Test suite errors in StructuralVariantAnnotation



My guess is that stricter input file header validation in the r-bioc-variantannotation/r-bioc-rsamtools package dependencies is making the test cases unhappy. It's a failure of the test harness, not the actual test case functionality. I'll look into it.

Cheers
Daniel

On Wed, Nov 23, 2022 at 7:39 PM Andreas Tille <andreas@an3as.eu> wrote:
Hi Daniel,

when upgrading the Debian packaged BioConductor packages I stumbled upon
test suite errors in StructuralVariantAnnotation which are starting
with:

[ FAIL 30 | WARN 4 | SKIP 0 | PASS 85 ]
══ Failed tests ════════════════════════════════════════════════════════════════
── Error ('test-extensions-VCF.R:28'): INFO column import ──────────────────────
Error in `h(simpleError(msg, call))`: error in evaluating the argument 'x' in selecting a method for function 'breakpointRanges': no 'header' line "#CHROM POS ID..."?
Backtrace:
     ▆
  1. ├─StructuralVariantAnnotation::breakpointRanges(.testrecord(c("chr10\t2991435\tINV\tN\t<INV>\t.\tLowQual\tSVTYPE=INV;CHR2=chr1;END=19357517;CT=3to5"))) at test-extensions-VCF.R:28:8
  2. ├─StructuralVariantAnnotation:::.testrecord(c("chr10\t2991435\tINV\tN\t<INV>\t.\tLowQual\tSVTYPE=INV;CHR2=chr1;END=19357517;CT=3to5"))
  3. │ ├─VariantAnnotation::readVcf(.testfile(filename), "")
  4. │ └─VariantAnnotation::readVcf(.testfile(filename), "")
  5. │   └─VariantAnnotation (local) .local(file, genome = genome, param = param, ...)
  6. │     └─VariantAnnotation:::.readVcf(...)
  7. │       └─VariantAnnotation:::.scanVcfToVCF(...)
  8. │         ├─VariantAnnotation::scanVcfHeader(file)
  9. │         └─VariantAnnotation::scanVcfHeader(file)
 10. │           ├─Rsamtools::scanBcfHeader(file[[1]], ...)
 11. │           └─Rsamtools::scanBcfHeader(file[[1]], ...)
 12. │             └─BiocGenerics::Map(...)
 13. │               ├─BiocGenerics (local) standardGeneric("Map")
 14. │               │ ├─BiocGenerics::eval(mc, env)
 15. │               │ └─base::eval(mc, env)
 16. │               │   └─base::eval(mc, env)
 17. │               └─base::Map(f = f, ...)
 18. │                 └─base::mapply(FUN = f, ..., SIMPLIFY = FALSE)
 19. │                   └─Rsamtools (local) `<fn>`(dots[[1L]][[1L]])
 20. │                     ├─Rsamtools::scanBcfHeader(bf)
 21. │                     └─Rsamtools::scanBcfHeader(bf)
 22. └─base::.handleSimpleError(...)
 23.   └─base (local) h(simpleError(msg, call))
── Error ('test-extensions-VCF.R:49'): Delly TRA ───────────────────────────────
Error in `h(simpleError(msg, call))`: error in evaluating the argument 'x' in selecting a method for function 'breakpointRanges': no 'header' line "#CHROM POS ID..."?
Backtrace:
     ▆
  1. ├─StructuralVariantAnnotation::breakpointRanges(.testrecord(c("chr10\t2991435\tTRA00000001\tN\t<TRA>\t.\tLowQual\tCIEND=0,100;CIPOS=0,50;SVTYPE=TRA;CHR2=chr1;END=19357517;CT=3to5"))) at test-extensions-VCF.R:49:4
  2. ├─StructuralVariantAnnotation:::.testrecord(c("chr10\t2991435\tTRA00000001\tN\t<TRA>\t.\tLowQual\tCIEND=0,100;CIPOS=0,50;SVTYPE=TRA;CHR2=chr1;END=19357517;CT=3to5"))
  3. │ ├─VariantAnnotation::readVcf(.testfile(filename), "")
  4. │ └─VariantAnnotation::readVcf(.testfile(filename), "")
  5. │   └─VariantAnnotation (local) .local(file, genome = genome, param = param, ...)
  6. │     └─VariantAnnotation:::.readVcf(...)
  7. │       └─VariantAnnotation:::.scanVcfToVCF(...)
  8. │         ├─VariantAnnotation::scanVcfHeader(file)
  9. │         └─VariantAnnotation::scanVcfHeader(file)
 10. │           ├─Rsamtools::scanBcfHeader(file[[1]], ...)
 11. │           └─Rsamtools::scanBcfHeader(file[[1]], ...)
 12. │             └─BiocGenerics::Map(...)
 13. │               ├─BiocGenerics (local) standardGeneric("Map")
 14. │               │ ├─BiocGenerics::eval(mc, env)
 15. │               │ └─base::eval(mc, env)
 16. │               │   └─base::eval(mc, env)
 17. │               └─base::Map(f = f, ...)
 18. │                 └─base::mapply(FUN = f, ..., SIMPLIFY = FALSE)
 19. │                   └─Rsamtools (local) `<fn>`(dots[[1L]][[1L]])
 20. │                     ├─Rsamtools::scanBcfHeader(bf)
 21. │                     └─Rsamtools::scanBcfHeader(bf)
 22. └─base::.handleSimpleError(...)
 23.   └─base (local) h(simpleError(msg, call))
── Error ('test-extensions-VCF.R:71'): empty VCF ───────────────────────────────
Error in `h(simpleError(msg, call))`: error in evaluating the argument 'x' in selecting a method for function 'breakpointRanges': no 'header' line "#CHROM POS ID..."?


You can find the full test log here:

   https://salsa.debian.org/r-pkg-team/r-bioc-structuralvariantannotation/-/jobs/3554988#875

Any hint what might went wrong here?

Kind regards

     Andreas.

--
http://fam-tille.de

Reply to: