statsample 0.15.1 → 0.16.0

data/lib/statsample/graph/gdchart.rb

@@ -1,45 +0,0 @@
-require 'GDChart'
-module Statsample
-	module Util # :nodoc:
-	class << self
-		def chart_gdchart(file,width,height,chart_type, labels, options,num_datasets,data)
-				require 'GDChart'
-				gdc=GDChart.new
-				gdc.title="Generic title"
-				gdc.bg_color=0xFFFFFF
-				gdc.image_type=GDChart::JPEG
-				options.each{|k,v|
-					gdc.send(k+"=",v)
-				}
-                
-				File.open(file,"w") {|f|
-					gdc.out_graph(width,height,f,chart_type, data.length/num_datasets,labels,num_datasets,data)
-				}
-		end
-		end
-	end
-	class Vector 
-		# Creates a barchart using ruby-gdchart
-		def gdchart_frequencies(file, width=300, height=150, chart_type=GDChart::BAR, options={}) # :nodoc:
-			labels,data=[],[]
-			self.frequencies.sort.each{|k,v|
-				labels.push(k.to_s)
-				data.push(v) 
-			}
-			options['ext_color']=[0xFF3399,0xFF9933,0xFFEE33,0x33FF33, 0x9966FF]
-			Statsample::Util.chart_gdchart(file,width,height,chart_type, labels,options,1,data)
-		end
-		def gdchart_histogram(bins,file, width=300, height=150, chart_type=GDChart::BAR, options={}) # :nodoc:
-      check_type :scale
-      labels=[]
-      h=histogram(bins)
-      data=[]
-      (0...bins).each{|bin|
-        data.push(h[bin])
-        range=h.get_range(bin)
-        labels.push(((range[0]+range[1]) / 2.to_f).to_s)
-      }
-      Statsample::Util.chart_gdchart(file, width, height, chart_type, labels,options, 1,data)
-		end
-	end
-end

data/lib/statsample/graph/svgboxplot.rb

@@ -1,108 +0,0 @@
-module Statsample
-  module Graph
-    class SvgBoxplot < SVG::Graph::Bar # :nodoc:
-        def initialize(config={})
-            config[:fields]=["dummy"]
-            super(config)
-        end
-        def get_x_labels
-            @data.collect{|d|
-                d[:title]
-            }
-        end
-        
-        def min_value
-            min = 0
-            if min_scale_value.nil? 
-                min = @data.collect{|x| x[:data].min}.min
-                if min > 0
-                    if min > 10
-                        min=min-2
-                    else
-                        min=0
-                    end
-                end
-            else
-                min = min_scale_value
-            end
-            return min
-        end
-        def draw_data
-        minvalue = min_value
-        fieldwidth = field_width
-        unit_size =  (@graph_height.to_f - font_size*2*top_font) / 
-                      (get_y_labels.max - get_y_labels.min)
-        bargap = bar_gap ? (fieldwidth < 10 ? fieldwidth / 2 : 10) : 0
-        
-        bar_width = fieldwidth - bargap
-        bar_width /= @data.length if stack == :side
-        x_mod = (@graph_width-bargap) / 2 - (stack==:side ? bar_width/2 : 0)
-        
-        bottom = @graph_height
-        
-        field_count = 0
-        for dataset in @data
-        
-        # cases (assume 0 = +ve):
-        #   value  min  length
-        #    +ve   +ve  value - min
-        #    +ve   -ve  value - 0
-        #    -ve   -ve  value.abs - 0
-        
-        min=dataset[:data].min
-        max=dataset[:data].max
-        median=dataset[:vector].median
-        q1=dataset[:vector].percentil(25)
-        q3=dataset[:vector].percentil(75)
-        iqr=q3-q1
-        left = (fieldwidth * field_count)
-        #length = (value.abs - (minvalue > 0 ? minvalue : 0)) * unit_size
-        # top is 0 if value is negative
-        top_wisk=(q3+iqr*1.5 < max) ? q3+iqr*1.5 : max
-        down_wisk= (q1-iqr*1.5 > min) ? q1-iqr*1.5 : min
-        
-        top=@graph_height-((top_wisk-minvalue)*unit_size)
-        down=@graph_height-((down_wisk-minvalue)*unit_size)
-        
-        median_bar=@graph_height-((median-minvalue)*unit_size)
-        middle= left+(bar_width / 2)
-        left_whis=left+(bar_width * 0.4)
-        rigth_whis=left+(bar_width*0.6)
-        left_rect= left+(bar_width * 0.25)
-        rigth_rect = left+ (bar_width * 0.75)
-        top_rect=@graph_height-((q3-minvalue)*unit_size)
-        height_rect=iqr*unit_size
-        path="M #{left_whis} #{top} H #{rigth_whis} M #{middle} #{top} V #{down} M #{left_whis} #{down} H #{rigth_whis} M #{left_rect} #{median_bar} H #{rigth_rect}"
-        
-        
-        # Marcamos Outliers
-        if top_wisk!=max or down_wisk!=min
-            dataset[:vector].valid_data.each{|d|
-                if(d<down_wisk) or (d>top_wisk)
-                y_out=(@graph_height - (d -minvalue)*unit_size).to_s
-                    @graph.add_element( "circle", {
-                     "cx" => (middle).to_s,
-                     "cy" => y_out,
-                    "r" => "3",
-                    "class" => "dataPoint#{field_count+1}"
-                    })
-                    @graph.add_element( "text", {
-                    "x" => (middle+20).to_s,
-                    "y" => y_out,
-                    "class" => "dataPointLabel",
-                    "style" => "#{style} stroke: #000;"
-                    }).text = d.to_s                    end
-            }
-        end
-        
-        @graph.add_element( "rect", { "x" => left_rect.to_s, "y" => top_rect.to_s, "width" => (bar_width / 2).to_s, "height" => (height_rect).to_s, "class" => "fill#{field_count+1}"})
-
-        @graph.add_element("path",{"d"=>path, "style"=>"stroke:black;stroke-width:2"})
-        
-        field_count += 1
-        end
-       
-        end
-    end
-  end
-end

