How Do I Convert A Log-Likelihood Obtained From Hmmer To An E-Value Or Vice-Versa?
Hi - I have created a hmm profile for some transport proteins using hmmer, and used the hmmsearch function to return the e-value of the query sequence. I have also built a hmm model using the bayesian network tool in matlab - where the output is in the form of a log-likelihood. My question is how do I convert a log-likelihood to an e-value or vice-versa.
Thanks very much,
Arron.
• 2,947 views
•
link
updated
by
Biostar
•
written
by
arronslacey
•
written
by
arronslacey
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
IQTree: Calculating the Mean Branch Length
written by VenGeno •
-
hmmsearch vs hmmscan - e-value, speed, & output differences?
written by jabaron.phd •
-
What's the point in estimating the background distribution in probabilstic models of motif finding …
written by guixien •
-
HMMER: hmmscan ---> hmmsearch
written by Metaxo •
-
Make a multiple sequence alignment with a HMMER hmmsearch output?
written by westin.kosater •
-
Exact match with query architecture option for local version of HMMER
written by Sej Modha
-
Does hmmsearch look at the reverse complement of the fasta?
written by uwe.menzel •
-
is there a tblastn-version of hmmsearch i.e. search protein profile against nucleotide database?
written by wl284 •
-
Speed of hmmsearch
written by kentnf •
-
Any Ideas To Increase Hmm Classification Accuracy Of Two Groups Of Proteins?
written by arronslacey
I am not expert in protein structure, but are you sure you don't have to use hmmsearch on the hmm model obtained in matlab? Looks like in the first step you built a protein profile and then searched it against a database (which returns E-value) and in the second you created a protein model (for which you obtained confidence in the form of log-likelihood), but to compare with the first profile I think you still need to perform a search against the very same database.