miga-base 1.2.2.2 → 1.2.2.3

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 4b1fb93676b3b4c48d2204f48a106575274ee50ab294bbf34a3790e7c9adb94b
-  data.tar.gz: 629d74fd9c861b9afc99c9764b1990e8ea8c9b89e5d10780201828fd68ab3c7c
+  metadata.gz: 3b3fd16e90667560b5cba41046eb92613187279fd7194d1825cc3dee9690273d
+  data.tar.gz: 84043b09c234b825ba0b89ca057b4301a97d5234a87282583516abffbc377e05
 SHA512:
-  metadata.gz: 46d9ad97e0033b0c99e06476594148baaf2bc630f59c4775d763a2ab322a2ad449229fa9f342509fdfe6ca9a0949405f50b17c719fc817eb426f32a9dd287f54
-  data.tar.gz: 4e6ae0abee96441ac65eb2e8e886f4c2f3d00cc9716205936563022d7858f0d966d28d110c3b18ecdaa8aa328ef58b3b53410e9c32560aa0368b88e7ea6604b8
+  metadata.gz: 2e0d76139cfd93f77a91ee266bd596d74e8881da94fb750d97b03bb6c491424e75af208177687525d284ae3fe96e77f1e373f0bdce0583445a4070f3ed09cac0
+  data.tar.gz: 5bd98ea9eb89bcb4af4b74128c2d84ac1f32047322bb2177f7d47543efcc4122a4f8c23f463e8ae8291c8e6dff753182a39fd11fbdf42679975e5b2dba0ca571

data/lib/miga/cli/action/doctor/base.rb

@@ -31,6 +31,7 @@ module MiGA::Cli::Action::Doctor::Base
         file = File.join(base, dir, file_db)
         blk[file, metric, result, rank] if File.exist? file
       end
+      # Query databases for reference databases refer to taxonomy runs
       base = File.join(base, '05.taxonomy')
       result = :taxonomy
     end
@@ -137,10 +138,11 @@ module MiGA::Cli::Action::Doctor::Base
   def save_bidirectional(a, dist)
     each_database_file(a) do |db_file, metric, result, rank|
       next if rank == :haai # No need for hAAI to be bidirectional
+      next if result == :taxonomy # Taxonomy is never bidirectional
 
       b2a = dist[rank].map { |b_name, v| b_name if v[a.name] }.compact
       a2b = dist[rank][a.name]&.keys || []
-      SQLite3::Database.new(db_file) do |db|
+      MiGA::SQLite.new(db_file).run do |db|
         sql = <<~SQL
           insert into #{metric}(seq1, seq2, #{metric}, sd, n, omega) \
           values(?, ?, ?, ?, ?, ?);

data/lib/miga/common/with_result.rb

@@ -26,9 +26,7 @@ module MiGA::Common::WithResult
     task = task.to_sym
     return nil if result_dirs[task].nil?
 
-    base = File.join(
-      project.path, "data/#{result_dirs[task]}/#{result_base}"
-    )
+    base = File.join(project.path, 'data', result_dirs[task], result_base)
     json = "#{base}.json"
     return MiGA::Result.load(json) unless save
 

data/lib/miga/sqlite.rb

@@ -31,16 +31,30 @@ class MiGA::SQLite < MiGA::MiGA
   def initialize(path, opts = {})
     @opts = MiGA::SQLite.default_opts(opts)
     @path = File.absolute_path(path)
+    MiGA::MiGA.DEBUG("Accessing database: #{path}")
   end
 
   ##
-  # Executes +cmd+ and returns the result
+  # Executes +cmd+ and returns the result. Alternatively, if a block is
+  # passed, the commands are ignored and the block is executed with a
+  # single parameter of the database connection
   def run(*cmd)
+    if block_given?
+      run_block { |conn| yield(conn) }
+    else
+      y = nil
+      run_block { |conn| y = conn.execute(*cmd) }
+      y
+    end
+  end
+
+  ##
+  # Executes +blk+ that accepts a single parameter for the database
+  # connection
+  def run_block(&blk)
     busy_attempts ||= 0
     io_attempts ||= 0
-    y = nil
-    SQLite3::Database.new(path) { |conn| y = conn.execute(*cmd) }
-    y
+    SQLite3::Database.new(path) { |conn| blk[conn] }
   rescue SQLite3::BusyException => e
     busy_attempts += 1
     raise "Database busy #{path}: #{e.message}" if busy_attempts >= 3

data/lib/miga/version.rb

@@ -12,7 +12,7 @@ module MiGA
   # - String indicating release status:
   #   - rc* release candidate, not released as gem
   #   - [0-9]+ stable release, released as gem
