num4regana 0.0.1-java

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA256:
+  metadata.gz: 616a5ce4aa14748f31750e7b1c2e69ca9af26295e7a173a8454ff9b3bc83cbfe
+  data.tar.gz: 74b2370bb501a69ef6546840d3ed0b61315e13ef35f584d5dbe3f2f38725cbb2
+SHA512:
+  metadata.gz: 7133ab852f61abbcb0e93bbc1be555bcfe16dfab0751ff1026be21dc8a46d829b3d03ba94313a7a196a0f385b7dbe01a6cb16d4678337bb063662e3faac88b43
+  data.tar.gz: a9851b753e8bf686a1a5bccc5dc392dc229822e51e47a3b1b82c8f1f76dbaf8629c04ecc825679aee250ff1225a7cf9ab6ed2f774f61e927f32fcb58ba83e7ee

data/CHANGELOG.md

@@ -0,0 +1,11 @@
+# Change Log
+
+## Unreleased
+
+## [0.0.1] - 2024-06-27
+
+### Fixed
+- fix first fixed.
+
+
+

data/Gemfile

@@ -0,0 +1,7 @@
+source "https://rubygems.org"
+
+gemspec
+platforms :jruby do
+    gem "rake-compiler",    ">= 1.2.5"
+end
+

data/LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2024 siranovel
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.

data/Rakefile

@@ -0,0 +1,10 @@
+require  'rake/javaextensiontask'
+
+jars = Dir.glob("lib/*.jar")
+Rake::JavaExtensionTask.new(name='num4regana')  do | ext |
+  ext.release          = '11'
+  ext.classpath        = jars.map { |x| File.expand_path x }.join ":"
+end
+task :default => [:compile]
+
+

data/ext/num4regana/OLSMultRegAna.java

@@ -0,0 +1,75 @@
+import org.apache.commons.math3.stat.regression.OLSMultipleLinearRegression;
+
+public class OLSMultRegAna {
+    private static OLSMultRegAna regana = new OLSMultRegAna();
+    public static OLSMultRegAna getInstance() {
+        return regana;
+    }
+    public LineReg lineRegAna(double[] yi, double xij[][]) {
+        LineRegAna line = new LineRegAna();
+
+        return line.lineRegAna(yi, xij);
+    }
+    public double getR2(double[] yi, double xij[][]) {
+        LineRegAna line = new LineRegAna();
+
+        return line.getR2(yi, xij);
+    }
+    public double getAdjR2(double[] yi, double xij[][]) {
+        LineRegAna line = new LineRegAna();
+
+        return line.getAdjR2(yi, xij);
+    }
+    /*********************************/
+    /* interface define              */
+    /*********************************/
+    /*********************************/
+    /* class define                  */
+    /*********************************/
+    public class LineReg {
+        private double a   = 0.0;
+        private double[] b = null;
+        public LineReg(double[] b) {
+            this.a = b[0];
+            this.b = new double[b.length - 1];
+            for (int i = 0; i < this.b.length; i++) {
+                this.b[i] = b[i + 1];
+            }
+        }
+        public double getIntercept() {
+            return a;
+        }
+        public double[] getSlope() {
+            return b;
+        }
+    }
+    private class LineRegAna {
+        // 最小２乗法
+        public LineReg lineRegAna(double[] yi, double xij[][]) {
+            OLSMultipleLinearRegression regression = new OLSMultipleLinearRegression();
+
+            regression.newSampleData(yi, xij);
+
+            double[] beta = regression.estimateRegressionParameters();
+
+            LineReg ret = new LineReg(beta);
+            return ret;
+        }
+        // 決定係数取得
+        public double getR2(double[] yi, double xij[][]) {
+            OLSMultipleLinearRegression regression = new OLSMultipleLinearRegression();
+
+            regression.newSampleData(yi, xij);
+            return regression.calculateRSquared();
+        }
+        // 自由度調整済み決定係数
+        public double getAdjR2(double[] yi, double xij[][]) {
+            OLSMultipleLinearRegression regression = new OLSMultipleLinearRegression();
+
+            regression.newSampleData(yi, xij);
+            return regression.calculateAdjustedRSquared();
+        }
+        
+    }
+}
+

data/ext/num4regana/SmplRegAna.java

