rumale-metric_learning 0.24.0 → 0.26.0

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: de8906421117ba7ceb261e21f9a469cdda66c4a582acddef8b7c25d9e81b5b2c
-  data.tar.gz: 1f89eca43ee6a34c4bfecb3779d6b7e2ba210833834e46105ab9d56cfbe3ef91
+  metadata.gz: 7f20a7e0bf72c2a9869a837c493121ce2ec70132d9818f5bfca3dd2b8dd7db94
+  data.tar.gz: af5c7339a4ef9164223c18b6e18c9453b081b90be787cc9a23ccb4fc9493566f
 SHA512:
-  metadata.gz: 6cd556eb2d5060f16f766fae6595b3cf302723f38788ab67680ad053306642a86ef8745e100d387059cc304ad1f07cc9e7e6493bddaa6bc9e066c2ebda5c0772
-  data.tar.gz: 77afbc98ad84566a932193ddb3c774b5edff63a7b7747958cf651ecbff41c1826e0571e9f6232e5261520ecf7545c9e3e8195680ca1f18aff5506233d0d5b682
+  metadata.gz: ff8a56b94c0834938085284cbc4290390554ebbe030a0a69dd19fc3efc0ff198c35686656e93e0b4b734c120ce88b4da4a7ef3180128b87bf0acff12c56eb2cb
+  data.tar.gz: 96eddc81041593518dd8139114e89b59634d658fc20bad5ebdde297a5517914b561182028fe18359da3df83b9c204f6aecd7262208c28fd1caba5423bdfa36cb

data/LICENSE.txt

@@ -1,4 +1,4 @@
-Copyright (c) 2022 Atsushi Tatsuma
+Copyright (c) 2022-2023 Atsushi Tatsuma
 All rights reserved.
 
 Redistribution and use in source and binary forms, with or without

data/lib/rumale/metric_learning/local_fisher_discriminant_analysis.rb

@@ -0,0 +1,115 @@
+# frozen_string_literal: true
+
+require 'rumale/base/estimator'
+require 'rumale/base/transformer'
+require 'rumale/validation'
+
+module Rumale
+  module MetricLearning
+    # LocalFisherDiscriminantAnalysis is a class that implements Local Fisher Discriminant Analysis.
+    #
+    # @example
+    #   require 'rumale/metric_learning/local_fisher_discriminant_analysis'
+    #
+    #   transformer = Rumale::MetricLearning::LocalFisherDiscriminantAnalysis.new
+    #   transformer.fit(training_samples, traininig_labels)
+    #   low_samples = transformer.transform(testing_samples)
+    #
+    # *Reference*
+    # - Sugiyama, M., "Local Fisher Discriminant Analysis for Supervised Dimensionality Reduction," Proc. ICML'06, pp. 905--912, 2006.
+    class LocalFisherDiscriminantAnalysis < ::Rumale::Base::Estimator
+      include ::Rumale::Base::Transformer
+
+      # Returns the transform matrix.
+      # @return [Numo::DFloat] (shape: [n_components, n_features])
+      attr_reader :components
+
+      # Return the class labels.
+      # @return [Numo::Int32] (shape: [n_classes])
+      attr_reader :classes
+
+      # Create a new transformer with LocalFisherDiscriminantAnalysis.
+      #
+      # @param n_components [Integer] The number of components.
+      # @param gamma [Float] The parameter of rbf kernel, if nil it is 1 / n_features.
+      def initialize(n_components: nil, gamma: nil)
+        super()
+        @params = {
+          n_components: n_components,
+          gamma: gamma
+        }
+      end
+
+      # Fit the model with given training data.
+      #
+      # @param x [Numo::DFloat] (shape: [n_samples, n_features]) The training data to be used for fitting the model.
+      # @param y [Numo::Int32] (shape: [n_samples]) The labels to be used for fitting the model.
+      # @return [LocalFisherDiscriminantAnalysis] The learned classifier itself.
+      def fit(x, y)
+        unless enable_linalg?(warning: false)
+          raise 'LocalFisherDiscriminatAnalysis#fit requires Numo::Linalg but that is not loaded.'
+        end
+
+        x = Rumale::Validation.check_convert_sample_array(x)
+        y = Rumale::Validation.check_convert_label_array(y)
+        Rumale::Validation.check_sample_size(x, y)
+
+        # initialize some variables.
+        n_samples, n_features = x.shape
+        @classes = Numo::Int32[*y.to_a.uniq.sort]
+        n_components = @params[:n_components] || n_features
+        @params[:gamma] ||= 1.fdiv(n_features)
+        affinity_mat = Rumale::PairwiseMetric.rbf_kernel(x, nil, @params[:gamma])
+        affinity_mat[affinity_mat.diag_indices] = 1.0
+
+        # calculate within and mixture scatter matricies.
+        class_mat = Numo::DFloat.zeros(n_samples, n_samples)
+        within_weight_mat = Numo::DFloat.zeros(n_samples, n_samples)
+        @classes.each do |label|
+          pos = y.eq(label)
+          n_class_samples = pos.count
+          pos_vec = Numo::DFloat.cast(pos)
+          pos_mat = pos_vec.outer(pos_vec)
+          class_mat += pos_mat
+          within_weight_mat += pos_mat * 1.fdiv(n_class_samples)
+        end
+
+        mixture_weight_mat = ((affinity_mat - 1) / n_samples) * class_mat + 1.fdiv(n_samples)
+        within_weight_mat *= affinity_mat
+        mixture_weight_mat = mixture_weight_mat.sum(axis: 1).diag - mixture_weight_mat
+        within_weight_mat = within_weight_mat.sum(axis: 1).diag - within_weight_mat
+
+        # calculate components.
+        mixture_mat = x.transpose.dot(mixture_weight_mat.dot(x))
+        within_mat = x.transpose.dot(within_weight_mat.dot(x))
+        _, evecs = Numo::Linalg.eigh(mixture_mat, within_mat, vals_range: (n_features - n_components)...n_features)
+        comps = evecs.reverse(1).transpose.dup
+        @components = n_components == 1 ? comps[0, true].dup : comps.dup
+        self
+      end
+
+      # Fit the model with training data, and then transform them with the learned model.
+      #
+      # @param x [Numo::DFloat] (shape: [n_samples, n_features]) The training data to be used for fitting the model.
+      # @param y [Numo::Int32] (shape: [n_samples]) The labels to be used for fitting the model.
+      # @return [Numo::DFloat] (shape: [n_samples, n_components]) The transformed data
+      def fit_transform(x, y)
+        x = Rumale::Validation.check_convert_sample_array(x)
+        y = Rumale::Validation.check_convert_label_array(y)
+        Rumale::Validation.check_sample_size(x, y)
+
+        fit(x, y).transform(x)
+      end
+
+      # Transform the given data with the learned model.
+      #
+      # @param x [Numo::DFloat] (shape: [n_samples, n_features]) The data to be transformed with the learned model.
+      # @return [Numo::DFloat] (shape: [n_samples, n_components]) The transformed data.
+      def transform(x)
+        x = Rumale::Validation.check_convert_sample_array(x)
+
+        x.dot(@components.transpose)
+      end
+    end
+  end
+end

