sapor 0.3.6

data/spec/unit/portugal_with_a_and_il_without_paf_spec.rb

@@ -0,0 +1,68 @@
+#
+# Statistical Analysis of Polling Results (SAPoR)
+# Copyright (C) 2020 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::PortugalWithAAndIlWithoutPaf, '#area_code' do
+  it 'returns PT∪{A,IL}\{PàF} as the area code' do
+    expect(Sapor::PortugalWithAAndIlWithoutPaf.instance.area_code).to \
+      eq('PT∪{A,IL}\{PàF}')
+  end
+end
+
+describe Sapor::PortugalWithAAndIlWithoutPaf, '#no_of_seats' do
+  it 'returns 230 as the number of seats' do
+    expect(Sapor::PortugalWithAAndIlWithoutPaf.instance.no_of_seats).to eq(230)
+  end
+end
+
+describe Sapor::PortugalWithAAndIlWithoutPaf, '#population_size' do
+  it 'returns a population size of 5,206,113' do
+    expect(Sapor::PortugalWithAAndIlWithoutPaf.instance.population_size).to eq(5_206_113)
+  end
+end
+
+describe Sapor::PortugalWithAAndIlWithoutPaf, '#seats' do
+  it 'calculates the number of seats for the election of 2015 adjusted for' \
+     ' the seat distribution for 2019' do
+    PortugalWithAAndIlWithoutPaf = Sapor::PortugalWithAAndIlWithoutPaf.instance
+    results = PortugalWithAAndIlWithoutPaf.overall_election_results_of_2015
+    seats = Sapor::PortugalWithAAndIlWithoutPaf.instance.seats(results)
+    expect(seats['Partido Social Democrata']).to eq(5 + 37)
+    expect(seats['Aliança']).to eq(0 + 37)
+    expect(seats['Iniciativa Liberal']).to eq(0 + 80)    
+    expect(seats['CDS–Partido Popular']).to eq(0 + 30)
+    expect(seats['Partido Socialista']).to eq(86 - 54)
+    expect(seats['Bloco de Esquerda']).to eq(19 - 15)
+    expect(seats['Coligação Democrática Unitária']).to eq(17 - 12)
+    expect(seats['Pessoas–Animais–Natureza']).to eq(1 - 1)
+    expect(seats['Partido Democrático Republicano']).to eq(0)
+    expect(seats['Partido Comunista dos Trabalhadores Portugueses']).to eq(0)
+    expect(seats['LIVRE']).to eq(0)
+    expect(seats['Partido Nacional Renovador']).to eq(0)
+    expect(seats['Partido da Terra']).to eq(0)
+    expect(seats['AGIR']).to eq(0)
+    expect(seats['Partido Trabalhista Português']).to eq(0)
+    expect(seats['Nós, Cidadãos!']).to eq(0)
+    expect(seats['Partido Popular Monárquico']).to eq(0)
+    expect(seats['Juntos pelo Povo']).to eq(0)
+    expect(seats['Partido Unido dos Reformados e Pensionistas']).to eq(0)
+    expect(seats['Aliança Açores']).to eq(0)
+    expect(seats['Partido Cidadania e Democracia Cristã']).to eq(0)
+  end
+end

data/spec/unit/portugal_with_a_ch_and_il_without_paf_spec.rb

