bio-ucsc-api 0.2.0 → 0.2.1

data/lib/bio-ucsc/pripac1/chaincaerem2link.rb

@@ -0,0 +1,81 @@
+# Copyright::
+#  Copyright (C) 2011 MISHIMA, Hiroyuki <missy at be.to / hmishima at nagasaki-u.ac.jp> 
+# License::     The Ruby licence (Ryby's / GPLv2 dual)
+#
+# In the hg18 database, this table is actually separated
+# into "chr1_*", "chr2_*", etc. This class dynamically
+# define *::Chr1_*, *::Chr2_*, etc. The
+# Rmsk.find_by_interval calls an appropreate class automatically.
+
+module Bio
+  module Ucsc
+    module PriPac1
+
+      class ChainCaeRem2Link
+        KLASS = "ChainCaeRem2Link"
+        KLASS_S = KLASS[0..0].downcase + KLASS[1..-1]
+
+        Bio::Ucsc::PriPac1::CHROMS.each do |chr|
+          class_eval %!
+            class #{chr[0..0].upcase + chr[1..-1]}_#{KLASS} < DBConnection
+              set_table_name "#{chr[0..0].downcase + chr[1..-1]}_#{KLASS_S}"
+              set_primary_key nil
+              set_inheritance_column nil
+
+              def self.find_by_interval(interval, opt = {:partial => true})
+                find_first_or_all_by_interval(interval, :first, opt)
+              end
+        
+              def self.find_all_by_interval(interval, opt = {:partial => true})
+                find_first_or_all_by_interval(interval, :all, opt)
+              end
+
+              def self.find_first_or_all_by_interval(interval, first_all, opt)
+                zstart = interval.zero_start
+                zend   = interval.zero_end
+                if opt[:partial] == true
+                  where = <<-SQL
+      tName = :chrom
+AND   bin in (:bins)
+AND ((tStart BETWEEN :zstart AND :zend)
+ OR  (tEnd BETWEEN :zstart AND :zend)
+ OR  (tStart <= :zstart AND tEnd >= :zend))
+                  SQL
+                else
+                  where = <<-SQL
+      tName = :chrom
+AND   bin in (:bins)
+AND ((tStart BETWEEN :zstart AND :zend)
+AND  (tEnd BETWEEN :zstart AND :zend))
+                  SQL
+                end
+                cond = {
+                  :chrom  => interval.chrom,
+                  :bins   => Bio::Ucsc::UcscBin.bin_all(zstart, zend),
+                  :zstart => zstart,
+                  :zend   => zend,
+                }
+                self.find(first_all,
+                          { :select => "*",
+                            :conditions => [where, cond], })
+              end
+            end
+          !
+        end # each chromosome
+
+        def self.find_by_interval(interval, opt = {:partial => true})
+          chrom = interval.chrom[0..0].upcase + interval.chrom[1..-1]
+          chr_klass = self.const_get("#{chrom}_#{KLASS}")
+          chr_klass.__send__(:find_by_interval, interval, opt)
+        end
+
+        def self.find_all_by_interval(interval, opt = {:partial => true})
+chrom = interval.chrom[0..0].upcase + interval.chrom[1..-1]
+          chr_klass = self.const_get("#{chrom}_#{KLASS}")
+          chr_klass.__send__(:find_all_by_interval, interval, opt)
+        end
+      end # class
+
+    end # module Hg18 
+  end # module Ucsc
+end # module Bio

data/lib/bio-ucsc/pripac1/chaincb3.rb

