public class MultipleAlignmentScorer extends Object
MultipleAlignments.| Modifier and Type | Field and Description |
|---|---|
static String |
AVGTM_SCORE |
static String |
CE_SCORE |
static String |
MC_SCORE |
static String |
PROBABILITY |
static String |
REF_RMSD |
static String |
REFTM_SCORE |
static String |
RMSD |
| Constructor and Description |
|---|
MultipleAlignmentScorer() |
| Modifier and Type | Method and Description |
|---|---|
static void |
calculateScores(MultipleAlignment alignment)
Calculates and puts the RMSD and the average TM-Score of the
MultipleAlignment.
|
static double |
getAvgTMScore(List<Atom[]> transformed,
List<Integer> lengths)
Calculates the average TMScore all the possible pairwise structure
comparisons of the given a set of superimposed Atoms and the original
structure lengths.
|
static double |
getAvgTMScore(MultipleAlignment alignment)
Calculates the average TMScore of all the possible pairwise structure
comparisons of the given alignment.
|
static double |
getMCScore(MultipleAlignment alignment,
double gapOpen,
double gapExtension,
double dCutoff)
Calculates the MC score, specific for the MultipleAlignment algorithm.
|
static double |
getRefRMSD(List<Atom[]> transformed,
int reference)
Calculates the average RMSD from all structures to a reference s
tructure, given a set of superimposed atoms.
|
static double |
getRefRMSD(MultipleAlignment alignment,
int ref)
/** Calculates the average RMSD from all structures to a reference s
tructure, given a set of superimposed atoms.
|
static double |
getRefTMScore(List<Atom[]> transformed,
List<Integer> lengths,
int reference)
Calculates the average TMScore from all structures to a reference
structure, given a set of superimposed atoms.
|
static double |
getRefTMScore(MultipleAlignment alignment,
int ref)
Calculates the average TMScore from all structures to a reference
structure, given a set of superimposed atoms.
|
static double |
getRMSD(List<Atom[]> transformed)
Calculates the RMSD of all-to-all structure comparisons (distances),
given a set of superimposed atoms.
|
static double |
getRMSD(MultipleAlignment alignment)
Calculates the RMSD of all-to-all structure comparisons (distances) of
the given MultipleAlignment.
|
public static final String PROBABILITY
public static final String CE_SCORE
public static final String RMSD
public static final String AVGTM_SCORE
public static final String MC_SCORE
public static final String REF_RMSD
public static final String REFTM_SCORE
public static void calculateScores(MultipleAlignment alignment) throws StructureException
alignment - StructureExceptiongetAvgTMScore(MultipleAlignment),
getRMSD(MultipleAlignment)public static double getRMSD(MultipleAlignment alignment)
Complexity: T(n,l) = O(l*n^2), if n=number of structures and l=alignment length.
The formula used is just the sqroot of the average distance of all possible pairs of atoms. Thus, for R structures aligned over C columns without gaps, we have
RMSD = \sqrt{1/(C*(R*(R-1)/2)) * \sum_{r1=1}^{R-1}
\sum_{r2=r1+1}^{R-1} \sum_{j=0}^{C-1} (atom[r1][c]-atom[r2][c])^2}
alignment - public static double getRMSD(List<Atom[]> transformed)
transformed - getRMSD(MultipleAlignment)public static double getRefRMSD(MultipleAlignment alignment, int ref)
Complexity: T(n,l) = O(l*n), if n=number of structures and l=alignment length.
For ungapped alignments, this is just the sqroot of the average distance from an atom to the aligned atom from the reference. Thus, for R structures aligned over C columns (with structure 0 as the reference), we have:
RefRMSD = \sqrt{ 1/(C*(R-1)) * \sum_{r=1}^{R-1} \sum_{j=0}^{C-1}
(atom[1][c]-atom[r][c])^2 }
For gapped alignments, null atoms are omitted from consideration, so that the RMSD is the average over all columns with non-null reference of the average RMSD within the non-null elements of the column.
alignment - MultipleAlignmentref - reference structure indexpublic static double getRefRMSD(List<Atom[]> transformed, int reference)
Complexity: T(n,l) = O(l*n), if n=number of structures and l=alignment length.
For ungapped alignments, this is just the sqroot of the average distance from an atom to the aligned atom from the reference. Thus, for R structures aligned over C columns (with structure 0 as the reference), we have:
RefRMSD = \sqrt{ 1/(C*(R-1)) * \sum_{r=1}^{R-1} \sum_{j=0}^{C-1}
(atom[1][c]-atom[r][c])^2 }
For gapped alignments, null atoms are omitted from consideration, so that the RMSD is the average over all columns with non-null reference of the average RMSD within the non-null elements of the column.
transformed - reference - public static double getAvgTMScore(MultipleAlignment alignment) throws StructureException
Complexity: T(n,l) = O(l*n^2), if n=number of structures and l=alignment length.
alignment - StructureExceptionpublic static double getAvgTMScore(List<Atom[]> transformed, List<Integer> lengths) throws StructureException
Complexity: T(n,l) = O(l*n^2), if n=number of structures and l=alignment length.
transformed - aligned Atoms transformedlengths - lengths of the structures in residue numberStructureExceptionpublic static double getRefTMScore(MultipleAlignment alignment, int ref) throws StructureException
Complexity: T(n,l) = O(l*n), if n=number of structures and l=alignment length.
alignment - reference - Index of the reference structureStructureExceptionpublic static double getRefTMScore(List<Atom[]> transformed, List<Integer> lengths, int reference) throws StructureException
Complexity: T(n,l) = O(l*n^2), if n=number of structures and l=alignment length.
transformed - Arrays of aligned atoms, after superpositionlengths - lengths of the full input structuresreference - Index of the reference structureStructureExceptionpublic static double getMCScore(MultipleAlignment alignment, double gapOpen, double gapExtension, double dCutoff) throws StructureException
The maximum score of a match is 20, and the penalties for gaps are part of the input. The half-score distance, d0, is chosen as in the TM-score.
Complexity: T(n,l) = O(l*n^2), if n=number of structures and l=alignment length.
alignment - gapOpen - penalty for gap opening (reasonable values are in the range
(1.0-20.0)gapExtension - penalty for extending a gap (reasonable values are in the
range (0.5-10.0)dCutoff - the distance cutoffStructureExceptionCopyright © 2000–2016 BioJava. All rights reserved.