spatial_stats 0.1.1 → 0.2.1

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: spatial_stats
 version: !ruby/object:Gem::Version
-  version: 0.1.1
+  version: 0.2.1
 platform: ruby
 authors:
 - Keith Doggett
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2020-03-29 00:00:00.000000000 Z
+date: 2020-04-03 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: numo-narray
@@ -129,8 +129,8 @@ files:
 - lib/spatial_stats/global/stat.rb
 - lib/spatial_stats/local.rb
 - lib/spatial_stats/local/bivariate_moran.rb
-- lib/spatial_stats/local/g.rb
 - lib/spatial_stats/local/geary.rb
+- lib/spatial_stats/local/getis_ord.rb
 - lib/spatial_stats/local/moran.rb
 - lib/spatial_stats/local/multivariate_geary.rb
 - lib/spatial_stats/local/stat.rb
@@ -150,7 +150,8 @@ files:
 homepage: https://www.github.com/keithdoggett/spatial_stats
 licenses:
 - MIT
-metadata: {}
+metadata:
+  documentation_uri: https://keithdoggett.github.io/spatial_stats/docs/index.html
 post_install_message: 
 rdoc_options: []
 require_paths:

data/lib/spatial_stats/local/g.rb

@@ -1,75 +0,0 @@
-# frozen_string_literal: true
-
-module SpatialStats
-  module Local
-    class G < Stat
-      def initialize(scope, field, weights, star = false)
-        super(scope, field, weights)
-        @star = star
-      end
-      attr_accessor :star
-      attr_writer :x, :z_lag
-
-      def i
-        x.each_with_index.map do |_x_val, idx|
-          i_i(idx)
-        end
-      end
-
-      def i_i(idx)
-        x_lag[idx] / denominators[idx]
-      end
-
-      def x
-        @x ||= SpatialStats::Queries::Variables.query_field(@scope, @field)
-      end
-      alias z x
-
-      def star?
-        @star ||= weights.full.trace.positive?
-      end
-
-      private
-
-      def w
-        @w ||= begin
-          if star?
-            # TODO: try to fix this because it will still likely be a
-            # bottleneck in mc testing
-            weights.full.windowed.row_standardized
-          else
-            weights.standardized
-          end
-        end
-      end
-
-      def z_lag
-        # window if star is true
-        @z_lag ||= begin
-          if star?
-            SpatialStats::Utils::Lag.window_sum(w, x)
-          else
-            SpatialStats::Utils::Lag.neighbor_sum(w, x)
-          end
-        end
-      end
-      alias x_lag z_lag
-
-      def denominators
-        @denominators ||= begin
-          n = w.shape[0]
-          if star?
-            [x.sum] * n
-          else
-            # add everything but i
-            (0..n - 1).each.map do |idx|
-              terms = x.dup
-              terms.delete_at(idx)
-              terms.sum
-            end
-          end
-        end
-      end
-    end
-  end
-end