pavement_condition_index 1.0.5

data/lib/tasks/regression.rake

@@ -0,0 +1,121 @@
+require 'pavement_condition_index/lookups/observed_deduct_values'
+require 'pavement_condition_index/lookups/observed_corrected_deduct_values'
+require 'pp'
+require 'matrix'
+
+namespace :regression do
+
+  task :generate_deduct_coefficients do
+    output = PavementConditionIndex::Lookups::ObservedDeductValues::OBSERVED_VALUES.map do |pavement_type_key,pavement_type_values|
+      remapped_pavement_type_values = pavement_type_values.map do |distress_type_key,distress_type_values|
+        coefficients = {low: nil, medium: nil, high: nil}.map do |severity_key,coefficient|
+          polynomial_degree = 0
+          @candidate_coefficients = []
+          @x_values = distress_type_values[:chart_type] == :log ? distress_type_values[:x_values].map {|x| Math.log10(x)} : distress_type_values[:x_values]
+          loop do 
+            polynomial_degree += 1
+            @candidate_coefficients = regress_polynomial(@x_values,distress_type_values[:expected_y_values][severity_key],polynomial_degree)
+            candidate_y_values = @x_values.map {|x| evaluate_polynomial(*@candidate_coefficients).call(x) }
+            sum_of_squares = 0.0
+            distress_type_values[:x_values].each_with_index do |x, index|
+              sum_of_squares += (candidate_y_values[index] - distress_type_values[:expected_y_values][severity_key][index])**2
+            end
+            break if sum_of_squares < 10.0
+          end
+          [severity_key, @candidate_coefficients]
+        end.to_h
+        remapped_distress_type_values = {
+          valid_min: distress_type_values[:valid_min],
+          valid_max: distress_type_values[:valid_max],
+          chart_type: distress_type_values[:chart_type],
+          coefficients: coefficients
+        }
+        [distress_type_key,remapped_distress_type_values]
+      end.to_h
+      [pavement_type_key,remapped_pavement_type_values]
+    end.to_h
+    pp output
+  end
+
+  task :generate_corrected_deduct_coefficients do
+    output = PavementConditionIndex::Lookups::ObservedCorrectedDeductValues::OBSERVED_VALUES.map do |pavement_type_key,pavement_type_values|
+      coefficients = pavement_type_values[:expected_y_values].map do |q_value_key,observed_values|
+        polynomial_degree = 0
+        @candidate_coefficients = []
+        @x_values = pavement_type_values[:x_values]
+        loop do 
+          polynomial_degree += 1
+          @candidate_coefficients = regress_polynomial(@x_values,observed_values,polynomial_degree)
+          candidate_y_values = @x_values.map {|x| evaluate_polynomial(*@candidate_coefficients).call(x) }
+          sum_of_squares = 0.0
+          pavement_type_values[:x_values].each_with_index do |x, index|
+            sum_of_squares += (candidate_y_values[index] - observed_values[index])**2
+          end
+          break if sum_of_squares < 10.0
+        end
+        [q_value_key, @candidate_coefficients]
+      end.to_h
+      remapped_pavement_type_values = {
+        coefficients: coefficients
+      }
+      [pavement_type_key,remapped_pavement_type_values]
+    end.to_h
+    pp output
+  end
+
+  task :generate_calculated_deduct_coefficients_file do
+    header =
+      "# generate this whole file by running the following command:\n"+
+      "# rake regression:generate_calculated_deduct_coefficients_file > lib/pavement_condition_index/lookups/calculated_deduct_coefficients.rb\n"+
+      "\n"+ 
+      "module PavementConditionIndex\n"+
+      "  module Lookups\n"+
+      "    class CalculatedDeductCoefficients\n"+
+      "\n"+
+      "# Output of `rake regression:generate_deduct_coefficients`\n"+
+      "COEFFICIENTS = "
+    puts header
+
+    Rake::Task["regression:generate_deduct_coefficients"].invoke
+
+    footer = 
+      "    end\n"+
+      "  end\n"+
+      "end\n"
+    puts footer
+  end
+
+  task :generate_calculated_corrected_deduct_coefficients_file do
+    header =
+      "# generate this whole file by running the following command:\n"+
+      "# rake regression:generate_calculated_corrected_deduct_coefficients_file > lib/pavement_condition_index/lookups/calculated_corrected_deduct_coefficients.rb\n"+
+      "\n"+ 
+      "module PavementConditionIndex\n"+
+      "  module Lookups\n"+
+      "    class CalculatedCorrectedDeductCoefficients\n"+
+      "\n"+
+      "# Output of `rake regression:generate_corrected_deduct_coefficients`\n"+
+      "COEFFICIENTS = "
+    puts header
+
+    Rake::Task["regression:generate_corrected_deduct_coefficients"].invoke
+
+    footer = 
+      "    end\n"+
+      "  end\n"+
+      "end\n"
+    puts footer
+  end  
+end
+
+def regress_polynomial x, y, degree
+  rows = x.map do |i|
+    (0..degree).map { |power| (i ** power).to_f }
+  end
+  mx, my = Matrix.rows(rows), Matrix.columns([y])
+  ((mx.transpose * mx).inv * mx.transpose * my).transpose.row(0).to_a
+end
+
+def evaluate_polynomial(*coefficients)
+  Proc.new {|x| coefficients.map.with_index{|c,i| c*(x**i)}.reduce(:+)}
+end

