SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Varscan-output interpretation bioman1 Bioinformatics 0 05-24-2012 12:50 PM
The VCF file interpretation ashkot Bioinformatics 11 03-13-2012 11:56 AM
Velvet Interpretation nxtgenkid10 De novo discovery 2 02-24-2012 01:43 AM
Help interpretation indel labmtl Illumina/Solexa 0 03-29-2011 07:01 AM
MAQ: SNPs interpretation, etc! nanelle Bioinformatics 0 08-25-2009 10:51 AM

Reply
 
Thread Tools
Old 07-05-2012, 05:18 AM   #1
Lien
Member
 
Location: Leuven

Join Date: Dec 2009
Posts: 47
Default Interpretation of Picard's MEDIAN_CV_COVERAGE

Dear all,

We just performed an RNA-seq experiment from human samples. To get some idea of the quality, I ran Picard's CollectRNASeqMetrics.
In the output, I find 'MEAN_CV_COVERAGE'. The explanation for this value is also on the website: the median CV of coverage of the 1000 most highly expressed transcripts.

I'm not really sure what CV means. And how can I interpret this value? The value I get is 0.48. If I think of this as 48x coverage, this seems really a lot to me. Especially since visualisation with IGV shows me a lower coverage.

Any help will be appreciated!

Thanks
Lien
Lien is offline   Reply With Quote
Old 07-05-2012, 09:20 AM   #2
Bukowski
Senior Member
 
Location: UK

Join Date: Jan 2010
Posts: 390
Default

Quote:
Originally Posted by Lien View Post
Dear all,

We just performed an RNA-seq experiment from human samples. To get some idea of the quality, I ran Picard's CollectRNASeqMetrics.
In the output, I find 'MEAN_CV_COVERAGE'. The explanation for this value is also on the website: the median CV of coverage of the 1000 most highly expressed transcripts.

I'm not really sure what CV means. And how can I interpret this value? The value I get is 0.48. If I think of this as 48x coverage, this seems really a lot to me. Especially since visualisation with IGV shows me a lower coverage.

Any help will be appreciated!

Thanks
Lien
I always assumed CV referred to the coefficient of variation:

http://en.wikipedia.org/wiki/Coefficient_of_variation

As in 'mean coefficient of variation of coverage'

Although happy to be proved wrong.
Bukowski is offline   Reply With Quote
Old 07-24-2013, 01:04 PM   #3
jstjohn
Member
 
Location: San Francisco, CA

Join Date: Jun 2010
Posts: 35
Default from the javadoc

According to the javadoc:
"The median CV of coverage of the 1000 most highly expressed transcripts. Ideal value = 0."

What is confusing me though is it looks like this involves the coefficient of variation (sample sd)/(sample mean). Since this number is reported for a single RNA-seq bam file, I am wondering where the population is coming from? It looks like the only way to get at an sd or mean from the above description is by looking at the coverage of the top 1000 transcripts. That would give you a single number, not something you could take the median of. Pretty sure I am missing something here.


http://picard.sourceforge.net/javadoc/net/sf/picard/analysis/RnaSeqMetrics.html#MEDIAN_CV_COVERAGE
jstjohn is offline   Reply With Quote
Old 04-28-2015, 01:50 AM   #4
pengchy
Senior Member
 
Location: China

Join Date: Feb 2009
Posts: 116
Default What use of MEDIAN_CV_COVERAGE

Quote:
Originally Posted by jstjohn View Post
According to the javadoc:
"The median CV of coverage of the 1000 most highly expressed transcripts. Ideal value = 0."

What is confusing me though is it looks like this involves the coefficient of variation (sample sd)/(sample mean). Since this number is reported for a single RNA-seq bam file, I am wondering where the population is coming from? It looks like the only way to get at an sd or mean from the above description is by looking at the coverage of the top 1000 transcripts. That would give you a single number, not something you could take the median of. Pretty sure I am missing something here.


http://picard.sourceforge.net/javadoc/net/sf/picard/analysis/RnaSeqMetrics.html#MEDIAN_CV_COVERAGE
According to the definition, MEDIAN_CV_COVERAGE is "The median CV of coverage of the 1000 most highly expressed transcripts. Ideal value = 0." I think the variation should be very large, because there will be several gene expressed at very high value. Then, why the ideal value is "0"? "0" means the variation is zero, it not consistent with the fact.
pengchy is offline   Reply With Quote
Reply

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 07:12 PM.


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