@@ -0,0 +1,69 @@
+#
+# Statistical Analysis of Polling Results (SAPoR)
+# Copyright (C) 2020 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::PortugalWithAChAndIlWithoutPaf, '#area_code' do
+  it 'returns PT∪{A,CH,IL}\{PàF} as the area code' do
+    expect(Sapor::PortugalWithAChAndIlWithoutPaf.instance.area_code).to \
+      eq('PT∪{A,CH,IL}\{PàF}')
+  end
+end
+
+describe Sapor::PortugalWithAChAndIlWithoutPaf, '#no_of_seats' do
+  it 'returns 230 as the number of seats' do
+    expect(Sapor::PortugalWithAChAndIlWithoutPaf.instance.no_of_seats).to eq(230)
+  end
+end
+
+describe Sapor::PortugalWithAChAndIlWithoutPaf, '#population_size' do
+  it 'returns a population size of 5,206,113' do
+    expect(Sapor::PortugalWithAChAndIlWithoutPaf.instance.population_size).to eq(5_206_113)
+  end
+end
+
+describe Sapor::PortugalWithAChAndIlWithoutPaf, '#seats' do
+  it 'calculates the number of seats for the election of 2015 adjusted for' \
+     ' the seat distribution for 2019' do
+    PortugalWithAChAndIlWithoutPaf = Sapor::PortugalWithAChAndIlWithoutPaf.instance
+    results = PortugalWithAChAndIlWithoutPaf.overall_election_results_of_2015
+    seats = Sapor::PortugalWithAChAndIlWithoutPaf.instance.seats(results)
+    expect(seats['Partido Social Democrata']).to eq(5 + 33)
+    expect(seats['Aliança']).to eq(0 + 33)
+    expect(seats['Chega']).to eq(0 + 28)
+    expect(seats['Iniciativa Liberal']).to eq(0 + 70)    
+    expect(seats['CDS–Partido Popular']).to eq(0 + 26)
+    expect(seats['Partido Socialista']).to eq(86 - 58)
+    expect(seats['Bloco de Esquerda']).to eq(19 - 15)
+    expect(seats['Coligação Democrática Unitária']).to eq(17 - 14)
+    expect(seats['Pessoas–Animais–Natureza']).to eq(1 - 1)
+    expect(seats['Partido Democrático Republicano']).to eq(0)
+    expect(seats['Partido Comunista dos Trabalhadores Portugueses']).to eq(0)
+    expect(seats['LIVRE']).to eq(0)
+    expect(seats['Partido Nacional Renovador']).to eq(0)
+    expect(seats['Partido da Terra']).to eq(0)
+    expect(seats['AGIR']).to eq(0)
+    expect(seats['Partido Trabalhista Português']).to eq(0)
+    expect(seats['Nós, Cidadãos!']).to eq(0)
+    expect(seats['Partido Popular Monárquico']).to eq(0)
+    expect(seats['Juntos pelo Povo']).to eq(0)
+    expect(seats['Partido Unido dos Reformados e Pensionistas']).to eq(0)
+    expect(seats['Aliança Açores']).to eq(0)
+    expect(seats['Partido Cidadania e Democracia Cristã']).to eq(0)
+  end
+end

data/spec/unit/portugal_with_a_without_paf_spec.rb

@@ -0,0 +1,67 @@
+#
+# Statistical Analysis of Polling Results (SAPoR)
+# Copyright (C) 2020 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::PortugalWithAWithoutPaf, '#area_code' do
+  it 'returns PT∪{A}\{PàF} as the area code' do
+    expect(Sapor::PortugalWithAWithoutPaf.instance.area_code).to \
+      eq('PT∪{A}\{PàF}')
+  end
+end
+
+describe Sapor::PortugalWithAWithoutPaf, '#no_of_seats' do
+  it 'returns 230 as the number of seats' do
+    expect(Sapor::PortugalWithAWithoutPaf.instance.no_of_seats).to eq(230)
+  end
+end
+
+describe Sapor::PortugalWithAWithoutPaf, '#population_size' do
+  it 'returns a population size of 5,206,113' do
+    expect(Sapor::PortugalWithAWithoutPaf.instance.population_size).to eq(5_206_113)
+  end
+end
+
+describe Sapor::PortugalWithAWithoutPaf, '#seats' do
+  it 'calculates the number of seats for the election of 2015 adjusted for' \
+     ' the seat distribution for 2019' do
+    PortugalWithAWithoutPaf = Sapor::PortugalWithAWithoutPaf.instance
+    results = PortugalWithAWithoutPaf.overall_election_results_of_2015
+    seats = Sapor::PortugalWithAWithoutPaf.instance.seats(results)
+    expect(seats['Partido Social Democrata']).to eq(5 + 60)
+    expect(seats['Aliança']).to eq(0 + 57)
+    expect(seats['CDS–Partido Popular']).to eq(0 + 47)
+    expect(seats['Partido Socialista']).to eq(86 - 41)
+    expect(seats['Bloco de Esquerda']).to eq(19 - 10)
+    expect(seats['Coligação Democrática Unitária']).to eq(17 - 10)
+    expect(seats['Pessoas–Animais–Natureza']).to eq(1 - 1)
+    expect(seats['Partido Democrático Republicano']).to eq(0)
+    expect(seats['Partido Comunista dos Trabalhadores Portugueses']).to eq(0)
+    expect(seats['LIVRE']).to eq(0)
+    expect(seats['Partido Nacional Renovador']).to eq(0)
+    expect(seats['Partido da Terra']).to eq(0)
+    expect(seats['AGIR']).to eq(0)
+    expect(seats['Partido Trabalhista Português']).to eq(0)
+    expect(seats['Nós, Cidadãos!']).to eq(0)
+    expect(seats['Partido Popular Monárquico']).to eq(0)
+    expect(seats['Juntos pelo Povo']).to eq(0)
+    expect(seats['Partido Unido dos Reformados e Pensionistas']).to eq(0)
+    expect(seats['Aliança Açores']).to eq(0)
+    expect(seats['Partido Cidadania e Democracia Cristã']).to eq(0)
+  end
+end

