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