This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Trinity error using Perl script

Hello all,

I am working on Trinity for de-novo assembly of my data but stuck at differential gene expression step (screenshot). I have 3 sample files and followed the steps as mentioned in the trinity tutorial but somehow getting this error in edgeR.

Error, cmd: R --vanilla -q < matrix.counts.matrix.align_estimate_rsem_P3_vs_align_estimate_rsem_unP3.align_estimate_rsem_unP3.vs.align_estimate_rsem_unP3. EdgeR.Rscript died with ret (32512) at /app/setup/trinityrnaseq-2.2.0/analysis/diffenetialexpression/run_DE_analysis.pl     line 712.

WARNING: This EdgeR comparison failed...

Plz, guide me.

Best Regards, Sakshi

trinity perl script
• 1,695 views • link
updated by ATpoint ATpoint • written by singh.bioinfo • singh.bioinfo

Your code does not allow any reproduction of the error. Where did you get the script from? If custom-made, please share it. Otherwise it is impossible to help.

• 0 views • link
written by ATpoint ATpoint

That script is part of trinity distribution.

@sakshi: What do you mean by you have 3 sample files? You just have three samples that you are trying to do DE analysis on?

• 0 views • link
written by GenoMax GenoMax

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • trouble abundance_estimates_to_matrix.pl (trinity2.4.0)
    written by 1353288408 • 1353288408
    /opt/software/trinityrnaseq-Trinity-v2.4.0/util/abundance_estimates_to_matrix.pl --est_method RSEM --out_prefix transyyy ../*.isoforms.results Output -reading file: ../a.isoforms.results -reading file: ../b.isoforms.results * Outputting combined matrix. /opt/software/trinityrnaseq-Trinity-v2.4.0/util/support_scripts/run_TMM_scale_matrix.pl --matrix transyyy.TPM.not_cross_norm > transyyy.TMM.EXPR.matrixCMD: R --vanilla -q …
  • problem in trinity differential gene expression protocol
    written by majeedaasim • majeedaasim
    After successful quantification of reads through trinity protocol., I tried to generate the count matrix through the provided perl script. as perl /home/aasim/anaconda3/pkgs/trinity-2.5.1-0/opt/trinity-2.5.1/util/abundance_estimates_to_matrix.pl --est_method RSEM …
  • problem in matrix generation after abundance estimation by RSEM?
    written by majeedaasim • majeedaasim
    I am performing a differential gene expression analysis for which I generated the count files through RSEM. Now to produce a count matrix from individual …
  • ERROR for Trinity abundance_estimates_to_matrix.pl
    written by nasreenbano91 • nasreenbano91
    Hello, Trying to run abundance_estimates_to_matrix.pl abundance_estimates_to_matrix.pl --est_method RSEM --out_prefix /home/nbri/nasreen/muku/matrix /home/nbri/nasreen/muku/RSEM_M1.genes.results /home/nbri/nasreen/muku/RSEM_M2.genes.results I get this: -reading file: /scratch/sbag/nasreen/mukuna/align11/RSEM_M11.genes.results -reading file: /scratch/sbag/nasreen/mukuna/align2/RSEM_M2.genes.results * Outputting combined matrix. …
  • Error:abundance_estimates_to_matrix.pl [.isAllZero(counts) : counts must be positive finite value…
    written by jieandze1314 • jieandze1314
    Hi guys, I'm doing the abundance_estimates_to_matrix with 6 RSEM.genes.results. But I met a problem, the error info as follows: -reading file: SS_1/RSEM.genes.results -reading file: SS_2/RSEM.genes.results …
  • problem with trinity
    written by mra8187 • mra8187
    dear all when i am using Trinity v 2.3.2 by this command .`/Trinity --seqType fq --max_memory 5G --single '/media/mra/4A3B44B64CA65A93/SNP/DATA/Queen-Larvae/SRR1571722/Trimmomatic/SE/SRR1571722-1.trimmed.fastq' --CPU 4` i see error ... …
  • ERROR for Trinity abundance_estimates_to_matrix.pl
    written by felix.kuo.1211 • felix.kuo.1211
    This is the ERROR. hadoop@Master:~/simple/DE$ /usr/local/soft/Trinity/util/abundance_estimates_to_matrix.pl --est_method RSEM --out_prefix Trinity_trans Ae.RSEM.isoforms.results AT.RSEM.isoforms.results Be.RSEM.isoforms.results BT.RSEM.isoforms.results Ce.RSEM.isoforms.results CT.RSEM.isoforms.results -reading file: Ae.RSEM.isoforms.results -reading file: AT.RSEM.isoforms.results -reading file: Be.RSEM.isoforms.results …
  • Trinity run failed -- Phase 1
    written by ddowlin ddowlin
    Hi all, I am using Trinity to de novo assemble a transcriptome. Unfortunately my run fails very early on. i.e. in Phase 1 :Clustering of …
  • Format convert error while using Trinity Galaxy
    written by utsafar utsafar
    I am using trinity galaxy to assemble [an RNA-Seq dataset][1] from saffron.I imported File1 and File2 directly to trinity galaxy. Now I am trying to …
  • Error at Inchworm during running Trinity
    written by snar86 • snar86
    Hello, I ran Trinity with this command Trinity --seqType fq --max_memory 20G --left 2025-bark-1b_R1.fastq --right 2025-bark-1b_R2.fastq --SS_lib_type RF --CPU 8 but i got this error: …