data/spec/unit/portugal_with_ch_and_il_without_paf_spec.rb

@@ -0,0 +1,68 @@
+#
+# Statistical Analysis of Polling Results (SAPoR)
+# Copyright (C) 2020 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::PortugalWithChAndIlWithoutPaf, '#area_code' do
+  it 'returns PT∪{CH,IL}\{PàF} as the area code' do
+    expect(Sapor::PortugalWithChAndIlWithoutPaf.instance.area_code).to \
+      eq('PT∪{CH,IL}\{PàF}')
+  end
+end
+
+describe Sapor::PortugalWithChAndIlWithoutPaf, '#no_of_seats' do
+  it 'returns 230 as the number of seats' do
+    expect(Sapor::PortugalWithChAndIlWithoutPaf.instance.no_of_seats).to eq(230)
+  end
+end
+
+describe Sapor::PortugalWithChAndIlWithoutPaf, '#population_size' do
+  it 'returns a population size of 5,206,113' do
+    expect(Sapor::PortugalWithChAndIlWithoutPaf.instance.population_size).to eq(5_206_113)
+  end
+end
+
+describe Sapor::PortugalWithChAndIlWithoutPaf, '#seats' do
+  it 'calculates the number of seats for the election of 2015 adjusted for' \
+     ' the seat distribution for 2019' do
+    PortugalWithChAndIlWithoutPaf = Sapor::PortugalWithChAndIlWithoutPaf.instance
+    results = PortugalWithChAndIlWithoutPaf.overall_election_results_of_2015
+    seats = Sapor::PortugalWithChAndIlWithoutPaf.instance.seats(results)
+    expect(seats['Partido Social Democrata']).to eq(5 + 37)
+    expect(seats['Chega']).to eq(0 + 37)
+    expect(seats['Iniciativa Liberal']).to eq(0 + 80)    
+    expect(seats['CDS–Partido Popular']).to eq(0 + 30)
+    expect(seats['Partido Socialista']).to eq(86 - 54)
+    expect(seats['Bloco de Esquerda']).to eq(19 - 15)
+    expect(seats['Coligação Democrática Unitária']).to eq(17 - 12)
+    expect(seats['Pessoas–Animais–Natureza']).to eq(1 - 1)
+    expect(seats['Partido Democrático Republicano']).to eq(0)
+    expect(seats['Partido Comunista dos Trabalhadores Portugueses']).to eq(0)
+    expect(seats['LIVRE']).to eq(0)
+    expect(seats['Partido Nacional Renovador']).to eq(0)
+    expect(seats['Partido da Terra']).to eq(0)
+    expect(seats['AGIR']).to eq(0)
+    expect(seats['Partido Trabalhista Português']).to eq(0)
+    expect(seats['Nós, Cidadãos!']).to eq(0)
+    expect(seats['Partido Popular Monárquico']).to eq(0)
+    expect(seats['Juntos pelo Povo']).to eq(0)
+    expect(seats['Partido Unido dos Reformados e Pensionistas']).to eq(0)
+    expect(seats['Aliança Açores']).to eq(0)
+    expect(seats['Partido Cidadania e Democracia Cristã']).to eq(0)
+  end
+end

