loganb-kder 0.1.2 → 0.1.3

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA1:
-  metadata.gz: 17148d0bed58d5e6a0c2c7d93e93322d4bb5497b
-  data.tar.gz: 5f4f7639de0cbd1d16526c46d4c083b904f9442e
+  metadata.gz: e5b8c12f1724491f1d7e1030dca865ca131aff67
+  data.tar.gz: dd269d0ef838ebe5d6c45a62288f2b688971a1fe
 SHA512:
-  metadata.gz: f2d329f6b49f4426b1cfb290801fc0845c10f93a45a97b1eaa87e9653e006513e403958ccc117b8ba350023d6ce94005baba5a078a67f84df596274d81006934
-  data.tar.gz: 83aa31ddc5861e47b54a9333c59594c591669a71c530a8d0986aa1a68e5545ff930e6b6bb1e0dd77d7f30427c49d7a6a768da5a69381c8128eb1a3a3e997c344
+  metadata.gz: 1964cd8f5e3446644fc5f26e623c953a94e56dd2f0efac20b2ec42164d458f71038eba37c521034166909b983875e82d6ca66b74f3a71dd625b2168d2d9bdc01
+  data.tar.gz: '07429b569c86eeef572489945412a478ef4329162fb0325b520ed6ac15bb8cabca511d10ee2fd994f72330cf804e1d3b443813bfcea28a0a425f8d003ba629fe'

data/lib/kder/kde.rb

@@ -1,4 +1,4 @@
-require 'set'
+require 'matrix'
 
 class Kder
   require_relative 'bandwidth'
@@ -83,5 +83,50 @@ class Kder
       output << [max,0]
       output.compact.transpose      
     end
+
+    #
+    # Smooths multiple data sets with a shared x-axis at the same time. 
+    #
+    # values in the form of: 
+    #
+    # [ 
+    #   [x1, y1, z1, k1],
+    #   [x2, y2, 0, z3], # "Hole" where there is no data for z
+    #   ...
+    # ]
+    #
+    #
+    def kdemultivec(values, bw, opts = {}.freeze)
+      opts = {sigmas: Sigmas, sampling_density: MeshCount, threshold: MinimumThresholdValue, minimum_delta: DifferenceThreshold}.merge(opts)
+
+      # Initialization steps
+      range = bw*opts[:sigmas]
+      min = values.first[0] - range
+      max = values.last[0]  + range
+      step_size = (max-min)/(opts[:sampling_density].to_f)
+      step_size = step_size < MinimumStepSize ? MinimumStepSize : step_size
+
+      # initialize the range variables
+      ranges = (min..max).step(step_size).to_a
+      output = [Array.new(values.first.length-1, 0).unshift(min)]
+      old_intensities = Vector::elements(Array.new(values.first.length-1, 0))
+      # Step through the range
+      ranges[1..-1].map do |mid|
+        high_end  = mid + range
+        lower_end = mid - range
+        selection_range = (lower_end..high_end)
+        included = values.select {|v| selection_range.include?(v[0])}
+
+        intensities = included.collect { |vec| Vector::elements(vec[1..-1],false) * Kder::Statistics.custom_pdf(vec[0] - mid, bw) }.
+          inject(&:+) || old_intensities 
+
+        unless intensities.all? { |v| v < opts[:threshold] } or (intensities - old_intensities).all? { |v| v.abs < opts[:minimum_delta] }
+          output << intensities.to_a.unshift(mid) 
+          old_intensities = intensities
+        end
+      end
+      output << Array.new(values.first.length - 1, 0).unshift(max)
+      output.compact.transpose
+    end
   end
 end

data/lib/kder/version.rb

@@ -1,3 +1,3 @@
 class Kder
-  VERSION = "0.1.2"
+  VERSION = "0.1.3"
 end

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: loganb-kder
 version: !ruby/object:Gem::Version
-  version: 0.1.2
+  version: 0.1.3
 platform: ruby
 authors:
 - Ryan Taylor
@@ -9,7 +9,7 @@ authors:
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2017-11-15 00:00:00.000000000 Z
+date: 2017-11-18 00:00:00.000000000 Z
 dependencies: []
 description: A KDE implementation in Ruby.
 email:
@@ -44,7 +44,7 @@ required_rubygems_version: !ruby/object:Gem::Requirement
       version: '0'
 requirements: []
 rubyforge_project: 
-rubygems_version: 2.5.2
+rubygems_version: 2.5.2.1
 signing_key: 
 specification_version: 4
 summary: Kernel Density Estimation, and an associated bandwidth estimator, rudimentary