estimating K and Lambda from an extreme value distribution

Gordon D. Pusch g_d_pusch_remove_underscores at
Mon Feb 23 00:13:46 EST 2004

ranjeeva_r at (Ranjeeva) writes:

> I'm trying to fit a set of scores I get from searching  a database of
> 1000 amino acid sequences with a HMM. I want to calculate a p-value
> for each matching score. My questions are
> a) How do you estimate the scalling factors K and Lambda to fit my
> scores (1000) to an extreme value distribution?

The obvious question would be: Why would you bother, since an HMM _directly_
yields a generative probability estimate?  Simply compare the HMM probability
estimate to that of a "fiducial model," e.g., the "random sequence" model. 

However, if you _insist_ on (ab)using extreme-value theory for this problem,
googling on the exact phrase "extreme value distribution" plus "fitting"
yields 1,740 hits.

-- Gordon D. Pusch   

perl -e '$_ = "gdpusch\\n"; s/NO\.//; s/SPAM\.//; print;'

More information about the Comp-bio mailing list