@@ -0,0 +1,81 @@
+# Copyright::
+#  Copyright (C) 2011 MISHIMA, Hiroyuki <missy at be.to / hmishima at nagasaki-u.ac.jp> 
+# License::     The Ruby licence (Ryby's / GPLv2 dual)
+#
+# In the hg18 database, this table is actually separated
+# into "chr1_*", "chr2_*", etc. This class dynamically
+# define *::Chr1_*, *::Chr2_*, etc. The
+# Rmsk.find_by_interval calls an appropreate class automatically.
+
+module Bio
+  module Ucsc
+    module PriPac1
+
+      class ChainCb3
+        KLASS = "ChainCb3"
+        KLASS_S = KLASS[0..0].downcase + KLASS[1..-1]
+
+        Bio::Ucsc::PriPac1::CHROMS.each do |chr|
+          class_eval %!
+            class #{chr[0..0].upcase + chr[1..-1]}_#{KLASS} < DBConnection
+              set_table_name "#{chr[0..0].downcase + chr[1..-1]}_#{KLASS_S}"
+              set_primary_key nil
+              set_inheritance_column nil
+
+              def self.find_by_interval(interval, opt = {:partial => true})
+                find_first_or_all_by_interval(interval, :first, opt)
+              end
+        
+              def self.find_all_by_interval(interval, opt = {:partial => true})
+                find_first_or_all_by_interval(interval, :all, opt)
+              end
+
+              def self.find_first_or_all_by_interval(interval, first_all, opt)
+                zstart = interval.zero_start
+                zend   = interval.zero_end
+                if opt[:partial] == true
+                  where = <<-SQL
+      tName = :chrom
+AND   bin in (:bins)
+AND ((tStart BETWEEN :zstart AND :zend)
+ OR  (tEnd BETWEEN :zstart AND :zend)
+ OR  (tStart <= :zstart AND tEnd >= :zend))
+                  SQL
+                else
+                  where = <<-SQL
+      tName = :chrom
+AND   bin in (:bins)
+AND ((tStart BETWEEN :zstart AND :zend)
+AND  (tEnd BETWEEN :zstart AND :zend))
+                  SQL
+                end
+                cond = {
+                  :chrom  => interval.chrom,
+                  :bins   => Bio::Ucsc::UcscBin.bin_all(zstart, zend),
+                  :zstart => zstart,
+                  :zend   => zend,
+                }
+                self.find(first_all,
+                          { :select => "*",
+                            :conditions => [where, cond], })
+              end
+            end
+          !
+        end # each chromosome
+
+        def self.find_by_interval(interval, opt = {:partial => true})
+          chrom = interval.chrom[0..0].upcase + interval.chrom[1..-1]
+          chr_klass = self.const_get("#{chrom}_#{KLASS}")
+          chr_klass.__send__(:find_by_interval, interval, opt)
+        end
+
+        def self.find_all_by_interval(interval, opt = {:partial => true})
+          chrom = interval.chrom[0..0].upcase + interval.chrom[1..-1]
+          chr_klass = self.const_get("#{chrom}_#{KLASS}")
+          chr_klass.__send__(:find_all_by_interval, interval, opt)
+        end
+      end # class
+
+    end # module Hg18 
+  end # module Ucsc
+end # module Bio

data/lib/bio-ucsc/pripac1/chaincb3link.rb

@@ -0,0 +1,81 @@
+# Copyright::
+#  Copyright (C) 2011 MISHIMA, Hiroyuki <missy at be.to / hmishima at nagasaki-u.ac.jp> 
+# License::     The Ruby licence (Ryby's / GPLv2 dual)
+#
+# In the hg18 database, this table is actually separated
+# into "chr1_*", "chr2_*", etc. This class dynamically
+# define *::Chr1_*, *::Chr2_*, etc. The
+# Rmsk.find_by_interval calls an appropreate class automatically.
+
+module Bio
+  module Ucsc
+    module PriPac1
+
+      class ChainCb3Link
+        KLASS = "ChainCb3Link"
+        KLASS_S = KLASS[0..0].downcase + KLASS[1..-1]
+
+        Bio::Ucsc::PriPac1::CHROMS.each do |chr|
+          class_eval %!
+            class #{chr[0..0].upcase + chr[1..-1]}_#{KLASS} < DBConnection
+              set_table_name "#{chr[0..0].downcase + chr[1..-1]}_#{KLASS_S}"
+              set_primary_key nil
+              set_inheritance_column nil
+
+              def self.find_by_interval(interval, opt = {:partial => true})
+                find_first_or_all_by_interval(interval, :first, opt)
+              end
+        
+              def self.find_all_by_interval(interval, opt = {:partial => true})
+                find_first_or_all_by_interval(interval, :all, opt)
+              end
+
+              def self.find_first_or_all_by_interval(interval, first_all, opt)
+                zstart = interval.zero_start
+                zend   = interval.zero_end
+                if opt[:partial] == true
+                  where = <<-SQL
+      tName = :chrom
+AND   bin in (:bins)
+AND ((tStart BETWEEN :zstart AND :zend)
+ OR  (tEnd BETWEEN :zstart AND :zend)
+ OR  (tStart <= :zstart AND tEnd >= :zend))
+                  SQL
+                else
+                  where = <<-SQL
+      tName = :chrom
+AND   bin in (:bins)
+AND ((tStart BETWEEN :zstart AND :zend)
+AND  (tEnd BETWEEN :zstart AND :zend))
+                  SQL
+                end
+                cond = {
+                  :chrom  => interval.chrom,
+                  :bins   => Bio::Ucsc::UcscBin.bin_all(zstart, zend),
+                  :zstart => zstart,
+                  :zend   => zend,
+                }
+                self.find(first_all,
+                          { :select => "*",
+                            :conditions => [where, cond], })
+              end
+            end
+          !
+        end # each chromosome
+
+        def self.find_by_interval(interval, opt = {:partial => true})
+          chrom = interval.chrom[0..0].upcase + interval.chrom[1..-1]
+          chr_klass = self.const_get("#{chrom}_#{KLASS}")
+          chr_klass.__send__(:find_by_interval, interval, opt)
+        end
+
+        def self.find_all_by_interval(interval, opt = {:partial => true})
+chrom = interval.chrom[0..0].upcase + interval.chrom[1..-1]
+          chr_klass = self.const_get("#{chrom}_#{KLASS}")
+          chr_klass.__send__(:find_all_by_interval, interval, opt)
+        end
+      end # class
+
+    end # module Hg18 
+  end # module Ucsc
+end # module Bio

