sapor 0.1b1

data/spec/unit/first_past_the_post_spec.rb

@@ -0,0 +1,53 @@
+# encoding: utf-8
+#
+# Statistical Analysis of Polling Results (SAPoR)
+# Copyright (C) 2014 Filip van Laenen <f.a.vanlaenen@ieee.org>
+#
+# This file is part of SAPoR.
+#
+# SAPoR is free software: you can redistribute it and/or modify it under the
+# terms of the GNU General Public License as published by the Free Software
+# Foundation, either version 3 of the License, or (at your option) any later
+# version.
+#
+# SAPoR is distributed in the hope that it will be useful, but WITHOUT ANY
+# WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
+# A PARTICULAR PURPOSE.  See the GNU General Public License for more details.
+#
+# You can find a copy of the GNU General Public License in /doc/gpl.txt
+#
+
+SAMPLE_RG_ELECTION_RESULT = { 'Red' => 91_811, 'Green' => 190_934 }
+
+SAMPLE_DETAILED_ELECTION_RESULT = { 'North' => { 'Red' => 50, 'Green' => 70
+                                             },
+                                    'South' => { 'Red' => 70, 'Green' => 50,
+                                                 'Blue' => 100 } }
+
+FPTP = Sapor::FirstPastThePost.new(SAMPLE_RG_ELECTION_RESULT,
+                                   SAMPLE_DETAILED_ELECTION_RESULT)
+
+describe Sapor::FirstPastThePost, '#project' do
+  it 'projects same result as last result if fed with last election result' do
+    projection = FPTP.project(SAMPLE_RG_ELECTION_RESULT)
+    expect(projection['Red']).to eq(0)
+    expect(projection['Green']).to eq(1)
+    expect(projection['Other']).to eq(1)
+  end
+
+  it 'it moves a seat to Red if Red is polling well enough' do
+    projection = FPTP.project('Red' => 91_811 * 70 / 49,
+                              'Green' => 190_934)
+    expect(projection['Red']).to eq(1)
+    expect(projection['Green']).to eq(0)
+    expect(projection['Other']).to eq(1)
+  end
+
+  it 'it moves a seat away from Other if Red is polling well enough' do
+    projection = FPTP.project('Red' => 91_811 * 101 / 70,
+                              'Green' => 190_934 * 70 / 101)
+    expect(projection['Red']).to eq(2)
+    expect(projection['Green']).to eq(0)
+    expect(projection['Other']).to eq(0)
+  end
+end

data/spec/unit/leveled_proportional_spec.rb

@@ -0,0 +1,51 @@
+# encoding: utf-8
+#
+# Statistical Analysis of Polling Results (SAPoR)
+# Copyright (C) 2014 Filip van Laenen <f.a.vanlaenen@ieee.org>
+#
+# This file is part of SAPoR.
+#
+# SAPoR is free software: you can redistribute it and/or modify it under the
+# terms of the GNU General Public License as published by the Free Software
+# Foundation, either version 3 of the License, or (at your option) any later
+# version.
+#
+# SAPoR is distributed in the hope that it will be useful, but WITHOUT ANY
+# WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
+# A PARTICULAR PURPOSE.  See the GNU General Public License for more details.
+#
+# You can find a copy of the GNU General Public License in /doc/gpl.txt
+#
+
+require 'spec_helper'
+
+SAMPLE_ELECTION_RESULT = { 'Red' => 120, 'Green' => 120, 'Blue' => 100 }
+
+SAMPLE_DETAILED_ELECTION_RESULT = { 'North' => { 'Red' => 50, 'Green' => 70
+                                             },
+                                    'South' => { 'Red' => 70, 'Green' => 50,
+                                                 'Blue' => 100 } }
+
+SAMPLE_SEAT_DISTRIBUTION = { 'North' => 3, 'South' => 5 }
+
+SAMPLE_LEVELING_SEATS = 2
+
+SAMPLE_LEVELING_THRESHOLD = 0.05
+
+SAMPLE_POLL_RESULT = { 'Red' => 120, 'Green' => 120 }
+
+LEVELED = Sapor::LeveledProportional.new(SAMPLE_ELECTION_RESULT,
+                                         SAMPLE_DETAILED_ELECTION_RESULT,
+                                         SAMPLE_SEAT_DISTRIBUTION,
+                                         SAMPLE_LEVELING_SEATS,
+                                         SAMPLE_LEVELING_THRESHOLD,
+                                         Sapor::SainteLague14Denominators)
+
+describe Sapor::LeveledProportional, '#project' do
+  it 'projects same result as last result if fed with last election result' do
+    projection = LEVELED.project(SAMPLE_POLL_RESULT)
+    expect(projection['Red']).to eq(4)
+    expect(projection['Green']).to eq(4)
+    expect(projection['Blue']).to eq(2)
+  end
+end

data/spec/unit/norway_spec.rb

@@ -0,0 +1,47 @@
+# encoding: utf-8
+#
+# Statistical Analysis of Polling Results (SAPoR)
+# Copyright (C) 2014 Filip van Laenen <f.a.vanlaenen@ieee.org>
+#
+# This file is part of SAPoR.
+#
+# SAPoR is free software: you can redistribute it and/or modify it under the
+# terms of the GNU General Public License as published by the Free Software
+# Foundation, either version 3 of the License, or (at your option) any later
+# version.
+#
+# SAPoR is distributed in the hope that it will be useful, but WITHOUT ANY
+# WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
+# A PARTICULAR PURPOSE.  See the GNU General Public License for more details.
+#
+# You can find a copy of the GNU General Public License in /doc/gpl.txt
+#
+
+require 'spec_helper'
+
+describe Sapor::Norway, '#seats' do
+  it 'calculates the number of seats for the election of 2013 correctly' do
+    seats = Sapor::Norway.instance.seats(Sapor::Norway::LAST_ELECTION_RESULT)
+    expect(seats['Arbeiderpartiet']).to eq(55)
+    expect(seats['Høyre']).to eq(48)
+    expect(seats['Fremskrittspartiet']).to eq(29)
+    expect(seats['Kristelig Folkeparti']).to eq(10)
+    expect(seats['Senterpartiet']).to eq(10)
+    expect(seats['Venstre']).to eq(9)
+    expect(seats['Sosialistisk Venstreparti']).to eq(7)
+    expect(seats['Miljøpartiet de Grønne']).to eq(1)
+    expect(seats['Rødt']).to eq(0)
+    expect(seats['De Kristne']).to eq(0)
+    expect(seats['Pensjonistpartiet']).to eq(0)
+    expect(seats['Piratpartiet']).to eq(0)
+    expect(seats['Kystpartiet']).to eq(0)
+    expect(seats['Demokratene i Norge']).to eq(0)
+    expect(seats['Kristent Samlingsparti']).to eq(0)
+    expect(seats['Det Liberale Folkepartiet']).to eq(0)
+    expect(seats['Norges Kommunistiske Parti']).to eq(0)
+    expect(seats['Samfunnspartiet']).to eq(0)
+  end
+end
+
+# TODO: Sanity check that LAST_DETAILED_ELECTION_RESULT sums up to
+#       LAST_ELECTION_RESULT

data/spec/unit/number_formatter_spec.rb

@@ -0,0 +1,173 @@
+# encoding: utf-8
+#
+# Statistical Analysis of Polling Results (SAPoR)
+# Copyright (C) 2014 Filip van Laenen <f.a.vanlaenen@ieee.org>
+#
+# This file is part of SAPoR.
+#
+# SAPoR is free software: you can redistribute it and/or modify it under the
+# terms of the GNU General Public License as published by the Free Software
+# Foundation, either version 3 of the License, or (at your option) any later
+# version.
+#
+# SAPoR is distributed in the hope that it will be useful, but WITHOUT ANY
+# WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
+# A PARTICULAR PURPOSE.  See the GNU General Public License for more details.
+#
+# You can find a copy of the GNU General Public License in /doc/gpl.txt
+#
+
+require 'spec_helper'
+
+#
+# Test class that mixes in the PercentageFormatter.
+#
+class Formatter
+  include Sapor::NumberFormatter
+end
+
+describe Sapor::NumberFormatter, '#three_digits_percentage' do
+  it 'formats 10 as 1000%' do
+    output = Formatter.new.three_digits_percentage(10)
+    expect(output).to eq('1000%')
+  end
+
+  it 'formats 1 as 100%' do
+    output = Formatter.new.three_digits_percentage(1)
+    expect(output).to eq('100%')
+  end
+
+  it 'formats 0.9995 as 100%' do
+    output = Formatter.new.three_digits_percentage(0.9995)
+    expect(output).to eq('100%')
+  end
+
+  it 'formats 0.99949 as 99.9%' do
+    output = Formatter.new.three_digits_percentage(0.99949)
+    expect(output).to eq('99.9%')
+  end
+
+  it 'formats 0.09995 as 10.0%' do
+    output = Formatter.new.three_digits_percentage(0.09995)
+    expect(output).to eq('10.0%')
+  end
+
+  it 'formats 0.099949 as 9.99%' do
+    output = Formatter.new.three_digits_percentage(0.099949)
+    expect(output).to eq('9.99%')
+  end
+
+  it 'formats 0.009995 as 1.00%' do
+    output = Formatter.new.three_digits_percentage(0.009995)
+    expect(output).to eq('1.00%')
+  end
+
+  it 'formats 0.0099949 as 0.999%' do
+    output = Formatter.new.three_digits_percentage(0.0099949)
+    expect(output).to eq('0.999%')
+  end
+
+  it 'formats 0.000005 as 0.001%' do
+    output = Formatter.new.three_digits_percentage(0.000005)
+    expect(output).to eq('0.001%')
+  end
+
+  it 'formats 0.0000049 as 0%' do
+    output = Formatter.new.three_digits_percentage(0.0000049)
+    expect(output).to eq('0%')
+  end
+
+  it 'formats 0 as 0%' do
+    output = Formatter.new.three_digits_percentage(0)
+    expect(output).to eq('0%')
+  end
+end
+
+describe Sapor::NumberFormatter, '#six_char_percentage' do
+  it 'formats 1 as 100.0%' do
+    output = Formatter.new.six_char_percentage(1)
+    expect(output).to eq('100.0%')
+  end
+
+  it 'formats 0.9995 as 100.0%' do
+    output = Formatter.new.six_char_percentage(0.9995)
+    expect(output).to eq('100.0%')
+  end
+
+  it 'formats 0.99949 as 99.9%' do
+    output = Formatter.new.six_char_percentage(0.99949)
+    expect(output).to eq(' 99.9%')
+  end
+
+  it 'formats 0.0995 as 10.0%' do
+    output = Formatter.new.six_char_percentage(0.0995)
+    expect(output).to eq(' 10.0%')
+  end
+
+  it 'formats 0.09949 as 9.9%' do
+    output = Formatter.new.six_char_percentage(0.09949)
+    expect(output).to eq('  9.9%')
+  end
+
+  it 'formats 0.00951 as 1.0%' do
+    output = Formatter.new.six_char_percentage(0.00951)
+    expect(output).to eq('  1.0%')
+  end
+
+  it 'formats 0.00949 as 0.9%' do
+    output = Formatter.new.six_char_percentage(0.00949)
+    expect(output).to eq('  0.9%')
+  end
+
+  it 'formats 0.0005 as 0.1%' do
+    output = Formatter.new.six_char_percentage(0.0005)
+    expect(output).to eq('  0.1%')
+  end
+
+  it 'formats 0.00049 as 0.0%' do
+    output = Formatter.new.six_char_percentage(0.00049)
+    expect(output).to eq('  0.0%')
+  end
+
+  it 'formats 0 as 0.0%' do
+    output = Formatter.new.six_char_percentage(0)
+    expect(output).to eq('  0.0%')
+  end
+end
+
+describe Sapor::NumberFormatter, '#with_thousands_separator' do
+  it 'formats 0 as 0' do
+    output = Formatter.new.with_thousands_separator(0)
+    expect(output).to eq('0')
+  end
+
+  it 'formats 1 as 1' do
+    output = Formatter.new.with_thousands_separator(1)
+    expect(output).to eq('1')
+  end
+
+  it 'formats 999 as 999' do
+    output = Formatter.new.with_thousands_separator(999)
+    expect(output).to eq('999')
+  end
+
+  it 'formats 1_000 as 1,000' do
+    output = Formatter.new.with_thousands_separator(1_000)
+    expect(output).to eq('1,000')
+  end
+
+  it 'formats 999_999 as 999,999' do
+    output = Formatter.new.with_thousands_separator(999_999)
+    expect(output).to eq('999,999')
+  end
+
+  it 'formats 1_000_000 as 1,000,000' do
+    output = Formatter.new.with_thousands_separator(1_000_000)
+    expect(output).to eq('1,000,000')
+  end
+
+  it 'formats 999_999_999 as 999,999,999' do
+    output = Formatter.new.with_thousands_separator(999_999_999)
+    expect(output).to eq('999,999,999')
+  end
+end

