tree_clusters 0.5.1 → 0.5.2

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA1:
-  metadata.gz: ad186617714194c3d20bb4bc0a6e79d9779c4e29
-  data.tar.gz: 30cbedb0728fffae00c185735e0ee3d583741e17
+  metadata.gz: 7d7f7982ea210cef7357e004ceb85d8b17389a32
+  data.tar.gz: a924d2bfd587ce03678c391f2a5f761848169726
 SHA512:
-  metadata.gz: a7406ddcd526a485e2721a82bd74ccd5295b7a4c07482d40d5127f2f9e6e50b97b59deb3215576f9420e73315d1f4e8394c8ebef0c278e482b9ce693dbcaf5fe
-  data.tar.gz: ceceec6c5f11037da5261e89754d02eeb463ea634f38df0ecfe150c78cf21ceb128c96a12e2ba662df4b808f625af5a97883d0638c9527ce9a977d148fa54c8b
+  metadata.gz: 7b4540a485a4635194a170c74c1bae7cd4eeb5b58df3d3a934933a161ef918903461ba75924f3baf65724050b2d3d122ac953a0822f5218c7ac0f8b37db72405
+  data.tar.gz: fced0e673bcacf6dc113d6112bd5d02b9d5fcf5756ae50d1b77cb41b1c8f97b76a3fe2ef86f8f6c7ea8ddeee96a5b5aa06600cdef310842276a486b829763eae

data/.gitignore

@@ -21,6 +21,8 @@ snazzy_clades.*
 
 TEST
 
+time.html
+
 # rspec failure tracking
 .rspec_status
 *.lock

data/exe/snazzy_clades_key_cols

@@ -19,9 +19,42 @@ opts = Trollop.options do
 
   banner <<-EOS
 
-  Note that if a clade's parent would be the root of the tree, no
-  columns will be subtracted when removing the parent columns as it
-  would be the entire alignment.
+
+  Checking IDs
+  ------------
+
+  IDs for the sequences must match between the three input files.
+
+  The tree file is allowed to have quoted taxa names, but the mapping
+  file and alignment file are not.
+
+  If your alignment file has spaces in the name, the ID part of the
+  header (i.e., the part up until the space) must match with the
+  sequence IDs in the tree and the mapping file.
+
+  Example: This would be okay.
+
+  tree file:
+    ('genome_A', 'genome_B');
+
+  aln file:
+    >genome_A apple pie
+    AAAAA
+    >genome_B brown sugar
+    AATTA
+
+  mapping file:
+    name	coolness
+    genome_A	cool
+    genome_B	notcool
+
+
+  Subtracting parent nodes
+  ------------------------
+
+  If a clade's parent would be the root of the tree, no columns will
+  be subtracted when removing the parent columns as it would be the
+  entire alignment.
 
   Options:
   EOS

data/lib/tree_clusters.rb

@@ -22,6 +22,10 @@ class NewickTree
     end
     return clades
   end