data/spec/unit/portugal_without_paf_spec.rb

@@ -0,0 +1,66 @@
+#
+# Statistical Analysis of Polling Results (SAPoR)
+# Copyright (C) 2020 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::PortugalWithoutPaf, '#area_code' do
+  it 'returns PT\{PàF} as the area code' do
+    expect(Sapor::PortugalWithoutPaf.instance.area_code).to \
+      eq('PT\{PàF}')
+  end
+end
+
+describe Sapor::PortugalWithoutPaf, '#no_of_seats' do
+  it 'returns 230 as the number of seats' do
+    expect(Sapor::PortugalWithoutPaf.instance.no_of_seats).to eq(230)
+  end
+end
+
+describe Sapor::PortugalWithoutPaf, '#population_size' do
+  it 'returns a population size of 5,206,113' do
+    expect(Sapor::PortugalWithoutPaf.instance.population_size).to eq(5_206_113)
+  end
+end
+
+describe Sapor::PortugalWithoutPaf, '#seats' do
+  it 'calculates the number of seats for the election of 2015 adjusted for' \
+     ' the seat distribution for 2019' do
+    PortugalWithoutPaf = Sapor::PortugalWithoutPaf.instance
+    results = PortugalWithoutPaf.overall_election_results_of_2015
+    seats = Sapor::PortugalWithoutPaf.instance.seats(results)
+    expect(seats['Partido Social Democrata']).to eq(5 + 73)
+    expect(seats['CDS–Partido Popular']).to eq(0 + 67)
+    expect(seats['Partido Socialista']).to eq(86 - 23)
+    expect(seats['Bloco de Esquerda']).to eq(19 - 6)
+    expect(seats['Coligação Democrática Unitária']).to eq(17 - 8)
+    expect(seats['Pessoas–Animais–Natureza']).to eq(1 - 1)
+    expect(seats['Partido Democrático Republicano']).to eq(0)
+    expect(seats['Partido Comunista dos Trabalhadores Portugueses']).to eq(0)
+    expect(seats['LIVRE']).to eq(0)
+    expect(seats['Partido Nacional Renovador']).to eq(0)
+    expect(seats['Partido da Terra']).to eq(0)
+    expect(seats['AGIR']).to eq(0)
+    expect(seats['Partido Trabalhista Português']).to eq(0)
+    expect(seats['Nós, Cidadãos!']).to eq(0)
+    expect(seats['Partido Popular Monárquico']).to eq(0)
+    expect(seats['Juntos pelo Povo']).to eq(0)
+    expect(seats['Partido Unido dos Reformados e Pensionistas']).to eq(0)
+    expect(seats['Aliança Açores']).to eq(0)
+    expect(seats['Partido Cidadania e Democracia Cristã']).to eq(0)
+  end
+end

data/spec/unit/pseudorandom_multirange_enumerator_spec.rb

