Description Usage Arguments Details Value Author(s) References See Also Examples
This function calculates cumulative relative entropy score according to: Hannenhalli and Russell (2000).
1 | CRE_conservativity(alignment, hmmbuild_path=NULL, pairwiseAlignemnt_scores=NULL)
|
alignment |
An alignment object read with |
hmmbuild_path |
(optional if running under UNIX) The aboslute path to the hmmbuild binary |
pairwiseAlignemnt_scores |
(optional) A matrix with pairwise alignment scores. For example created by |
PSEUDO-ALGORITHM (According to Hannenhalli and Russell (2000)):
(If score matrix is not provided) Run pairwise alignments for all available sequences in the input MSA and save scores to a matrix
(If score matrix is not provided) Calculate a distance matrix based off of the alignment scores one
Perform hierarchical clustering on the distance matrix (UPGMA method)
Get the sequence clusters
Divide the alignment into sub_groups
which are the clusters
Run hmmbuild for whole_alignment
without sub-group
and sub_group
Calculate relative entropy using these two as indicated in the Reference and repeat for each sub_group
Calculate the cumulative relative entropy
hmmbuild program:
This function uses hmmbuild program of HMMER suite for HMM profile generation for MSA.
We recommend downloading and installing HMMER by following the instructions and steps in the HMMER installation website .
score |
A vector of length equal to the length of aligned sequences |
Michal Stolarczyk & Alicja Pluciennik
Hannenhalli, S. S. & Russell, R. B. Analysis and prediction of functional sub-types from protein sequence alignments11Edited by J. Thornton. Journal of Molecular Biology 303, 61–76 (2000).
consensus
, cons2seqs_ident
, read.alignment
1 | #No example due to external software requirements
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.