data/lib/statsample/graph/svggraph.rb

@@ -1,184 +0,0 @@
-require 'SVG/Graph/Bar'
-require 'SVG/Graph/BarHorizontal'
-require 'SVG/Graph/Pie'
-require 'SVG/Graph/Line'
-require 'SVG/Graph/Plot'
-
-module Statsample
-  class Vector
-    # Creates a barchart using ruby-gdchart
-    def svggraph_frequencies(file, width=600, height=300, chart_type=SVG::Graph::BarNoOp, options={})
-      labels, data1=[],[]
-      self.frequencies.sort.each{|k,v|
-        labels.push(k.to_s)
-        data1.push(v)
-      }
-      options[:height]=height
-      options[:width]=width
-      options[:fields]=labels
-      graph = chart_type.new(options)
-      graph.add_data(
-      :data => data1,
-      :title => "Frequencies"
-      )
-      if file.respond_to? :write
-        file.write(graph.burn)
-      else
-        File.open(file.to_s,"wb") {|f|
-          f.puts(graph.burn)
-        }
-      end
-      
-    end
-    
-    def svggraph_histogram(bins, options={})
-      check_type :scale
-      options={:graph_title=>"Histogram", :show_graph_title=>true,:show_normal=>true, :mean=>self.mean, :sigma=>sdp }.merge! options
-      graph = Statsample::Graph::SvgHistogram.new(options)
-      graph.histogram=histogram(bins)
-      graph
-    end
-    # Returns a Run-Sequence Plot
-    # Reference: http://www.itl.nist.gov/div898/handbook/eda/section3/runseqpl.htm
-    def svggraph_runsequence_plot(options={})
-      check_type :scale
-      options={:graph_title=>"Run-Sequence Plot", :show_graph_title=>true, :scale_x_integers => true, :add_popups=>true }.merge! options
-      vx=(1..@data.size).to_a.to_vector(:scale)
-      vy=@data.to_vector(:scale)
-      ds={'index'=>vx,'value'=>vy}.to_dataset
-      graph = Statsample::Graph::SvgScatterplot.new(ds,options)
-      graph.set_x('index')
-      graph.parse
-      graph
-    end
-    def svggraph_boxplot(options={})
-      check_type :scale
-      options={:graph_title=>"Boxplot", :fields=>['vector'], :show_graph_title=>true}.merge! options
-      vx=@valid_data.to_a.to_vector(:scale)
-      graph = Statsample::Graph::SvgBoxplot.new(options)
-      graph.add_data(:title=>"vector", :data=>@data.to_a)
-      graph
-    end
-
-    def svggraph_lag_plot(options={})
-      check_type :scale
-      options={:graph_title=>"Lag Plot", :show_graph_title=>true}.merge! options
-      vx=@valid_data[0...(@valid_data.size-1)].to_vector(:scale)
-      vy=@valid_data[1...@valid_data.size].to_vector(:scale)
-      ds={'x_minus_1'=>vx,'x'=>vy}.to_dataset
-      graph = Statsample::Graph::SvgScatterplot.new(ds,options)
-      graph.set_x('x_minus_1')
-      graph.parse
-      graph
-    end
-    # Returns a Normal Probability Plot
-    # Reference: http://www.itl.nist.gov/div898/handbook/eda/section3/normprpl.htm
-    def svggraph_normalprobability_plot(options={})
-      extend Statsample::Util
-      check_type :scale
-      options={:graph_title=>"Normal Probability Plot", :show_graph_title=>true}.merge! options
-      n=@valid_data.size
-      vx=(1..@valid_data.size).to_a.collect{|i|
-        Distribution::Normal.p_value(normal_order_statistic_medians(i,n))
-      }.to_vector(:scale)
-      vy=@valid_data.sort.to_vector(:scale)
-      ds={'normal_order_statistics_medians'=>vx, 'ordered_response'=>vy}.to_dataset
-      graph = Statsample::Graph::SvgScatterplot.new(ds,options)
-      graph.set_x('normal_order_statistics_medians')
-      graph.parse
-      graph
-    end
-  end
-end
-
-# replaces all key and fill classes with similar ones, without opacity
-# this allows rendering of svg and png on rox and gqview without problems
-module SVG #:nodoc:
-  module Graph
-    class BarNoOp < Bar # :nodoc:
-      def get_css; SVG::Graph.get_css_standard; end
-    end
-    class BarHorizontalNoOp < BarHorizontal # :nodoc:
-      def get_css; SVG::Graph.get_css_standard; end
-    end
-
-    class LineNoOp < Line # :nodoc:
-      def get_css; SVG::Graph.get_css_standard; end
-
-    end
-    class PlotNoOp < Plot # :nodoc:
-      def get_css; SVG::Graph.get_css_standard; end
-    end
-    class PieNoOp < Pie # :nodoc:
-      def get_css; SVG::Graph.get_css_standard; end
-
-    end
-    class << self
-      def get_css_standard
-        return <<-EOL
-        /* default fill styles for multiple datasets (probably only use a single dataset on this graph though) */
-        .key1,.fill1{
-          fill: #ff0000;
-          stroke: none;
-          stroke-width: 0.5px;
-        }
-        .key2,.fill2{
-          fill: #0000ff;
-          stroke: none;
-          stroke-width: 1px;
-        }
-        .key3,.fill3{
-          fill: #00ff00;
-          stroke: none;
-          stroke-width: 1px;
-        }
-        .key4,.fill4{
-          fill: #ffcc00;
-          stroke: none;
-          stroke-width: 1px;
-        }
-        .key5,.fill5{
-          fill: #00ccff;
-          stroke: none;
-          stroke-width: 1px;
-        }
-        .key6,.fill6{
-          fill: #ff00ff;
-          stroke: none;
-          stroke-width: 1px;
-        }
-        .key7,.fill7{
-          fill: #00ffff;
-          stroke: none;
-          stroke-width: 1px;
-        }
-        .key8,.fill8{
-          fill: #ffff00;
-          stroke: none;
-          stroke-width: 1px;
-        }
-        .key9,.fill9{
-          fill: #cc6666;
-          stroke: none;
-          stroke-width: 1px;
-        }
-        .key10,.fill10{
-          fill: #663399;
-          stroke: none;
-          stroke-width: 1px;
-        }
-        .key11,.fill11{
-          fill: #339900;
-          stroke: none;
-          stroke-width: 1px;
-        }
-        .key12,.fill12{
-          fill: #9966FF;
-          stroke: none;
-          stroke-width: 1px;
-        }
-        EOL
-      end
-    end
-  end
-end

