Skip to content

No transcripts with expression above threshold #5

@dfermin

Description

@dfermin

Hello.
I'm trying to analyze my data with sqtlseeker2-nf.
Whenever I run the tool with my data I get this error message:


Caused by:
  Process `prepare (all)` terminated with an error exit status (1)

Command executed:

  prepare_trexp.R --group "all" -t transExpression.v2.tsv -m metadata.tsv --gene_location genes.bed --output_tre tre.df.RData --output_gene genes.ss.bed --output_cov covariates.df.RData

Command exit status:
  1

Command output:
  (empty)

Command error:
  Error in prepare.trans.exp(te.df, min.gene.exp = opt$min_gene_expr, min.transcript.exp = opt$min_transcript_expr,  :
    No transcripts with expression above threshold
  Execution halted

My dataset is human RNAseq data for 230 samples. The transcript data is in FPKM (not scaled or transformed).
I've filtered the data to include only genes with > 1 transcript. In addition, I am only keeping transcripts with a relative expression > 0.01 (relative to the total expression for the gene).

Do you have any suggestions for what might be wrong?

Thanks,
Damian

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions