[Bioperl-l] $hit_object->frac_aligned_hit/$hit_object->frac_aligned_query

Andreo Beck andreo_beck at yahoo.com
Thu Aug 3 14:14:34 EDT 2006


It still gives more than 1. Can you tell me how good approximation it does? I mean if its 1.3 can I assume it 1 or something like that?

Sendu Bala <bix at sendu.me.uk> wrote: Andreo Beck wrote:
> Can $hit_object->frac_aligned_hit or $hit_object->frac_aligned_query
> give outputs > 1 ? I get some > 1 values.

This should now be fixed in CVS. You should be able to grab just 
Bio/Search/SearchUtils.pm if you don't want to install all of bioperl-live.

( http://code.open-bio.org/cgi/viewcvs.cgi/bioperl-live/ )

It should be noted that all frac_* statistics and probably others from 
hit objects have had a high chance of being wrong in the past, and 
frac_identical and frac_conserved can still be very wrong*. It would be 
a good idea if someone were to make these methods return slightly more 
sane numbers.


[*] which is to say, more wrong than you might reasonably expect, given 
the limitations with gapped or alignment-free blasts


 __________________________________________________
Do You Yahoo!?
Tired of spam?  Yahoo! Mail has the best spam protection around 
http://mail.yahoo.com 


More information about the Bioperl-l mailing list