+
+  def unquoted_taxa
+    self.taxa.map { |str| str.tr %q{"'}, "" }
+  end
 end
 
 # Top level namespace of the Gem.
@@ -82,8 +86,10 @@ module TreeClusters
     Set.new low_ent_cols
   end
 
+  # @note If there are quoted names in the tree file, they are
+  #   unquoted first.
   def check_ids tree, mapping, aln
-    tree_ids = Set.new(NewickTree.fromFile(tree).taxa)
+    tree_ids = Set.new(NewickTree.fromFile(tree).unquoted_taxa)
 
     mapping_ids = Set.new
     File.open(mapping, "rt").each_line.with_index do |line, idx|
@@ -143,8 +149,11 @@ module TreeClusters
     metadata.each do |md_cat, leaf2mdtag|
       already_checked = Set.new
       single_tag_clades = {}
+      p [md_cat, leaf2mdtag]
 
       clades.each do |clade|
+        p [clade.name, clade.all_leaves]
+
         assert clade.all_leaves.count > 1,
                "A clade cannot also be a leaf"
 
@@ -173,7 +182,7 @@ module TreeClusters
       end
 
       single_tag_clades.each do |clade, md_tag|
-        non_clade_leaves = tree.taxa - clade.all_leaves
+        non_clade_leaves = tree.unquoted_taxa - clade.all_leaves
 
         non_clade_leaves_with_this_md_tag = non_clade_leaves.map do |leaf|
           [leaf, leaf2mdtag[leaf]]
@@ -288,10 +297,15 @@ module TreeClusters
                   :single_tag_info,
                   :all_tags
 
+    # @note If a node name is quoted, then those quotes are removed
+    #   first.
+    #
     # @param node [NewickNode] a NewickNode from a NewickTree
     # @param tree [NewickTree] a NewickTree
     def initialize node, tree, metadata=nil
-      @name = node.name
+      tree_taxa = tree.unquoted_taxa
+
+      @name = unquote node.name
       @all_leaves = descendant_leaves node
 
       if (children = node.children).count == 2
@@ -317,10 +331,10 @@ module TreeClusters
       @parent_leaves = descendant_leaves parent
 
       @other_leaves =
-        Object::Set.new(tree.taxa) - Object::Set.new(all_leaves)
+        Object::Set.new(tree_taxa) - Object::Set.new(all_leaves)
 
       @non_parent_leaves =
-        Object::Set.new(tree.taxa) - Object::Set.new(parent_leaves)
+        Object::Set.new(tree_taxa) - Object::Set.new(parent_leaves)
 
       if metadata
         @metadata = metadata
@@ -345,7 +359,8 @@ module TreeClusters
         self.each_sibling_leaf_set == clade.each_sibling_leaf_set &&
         self.parent_leaves == clade.parent_leaves &&
         self.other_leaves == clade.other_leaves &&
-        self.single_tag_info == clade.single_tag_info
+        self.single_tag_info == clade.single_tag_info &&
+        self.all_tags == clade.all_tags
       )
     end
 
@@ -379,14 +394,19 @@ module TreeClusters
 
     def descendant_leaves node
       if node.leaf?
-        [node.name]
+        [unquote(node.name)]
       else
         node.
           descendants.
           flatten.
           uniq.
-          select { |node| node.leaf? }.map(&:name)
+          select { |node| node.leaf? }.
+          map { |node| unquote(node.name) }
       end
     end
+
+    def unquote str
+      str.tr %q{"'}, ""
+    end
   end
 end

data/lib/tree_clusters/version.rb

@@ -1,3 +1,3 @@
 module TreeClusters
-  VERSION = "0.5.1"
+  VERSION = "0.5.2"
 end

data/test_files/small.tre

@@ -1 +1 @@
-((a-1, a-2)cluster_A, ((b-1, b-2)cluster_B1, (bb-1, (bbb-1, bbb-2)cluster_B3)cluster_B2)cluster_B)cluster_C;
+(('a-1', 'a-2')cluster_A, (("b-1", b-2)cluster_B1, (bb-1, (bbb-1, bbb-2)cluster_B3)cluster_B2)cluster_B)cluster_C;

data/test_files/test.tre

@@ -1 +1 @@
-((((((g1:1, g2:1)cluster1:1, g3:2)cluster4:1, (g4a:1, g4b:1)cluster7:2)cluster6:1, ((g5:2, g6:2)cluster11:1, g7:3)cluster14:1)cluster10:1, g8:5)cluster16:1, ((g9:1, g10:1)cluster19:1, g11:2)cluster22:4)cluster18:1;
+((((((g1:1, 'g2':1)cluster1:1, g3:2)cluster4:1, (g4a:1, g4b:1)cluster7:2)cluster6:1, ((g5:2, "g6":2)cluster11:1, g7:3)cluster14:1)cluster10:1, g8:5)cluster16:1, ((g9:1, g10:1)cluster19:1, g11:2)cluster22:4)cluster18:1;

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: tree_clusters
 version: !ruby/object:Gem::Version
-  version: 0.5.1
+  version: 0.5.2
 platform: ruby
 authors:
 - Ryan Moore