data/lib/bio-ucsc/pripac1/chaince6.rb

@@ -0,0 +1,81 @@
+# Copyright::
+#  Copyright (C) 2011 MISHIMA, Hiroyuki <missy at be.to / hmishima at nagasaki-u.ac.jp> 
+# License::     The Ruby licence (Ryby's / GPLv2 dual)
+#
+# In the hg18 database, this table is actually separated
+# into "chr1_*", "chr2_*", etc. This class dynamically
+# define *::Chr1_*, *::Chr2_*, etc. The
+# Rmsk.find_by_interval calls an appropreate class automatically.
+
+module Bio
+  module Ucsc
+    module PriPac1
+
+      class ChainCe6
+        KLASS = "ChainCe6"
+        KLASS_S = KLASS[0..0].downcase + KLASS[1..-1]
+
+        Bio::Ucsc::PriPac1::CHROMS.each do |chr|
+          class_eval %!
+            class #{chr[0..0].upcase + chr[1..-1]}_#{KLASS} < DBConnection
+              set_table_name "#{chr[0..0].downcase + chr[1..-1]}_#{KLASS_S}"
+              set_primary_key nil
+              set_inheritance_column nil
+
+              def self.find_by_interval(interval, opt = {:partial => true})
+                find_first_or_all_by_interval(interval, :first, opt)
+              end
+        
+              def self.find_all_by_interval(interval, opt = {:partial => true})
+                find_first_or_all_by_interval(interval, :all, opt)
+              end
+
+              def self.find_first_or_all_by_interval(interval, first_all, opt)
+                zstart = interval.zero_start
+                zend   = interval.zero_end
+                if opt[:partial] == true
+                  where = <<-SQL
+      tName = :chrom
+AND   bin in (:bins)
+AND ((tStart BETWEEN :zstart AND :zend)
+ OR  (tEnd BETWEEN :zstart AND :zend)
+ OR  (tStart <= :zstart AND tEnd >= :zend))
+                  SQL
+                else
+                  where = <<-SQL
+      tName = :chrom
+AND   bin in (:bins)
+AND ((tStart BETWEEN :zstart AND :zend)
+AND  (tEnd BETWEEN :zstart AND :zend))
+                  SQL
+                end
+                cond = {
+                  :chrom  => interval.chrom,
+                  :bins   => Bio::Ucsc::UcscBin.bin_all(zstart, zend),
+                  :zstart => zstart,
+                  :zend   => zend,
+                }
+                self.find(first_all,
+                          { :select => "*",
+                            :conditions => [where, cond], })
+              end
+            end
+          !
+        end # each chromosome
+
+        def self.find_by_interval(interval, opt = {:partial => true})
+          chrom = interval.chrom[0..0].upcase + interval.chrom[1..-1]
+          chr_klass = self.const_get("#{chrom}_#{KLASS}")
+          chr_klass.__send__(:find_by_interval, interval, opt)
+        end
+
+        def self.find_all_by_interval(interval, opt = {:partial => true})
+          chrom = interval.chrom[0..0].upcase + interval.chrom[1..-1]
+          chr_klass = self.const_get("#{chrom}_#{KLASS}")
+          chr_klass.__send__(:find_all_by_interval, interval, opt)
+        end
+      end # class
+
+    end # module Hg18 
+  end # module Ucsc
+end # module Bio

data/lib/bio-ucsc/pripac1/chaince6link.rb