@@ -0,0 +1,77 @@
+import org.apache.commons.math3.stat.regression.SimpleRegression;
+
+public class SmplRegAna {
+    private static SmplRegAna regana = new SmplRegAna();
+    public static SmplRegAna getInstance() {
+        return regana;
+    }
+    public LineReg lineRegAna(double[] yi, double xi[]) {
+        LineRegAna line = new LineRegAna();
+
+        return line.lineRegAna(yi, xi);
+    }
+    public double getR2(double[] yi, double xi[]) {
+        LineRegAna line = new LineRegAna();
+
+        return line.getR2(yi, xi);
+    }
+    public double getR(double[] yi, double xi[]) {
+        LineRegAna line = new LineRegAna();
+
+        return line.getR(yi, xi);
+    }
+    /*********************************/
+    /* interface define              */
+    /*********************************/
+    /*********************************/
+    /* class define                  */
+    /*********************************/
+    public class LineReg {
+        private double a = 0.0;
+        private double b = 0.0;
+        public LineReg(double a, double b) {
+            this.a = a;
+            this.b = b;
+        }
+        public double getIntercept() {
+            return a;
+        }
+        public double getSlope() {
+            return b;
+        }
+    }
+    private class LineRegAna {
+        // 単回帰の分析
+        public LineReg lineRegAna(double[] yi, double xi[]) {
+            SimpleRegression simpleReg = new SimpleRegression(true);
+
+            for(int i = 0; i < yi.length; i++) {
+                simpleReg.addData(xi[i], yi[i]);
+            }
+            LineReg ret = new LineReg(simpleReg.getIntercept(), simpleReg.getSlope());
+
+            return ret;
+        }
+        // 決定係数取得
+        public double getR2(double[] yi, double xi[]) {
+            SimpleRegression simpleReg = new SimpleRegression(true);
+
+            for(int i = 0; i < yi.length; i++) {
+                simpleReg.addData(xi[i], yi[i]);
+            }
+
+            return simpleReg.getRSquare();
+        }
+        // 相関係数取得
+        public double getR(double[] yi, double xi[]) {
+            SimpleRegression simpleReg = new SimpleRegression(true);
+
+            for(int i = 0; i < yi.length; i++) {
+                simpleReg.addData(xi[i], yi[i]);
+            }
+
+            return simpleReg.getR();
+        }
+    }
+}
+

data/lib/num4regana.rb