data/lib/statsample/graph/svghistogram.rb

@@ -1,206 +0,0 @@
-module Statsample
-  module Graph
-    class SvgHistogram < SVG::Graph::BarBase # :nodoc:
-      attr_accessor :inner_margin, :mean, :sigma, :show_normal
-      def initialize(config)
-        config[:fields]=[:dummy]
-        super(config)
-        
-        @histogram=nil
-      end
-      include REXML
-      
-      # In addition to the defaults set in Graph::initialize, sets
-      # [inner_margin] 14
-      # [key] false
-      def set_defaults 
-      super
-      self.top_align = self.top_font = 0
-      init_with({
-        :inner_margin=>16,
-        :key=>false,
-        :show_normal=>false
-        })
-      
-      end
-      
-      def histogram=(h)
-        @histogram=h
-        @data=[{:data=>(0...@histogram.bins).to_a.collect {|i|
-            @histogram[i]
-        }}]
-      end
-      def get_x_labels
-        [""]
-      end
-      
-      def get_y_labels
-        maxvalue = max_value
-        minvalue = min_value
-        range = maxvalue - minvalue
-        
-        top_pad = range == 0 ? 10 : range / 20.0
-        scale_range = (maxvalue + top_pad) - minvalue
-        
-        scale_division = scale_divisions || (scale_range / 10.0)
-        
-        if scale_integers
-        scale_division = scale_division < 1 ? 1 : scale_division.round
-        end
-        
-        rv = []
-        maxvalue = maxvalue%scale_division == 0 ? 
-        maxvalue : maxvalue + scale_division
-        minvalue.step( maxvalue, scale_division ) {|v| rv << v}
-        rv
-      end
-      
-      def unit_width
-        (@graph_width-(@inner_margin*2)) / (@histogram.max-@histogram.min).to_f
-      end
-      def draw_x_label(v)
-        left = (v - @histogram.min)*unit_width                    
-        x=@inner_margin+left
-        text = @graph.add_element( "text" )
-        text.attributes["class"] = "xAxisLabels"
-        text.text = sprintf("%0.2f",v)
-        y = @graph_height + x_label_font_size + 3
-        text.attributes["x"] = x.to_s
-        text.attributes["y"] = y.to_s
-      end
-      
-      def draw_x_labels
-        if show_x_labels
-          (0...@histogram.bins).each do |i|
-            value = @histogram[i]
-            range = @histogram.get_range(i)
-            draw_x_label(range[0])
-            if(i==(@histogram.bins)-1)
-              draw_x_label(range[1])
-            end
-          end
-        end
-      end
-      def draw_data
-        minvalue = min_value
-        fieldwidth = field_width
-        unit_size =  (@graph_height.to_f - font_size*2*top_font) / 
-                    (get_y_labels.max - get_y_labels.min)		
-        bottom = @graph_height
-        field_count = 0
-        hist_min=@histogram.min
-        hist_max=@histogram.max
-        range_hist=hist_max-hist_min
-        total=0
-        
-        (0...@histogram.bins).each do |i|
-          dataset_count = 0
-          value = @histogram[i]
-          total=total+value
-          range = @histogram.get_range(i)
-          left = (range[0] - hist_min)*unit_width
-          bar_width = (range[1] - hist_min)*unit_width - left
-          length = (value.abs - (minvalue > 0 ? minvalue : 0)) * unit_size
-          # top is 0 if value is negative
-          top = bottom - (((value < 0 ? 0 : value) - minvalue) * unit_size)
-
-          @graph.add_element( "rect", {
-            "x" => (@inner_margin+left).to_s,
-            "y" => top.to_s,
-            "width" => bar_width.to_s,
-            "height" => length.to_s,
-            "class" => "fill#{dataset_count+1}"
-          })
-
-          make_datapoint_text(left + @inner_margin+ (bar_width/2), top - 6, value.to_s)
-          field_count += 1
-        end
-        if(show_normal)
-          divs=30
-          path=""
-          0.upto(divs) do |i|
-            x_abs=hist_min+(range_hist/divs)*i
-            y=GSL::Ran::gaussian_pdf((x_abs-mean) / sigma)*total
-            xg=@inner_margin+((x_abs-hist_min)*unit_width)
-            yg=bottom-(y-minvalue)*unit_size
-            if i==0
-              path="M#{xg} #{yg} "
-            else
-              path+="L#{xg} #{yg} "
-            end
-          end
-          @graph.add_element("path", 
-            { "d"=>path, "style"=>"stroke:black;fill:none" })
-        end
-      end
-      
-      
-      def get_css
-      return <<EOL
-/* default fill styles for multiple datasets (probably only use a single dataset on this graph though) */
-
-.key1,.fill1{
-fill: #ff0000;
-stroke: black;
-stroke-width: 1px;	
-}
-.key2,.fill2{
-fill: #0000ff;
-stroke: black;
-stroke-width: 1px;	
-}
-.key3,.fill3{
-fill: #00ff00;
-stroke: none;
-stroke-width: 1px;	
-}
-.key4,.fill4{
-fill: #ffcc00;
-stroke: none;
-stroke-width: 1px;	
-}
-.key5,.fill5{
-fill: #00ccff;
-stroke: none;
-stroke-width: 1px;	
-}
-.key6,.fill6{
-fill: #ff00ff;
-stroke: none;
-stroke-width: 1px;	
-}
-.key7,.fill7{
-fill: #00ffff;
-stroke: none;
-stroke-width: 1px;	
-}
-.key8,.fill8{
-fill: #ffff00;
-stroke: none;
-stroke-width: 1px;	
-}
-.key9,.fill9{
-fill: #cc6666;
-stroke: none;
-stroke-width: 1px;	
-}
-.key10,.fill10{
-fill: #663399;
-stroke: none;
-stroke-width: 1px;	
-}
-.key11,.fill11{
-fill: #339900;
-stroke: none;
-stroke-width: 1px;	
-}
-.key12,.fill12{
-fill: #9966FF;
-stroke: none;
-stroke-width: 1px;	
-}
-EOL
-      end
-    end
-  end
-end