drop icon indicating copy to clipboard operation
drop copied to clipboard

Rows (gene IDs) in pipeline's Rds expression counts and rows in the recommended external geneCounts.tsv.gz are in different order

Open sitems opened this issue 5 months ago • 0 comments

read 19 files Error: The rows (genes) of the count matrices to be merged are not the same. Execution halted Error in rule AberrantExpression_pipeline_Counting_mergeCounts_R:

I am using drop 1.5.0 with my own samples (bams obtained from rnaseq nf core pipeline), mapped to gencode29, together with your external samples from https://zenodo.org/records/6078397. I inspected the gene IDs (i.e. rownames) in external as well as in my data and values are exactly the same, just the order of rows is different (and that probably causes failing the pipeline). External rownames (gtex/v8/exported_counts/Whole_Blood--GRCh38--gencode29/geneCounts.tsv.gz) are sorted, but rownames from your RDS files (Output/processed_data/aberrant_expression/v29/counts/somename.Rds) are not.

sitems avatar Nov 10 '25 20:11 sitems