statsample 0.16.0 → 0.17.0

data/History.txt

@@ -1,3 +1,9 @@
+=== 0.17.0 / 2010-12-09
+* Added Statsample::Graph::Histogram and Statsample::Graph::Boxplot
+* Added Statsample::Reliability::SkillScaleAnalysis for analysis of skill based scales.
+* Delete combination and permutation clases. Backport for ruby 1.8.7 widely available
+* Deleted unused variables (thanks, ruby-head)
+
 === 0.16.0 / 2010-11-13
 * Works on ruby 1.9.2 and HEAD. Updated Rakefile and manifest
 * Removed all graph based on Svg::Graph. 
@@ -6,6 +12,7 @@
 * Added reference on references.txt
 * Ruby-based random gaussian distribution generator when gsl not available
 * Added population average deviation [Al Chou]
+
 === 0.15.1 / 2010-10-20
 * Statsample::Excel and Statsample::PlainText add name to vectors equal to field name 
 * Statsample::Dataset.delete_vector accept multiple fields.

data/Manifest.txt

@@ -12,10 +12,12 @@ data/test_binomial.csv
 data/tetmat_matrix.txt
 data/tetmat_test.txt
 doc_latex/manual/equations.tex
+examples/boxplot.rb
 examples/correlation_matrix.rb
 examples/dataset.rb
 examples/dominance_analysis.rb
 examples/dominance_analysis_bootstrap.rb
+examples/histogram.rb
 examples/icc.rb
 examples/levene.rb
 examples/multiple_regression.rb
@@ -47,7 +49,6 @@ lib/statsample/anova/twoway.rb
 lib/statsample/bivariate.rb
 lib/statsample/bivariate/pearson.rb
 lib/statsample/codification.rb
-lib/statsample/combination.rb
 lib/statsample/converter/csv.rb
 lib/statsample/converter/spss.rb
 lib/statsample/converters.rb
@@ -62,6 +63,8 @@ lib/statsample/factor/pca.rb
 lib/statsample/factor/principalaxis.rb
 lib/statsample/factor/rotation.rb
 lib/statsample/graph.rb
+lib/statsample/graph/boxplot.rb
+lib/statsample/graph/histogram.rb
 lib/statsample/graph/scatterplot.rb
 lib/statsample/histogram.rb
 lib/statsample/matrix.rb
@@ -70,7 +73,6 @@ lib/statsample/mle/logit.rb
 lib/statsample/mle/normal.rb
 lib/statsample/mle/probit.rb
 lib/statsample/multiset.rb
-lib/statsample/permutation.rb
 lib/statsample/regression.rb
 lib/statsample/regression/binomial.rb
 lib/statsample/regression/binomial/logit.rb
@@ -86,6 +88,7 @@ lib/statsample/reliability.rb
 lib/statsample/reliability/icc.rb
 lib/statsample/reliability/multiscaleanalysis.rb
 lib/statsample/reliability/scaleanalysis.rb
+lib/statsample/reliability/skillscaleanalysis.rb
 lib/statsample/resample.rb
 lib/statsample/rserve_extension.rb
 lib/statsample/srs.rb
@@ -111,7 +114,6 @@ test/test_anovawithvectors.rb
 test/test_bartlettsphericity.rb
 test/test_bivariate.rb
 test/test_codification.rb
-test/test_combination.rb
 test/test_crosstab.rb
 test/test_csv.csv
 test/test_csv.rb
@@ -126,10 +128,10 @@ test/test_logit.rb
 test/test_matrix.rb
 test/test_mle.rb
 test/test_multiset.rb
-test/test_permutation.rb
 test/test_regression.rb
 test/test_reliability.rb
 test/test_reliability_icc.rb
+test/test_reliability_skillscale.rb
 test/test_resample.rb
 test/test_rserve_extension.rb
 test/test_srs.rb

data/README.txt

@@ -21,6 +21,7 @@ Include:
 * Sample calculation related formulas
 * Structural Equation Modeling (SEM), using R libraries +sem+ and +OpenMx+
 * Creates reports on text, html and rtf, using ReportBuilder gem
+* Graphics: Histogram, Boxplot and Scatterplot
 
 == FEATURES:
 
@@ -69,8 +70,11 @@ Include:
   * Statsample::Test::UMannWhitney
   * Statsample::Test::T
   * Statsample::Test::F  