data/pavement_condition_index.gemspec

@@ -0,0 +1,24 @@
+# coding: utf-8
+lib = File.expand_path("../lib", __FILE__)
+$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+require "pavement_condition_index/version"
+
+Gem::Specification.new do |spec|
+  spec.name          = "pavement_condition_index"
+  spec.version       = PavementConditionIndex::VERSION
+  spec.authors       = ["Andy Monroe", "Nate Clark"]
+  spec.email         = ["hello@brandnewbox.com"]
+
+  spec.summary       = %q{Summary}
+  spec.homepage      = "https://github.com/brandnewbox/pavement_condition_index"
+  spec.license       = "MIT"
+
+  spec.files         = `git ls-files -z`.split("\x0").reject { |f| f.match(%r{^(test|spec|features)/}) }
+  spec.bindir        = "exe"
+  spec.executables   = spec.files.grep(%r{^exe/}) { |f| File.basename(f) }
+  spec.require_paths = ["lib"]
+
+  spec.add_development_dependency "bundler", "~> 1.15"
+  spec.add_development_dependency "rake", "~> 10.0"
+  spec.add_development_dependency "rspec", "~> 3.0"
+end

metadata

@@ -0,0 +1,114 @@
+--- !ruby/object:Gem::Specification
+name: pavement_condition_index
+version: !ruby/object:Gem::Version
+  version: 1.0.5
+platform: ruby
+authors:
+- Andy Monroe
+- Nate Clark
+autorequire: 
+bindir: exe
+cert_chain: []
+date: 2017-08-01 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: bundler
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.15'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.15'
+- !ruby/object:Gem::Dependency
+  name: rake
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '10.0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '10.0'
+- !ruby/object:Gem::Dependency
+  name: rspec
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '3.0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '3.0'
+description: 
+email:
+- hello@brandnewbox.com
+executables: []
+extensions: []
+extra_rdoc_files: []
+files:
+- ".gitignore"
+- ".rspec"
+- ".travis.yml"
+- Gemfile
+- LICENSE.txt
+- README.md
+- Rakefile
+- bin/console
+- bin/setup
+- lib/pavement_condition_index.rb
+- lib/pavement_condition_index/cdv_iteration.rb
+- lib/pavement_condition_index/core_ext/numeric.rb
+- lib/pavement_condition_index/distress_group.rb
+- lib/pavement_condition_index/lookups/calculated_corrected_deduct_coefficients.rb
+- lib/pavement_condition_index/lookups/calculated_deduct_coefficients.rb
+- lib/pavement_condition_index/lookups/corrected_deduct_values.rb
+- lib/pavement_condition_index/lookups/deduct_values.rb
+- lib/pavement_condition_index/lookups/observed_corrected_deduct_values.rb
+- lib/pavement_condition_index/lookups/observed_deduct_values.rb
+- lib/pavement_condition_index/pci.rb
+- lib/pavement_condition_index/sample_unit_condition_survey/asphalt_survey.rb
+- lib/pavement_condition_index/sample_unit_condition_survey/base_survey.rb
+- lib/pavement_condition_index/sample_unit_condition_survey/concrete_survey.rb
+- lib/pavement_condition_index/section_condition_survey.rb
+- lib/pavement_condition_index/version.rb
+- lib/tasks/regression.rake
+- pavement_condition_index.gemspec
+homepage: https://github.com/brandnewbox/pavement_condition_index
+licenses:
+- MIT
+metadata: {}
+post_install_message: 
+rdoc_options: []
+require_paths:
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+required_rubygems_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+requirements: []
+rubyforge_project: 
+rubygems_version: 2.5.1
+signing_key: 
+specification_version: 4
+summary: Summary
+test_files: []