# Copyright:: # Copyright (C) 2011 MISHIMA, Hiroyuki # 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 DroSim1 class ChainDm3 KLASS = "ChainDm3" KLASS_S = KLASS[0..0].downcase + KLASS[1..-1] Bio::Ucsc::DroSim1::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