SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
EDGE-pro into DESeq, DESeq run error? epistatic Bioinformatics 5 09-02-2015 04:32 AM
RNA-Seq: EDGE-pro: Estimated Degree of Gene Expression in Prokaryotic Genomes. Newsbot! Literature Watch 0 03-28-2013 02:00 AM
bowtie: negative values rejected? Nicolas902 Bioinformatics 0 12-11-2012 05:36 AM
SAM format: negative values in tags viola1 Bioinformatics 3 03-12-2010 07:29 AM
TopHat results-- negative coverage values seqseqseq Bioinformatics 0 01-21-2010 12:54 PM

Reply
 
Thread Tools
Old 08-28-2014, 09:37 AM   #1
hi-koike
Member
 
Location: Japan

Join Date: Jul 2013
Posts: 13
Default Negative RPKM values by EDGE-pro

Hello,

I am using your EDGE-pro programs for RNA-seq analysis
of bacterial species.
I could run the program, but I got negative values for some of
the samples.

I found rRNA.numberReads was more than numberReads or
numberUniqueReads, thus I got negative value for
rpkm.numberReads.

I tried to find the reason by reading the perl script, but I found
counting of the number of reads was done by a binary code.

I would appreciate it if you could tell me why and how to
do to run the program.

Thanks in advance,
Hideaki
hi-koike is offline   Reply With Quote
Old 08-28-2014, 11:15 AM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

Why don't you send an email to the address listed on the EDGE-pro webpage?
dpryan is offline   Reply With Quote
Old 08-28-2014, 11:22 AM   #3
hi-koike
Member
 
Location: Japan

Join Date: Jul 2013
Posts: 13
Default

Hi dpryan,

Thank you for your advice. Actually I did, but it didn't work.
I would appreciate any idea.

Thank you,
Hideaki
hi-koike is offline   Reply With Quote
Old 08-28-2014, 11:40 AM   #4
blakeoft
Member
 
Location: Connecticut

Join Date: Oct 2013
Posts: 79
Default

The only thing I can think of is that the RPKM values are log transformed, but the manual doesn't seem to mention that they're log transformed. When did you try emailing them? It's possible that they haven't had a chance to reply yet.
blakeoft is offline   Reply With Quote
Old 08-28-2014, 11:51 AM   #5
hi-koike
Member
 
Location: Japan

Join Date: Jul 2013
Posts: 13
Default

Hi blakeoft,

Thank you. As you pointed, I will wait the answer for my email.

I tried to read the perl script, but the count was done by
a binary code "count". So I need a help from a well experienced person.

I have waited only for 3 days, so I will wait.

Thank you,
Hideaki
hi-koike is offline   Reply With Quote
Old 08-28-2014, 11:54 PM   #6
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

You might have to give it a few weeks. Many authors are pretty bad at getting back to people in a timely manner.
dpryan is offline   Reply With Quote
Old 06-05-2015, 03:23 PM   #7
freestile
Member
 
Location: Chile

Join Date: Aug 2014
Posts: 11
Default

Quote:
Originally Posted by hi-koike View Post
Hi blakeoft,

Thank you. As you pointed, I will wait the answer for my email.

I tried to read the perl script, but the count was done by
a binary code "count". So I need a help from a well experienced person.

I have waited only for 3 days, so I will wait.

Thank you,
Hideaki
Finally you solved ?? I have the same problem.
Regards
freestile is offline   Reply With Quote
Old 03-12-2017, 05:11 PM   #8
dbouzo
Junior Member
 
Location: Sydney, Australia

Join Date: Mar 2017
Posts: 1
Default

I am also experiencing this same problem, however my .numberReads file (total number of reads) is a value smaller than that of the rRNA.numberReads file (number of reads mapped to rRNAs). I don't understand how this is possible as I would interpret 'total number of mapped reads' as including the reads from rRNA.

Has anyone else had this issue at all, or an update with a way to resolve this? I have contacted the EDGE-pro developers and hope to hear back soon.
dbouzo is offline   Reply With Quote
Old 03-13-2017, 03:16 AM   #9
Markiyan
Senior Member
 
Location: Cambridge

Join Date: Sep 2010
Posts: 115
Lightbulb Try removing rRNA, tRNA, transposon reads.

In order to minimise the impact of the bug, try removing the reads matching rRNA, tRNA and transposase genes present in a genome of the analysed organism and rerunning the mapping/analysis.

Also check the definition of the rRNA regions (how the program gets those reads numbers).

So first make sure you have a multifasta file with rRNA (full rrn loci) and tRNA seqs for your organism, than you can use my tool as a simple contaminant remover:
(assuming adapters had been clipped already):

fastq_miseq_trimmer [R1.fastq.gz] [R2.fastq.gz] -base=[output_files_basename] -vs=[rtRNAs.fasta] -nohproc -noclip

Other processing options:
To remove illumina adapters: replace -noclip with -ts for Illumina TruSeq adapters clipping. To set quality clipping threshold add: -epc=0.5 (0.5=cumulative error probability in 10bp sliding window).

To set the ends to /1 /2, If your tool wants /1 or /2 as ends specifier, and replace ":" in a read_name with "_"
replace -nohproc with -ends=3

http://seqanswers.com/forums/attachm...3&d=1463489447
Markiyan is offline   Reply With Quote
Old 08-05-2019, 11:20 AM   #10
dmking
Junior Member
 
Location: Michigan

Join Date: Aug 2019
Posts: 1
Default

hi-koike was there a resolution of this issue? Or did you use the solution suggested by Markiyan?

I'm having the same problem and have emailed the authors but would like to know what, if any, approach worked for you.
dmking is offline   Reply With Quote
Old 08-05-2019, 01:37 PM   #11
tdcollingsworth
Junior Member
 
Location: Iowa

Join Date: Aug 2019
Posts: 1
Default Checking In

Hey, guys -

Pretty stumped at this as well: am in the process of trying to debug 'edge.pl' as well as 'count.cpp'.

The code is good but is pretty convoluted and not well documented: pretty standard scientific code, I suppose.

Interested in hearing about anyone else's success in figuring out where the bug is!

Cheers,

Timothy
tdcollingsworth is offline   Reply With Quote
Reply

Tags
bacterial rna seq, edge-pro, rna-seq, rpkm

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 11:05 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO