# = chainbraflo1.rb # 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 Mm9 class ChainBraFlo1 KLASS = "ChainBraFlo1" KLASS_S = "chainBraFlo1" Bio::Ucsc::Mm9::CHROMS.each do |chr| class_eval %! class #{chr[0..0].upcase + chr[1..-1]}_#{KLASS} < DBConnection self.table_name = "#{chr[0..0].downcase + chr[1..-1]}_#{KLASS_S}" self.primary_key = nil self.inheritance_column = nil def self.find_by_interval(interval, opt = {:partial => true}); interval = Bio::Ucsc::Gi.wrap(interval) find_first_or_all_by_interval(interval, :first, opt) end def self.find_all_by_interval(interval, opt = {:partial => true}); interval = Bio::Ucsc::Gi.wrap(interval) find_first_or_all_by_interval(interval, :all, opt) end def self.find_first_or_all_by_interval(interval, first_all, opt); interval = Bio::Ucsc::Gi.wrap(interval) 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}); interval = Bio::Ucsc::Gi.wrap(interval) 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}); interval = Bio::Ucsc::Gi.wrap(interval) 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