[Biopython] alignment scores -fitting to Poisson distribution

[Iddo Friedberg]

It's an extreme value distribution.
http://en.wikipedia.org/wiki/Generalized_extreme_value_distribution which
is what you generally get from selecting the maximum function from a set of
random variables. (You get a normal distribution from a sum of random
variables).  Aligning proteins provides a maximum function, so multiple
alignment scores are EVD distributed.

http://www.ncbi.nlm.nih.gov/BLAST/tutorial/Altschul-1.html

./I

On Wed, Nov 16, 2011 at 3:39 PM, George Devaniranjan <devaniranjan at gmail.com
> wrote:

> Hello,
>
> I did an alignment of a protein with hundreds of other proteins--I got
> a distribution that is NOT normal (see attached picture).
> I was wondering if biopython had a way of fitting this to one of the
> distribution  methods?
> I don't have mathlab/maple software, so if biopython cannot do it but
> if someone knows of an alternate software that I can use to fit this
> data--I would appreciate the help.
>
> What I want is to say......the target PDB was 'X' sigma away from the
> mean of the distribution.
> If it was a normal distribution, it would be easy to do, since the
> distribution is skewed I am wondering the best way to do it.
>
> Thank you,
> George
>
> _______________________________________________
> Biopython mailing list  -  Biopython at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/biopython
>
>


-- 
Iddo Friedberg
http://iddo-friedberg.net/contact.html
++++++++++[>+++>++++++>++++++++>++++++++++>+++++++++++<<<<<-]>>>>++++.>
++++++..----.<<<<++++++++++++++++++++++++++++.-----------..>>>+.-----.
.>-.<<<<--.>>>++.>+++.<+++.----.-.<++++++++++++++++++.>+.>.<++.<<<+.>>
>>----.<--.>++++++.<<<<------------------------------------.