+* Module Graph provides several classes to create beautiful graphs using rubyvis
+  * Statsample::Graph::Boxplot
+  * Statsample::Graph::Histogram
+  * Statsample::Graph::Scatterplot
 * Gem +statsample-sem+ provides a DSL to R libraries +sem+ and +OpenMx+
-* Interfaces to gdchart, gnuplot and SVG::Graph (experimental)
 * Close integration with gem <tt>reportbuilder</tt>, to easily create reports on text, html and rtf formats.
 
 == Examples of use:

data/Rakefile

@@ -41,7 +41,7 @@ h=Hoe.spec('statsample') do
   #self.testlib=:minitest
 	self.rubyforge_name = "ruby-statsample"
 	self.developer('Claudio Bustos', 'clbustos@gmail.com')
-	self.extra_deps << ["spreadsheet","~>0.6.0"] <<  ["reportbuilder", "~>1.0"] << ["minimization", "~>0.2.0"] << ["fastercsv"] << ["dirty-memoize", "~>0.0"] << ["extendmatrix","~>0.3.1"] << ["statsample-bivariate-extension", ">0"] << ["rserve-client", "~>0.2.5"] << ["rubyvis", "~>0.2.2"]
+	self.extra_deps << ["spreadsheet","~>0.6.0"] <<  ["reportbuilder", "~>1.4"] << ["minimization", "~>0.2.0"] << ["fastercsv", ">0"] << ["dirty-memoize", "~>0.0"] << ["extendmatrix","~>0.3.1"] << ["statsample-bivariate-extension", ">0"] << ["rserve-client", "~>0.2.5"] << ["rubyvis", "~>0.3.3"]
   
 	self.extra_dev_deps << ["shoulda"] << ["minitest", "~>2.0"]
   self.clean_globs << "test/images/*" << "demo/item_analysis/*" << "demo/Regression"

data/examples/boxplot.rb

@@ -0,0 +1,17 @@
+#!/usr/bin/ruby
+$:.unshift(File.dirname(__FILE__)+'/../lib/')
+$:.unshift('/home/cdx/dev/reportbuilder/lib/')
+
+require 'benchmark'
+require 'statsample'
+n=100
+a=(n-1).times.map {|i| rand()*20+50}
+b=n.times.map {|i| rand()*10+50}.to_scale
+c=n.times.map {|i| rand()*5+50}.to_scale
+ 
+a.push(30)
+a=a.to_scale
+sp=Statsample::Graph::Boxplot.new(:vectors=>[a,b,c],:width=>300, :height=>300, :groups=>%w{first first second}, :minimum=>0)
+rb=ReportBuilder.new
+rb.add(sp)
+puts rb.to_text

data/examples/dominance_analysis_bootstrap.rb

@@ -8,6 +8,11 @@ b=100.times.collect {rand}.to_scale
 c=100.times.collect {rand}.to_scale
 d=100.times.collect {rand}.to_scale
 
+a.name="a"
+b.name="b"
+c.name="c"
+d.name="d"
+
 ds={'a'=>a,'b'=>b,'c'=>c,'d'=>d}.to_dataset
 
 ds['y1']=ds.collect{|row| row['a']*5+row['b']*2+row['c']*2+row['d']*2+10*rand()}

data/examples/histogram.rb

@@ -0,0 +1,14 @@
+#!/usr/bin/ruby
+$:.unshift(File.dirname(__FILE__)+'/../lib/')
+$:.unshift('/home/cdx/dev/reportbuilder/lib/')
+
+require 'benchmark'
+require 'statsample'
+n=1000
+a=n.times.map {|i| rand()*20}.to_scale
+hg=Statsample::Graph::Histogram.new(a, :bins=>15)
+
+rb=ReportBuilder.new
+rb.add(a.histogram)
+rb.add(hg)
+puts rb.to_text

data/examples/scatterplot.rb

@@ -1,12 +1,13 @@
 #!/usr/bin/ruby
 $:.unshift(File.dirname(__FILE__)+'/../lib/')
+$:.unshift('/home/cdx/dev/reportbuilder/lib/')
+
 require 'benchmark'
 require 'statsample'
 n=100
 a=n.times.map {|i| rand(10)+i}.to_scale
 b=n.times.map {|i| rand(10)+i}.to_scale
 sp=Statsample::Graph::Scatterplot.new(a,b, :width=>200, :height=>200)
-rb=ReportBuilder.new do |b|
-  b.parse_element(sp)
-end  
+rb=ReportBuilder.new
+rb.add(sp)
 puts rb.to_text

