gene-matcher 0.1.5 → 0.1.6

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: f451b24184e517be0bfdf398c1ad227be0c956a7e3173f7a31a45530af64291d
-  data.tar.gz: ee44f41ce12e8ac4a311c58a59339621230513bbe6202083bae3f044f830f4c6
+  metadata.gz: 16d36396bec9b08dc0cd41f911274feef3b444c7c5408e406b528d627f74e92b
+  data.tar.gz: 71f75988056e3de2a0e022ed7b501e88209866f373e753e471298462ff7f428e
 SHA512:
-  metadata.gz: 7237131d90e3fed9a1f935c818bd75c8179a3de6a0bbcf4b7073cced5c3bc314a769680e73cfe46f3f228328427e447e0ee7fb22ec632920eea0b14c29510c98
-  data.tar.gz: 27614deb05f5649814174f3c5e4f8b771d3610c4dd9d8034409326f7dff77a32140e1fd038e6f055b1aa1661f0d3ea7ffe18cd0856617bd30c524ef5a28ee9b2
+  metadata.gz: 47f79bf0bf8be63ada87103503d2ec0ef7717cf2d008ff836dbf9d5582c035f082e6049bb3add93c6130c15024a015bb7f52aa8f04300963cafed1119a0f1d4c
+  data.tar.gz: 38e05772aacdd96a57d43ff0cb16f19c00adf300c083d85b6fbb93cbb0449ab51212981d94ef1cbe5a38c542fa707d9bd9cb39dc0b58b3a3782ed9f5ab5eeda5

data/lib/alignment.rb

@@ -1,12 +1,10 @@
 class Alignment
-  attr_accessor :score, :alignmentI, :alignmentJ, :startI, :startJ, :endI, :endJ, :reversed, :aside, :source
+  attr_accessor :alignmentI, :alignmentJ, :startI, :startJ, :endI, :endJ, :reversed, :aside, :source
 
   BLANK = "-"
 
   def initialize(h=nil)
     return init_from_hash(h) if h
-    # スコア
-    @score = 0
     # 検索対象（データベースに入っていた）配列
     @alignmentI = ""
     # 検索した（クエリ=入力された）配列
@@ -27,7 +25,6 @@ class Alignment
 
   def init_from_hash(h)
     # 引数のハッシュから初期化
-    @score = h["score"]
     @alignmentI = h["alignmentI"]
     @alignmentJ = h["alignmentJ"]
     @startI = h["startI"]
@@ -42,7 +39,6 @@ class Alignment
   def to_h
     # ハッシュに変換
     h = {}
-    h["score"] = @score
     h["alignmentI"] = @alignmentI
     h["alignmentJ"] = @alignmentJ
     h["startI"] = @startI
@@ -87,7 +83,17 @@ class Alignment
   def number_of_blank
   end
 
+  def score_a
+    return alignment_count / @alignmentI.length.to_f
+  end
+
+  def score(thresh=20)
+    # 長さがthresh以下なら0を返す
+    return 0 if @alignmentI.length < thresh
+    return alignment_count * score_a
+  end
+
   def to_s
-    "score=#{@score} I=" + @alignmentI + " J=" + @alignmentJ;
+    "score=#{score} I=" + @alignmentI + " J=" + @alignmentJ;
   end
 end

data/lib/gene-matcher.rb

@@ -14,8 +14,9 @@ class Matcher
   def scan(target_sequence,source = {})
     sw = SmithWaterman.instance
     a = sw.alignment(target_sequence, @input_sequence)
-    a.source = source
-    @alignments += [a] if a.score >= @limit
+    a.source = {target_sequence: target_sequence}.merge(source)
+    puts "#{a.score} / #{a.alignment_count} = #{a.score / a.alignment_count}" if ENV["DEBUG"]
+    @alignments += [a] if a.score / a.alignment_count >= @limit
   end
 end
 

data/lib/smith-waterman.rb

@@ -88,11 +88,7 @@ class SmithWaterman
     a.startI = i
     a.startJ = j
 
-    if a.alignmentI.length <= 20
-      a.score = 0
-    else
-      a.score = a.alignment_count ** 2 / a.alignmentI.length.to_f
-    end
+    # スコアはここでは求めない(Java版ではここで決定していた)
     return a
   end
 

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: gene-matcher
 version: !ruby/object:Gem::Version
-  version: 0.1.5
+  version: 0.1.6
 platform: ruby
 authors:
 - ITO Yosei
 autorequire:
 bindir: bin
 cert_chain: []
-date: 2023-08-31 00:00:00.000000000 Z
+date: 2023-09-11 00:00:00.000000000 Z
 dependencies: []
 description: Algorithm for determining similar regions between nucleic acid sequences.
 email: y-itou@lumber-mill.co.jp