-  VERSION = [1.2, 2, 2].freeze
+  VERSION = [1.2, 2, 3].freeze
 
   ##
   # Nickname for the current major.minor version.
@@ -20,7 +20,7 @@ module MiGA
 
   ##
   # Date of the current gem relese.
-  VERSION_DATE = Date.new(2022, 2, 2)
+  VERSION_DATE = Date.new(2022, 2, 9)
 
   ##
   # References of MiGA

data/utils/core-pan-plot.R

@@ -1,29 +1,48 @@
 #!/usr/bin/env Rscript
 
-argv <- commandArgs(trailingOnly=T)
+argv <- commandArgs(trailingOnly = TRUE)
 
-plot_core_pan <- function(core_pan, pdf){
-  a <- read.table(core_pan, sep="\t", header=TRUE)
+plot_core_pan <- function (core_pan, pdf) {
+  a <- read.table(core_pan, sep = "\t", header = TRUE)
+  col <- c(
+    rgb(96, 11, 64, max = 255),
+    rgb(0, 121, 166, max = 255),
+    rgb(0.5, 0.5, 0.5, 1/3),
+    rep(rgb(0.5, 0.5, 0.5), 2)
+  )
   pdf(pdf, 7, 5)
-  plot(1, type="n", xlim=c(0, max(a$genomes)*1.05), xaxs="i", yaxs="i",
-    ylim=c(0, max(a$pan_q3)*1.05), xlab="Genomes", ylab="Orthologous Groups")
+  plot(
+    1, type = "n", xlim = c(0, max(a$genomes) * 1.05), xaxs = "i", yaxs = "i",
+    ylim = c(0, max(a$pan_q3) * 1.05)/1e3, xlab = "Genomes",
+    ylab = "Orthologous Groups (thousands)"
+  )
+
   # Core
-  polygon(c(a$genomes, rev(a$genomes)), c(a$core_q1, rev(a$core_q3)),
-    border=NA, col=rgb(0, 121, 166, 128/2, max=255))
-  lines(a$genomes, a$core_avg, col=rgb(0,121,166,max=255), lty=2)
-  lines(a$genomes, a$core_q2, col=rgb(0,121,166,max=255), lty=1)
+  polygon(
+    c(a$genomes, rev(a$genomes)), c(a$core_q1, rev(a$core_q3))/1e3,
+    border = NA, col = rgb(0, 121, 166, 256/4, max = 255)
+  )
+  lines(a$genomes, a$core_avg/1e3, col = col[2], lty = 2)
+  lines(a$genomes, a$core_q2/1e3, col = col[2], lty = 1)
+
   # Pan
-  polygon(c(a$genomes, rev(a$genomes)), c(a$pan_q1, rev(a$pan_q3)),
-    border=NA, col=rgb(96, 11, 64, 128/2, max=255))
-  lines(a$genomes, a$pan_avg, col=rgb(96,11,64,max=255), lty=2)
-  lines(a$genomes, a$pan_q2, col=rgb(96,11,64,max=255), lty=1)
+  polygon(
+    c(a$genomes, rev(a$genomes)), c(a$pan_q1, rev(a$pan_q3))/1e3,
+    border = NA, col = rgb(96, 11, 64, 256/4, max = 255)
+  )
+  lines(a$genomes, a$pan_avg/1e3, col = col[1], lty = 2)
+  lines(a$genomes, a$pan_q2/1e3, col = col[1], lty = 1)
+
   # Legend
-  legend("topleft",
-    legend=c("pangenome","core genome","Inter-Quartile","Median","Average"),
-    pch=c(16,16,15,NA,NA),lty=c(NA,NA,NA,1,2), pt.cex=c(1,1,2,NA,NA),
-    col=c(rgb(96,11,64,max=255), rgb(0,121,166,max=255),
-      rgb(0.5,0.5,0.5,166/255), rep(rgb(0.5,0.5,0.5),2)), bty="n")
+  legend <- c("pangenome", "core genome", "Inter-Quartile", "Median", "Average")
+  legend(
+    "topleft", legend = legend, col = col, bty = "n",
+    pch = c(16, 16, 15, NA, NA),
+    lty = c(NA, NA, NA, 1, 2),
+    pt.cex = c(1, 1, 2, NA, NA),
+  )
   dev.off()
 }
 
 plot_core_pan(argv[1], argv[2])
+

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: miga-base
 version: !ruby/object:Gem::Version
-  version: 1.2.2.2
+  version: 1.2.2.3
 platform: ruby
 authors:
 - Luis M. Rodriguez-R
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2022-02-02 00:00:00.000000000 Z
+date: 2022-02-09 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: daemons