data/lib/distribution/normalbivariate.rb

@@ -220,7 +220,7 @@ module Distribution
             asr = Math::asin(r)
             (1..lg).each do |i|
               [-1,1].each do |is|
-                sn = Math::sin( asr *(  is * x[i][ng] + 1 ).quo(2) )
+                sn = Math::sin(asr*(is* x[i][ng]+1).quo(2) )
                 bvn = bvn + w[i][ng] * Math::exp( ( sn*hk-hs ).quo( 1-sn*sn ) )
               end # do
             end # do

data/lib/statsample.rb

@@ -118,12 +118,10 @@ module Statsample
     @@has_gsl
   end
   
-  VERSION = '0.16.0'
+  VERSION = '0.17.0'
   SPLIT_TOKEN = ","
   autoload(:Database, 'statsample/converters')
   autoload(:Anova, 'statsample/anova')
-  autoload(:Combination, 'statsample/combination')
-  autoload(:Permutation, 'statsample/permutation')
   autoload(:CSV, 'statsample/converters')
   autoload(:PlainText, 'statsample/converters')
   autoload(:Excel, 'statsample/converters')
@@ -219,6 +217,21 @@ module Statsample
       end
       u
     end
+    
+    def self.nice(s,e) # :nodoc:
+      reverse = e<s
+      min = reverse ? e : s
+      max = reverse ? s : e
+      span=max-min
+      return [s, e] if (!span or (span.respond_to? :infinite? and span.infinite?))
+      
+      step=10**((Math::log(span).quo(Math::log(10))).round - 1).to_f  
+      out=[(min.quo(step)).floor * step, (max.quo(step)).ceil * step]
+      out.reverse! if reverse
+      out
+    end
+    
+    
   end
   
   

data/lib/statsample/bivariate.rb

@@ -1,5 +1,7 @@
 require 'statsample/bivariate/pearson'
 
+
+
 module Statsample
   # Diverse methods and classes to calculate bivariate relations
   # Specific classes: 
@@ -7,6 +9,8 @@ module Statsample
   # * Statsample::Bivariate::Tetrachoric : Tetrachoric correlation
   # * Statsample::Bivariate::Polychoric  : Polychoric correlation (using joint, two-step and polychoric series)
   module Bivariate
+    autoload(:Polychoric, 'statsample/bivariate/polychoric')
+    autoload(:Tetrachoric, 'statsample/bivariate/tetrachoric')
     
     class << self
       # Covariance between two vectors
@@ -335,6 +339,4 @@ module Statsample
   end
 end
 
-require 'statsample/bivariate/polychoric'
-require 'statsample/bivariate/tetrachoric'
 

data/lib/statsample/converter/csv.rb

@@ -1,5 +1,3 @@
-
-
 module Statsample
   class CSV < SpreadsheetBase
     if RUBY_VERSION<"1.9"

data/lib/statsample/converters.rb

@@ -181,8 +181,20 @@ module Statsample
       # USE:
       #     ds = Statsample::Excel.read("test.xls")
       #
-      def read(filename, worksheet_id=0, ignore_lines=0, empty=[''])
+      def read(filename, opts=Hash.new)
         require 'spreadsheet'
+        opts_default={
+          :worksheet_id=>0, 
+          :ignore_lines=>0, 
+          :empty=>['']
+        }
+        
+        opts=opts_default.merge opts
+        
+        worksheet_id=opts[:worksheet_id]
+        ignore_lines=opts[:ignore_lines]
+        empty=opts[:empty]
+        
         first_row=true
         fields=[]
         fields_data={}

data/lib/statsample/dataset.rb

@@ -331,7 +331,7 @@ module Statsample
     def bootstrap(n=nil)
       n||=@cases
       ds_boot=dup_empty
-      for i in 1..n
+      n.times do
         ds_boot.add_case_array(case_as_array(rand(n)))
       end
       ds_boot.update_valid_data
@@ -418,7 +418,6 @@ module Statsample
     # Returns a vector with sumatory of fields
     # if fields parameter is empty, sum all fields 
     def vector_sum(fields=nil)
-      a=[]
       fields||=@fields
       collect_with_index do |row, i|
         if(fields.find{|f| !@vectors[f].data_with_nils[i]})
@@ -557,6 +556,7 @@ module Statsample
         raise DatasetException.new(self, e)
       end
     end
+    
     # Returns each case as an array, coding missing values as nils
     def each_array_with_nils
       m=fields.size
