[Bioperl-l] BLAST scores
Josh Lauricha
laurichj at bioinfo.ucr.edu
Thu Jun 30 17:55:53 EDT 2005
On Thu 06/30/05 16:48, Joseph Bedell wrote:
> You can calculate the score given the bit score (from the tabular
> output) and Lambda (calculated from the matrix). The equation is Score =
> (Bits)/(Lambda in bits).
>
> Lambda is only dependent upon the matrix. Did you use NCBI-blast or
> WU-BLAST? Which flavor of blast (blastn, blastp, etc)? In any case, you
> can just run a single blast and look at the stats at the bottom of the
> report to get the value of lambda. For example, a default NCBI-blastn
> (+1/-3) search has a lambda of 1.37
>
> ============================
> Lambda K H
> 1.37 0.711 1.31
>
> Gapped
> Lambda K H
> 1.37 0.711 1.31
> ===============================
>
> But, what is difficult to discover is this lambda is in NATS. To convert
> it to bits, divide it by the natural log of 2, or in perl:
>
> perl -e 'print 1.37/log(2),"\n"'
> 1.97649220601788
>
> So, now you can take all of your bit scores divided by 1.97649220601788
> to get the Score.
>
> HTH,
> Joey
Cool, thanks. That'll save me a bunch of time ;) This was NCBI blastp,
so I've already got it calculated ;)
Thanks.
--
------------------------------------------------------
| Josh Lauricha | Ford, you're turning |
| laurichj at bioinfo.ucr.edu | into a penguin. Stop |
| Bioinformatics, UCR | it |
|----------------------------------------------------|
| OpenPG: |
| 4E7D 0FC0 DB6C E91D 4D7B C7F3 9BE9 8740 E4DC 6184 |
|----------------------------------------------------|
| Geek Code: Version 3.12 |
| GAT/CS$/IT$ d+ s-: a-->--- C++++$ UL++++$ P++ L++++|
| $E--- W+ N o? K? w--(---) O? M+(++) V? PS++ PE-(--)|
| Y+ PGP+++ t--- 5+++ X+ R tv DI++ D--- G++ |
| e++ h- r++ z? |
|----------------------------------------------------|
More information about the Bioperl-l
mailing list