@@ -0,0 +1,165 @@
+require 'java'
+require 'num4regana.jar'
+require 'commons-math3-3.6.1.jar'
+
+java_import 'SmplRegAna'
+java_import 'OLSMultRegAna'
+# 回帰分析
+#  (Apache commoms math3使用)
+module Num4RegAnaLib
+    # 単回帰分析
+    class SmplRegAnaLib
+        def initialize
+            @regana = SmplRegAna.getInstance()
+        end
+        # 単回帰分析
+        #
+        # @overload line_reg_ana(yi, xi)
+        #   @param [Array] yi yの値(double[])
+        #   @param [Array] xi xの値(double[])
+        #   @return [Hash] (intercept:定数項 slope:回帰係数)
+        # @example
+        #   yi = [286, 851, 589, 389, 158, 1037, 463, 563, 372, 1020]
+        #   xi = [107, 336, 233,  82,  61,  378, 129, 313, 142,  428]
+        #   regana = Num4RegAnaLib::SmplRegAnaLib.new
+        #   regana.line_reg_ana(yi, xi)
+        #   => 
+        #     {
+        #        "intercept":  99.075, # 定数項
+        #        "slope":      2.145,  # 回帰係数
+        #     }
+        def line_reg_ana(yi, xi)
+            ret = @regana.lineRegAna(yi.to_java(Java::double), xi.to_java(Java::double))
+            retRb = {
+                "intercept":  ret.getIntercept(), # 定数項
+                "slope":      ret.getSlope(),     # 回帰係数
+            }
+            return retRb
+        end
+        # 決定係数
+        #
+        # @overload getr2(yi, xi)
+        #   @param [Array] yi yの値(double[])
+        #   @param [Array] xi xの値(double[])
+        #   @return [double] 決定係数
+        # @example
+        #   yi = [286, 851, 589, 389, 158, 1037, 463, 563, 372, 1020]
+        #   xi = [107, 336, 233,  82,  61,  378, 129, 313, 142,  428]
+        #   regana = Num4RegAnaLib::SmplRegAnaLib.new
+        #   regana.getr2(yi, xi)
+        #   => 0.893
+        def getr2(yi, xi)
+            return @regana.getR2(yi.to_java(Java::double), xi.to_java(Java::double))
+        end
+        # 相関係数
+        #
+        # @overload getr(yi, xi)
+        #   @param [Array] yi yの値(double[])
+        #   @param [Array] xi xの値(double[])
+        #   @return [double] 決定係数
+        # @example
+        #   yi = [286, 851, 589, 389, 158, 1037, 463, 563, 372, 1020]
+        #   xi = [107, 336, 233,  82,  61,  378, 129, 313, 142,  428]
+        #   regana = Num4RegAnaLib::SmplRegAnaLib.new
+        #   regana.getr(yi, xi)
+        #   => 0.945
+        def getr(yi, xi)
+            return @regana.getR(yi.to_java(Java::double), xi.to_java(Java::double))
+        end
+    end
+    # 重回帰分析(最小２乗法)
+    class OLSMultRegAnaLib
+        def initialize
+            @regana = OLSMultRegAna.getInstance()
+        end
+        # 重回帰分析
+        #
+        # @overload line_reg_ana(yi, xij)
+        #   @param [Array] yi yの値(double[])
+        #   @param [Array] xij xの値(double[][])
+        #   @return [Hash] (intercept:定数項 slope:回帰係数)
+        # @example
+        #   olsyi = [45, 38, 41, 34, 59, 47, 35, 43, 54, 52]
+        #   olsxij = [
+        #       [17.5, 30],
+        #       [17.0, 25],
+        #       [18.5, 20],
+        #       [16.0, 30],
+        #       [19.0, 45],
+        #       [19.5, 35],
+        #       [16.0, 25],
+        #       [18.0, 35],
+        #       [19.0, 35],
+        #       [19.5, 40],
+        #   ]
+        #   regana = Num4RegAnaLib::OLSMultRegAnaLib.new
+        #   regana.line_reg_ana(olsyi, olsxij)
+        #   => 
+        #     {
+        #        "intercept":  -34.71,           # 定数項
+        #        "slope":      [3.47, 0.53],     # 回帰係数
+        #     }
+        def line_reg_ana(yi, xij)
+            ret = @regana.lineRegAna(yi.to_java(Java::double), xij.to_java(Java::double[]))
+            retRb = {
+                "intercept":  ret.getIntercept(), # 定数項
+                "slope":      ret.getSlope().to_a,     # 回帰係数
+            }
+            return retRb
+        end
+        # 決定係数
+        #
+        # @overload getr2(yi, xij)
+        #   @param [Array] yi yの値(double[])
+        #   @param [Array] xij xの値(double[][])
+        #   @return [double] 決定係数
+        # @example
+        #   olsyi = [45, 38, 41, 34, 59, 47, 35, 43, 54, 52]
+        #   olsxij = [
+        #       [17.5, 30],
+        #       [17.0, 25],
+        #       [18.5, 20],
+        #       [16.0, 30],
+        #       [19.0, 45],
+        #       [19.5, 35],
+        #       [16.0, 25],
+        #       [18.0, 35],
+        #       [19.0, 35],
+        #       [19.5, 40],
+        #   ]
+        #   regana = Num4RegAnaLib::OLSMultRegAnaLib.new
+        #   regana.getr2(yi, xi)
+        #   => 0.858
+        def getr2(yi, xij)
+            return @regana.getR2(yi.to_java(Java::double), xij.to_java(Java::double[]))
+        end
+        # 自由度調整済み決定係数
+        #
+        # @overload getadjr2(yi, xij)
+        #   @param [Array] yi yの値(double[])
+        #   @param [Array] xij xの値(double[][])
+        #   @return [double] 決定係数
+        # @example
+        #   olsyi = [45, 38, 41, 34, 59, 47, 35, 43, 54, 52]
+        #   olsxij = [
+        #       [17.5, 30],
+        #       [17.0, 25],
+        #       [18.5, 20],
+        #       [16.0, 30],
+        #       [19.0, 45],
+        #       [19.5, 35],
+        #       [16.0, 25],
+        #       [18.0, 35],
+        #       [19.0, 35],
+        #       [19.5, 40],
+        #   ]
+        #   regana = Num4RegAnaLib::OLSMultRegAnaLib.new
+        #   regana.getadjr2(yi, xij)
+        #   => 0.8176
+        def getadjr2(yi, xij)
+            return @regana.getAdjR2(yi.to_java(Java::double), xij.to_java(Java::double[]))
+        end
+
+    end
+end
+

metadata

@@ -0,0 +1,91 @@
+--- !ruby/object:Gem::Specification
+name: num4regana
+version: !ruby/object:Gem::Version
+  version: 0.0.1
+platform: java
+authors:
+- siranovel
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2024-06-27 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: rake
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '12.3'
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 12.3.3
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '12.3'
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 12.3.3
+- !ruby/object:Gem::Dependency
+  name: rake-compiler
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.2'
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 1.2.5
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.2'
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 1.2.5
+description: numerical solution for regression analysis.
+email: siranovel@gmail.com
+executables: []
+extensions:
+- Rakefile
+extra_rdoc_files: []
+files:
+- CHANGELOG.md
+- Gemfile
+- LICENSE
+- Rakefile
+- ext/num4regana/OLSMultRegAna.java
+- ext/num4regana/SmplRegAna.java
+- lib/commons-math3-3.6.1.jar
+- lib/num4regana.rb
+homepage: http://github.com/siranovel/num4regana
+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: []
+rubygems_version: 3.3.7
+signing_key: 
+specification_version: 4
+summary: num for regression analysis
+test_files: []