data/spec/unit/poll_spec.rb

@@ -0,0 +1,105 @@
+# encoding: utf-8
+#
+# Statistical Analysis of Polling Results (SAPoR)
+# Copyright (C) 2014 Filip van Laenen <f.a.vanlaenen@ieee.org>
+#
+# This file is part of SAPoR.
+#
+# SAPoR is free software: you can redistribute it and/or modify it under the
+# terms of the GNU General Public License as published by the Free Software
+# Foundation, either version 3 of the License, or (at your option) any later
+# version.
+#
+# SAPoR is distributed in the hope that it will be useful, but WITHOUT ANY
+# WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
+# A PARTICULAR PURPOSE.  See the GNU General Public License for more details.
+#
+# You can find a copy of the GNU General Public License in /doc/gpl.txt
+#
+
+require 'spec_helper'
+
+SAMPLE_POLL_ARRAY = ['foo=bar', 'baz=qux', '==', 'a=1', 'b=2']
+SAMPLE_METADATA = { 'Area' => 'UT' }
+SAMPLE_RESULTS = { 'Red' => 1, 'Green' => 2, 'Blue' => 3, 'Other' => 1 }
+MAX_ERROR = 0.05
+MAX_VALUE_ERROR = 50_000
+
+def sample_poll
+  Sapor::Poll.new(SAMPLE_METADATA.dup, SAMPLE_RESULTS.dup)
+end
+
+describe Sapor::Poll, '#as_hashes' do
+  it 'converts an array of lines into an array with two hashes' do
+    hashes = Sapor::Poll.as_hashes(SAMPLE_POLL_ARRAY)
+    expect(hashes.size).to eq(2)
+  end
+
+  it 'extracts the metadata into the first array' do
+    hashes = Sapor::Poll.as_hashes(SAMPLE_POLL_ARRAY)
+    expect(hashes[0]).to eq('foo' => 'bar', 'baz' => 'qux')
+  end
+
+  it 'extracts the results into the second array' do
+    hashes = Sapor::Poll.as_hashes(SAMPLE_POLL_ARRAY)
+    expect(hashes[1]).to eq('a' => '1', 'b' => '2')
+  end
+end
+
+describe Sapor::Poll, '#new' do
+  it 'sets the area' do
+    expect(sample_poll.area.area_code).to eq('UT')
+  end
+
+  it 'sets the results' do
+    expect(sample_poll.result('Red')).to eq(1)
+    expect(sample_poll.result('Green')).to eq(2)
+    expect(sample_poll.result('Blue')).to eq(3)
+  end
+end
+
+describe Sapor::Poll, '#most_probable_value' do
+  it 'returns nil if no analysis has been run' do
+    expect(sample_poll.most_probable_value('Blue')).to be_nil
+  end
+
+  it 'gets the most probable value after analysis (50%)' do
+    poll = sample_poll
+    poll.analyze(MAX_ERROR)
+    expect(poll.most_probable_value('Blue')).to \
+      be_within(MAX_VALUE_ERROR).of(428_571)
+  end
+
+  it 'gets the most probable value after analysis (33%)' do
+    poll = sample_poll
+    poll.analyze(MAX_ERROR)
+    expect(poll.most_probable_value('Green')).to \
+      be_within(MAX_VALUE_ERROR).of(285_714)
+  end
+end
+
+describe Sapor::Poll, '#most_probable_fraction' do
+  it 'returns nil if no analysis has been run' do
+    expect(sample_poll.most_probable_fraction('Blue')).to be_nil
+  end
+
+  it 'gets the most probable fraction after analysis (50%)' do
+    poll = sample_poll
+    poll.analyze(MAX_ERROR)
+    expect(poll.most_probable_fraction('Blue')).to \
+      be_within(MAX_ERROR).of(3.to_f / 7.to_f)
+  end
+
+  it 'gets the most probable fraction after analysis (33%)' do
+    poll = sample_poll
+    poll.analyze(MAX_ERROR)
+    expect(poll.most_probable_fraction('Green')).to \
+      be_within(MAX_ERROR).of(2.to_f / 7.to_f)
+  end
+end
+
+describe Sapor::Poll, '#confidence_interval' do
+  it 'returns nil if no analysis has been run' do
+    expect(sample_poll.confidence_interval('Blue')).to be_nil
+  end
+end