@@ -0,0 +1,81 @@
+# Copyright::
+#  Copyright (C) 2011 MISHIMA, Hiroyuki <missy at be.to / hmishima at nagasaki-u.ac.jp> 
+# License::     The Ruby licence (Ryby's / GPLv2 dual)
+#
+# In the hg18 database, this table is actually separated
+# into "chr1_*", "chr2_*", etc. This class dynamically
+# define *::Chr1_*, *::Chr2_*, etc. The
+# Rmsk.find_by_interval calls an appropreate class automatically.
+
+module Bio
+  module Ucsc
+    module PriPac1
+
+      class ChainCe6Link
+        KLASS = "ChainCe6Link"
+        KLASS_S = KLASS[0..0].downcase + KLASS[1..-1]
+
+        Bio::Ucsc::PriPac1::CHROMS.each do |chr|
+          class_eval %!
+            class #{chr[0..0].upcase + chr[1..-1]}_#{KLASS} < DBConnection
+              set_table_name "#{chr[0..0].downcase + chr[1..-1]}_#{KLASS_S}"
+              set_primary_key nil
+              set_inheritance_column nil
+
+              def self.find_by_interval(interval, opt = {:partial => true})
+                find_first_or_all_by_interval(interval, :first, opt)
+              end
+        
+              def self.find_all_by_interval(interval, opt = {:partial => true})
+                find_first_or_all_by_interval(interval, :all, opt)
+              end
+
+              def self.find_first_or_all_by_interval(interval, first_all, opt)
+                zstart = interval.zero_start
+                zend   = interval.zero_end
+                if opt[:partial] == true
+                  where = <<-SQL
+      tName = :chrom
+AND   bin in (:bins)
+AND ((tStart BETWEEN :zstart AND :zend)
+ OR  (tEnd BETWEEN :zstart AND :zend)
+ OR  (tStart <= :zstart AND tEnd >= :zend))
+                  SQL
+                else
+                  where = <<-SQL
+      tName = :chrom
+AND   bin in (:bins)
+AND ((tStart BETWEEN :zstart AND :zend)
+AND  (tEnd BETWEEN :zstart AND :zend))
+                  SQL
+                end
+                cond = {
+                  :chrom  => interval.chrom,
+                  :bins   => Bio::Ucsc::UcscBin.bin_all(zstart, zend),
+                  :zstart => zstart,
+                  :zend   => zend,
+                }
+                self.find(first_all,
+                          { :select => "*",
+                            :conditions => [where, cond], })
+              end
+            end
+          !
+        end # each chromosome
+
+        def self.find_by_interval(interval, opt = {:partial => true})
+          chrom = interval.chrom[0..0].upcase + interval.chrom[1..-1]
+          chr_klass = self.const_get("#{chrom}_#{KLASS}")
+          chr_klass.__send__(:find_by_interval, interval, opt)
+        end
+
+        def self.find_all_by_interval(interval, opt = {:partial => true})
+chrom = interval.chrom[0..0].upcase + interval.chrom[1..-1]
+          chr_klass = self.const_get("#{chrom}_#{KLASS}")
+          chr_klass.__send__(:find_all_by_interval, interval, opt)
+        end
+      end # class
+
+    end # module Hg18 
+  end # module Ucsc
+end # module Bio

data/lib/bio-ucsc/pripac1/db_connection.rb

@@ -0,0 +1,57 @@
+#
+# = db_connection.rb - UCSC DB connection
+#
+# Copyright::   Cioyrught (C) 2011
+#               MISHIMA, Hiroyuki <missy at be.to / hmishima at nagasaki-u.ac.jp> 
+# Copyright::   Copyright (C) 2010
+#               Jan Aerts <jan.aerts@gmail.com>
+# License::     Ruby licence (Ryby's / GPLv2 dual)
+
+require 'rubygems'
+require 'active_record'
+
+module Bio
+  module Ucsc
+
+    module PriPac1
+      DATABASE_NAME = "priPac1"
+
+      # = DESCRIPTION
+      # Actual connection established to the UCSC mysql server.
+      class DBConnection < ActiveRecord::Base
+        # = DESCRIPTION
+        # The #connect method make the connection
+        #
+        # = USAGE
+        # Bin::Ucsc::DB_NAME::DBConnection.connect
+        #
+        # ---
+        # *Arguments*: none
+        @@db_adapter  ||= 'mysql'
+        @@db_host     ||= 'genome-mysql.cse.ucsc.edu'
+        @@db_username ||= 'genome'
+        @@db_password ||= ''
+
+        cattr_accessor :db_adapter, :db_host, :db_username, :db_password
+
+        self.abstract_class = true
+
+        def self.default
+          @@db_adapter  = 'mysql'
+          @@db_host     = 'genome-mysql.cse.ucsc.edu'
+          @@db_username = 'genome'
+          @@db_password = ''
+        end
+
+        def self.connect
+          establish_connection({ :adapter =>  @@db_adapter,
+                                 :host =>     @@db_host,
+                                 :database => DATABASE_NAME,
+                                 :username => @@db_username,
+                                 :password => @@db_password, })
+        end
+      end
+    end # module
+    
+  end # module Ucsc
+end # module Bio