@@ -586,8 +586,9 @@ module Statsample
       @fields=f
       check_order
     end
-    
-    def check_order
+    # Check congruence between +fields+ attribute
+    # and keys on +vectors
+    def check_order #:nodoc:
       if(@vectors.keys.sort!=@fields.sort)
         @fields=@fields&@vectors.keys
         @fields+=@vectors.keys.sort-@fields
@@ -598,7 +599,7 @@ module Statsample
       if i.is_a? Range
         fields=from_to(i.begin,i.end)
         vectors=fields.inject({}) {|a,v| a[v]=@vectors[v];a}
-        ds=Dataset.new(vectors,fields)
+        Dataset.new(vectors,fields)
       else
         raise Exception,"Vector '#{i}' doesn't exists on dataset" unless @vectors.has_key?(i)
         @vectors[i]
@@ -613,7 +614,7 @@ module Statsample
       }
       Statsample::Vector.new(data,type)
     end
-    # Same as #collect, but giving case index as second parameter on yield.
+    # Same as Statsample::Vector.collect, but giving case index as second parameter on yield.
     def collect_with_index(type=:scale)
       data=[]
       each_with_index {|row, i|
@@ -661,14 +662,7 @@ module Statsample
       end
     end
 		
-    def to_multiset_by_split(*fields)
-			require 'statsample/multiset'
-			if fields.size==1
-				to_multiset_by_split_one_field(fields[0])
-			else
-				to_multiset_by_split_multiple_fields(*fields)
-			end
-    end
+   
     
     # Create a new dataset with all cases which the block returns true
     def filter
@@ -689,6 +683,20 @@ module Statsample
       a.to_vector(@vectors[field].type)
     end
     
+    # Creates a Stastample::Multiset, using one or more fields
+    # to split the dataset.
+    
+   
+    def to_multiset_by_split(*fields)
+			require 'statsample/multiset'
+			if fields.size==1
+				to_multiset_by_split_one_field(fields[0])
+			else
+				to_multiset_by_split_multiple_fields(*fields)
+			end
+    end
+    # Creates a Statsample::Multiset, using one field
+    
     def to_multiset_by_split_one_field(field)
       raise ArgumentError,"Should use a correct field name" if !@fields.include? field
       factors=@vectors[field].factors
@@ -831,7 +839,7 @@ module Statsample
     #      ]
     # 
     def one_to_many(parent_fields, pattern)
-      base_pattern=pattern.gsub(/%v|%n/,"")
+      #base_pattern=pattern.gsub(/%v|%n/,"")
       re=Regexp.new pattern.gsub("%v","(.+?)").gsub("%n","(\\d+?)")
       ds_vars=parent_fields
       vars=[]

data/lib/statsample/dominanceanalysis.rb

@@ -156,8 +156,9 @@ module Statsample
       @models=[]
       @models_data={}
       for i in 1..@predictors.size
-        c=Statsample::Combination.new(i,@predictors.size)
+        c=(0...@predictors.size).to_a.combination(i)
         c.each  do |data|
+          
           independent=data.collect {|i1| @predictors[i1] }
           @models.push(independent)
           if (@build_from_dataset)
@@ -268,7 +269,7 @@ module Statsample
     # Get all model of size k
     def md_k(k)
       out=[]
-      models=@models.each{|m| out.push(md(m)) if m.size==k }
+      @models.each{|m| out.push(md(m)) if m.size==k }
       out
     end
     

data/lib/statsample/dominanceanalysis/bootstrap.rb

@@ -158,8 +158,9 @@ module Statsample
         @samples_cd={}
         @samples_gd={}
         @pairs=[]
-        c=Statsample::Combination.new(2,@fields.size)
+        c=(0...@fields.size).to_a.combination(2)
         c.each do |data|
+          p data
           convert=data.collect {|i| @fields[i] }
           @pairs.push(convert)
           [@samples_td, @samples_cd, @samples_gd].each{|s|

data/lib/statsample/factor/parallelanalysis.rb

@@ -27,7 +27,7 @@ module Statsample
         ds.fields=vars.times.map {|i| "v#{i+1}"}
         ds.cases=cases
         opts=opts.merge({:bootstrap_method=> :random, :no_data=>true})
-        pa=new(ds, opts)
+        new(ds, opts)
       end
       include DirtyMemoize
       include Summarizable

data/lib/statsample/factor/principalaxis.rb

@@ -132,7 +132,7 @@ module Factor
         @communalities=pca.communalities(m)
         @eigenvalues=pca.eigenvalues
         com_sum = @communalities.inject(0) {|ac,v| ac+v}
-        jump=true
+        #jump=true
         
         break if (com_sum-prev_sum).abs < @delta
         @communalities.each_with_index do |v2,i2|

data/lib/statsample/graph.rb

@@ -1,4 +1,6 @@
 require 'statsample/graph/scatterplot'
+require 'statsample/graph/boxplot'
+require 'statsample/graph/histogram'
 module Statsample
   module Graph
   end

data/lib/statsample/graph/boxplot.rb

@@ -0,0 +1,234 @@
+require 'rubyvis'
+module Statsample
+  module Graph
+    # = Boxplot
+    # 
+    # From Wikipedia:
+    # In descriptive statistics, a box plot or boxplot (also known as a box-and-whisker diagram or plot) is a convenient way of graphically depicting groups of numerical data through their five-number summaries: the smallest observation (sample minimum), lower quartile (Q1), median (Q2), upper quartile (Q3), and largest observation (sample maximum). A boxplot may also indicate which observations, if any, might be considered outliers.
+    # 
+    # == Usage
+    # === Svg output
+    #  a=[1,2,3,4].to_scale
+    #  b=[3,4,5,6].to_scale
+    #  puts Statsample::Graph::Boxplot.new(:vectors=>[a,b]).to_svg
+    # === Using ReportBuilder
+    #  a=[1,2,3,4].to_scale
+    #  b=[3,4,5,6].to_scale
+    #  rb=ReportBuilder.new
+    #  rb.add(Statsample::Graph::Boxplot.new(:vectors=>[a,b]))
+    #  rb.save_html('boxplot.html')
+    
+    class Boxplot
+      include Summarizable
+      attr_accessor :name
+      # Total width of Boxplot
+      attr_accessor :width
+      # Total height of Boxplot
+      attr_accessor :height
+      # Top margin
+      attr_accessor :margin_top
+      # Bottom margin
+      attr_accessor :margin_bottom
+      # Left margin
+      attr_accessor :margin_left
+      # Right margin
+      attr_accessor :margin_right
+      # Array with assignation to groups of bars
+      # For example, for four vectors, 
+      #   boxplot.groups=[1,2,1,3]
+      # Assign same color to first and third element, and different to
+      # second and fourth
+      attr_accessor :groups
+      # Minimum value on y-axis. Automaticly defined from data
+      attr_accessor :minimum
+      # Maximum value on y-axis. Automaticly defined from data
+      attr_accessor :maximum
+      # Vectors to box-ploting
+      attr_accessor :vectors
+      
+      attr_reader :x_scale, :y_scale
+      # Create a new Boxplot.
+      # Parameters: Hash of options
+      # * :vectors: Array of vectors
+      # * :groups: Array of same size as :vectors:, with name of groups
+      #           to colorize vectors
+      def initialize(opts=Hash.new)
+        @vectors=opts.delete :vectors
+        raise "You should define vectors" if @vectors.nil?
+        
+        opts_default={
+          :name=>_("Boxplot"),
+          :groups=>nil,
+          :width=>400,
+          :height=>300,
+          :margin_top=>10,
+          :margin_bottom=>20,
+          :margin_left=>20,
+          :margin_right=>20,
+          :minimum=>nil,
+          :maximum=>nil
+        }
+        @opts=opts_default.merge(opts)
+        opts_default.keys.each {|k| send("#{k}=", @opts[k]) }
+      end
+      
+      # Returns a Rubyvis panel with scatterplot
+      def rubyvis_panel # :nodoc:
+        that=self
+        
+        min,max=@minimum, @maximum
+        
+        min||=@vectors.map {|v| v.min}.min
+        max||=@vectors.map {|v| v.max}.max
+        
+        
+        
+        margin_hor=margin_left + margin_right
+        margin_vert=margin_top  + margin_bottom
+        x_scale = pv.Scale.ordinal(@vectors.size.times.map.to_a).split_banded(0, width-margin_hor, 4.0/5)
+        y_scale=Rubyvis::Scale.linear(min,max).range(0,height-margin_vert)
+        y_scale.nice
+        # cache data
+        
+        colors=Rubyvis::Colors.category10
+        
+        data=@vectors.map {|v|
+          out={:percentil_25=>v.percentil(25), :median=>v.median, :percentil_75=>v.percentil(75), :name=>v.name}
+          out[:iqr]=out[:percentil_75]-out[:percentil_25]
+          
+          irq_max=out[:percentil_75]+out[:iqr]
+          irq_min=out[:percentil_25]-out[:iqr]
+          
+          # Find the last data inside the margin
+          min=out[:percentil_25]
+          max=out[:percentil_75]
+          
+          v.each {|d|
+            min=d if d<min and d>irq_min
+            max=d if d>max and d<irq_max
+          }
+          # Whiskers!
+          out[:low_whisker]=min
+          out[:high_whisker]=max
+          # And now, data outside whiskers
+          out[:outliers]=v.data_with_nils.find_all {|d| 
+            d<min or d>max
+          }
+          out
+        }
+        
+
+        
+        vis=Rubyvis::Panel.new do |pan| 
+          pan.width  width  - margin_hor
+          pan.height height - margin_vert
+          pan.bottom margin_bottom
+          pan.left   margin_left
+          pan.right  margin_right
+          pan.top    margin_top
+           # Y axis
+          pan.rule do
+            data y_scale.ticks
+            bottom y_scale
+            stroke_style {|d| d!=0 ? "#eee" : "#000"}
+            label(:anchor=>'left') do
+              visible {|d| true}
+              text y_scale.tick_format
+            end
+          end
+          pan.rule do
+            bottom 0
+            stroke_style 'black'
+          end
+          pan.label  do |l|
+            l.data data
+            l.left  {|v|  x_scale.scale(index)}
+            l.bottom -15
+            l.text {|v,x| v[:name]}
+          end
+          
+          pan.panel do |bp|
+            bp.data data
+            bp.left {|v|  x_scale.scale(index)}
+            bp.width x_scale.range_band
+            
+            
+            # Bar
+            bp.bar do |b|
+              b.bottom {|v| y_scale.scale(v[:percentil_25])}
+              b.height {|v| y_scale.scale(v[:percentil_75]) - y_scale.scale(v[:percentil_25]) }
+              b.line_width 1
+              b.stroke_style  {|v| 
+                if that.groups
+                  colors.scale(that.groups[parent.index]).darker
+                else
+                  colors.scale(index).darker
+                end
+              
+              
+              }
+              b.fill_style {|v| 
+                if that.groups
+                  colors.scale(that.groups[parent.index])
+                else
+                  colors.scale(index)
+                end
+              }
+            end
+            # Median
+            bp.rule do |r|
+              r.bottom {|v| y_scale.scale(v[:median])}
+              r.width x_scale.range_band
+              r.line_width 2
+            end
+            
+            # Whiskeys
+            bp.rule do |r|
+              r.visible {|v| v[:percentil_25]>v[:low_whisker]}
+              r.bottom {|v| y_scale.scale(v[:low_whisker])}              
+            end
+            bp.rule do |r|
+              r.visible {|v| v[:percentil_25]>v[:low_whisker]}
+              r.bottom {|v| y_scale.scale(v[:low_whisker])}              
+              r.left {|v| x_scale.range_band / 2.0}
+              r.height {|v| y_scale.scale(v[:percentil_25])-y_scale.scale(v[:low_whisker])}
+            end
+            bp.rule do |r|
+              r.visible {|v| v[:percentil_75]<v[:high_whisker]}
+              r.bottom {|v| y_scale.scale(v[:high_whisker])}              
+            end
+            
+             bp.rule do |r|
+              r.visible {|v| v[:percentil_75]<v[:high_whisker]}
+              r.bottom {|v| y_scale.scale(v[:percentil_75])}              
+              r.left {|v| x_scale.range_band / 2.0}
+              r.height {|v| y_scale.scale(v[:high_whisker])-y_scale.scale(v[:percentil_75])}
+            end
+            
+            bp.dot do |dot|
+              dot.shape_size 4
+              dot.data {|v| v[:outliers]}
+              dot.left {|v| x_scale.range_band / 2.0}
+              dot.bottom {|v| y_scale.scale(v)}
+              dot.title {|v| v}
+            end
+            
+            
+          end          
+        end
+      end
+      # Returns SVG with scatterplot
+      def to_svg
+        rp=rubyvis_panel
+        rp.render
+        rp.to_svg
+      end
+      def report_building(builder) # :nodoc:
+        builder.section(:name=>name) do |b|
+          b.image(to_svg, :type=>'svg', :width=>width, :height=>height)
+        end
+        
+      end
+    end
+  end
+end