PapenfussLab / gridss

GRIDSS: the Genomic Rearrangement IDentification Software Suite
Other
258 stars 71 forks source link

Running into errors with conda version while docker/singularity work fine #676

Open famosab opened 2 months ago

famosab commented 2 months ago

I am working on a module for nf-core (https://github.com/nf-core/modules/pull/5932/) where we are now porting to a new testing framework. Running these tests means giving the same input to the conda, docker and singularity version of the same tool. While both docker and singularity complete without errors, I always run into the following error when running gridss somaticfilter:

Attaching package: ‘Matrix’

The following object is masked from ‘package:S4Vectors’:

    expand

The following objects are masked from ‘package:tidyr’:

    expand, pack, unpack

Attaching package: ‘MatrixGenerics’

The following objects are masked from ‘package:matrixStats’:

    colAlls, colAnyNAs, colAnys, colAvgsPerRowSet, colCollapse,
    colCounts, colCummaxs, colCummins, colCumprods, colCumsums,
    colDiffs, colIQRDiffs, colIQRs, colLogSumExps, colMadDiffs,
    colMads, colMaxs, colMeans2, colMedians, colMins, colOrderStats,
    colProds, colQuantiles, colRanges, colRanks, colSdDiffs, colSds,
    colSums2, colTabulates, colVarDiffs, colVars, colWeightedMads,
    colWeightedMeans, colWeightedMedians, colWeightedSds,
    colWeightedVars, rowAlls, rowAnyNAs, rowAnys, rowAvgsPerColSet,
    rowCollapse, rowCounts, rowCummaxs, rowCummins, rowCumprods,
    rowCumsums, rowDiffs, rowIQRDiffs, rowIQRs, rowLogSumExps,
    rowMadDiffs, rowMads, rowMaxs, rowMeans2, rowMedians, rowMins,
    rowOrderStats, rowProds, rowQuantiles, rowRanges, rowRanks,
    rowSdDiffs, rowSds, rowSums2, rowTabulates, rowVarDiffs, rowVars,
    rowWeightedMads, rowWeightedMeans, rowWeightedMedians,
    rowWeightedSds, rowWeightedVars

Attaching package: ‘S4Arrays’

The following object is masked from ‘package:abind’:

    abind

The following object is masked from ‘package:base’:

    rowsum

No reference genome supplied using --ref. Not performing variant equivalence checks.
2024-09-03 17:14:56.3009 Reading test.vcf.gz
Tumour samples:
Matched normals:
Error: array-like subscript has more than one effective dimension
Execution halted

I am thankful for any hints what could cause this problem.

Command used:

#!/bin/bash -ue
gridss_somatic_filter \
    --input test.vcf.gz \
     \
    --output test.high_confidence_somatic.vcf \
    --fulloutput test.all_somatic.vcf \
    --normalordinal 0

cat <<-END_VERSIONS > versions.yml
"GRIDSS_SOMATICFILTER":
    gridss: 2.13.2
END_VERSIONS
lwtan90 commented 2 months ago

I am trying to understand normalordinal parameter. The default was 1. Is there a reason why your command involves --normalordinal 0 ?

famosab commented 2 months ago

I am trying to understand normalordinal parameter. The default was 1. Is there a reason why your command involves --normalordinal 0 ?

I'm afraid I cannot help you there. I did not configure it that way. I think 0 just means that we do not have a matched normal in that case.