data/lib/rumale/metric_learning/version.rb

@@ -5,6 +5,6 @@ module Rumale
   # Module for metric learning algorithms.
   module MetricLearning
     # @!visibility private
-    VERSION = '0.24.0'
+    VERSION = '0.26.0'
   end
 end

data/lib/rumale/metric_learning.rb

@@ -3,6 +3,7 @@
 require 'numo/narray'
 
 require_relative 'metric_learning/fisher_discriminant_analysis'
+require_relative 'metric_learning/local_fisher_discriminant_analysis'
 require_relative 'metric_learning/mlkr'
 require_relative 'metric_learning/neighbourhood_component_analysis'
 require_relative 'metric_learning/version'

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: rumale-metric_learning
 version: !ruby/object:Gem::Version
-  version: 0.24.0
+  version: 0.26.0
 platform: ruby
 authors:
 - yoshoku
 autorequire:
 bindir: exe
 cert_chain: []
-date: 2022-12-31 00:00:00.000000000 Z
+date: 2023-02-19 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: lbfgsb
@@ -44,28 +44,28 @@ dependencies:
     requirements:
     - - "~>"
       - !ruby/object:Gem::Version
-        version: 0.24.0
+        version: 0.26.0
   type: :runtime
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
     - - "~>"
       - !ruby/object:Gem::Version
-        version: 0.24.0
+        version: 0.26.0
 - !ruby/object:Gem::Dependency
   name: rumale-decomposition
   requirement: !ruby/object:Gem::Requirement
     requirements:
     - - "~>"
       - !ruby/object:Gem::Version
-        version: 0.24.0
+        version: 0.26.0
   type: :runtime
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
     - - "~>"
       - !ruby/object:Gem::Version
-        version: 0.24.0
+        version: 0.26.0
 description: |
   Rumale::MetricLearning provides metric learning algorithms,
   such as Fisher Discriminant Analysis and Neighboourhood Component Analysis
@@ -80,6 +80,7 @@ files:
 - README.md
 - lib/rumale/metric_learning.rb
 - lib/rumale/metric_learning/fisher_discriminant_analysis.rb
+- lib/rumale/metric_learning/local_fisher_discriminant_analysis.rb
 - lib/rumale/metric_learning/mlkr.rb
 - lib/rumale/metric_learning/neighbourhood_component_analysis.rb
 - lib/rumale/metric_learning/version.rb