@@ -0,0 +1,82 @@
+# encoding: utf-8
+#
+# Statistical Analysis of Polling Results (SAPoR)
+# Copyright (C) 2020 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::PseudoRandomMultiRangeEnumerator, '#each' do
+  it 'sets the size correctly for one dimension' do
+    enum = Sapor::PseudoRandomMultiRangeEnumerator.new([5]).each
+    expect(enum.size).to eq(5)
+  end
+
+  it 'sets the size correctly for two dimensions' do
+    enum = Sapor::PseudoRandomMultiRangeEnumerator.new([5, 3]).each
+    expect(enum.size).to eq(15)
+  end
+
+  it 'returns a simple enumeration for one dimension' do
+    enum = Sapor::PseudoRandomMultiRangeEnumerator.new([5]).each
+    expect(enum.to_a).to eq([[0], [1], [2], [3], [4]])
+  end
+
+  it 'returns [0] as the first element for one dimension' do
+    enum = Sapor::PseudoRandomMultiRangeEnumerator.new([5]).each
+    expect(enum.next).to eq([0])
+  end
+
+  it 'returns [0, 0] as the first element for two dimensions' do
+    enum = Sapor::PseudoRandomMultiRangeEnumerator.new([5, 3]).each
+    expect(enum.next).to eq([0, 0])
+  end
+
+  it 'returns [1] as the second element for one dimension' do
+    enum = Sapor::PseudoRandomMultiRangeEnumerator.new([5]).each
+    enum.next
+    expect(enum.next).to eq([1])
+  end
+
+  it 'returns [1, 2] as the second element of a 3×5 hyper-rectangle' do
+    enum = Sapor::PseudoRandomMultiRangeEnumerator.new([3, 5]).each
+    enum.next
+    expect(enum.next).to eq([1, 2])
+  end
+
+  it 'returns [2, 1] as the second element of a 5×3 hyper-rectangle' do
+    enum = Sapor::PseudoRandomMultiRangeEnumerator.new([5, 3]).each
+    enum.next
+    expect(enum.next).to eq([2, 1])
+  end
+
+  it 'returns [1, 3] as the second element of a 2×5 hyper-rectangle' do
+    enum = Sapor::PseudoRandomMultiRangeEnumerator.new([2, 5]).each
+    enum.next
+    expect(enum.next).to eq([1, 3])
+  end
+
+  it 'returns [1, 2, 5] as the second element of a 9×9×9 hyper-rectangle' do
+    enum = Sapor::PseudoRandomMultiRangeEnumerator.new([9, 9, 9]).each
+    enum.next
+    expect(enum.next).to eq([1, 2, 5])
+  end
+
+  it "raises an ArgumentError when it can't construct incrementers" do
+    expect { Sapor::PseudoRandomMultiRangeEnumerator.new([3, 2]) }.to \
+      raise_error(ArgumentError, 'Could not construct suitable incrementers.')
+  end
+end

data/spec/unit/referendum_polychotomy_spec.rb

@@ -0,0 +1,289 @@
+# encoding: utf-8
+#
+# Statistical Analysis of Polling Results (SAPoR)
+# Copyright (C) 2020 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 ReferendumTestArea
+  include Singleton
+
+  def area_code
+    'RA'
+  end
+
+  def population_size
+    1_000
+  end
+end
+
+REFERENDUM_TEST_AREA = ReferendumTestArea.instance
+
+def referendum_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::ReferendumPolychotomy.new(results, REFERENDUM_TEST_AREA, dichotomies, max_error)
+end
+
+describe Sapor::ReferendumPolychotomy, '#new' do
+  it 'extracts the 99.99% confidence intervals as ranges for max 1% error' do
+    expect(referendum_pentachotomy.range('Red').size).to eq(17)
+    expect(referendum_pentachotomy.range('Green').size).to eq(19)
+    expect(referendum_pentachotomy.range('Blue').size).to eq(20)
+    expect(referendum_pentachotomy.range('Yellow').size).to eq(22)
+  end
+
+  it 'extracts the 99.9999% confidence intervals as ranges for max 0.1%' \
+     ' error' do
+    expect(referendum_pentachotomy(0.001).range('Red').size).to eq(20)
+    expect(referendum_pentachotomy(0.001).range('Green').size).to eq(22)
+    expect(referendum_pentachotomy(0.001).range('Blue').size).to eq(23)
+    expect(referendum_pentachotomy(0.001).range('Yellow').size).to eq(25)
+  end
+
+  # TODO: A test case where the construction of the enumerator fails
+end
+
+describe Sapor::ReferendumPolychotomy, '#error_estimate' do
+  it 'is 1 by default (no refinement)' do
+    polychotomy = referendum_pentachotomy
+    expect(polychotomy.error_estimate).to eq(1.0)
+  end
+
+  it 'is 1 after the first refinement' do
+    polychotomy = referendum_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 = referendum_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 = referendum_pentachotomy
+    polychotomy.refine
+    polychotomy.refine
+    polychotomy.refine
+    expect(polychotomy.error_estimate.round(3)).to eq(0.259)
+  end
+end
+
+describe Sapor::ReferendumPolychotomy, '#most_probable_fraction' do
+  it 'is nil by default (no refinement)' do
+    polychotomy = referendum_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 = referendum_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 = referendum_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 = referendum_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::ReferendumPolychotomy, '#no_of_data_points' do
+  it 'is 0 by default (no refinement)' do
+    polychotomy = referendum_pentachotomy
+    expect(polychotomy.no_of_data_points).to eq(0)
+  end
+
+  it 'is 1 after 1 refinement' do
+    polychotomy = referendum_pentachotomy
+    polychotomy.refine
+    expect(polychotomy.no_of_data_points).to eq(1)
+  end
+
+  it 'is 19 after 2 refinements' do
+    polychotomy = referendum_pentachotomy
+    polychotomy.refine
+    polychotomy.refine
+    expect(polychotomy.no_of_data_points).to eq(8)
+  end
+
+  it 'is 38 after 3 refinements' do
+    polychotomy = referendum_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::ReferendumPolychotomy, '#no_of_simulations' do
+  it 'is 0 by default (no refinement)' do
+    polychotomy = referendum_pentachotomy
+    expect(polychotomy.no_of_simulations).to eq(0)
+  end
+
+  it 'is 1 after 1 refinement' do
+    polychotomy = referendum_pentachotomy
+    polychotomy.refine
+    expect(polychotomy.no_of_simulations).to eq(1)
+  end
+
+  it 'is 2 after 2 refinements' do
+    polychotomy = referendum_pentachotomy
+    polychotomy.refine
+    polychotomy.refine
+    expect(polychotomy.no_of_simulations).to eq(2)
+  end
+
+  it 'is 4 after 3 refinements' do
+    polychotomy = referendum_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::ReferendumPolychotomy, '#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 = referendum_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 = referendum_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.3' \
+                      ' of search space size (142,120).'
+    polychotomy = referendum_pentachotomy
+    12.times { polychotomy.refine }
+    expect(polychotomy.progress_report).to eq(expected_report)
+  end
+end
+
+describe Sapor::ReferendumPolychotomy, '#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(>↓)   P(R↑)    P(R)   P(R↓)\n" \
+                      "Yellow   46.2%    5.6%    9.3%    9.3%– 11.1%  100.0%            0.0%  100.0%\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%    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%    0.0%  100.0%\n" + # TODO: Shouldn't CI be 1.9%–1.9%?
+                      "Red       7.7%    1.9%    1.9%    1.9%–  3.7%          100.0%    0.0%" # TODO: Shouldn't CI be 1.9%–1.9%?
+    polychotomy = referendum_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::ReferendumPolychotomy.new(results, REFERENDUM_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(>↓)   P(R↑)    P(R)   P(R↓)\n" \
+                      "Medium Blue   42.9%    1.9%    1.9%    1.9%–  3.7%    0.0%            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%    0.0%  100.0%\n" + # TODO: Shouldn't CI be 1.9%–1.9%?
+                      "Dark Red      14.3%    1.9%    1.9%    1.9%–  3.7%          100.0%    0.0%" # TODO: Shouldn't CI be 1.9%–1.9%?
+    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::ReferendumPolychotomy.new(results, REFERENDUM_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(>↓)   P(R↑)    P(R)   P(R↓)\n" \
+                      "Blue     42.9%    1.9%    1.9%    1.9%–  3.7%    0.0%            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%    0.0%  100.0%\n" + # TODO: Shouldn't CI be 1.9%–1.9%?
+                      "Red      14.3%    1.9%    1.9%    1.9%–  3.7%          100.0%    0.0%" # TODO: Shouldn't CI be 1.9%–1.9%?
+    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::ReferendumPolychotomy.new(results, REFERENDUM_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(>↓)   P(R↑)    P(R)   P(R↓)\n" \
+                      "Blue     25.0%    1.9%    1.9%    1.9%–  3.7%    0.0%            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%    0.0%  100.0%\n" + # TODO: Shouldn't CI be 1.9%–1.9%?
+                      "Red      25.0%    1.9%    1.9%    1.9%–  3.7%          100.0%    0.0%" # TODO: Shouldn't CI be 1.9%–1.9%?
+    expect(polychotomy.report).to eq(expected_report)
+  end
+end