SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
PubMed: Identification of novel and known miRNAs in deep-sequencing data with miRDeep Newsbot! Literature Watch 2 07-12-2013 06:06 AM
Parameters of PaciBio SOLiDance Pacific Biosciences 5 11-22-2011 04:00 PM
Novoalign parameters madsaan Bioinformatics 1 10-03-2011 05:34 PM
CAP3 parameters BOYD Bioinformatics 7 07-29-2011 02:43 AM
Velvet Parameters Mona De novo discovery 5 02-14-2010 11:50 PM

Reply
 
Thread Tools
Old 09-29-2010, 11:22 PM   #1
NicoBxl
not just another member
 
Location: Belgium

Join Date: Aug 2010
Posts: 264
Default miRDeep -v and -z parameters

Hi,

I'm using miRDeep to discover new miRNA in Illumina HTS data. I don't understand the -v and the -z parameters.

here's the definition in the readme of the script

-v (for instance, -v 5) varies the cut-off. This is the single most important and difficult option. In
general, if you are missing a large fraction of known miRNAs that you know are in the dataset, you should
lower the cut-off. If you make the permuted controls (see below) and get many false positives, you should
raise the cut-off. As a rule of thumb, the cut-off should probably not be lowered below 0 or 1 unless the
dataset is _highly_ enriched in miRNA sequences.


the definiton is ok, but how to set a good thresold ? on what is it filtering ?


-z scores the stems below the potential precursor sequence for potential Drosha recognition. In general,
more base pairings between the 10 nt flanking the potential precursor suggest a tighter stem and gives a
higher score contribution. For each number of base pairings b in the stem, the probability p(b) of a
miRNA precursor hairpin or a background hairpin having p(b) was estimated, and log-odd score contributions
were calculated. When this option is used, we recommend that a slightly higher cut-off is used (as default
1.5).



can someone explain me that in other words, I'm a little lost there.


Thanks
NicoBxl is offline   Reply With Quote
Old 10-04-2010, 06:37 AM   #2
NicoBxl
not just another member
 
Location: Belgium

Join Date: Aug 2010
Posts: 264
Default

another question about the permutation script : how can i calculate the false positive rate from the result of the permutations ?
NicoBxl 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 06:08 PM.


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