| -rw-r--r-- | test/entropy/Align2Ref.m | 24 | ||||
| -rw-r--r-- | test/entropy/CalculateEntropy.m | 15 | ||||
| -rw-r--r-- | test/entropy/CalculateProteinEntropy.m | 25 | ||||
| -rw-r--r-- | test/entropy/FastNWalign2.c | 94 | ||||
| -rw-r--r-- | test/entropy/GenomeAlignments.m | 31 | ||||
| -rw-r--r-- | test/entropy/GenomePairwiseDist.m | 98 | ||||
| -rw-r--r-- | test/entropy/RefineAlignments.m | 276 | ||||
| -rw-r--r-- | test/entropy/don_anal.m | 40 | ||||
| -rw-r--r-- | test/entropy/nwalign_mod.m | 637 |
9 files changed, 1240 insertions, 0 deletions
diff --git a/test/entropy/CalculateEntropy.m b/test/entropy/CalculateEntropy.m new file mode 100644 index 0000000..18bfda5 --- a/dev/null +++ b/test/entropy/CalculateEntropy.m @@ -0,0 +1,15 @@ +function OUT_VEC = CalculateEntropy(ARRAY)
+% CalculateEntropy
+% Calculate the Informationational Entropy of each position in a
+% multialignment. In this function I remove all gap characters from
+% each column of the alignment before doing the calculation.
+%
+%
+%
+
+num_array = aa2int(ARRAY);
+bins = histc(num_array, 0:21, 1);
+
+probs=bsxfun(@rdivide, bins(2:end-1,:),sum(bins(2:end-1,:)));
+
+OUT_VEC = -nansum(log(probs).*double(probs~=0).*probs);
\ No newline at end of file |