data/spec/unit/polychotomy_spec.rb

@@ -0,0 +1,332 @@
+# encoding: utf-8
+#
+# Statistical Analysis of Polling Results (SAPoR)
+# Copyright (C) 2014 Filip van Laenen <f.a.vanlaenen@ieee.org>
+#
+# This file is part of SAPoR.
+#
+# SAPoR is free software: you can redistribute it and/or modify it under the
+# terms of the GNU General Public License as published by the Free Software
+# Foundation, either version 3 of the License, or (at your option) any later
+# version.
+#
+# SAPoR is distributed in the hope that it will be useful, but WITHOUT ANY
+# WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
+# A PARTICULAR PURPOSE.  See the GNU General Public License for more details.
+#
+# You can find a copy of the GNU General Public License in /doc/gpl.txt
+#
+
+require 'spec_helper'
+
+class TestArea
+  include Singleton
+
+  LAST_ELECTION_RESULT = { 'Red' => 91_811, 'Green' => 190_934,
+                           'Blue' => 290_647, 'Yellow' => 356_473 }
+
+  LAST_DETAILED_ELECTION_RESULT = { 'North' => { 'Red' => 50, 'Green' => 70
+                                               },
+                                    'South' => { 'Red' => 70, 'Green' => 50,
+                                                 'Blue' => 100 },
+                                    'East' => { 'Red' => 90, 'Green' => 70,
+                                                'Blue' => 90 },
+                                    'West' => { 'Red' => 110, 'Green' => 50,
+                                                'Yellow' => 120 } }
+
+  ELECTORAL_SYSTEM = Sapor::FirstPastThePost.new(LAST_ELECTION_RESULT,
+                                                 LAST_DETAILED_ELECTION_RESULT)
+
+  def area_code
+    'TA'
+  end
+
+  def coalitions
+    [['Red', 'Green'], ['Red', 'Blue']]
+  end
+
+  def population_size
+    1_000
+  end
+
+  def no_of_seats
+    LAST_DETAILED_ELECTION_RESULT.size
+  end
+
+  def seats(simulation)
+    ELECTORAL_SYSTEM.project(simulation)
+  end
+end
+
+TEST_AREA = TestArea.instance
+
+def pentachotomy(max_error = 0.01)
+  results = { 'Red' => 1, 'Green' => 2, 'Blue' => 3, 'Yellow' => 6,
+              'Other' => 1 }
+  dichotomies = Sapor::Dichotomies.new(results, 1000)
+  dichotomies.refine
+  dichotomies.refine
+  dichotomies.refine
+  Sapor::Polychotomy.new(results, TEST_AREA, dichotomies, max_error)
+end
+
+describe Sapor::Polychotomy, '#new' do
+  it 'extracts the 99.99% confidence intervals as ranges for max 1% error' do
+    expect(pentachotomy.range('Red').size).to eq(17)
+    expect(pentachotomy.range('Green').size).to eq(19)
+    expect(pentachotomy.range('Blue').size).to eq(20)
+    expect(pentachotomy.range('Yellow').size).to eq(22)
+  end
+
+  it 'extracts the 99.9999% confidence intervals as ranges for max 0.1%' \
+     ' error' do
+    expect(pentachotomy(0.001).range('Red').size).to eq(20)
+    expect(pentachotomy(0.001).range('Green').size).to eq(22)
+    expect(pentachotomy(0.001).range('Blue').size).to eq(23)
+    expect(pentachotomy(0.001).range('Yellow').size).to eq(25)
+  end
+
+  # TODO: A test case where the construction of the enumerator fails
+end
+
+describe Sapor::Polychotomy, '#error_estimate' do
+  it 'is 1 by default (no refinement)' do
+    polychotomy = pentachotomy
+    expect(polychotomy.error_estimate).to eq(1.0)
+  end
+
+  it 'is 1 after the first refinement' do
+    polychotomy = pentachotomy
+    polychotomy.refine
+    expect(polychotomy.error_estimate).to eq(1.0)
+  end
+
+  it 'is at least the difference between the most probable fractions (two' \
+     ' refinements)' do
+    polychotomy = pentachotomy
+    polychotomy.refine
+    polychotomy.refine
+    expect(polychotomy.error_estimate.round(3)).to eq(0.333)
+  end
+
+  it 'is at least the difference between the most probable fractions (three' \
+     ' refinements)' do
+    polychotomy = pentachotomy
+    polychotomy.refine
+    polychotomy.refine
+    polychotomy.refine
+    expect(polychotomy.error_estimate.round(3)).to eq(0.259)
+  end
+end
+
+describe Sapor::Polychotomy, '#most_probable_fraction' do
+  it 'is nil by default (no refinement)' do
+    polychotomy = pentachotomy
+    expect(polychotomy.most_probable_fraction('Red')).to be_nil
+    expect(polychotomy.most_probable_fraction('Green')).to be_nil
+    expect(polychotomy.most_probable_fraction('Blue')).to be_nil
+    expect(polychotomy.most_probable_fraction('Yellow')).to be_nil
+  end
+
+  it 'returns a value after the first refinement' do
+    polychotomy = pentachotomy
+    polychotomy.refine
+    expect(polychotomy.most_probable_fraction('Red').round(2)).to eq(0.02)
+    expect(polychotomy.most_probable_fraction('Green').round(2)).to eq(0.02)
+    expect(polychotomy.most_probable_fraction('Blue').round(2)).to eq(0.02)
+    expect(polychotomy.most_probable_fraction('Yellow').round(2)).to eq(0.09)
+  end
+
+  it 'returns an updated value after the second refinement' do
+    polychotomy = pentachotomy
+    polychotomy.refine
+    polychotomy.refine
+    expect(polychotomy.most_probable_fraction('Red').round(2)).to eq(0.28)
+    expect(polychotomy.most_probable_fraction('Green').round(2)).to eq(0.09)
+    expect(polychotomy.most_probable_fraction('Blue').round(2)).to eq(0.35)
+    expect(polychotomy.most_probable_fraction('Yellow').round(2)).to eq(0.20)
+  end
+
+  it 'returns an updated value after the third refinement' do
+    polychotomy = pentachotomy
+    polychotomy.refine
+    polychotomy.refine
+    polychotomy.refine
+    expect(polychotomy.most_probable_fraction('Red').round(2)).to eq(0.35)
+    expect(polychotomy.most_probable_fraction('Green').round(2)).to eq(0.09)
+    expect(polychotomy.most_probable_fraction('Blue').round(2)).to eq(0.09)
+    expect(polychotomy.most_probable_fraction('Yellow').round(2)).to eq(0.39)
+  end
+end
+
+describe Sapor::Polychotomy, '#no_of_data_points' do
+  it 'is 0 by default (no refinement)' do
+    polychotomy = pentachotomy
+    expect(polychotomy.no_of_data_points).to eq(0)
+  end
+
+  it 'is 1 after 1 refinement' do
+    polychotomy = pentachotomy
+    polychotomy.refine
+    expect(polychotomy.no_of_data_points).to eq(1)
+  end
+
+  it 'is 19 after 2 refinements' do
+    polychotomy = pentachotomy
+    polychotomy.refine
+    polychotomy.refine
+    expect(polychotomy.no_of_data_points).to eq(8)
+  end
+
+  it 'is 38 after 3 refinements' do
+    polychotomy = pentachotomy
+    polychotomy.refine
+    polychotomy.refine
+    polychotomy.refine
+    expect(polychotomy.no_of_data_points).to eq(27)
+  end
+
+  # TODO: A test case that runs all the simulations
+end
+
+describe Sapor::Polychotomy, '#no_of_simulations' do
+  it 'is 0 by default (no refinement)' do
+    polychotomy = pentachotomy
+    expect(polychotomy.no_of_simulations).to eq(0)
+  end
+
+  it 'is 1 after 1 refinement' do
+    polychotomy = pentachotomy
+    polychotomy.refine
+    expect(polychotomy.no_of_simulations).to eq(1)
+  end
+
+  it 'is 2 after 2 refinements' do
+    polychotomy = pentachotomy
+    polychotomy.refine
+    polychotomy.refine
+    expect(polychotomy.no_of_simulations).to eq(2)
+  end
+
+  it 'is 4 after 3 refinements' do
+    polychotomy = pentachotomy
+    polychotomy.refine
+    polychotomy.refine
+    polychotomy.refine
+    expect(polychotomy.no_of_simulations).to eq(4)
+  end
+
+  # TODO: A test case that runs until the end of the enumerator
+end
+
+describe Sapor::Polychotomy, '#progress_report' do
+  it 'reports progress with the number of simulations and data points,' \
+     ' together with the search space size and the fraction that already has' \
+     ' been searched' do
+    expected_report = '1 simulations out of 1 data points, 1 / 142,120 of' \
+                      ' search space size (142,120).'
+    polychotomy = pentachotomy
+    polychotomy.refine
+    expect(polychotomy.progress_report).to eq(expected_report)
+  end
+
+  it 'reports no of simulations and data points with thousands separator' do
+    expected_report = '1,024 simulations out of 9,287 data points, 1 / 15' \
+                      ' of search space size (142,120).'
+    polychotomy = pentachotomy
+    11.times { polychotomy.refine }
+    expect(polychotomy.progress_report).to eq(expected_report)
+  end
+
+  it 'rounds the fraction of the space size searched with 1 decimal if' \
+     ' larger than one tenth' do
+    expected_report = '2,048 simulations out of 19,433 data points, 1 / 7.0' \
+                      ' of search space size (142,120).'
+    polychotomy = pentachotomy
+    12.times { polychotomy.refine }
+    expect(polychotomy.progress_report).to eq(expected_report)
+  end
+end
+
+describe Sapor::Polychotomy, '#report' do
+  it 'produces a report after first refinement for short choice labels' do
+    expected_report = 'Most probable rounded fractions, fractions and 95%' \
+                      " confidence intervals:\n" \
+                      "Choice  Result    MPRF    MPF      CI(95%)      P(>↓)  Seats\n" \
+                      "Yellow   46.2%    5.6%    9.3%    9.3%– 11.1%  100.0%   1–1\n" + # TODO: Shouldn't MPRF be 5.6% and CI 9.3%–9.3%?
+                      "Blue     23.1%    1.9%    1.9%    1.9%–  3.7%    0.0%   0–0\n" + # TODO: Shouldn't CI be 1.9%–1.9%?
+                      "Green    15.4%    1.9%    1.9%    1.9%–  3.7%    0.0%   0–0\n" + # TODO: Shouldn't CI be 1.9%–1.9%?
+                      "Red       7.7%    1.9%    1.9%    1.9%–  3.7%           3–3\n" + # TODO: Shouldn't CI be 1.9%–1.9%?
+                      "Coalition    Result    MPRF    MPF      CI(95%)     P(>50%)  Seats  P(>50%)\n" +
+                      "Blue + Red    30.8%   50.0%    3.8%    3.8%–  3.8%     0.0%   3–3    100.0%\n" +
+                      'Green + Red   23.1%   50.0%    3.8%    3.8%–  3.8%     0.0%   3–3    100.0%'
+    polychotomy = pentachotomy
+    polychotomy.refine
+    expect(polychotomy.report).to eq(expected_report)
+  end
+
+  it 'produces a report by default for long choice labels' do
+    results = { 'Dark Red' => 1, 'Light Green' => 2, 'Medium Blue' => 3,
+                'Other' => 1 }
+    dichotomies = Sapor::Dichotomies.new(results, 1000)
+    dichotomies.refine
+    dichotomies.refine
+    dichotomies.refine
+    polychotomy = Sapor::Polychotomy.new(results, TEST_AREA, dichotomies, 0.01)
+    polychotomy.refine
+    expected_report = 'Most probable rounded fractions, fractions and 95%' \
+                      " confidence intervals:\n" \
+                      "Choice       Result    MPRF    MPF      CI(95%)      P(>↓)  Seats\n" \
+                      "Medium Blue   42.9%    1.9%    1.9%    1.9%–  3.7%    0.0%   0–0\n" + # TODO: Shouldn't CI be 1.9%–1.9%?
+                      "Light Green   28.6%    1.9%    1.9%    1.9%–  3.7%    0.0%   0–0\n" + # TODO: Shouldn't CI be 1.9%–1.9%?
+                      "Dark Red      14.3%    1.9%    1.9%    1.9%–  3.7%           0–0\n" + # TODO: Shouldn't CI be 1.9%–1.9%?
+# TODO:                      "Red                                                  3–3\n" \
+# TODO:                      'Yellow                                               1–1'
+                      "Coalition    Result    MPRF    MPF      CI(95%)     P(>50%)  Seats  P(>50%)\n" +
+                      "Blue + Red     0.0%   50.0%    0.0%    0.0%–  0.0%     0.0%   0–0      0.0%\n" +
+                      'Green + Red    0.0%   50.0%    0.0%    0.0%–  0.0%     0.0%   0–0      0.0%'
+    expect(polychotomy.report).to eq(expected_report)
+  end
+
+  it 'produces a report by default for short choice labels' do
+    results = { 'Red' => 1, 'Green' => 2, 'Blue' => 3, 'Other' => 1 }
+    dichotomies = Sapor::Dichotomies.new(results, 1000)
+    dichotomies.refine
+    dichotomies.refine
+    dichotomies.refine
+    polychotomy = Sapor::Polychotomy.new(results, TEST_AREA, dichotomies, 0.01)
+    polychotomy.refine
+    expected_report = 'Most probable rounded fractions, fractions and 95%' \
+                      " confidence intervals:\n" \
+                      "Choice  Result    MPRF    MPF      CI(95%)      P(>↓)  Seats\n" \
+                      "Blue     42.9%    1.9%    1.9%    1.9%–  3.7%    0.0%   0–0\n" + # TODO: Shouldn't CI be 1.9%–1.9%?
+                      "Green    28.6%    1.9%    1.9%    1.9%–  3.7%    0.0%   0–0\n" + # TODO: Shouldn't CI be 1.9%–1.9%?
+                      "Red      14.3%    1.9%    1.9%    1.9%–  3.7%           3–3\n" + # TODO: Shouldn't CI be 1.9%–1.9%?
+# TODO:                      'Yellow                                          1–1'
+                      "Coalition    Result    MPRF    MPF      CI(95%)     P(>50%)  Seats  P(>50%)\n" +
+                      "Blue + Red    57.1%   50.0%    3.8%    3.8%–  3.8%     0.0%   3–3    100.0%\n" +
+                      'Green + Red   42.9%   50.0%    3.8%    3.8%–  3.8%     0.0%   3–3    100.0%'
+    expect(polychotomy.report).to eq(expected_report)
+  end
+
+  it 'sorts line with the same results alphabetically' do
+    results = { 'Red' => 1, 'Green' => 1, 'Blue' => 1, 'Other' => 1 }
+    dichotomies = Sapor::Dichotomies.new(results, 1000)
+    dichotomies.refine
+    dichotomies.refine
+    dichotomies.refine
+    polychotomy = Sapor::Polychotomy.new(results, TEST_AREA, dichotomies, 0.01)
+    polychotomy.refine
+    expected_report = 'Most probable rounded fractions, fractions and 95%' \
+                      " confidence intervals:\n" \
+                      "Choice  Result    MPRF    MPF      CI(95%)      P(>↓)  Seats\n" \
+                      "Blue     25.0%    1.9%    1.9%    1.9%–  3.7%    0.0%   0–0\n" + # TODO: Shouldn't CI be 1.9%–1.9%?
+                      "Green    25.0%    1.9%    1.9%    1.9%–  3.7%    0.0%   0–0\n" + # TODO: Shouldn't CI be 1.9%–1.9%?
+                      "Red      25.0%    1.9%    1.9%    1.9%–  3.7%           3–3\n" + # TODO: Shouldn't CI be 1.9%–1.9%?
+# TODO:                      'Yellow                                          1–1'
+                      "Coalition    Result    MPRF    MPF      CI(95%)     P(>50%)  Seats  P(>50%)\n" +
+                      "Blue + Red    50.0%   50.0%    3.8%    3.8%–  3.8%     0.0%   3–3    100.0%\n" +
+                      'Green + Red   50.0%   50.0%    3.8%    3.8%–  3.8%     0.0%   3–3    100.0%'
+    expect(polychotomy.report).to eq(expected_report)
+  end
+end