| -rw-r--r-- | src/Makefile.am | 2 | ||||
| -rw-r--r-- | src/controller/callbacks/reshape.c | 2 | ||||
| -rw-r--r-- | src/db/dbconnect.sqc | 12 | ||||
| -rw-r--r-- | src/model/data/base.sqc | 2 | ||||
| -rw-r--r-- | src/model/state/state.h | 3 | ||||
| -rw-r--r-- | src/util/check_error.c | 6 | ||||
| -rw-r--r-- | src/util/check_error_db.c | 13 | ||||
| -rw-r--r-- | src/util/check_error_db.h | 9 | ||||
| -rw-r--r-- | test/distance_sanity_check/gi_227977170_78032581.png | bin | 0 -> 51161 bytes | |||
| -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 |
18 files changed, 1280 insertions, 9 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 |
