num4regana 0.0.2-java → 0.0.4-java

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: ad6fa72619eb7a02bdc561494b9aae197c37b6a8558fd21106bdf15210a9af81
-  data.tar.gz: c7a55ac44970c164dc7379a11955f29416392ceba4605118f8bcbba459724e94
+  metadata.gz: d2c06b0a7caeac79ec2e1d451dccabd387449af22574e60520125374bd4296d7
+  data.tar.gz: 5744160fa6bec668de85a0ae0fd4b71625b6f4976a9b01e831ca71b4d654ec8c
 SHA512:
-  metadata.gz: dd9e09e45d35384a3bfa9eb50a0132f4b3b3e970c2f97d241264a18a08986c3a46753d2292caf10c0ae1fcab8ec4f45f6456c97fe14df2ee9334168cd429a60a
-  data.tar.gz: 1cf8a55a5444c5e1e1a151adbb4d142bceaa9344cde28a859848b5642ce36660b6b919410907f8964105ec35df991591670ba8722658a708822e0840d790ce23
+  metadata.gz: 7030429dc48f2807211cf07446cd6fbcee4094fe55410d3d54ef01c9151ead39123ef869da446e6497db1f7e95ea59810281111a0f2a05922cdd3396fc9da6ba
+  data.tar.gz: 3967a10a750896e1b1a8f6bba7de42ae6a88862db18af456ba49f0018fc9e9e9d931aa792aceaa3e37d1c281dabfb52e461de484b0f349e59570e8e8ae73f3fe

data/CHANGELOG.md

@@ -2,6 +2,18 @@
 
 ## Unreleased
 
+## [0.0.4] - 2024-09-13
+
+### add
+- add function of getAIC
+
+## [0.0.3] - 2024-09-05
+
+### add
+- add GLM
+  logistic regression analystis
+  poisson regression analystis
+
 ## [0.0.2] - 2024-08-08
 
 ### fix

data/ext/num4regana/AbstratGLM.java

@@ -0,0 +1,62 @@
+import java.util.Arrays;
+
+abstract class AbstratGLM {
+    private final double eta = 0.005;
+    abstract double regression(double[] b, double[] xi);
+    abstract double linkFunc(double q);
+    // 勾配降下法
+    protected double[] grand_metod(double[] yi, double[] b, double[][] xij) {
+        // 交差エントロピー計算
+        double[] ei = calcE(yi, b, xij);
+
+        // パラメータ更新
+        for(int i = 0; i < ei.length; i++) {
+            b[i] -= eta * ei[i];
+        }
+        return b;
+    }
+    // AIC
+    protected double calcAIC(double[] b, double[][] xij) {
+        // 尤度計算
+        double maxL = calcL(b,xij); 
+        int k = 1 + xij[0].length;
+
+        return -2 * (maxL - k);
+    }
+    // 交差エントロピー計算
+    private double[] calcE(double[] yi, double[] b, double[][] xij) {
+        double[] xi = new double[1 + xij[0].length];
+        double[] ei = new double[1 + xij[0].length];
+
+        Arrays.fill(ei, 0.0);
+        for(int i = 0; i < yi.length; i++) {
+            xi[0] = 1.0;
+            System.arraycopy(xij[i], 0, xi, 1, xij[0].length);
+
+            double q = regression(b, xi);
+            double p = linkFunc(q);
+
+            for(int j = 0; j < xi.length; j++) {                
+                ei[j] += (p - yi[i]) * xi[j];
+            }
+        }
+
+        return ei;
+    }
+    // 尤度計算(パラメータ)
+    private double calcL(double[] b, double[][] xij) {
+        double l = 0.0;
+        double[] xi = new double[1 + xij[0].length];
+
+        for(int i = 0; i < xij.length; i++) {
+            xi[0] = 1.0;
+            System.arraycopy(xij[i], 0, xi, 1, xij[0].length);
+            double q = regression(b, xi);
+            double p = linkFunc(q);
+
+            l += Math.log(p);
+        }
+        return l;
+    }   
+}
+

data/ext/num4regana/AbstratGLMM.java

@@ -0,0 +1,18 @@
+import java.util.Arrays;
+import org.apache.commons.math3.distribution.BetaDistribution;
+
+abstract class AbstratGLMM {
+    abstract double rereion(double[] b, double[] xi, double r);
+    abstract double linkFunc(double q);
+    protected double[] mcmc(double[] yi, double[] b, double[][] xij) {
+        double[] bnew = new double[1 + xij[0].length];
+        BetaDistribution beDist = new BetaDistribution(1,1);
+
+        for(int i= 0; i < bnew.length; i++) {
+System.out.printf("%f ", beDist.sample());
+        }
+        System.out.println();
+        return null;
+    }    
+}
+

data/ext/num4regana/LogitRegAna.java

@@ -0,0 +1,68 @@
+import java.util.Arrays;
+import java.util.Map;
+
+public class LogitRegAna extends AbstratGLM {
+    private final int NUM = 1000;
+    private static LogitRegAna regana = new LogitRegAna();
+    public static LogitRegAna getInstance() {
+        return regana;
+    }
+    public LineReg nonLineRegAna(double[] yi, double xij[][]) {
+        double[] b = initB(xij[0].length);
+         
+        for  (int i = 0; i < NUM; i++) {
+            b = grand_metod(yi, b, xij);
+        }
+        return new LineReg(b);
+    }
+    public double getAIC(Map<String, Object> regCoe, double[][] xij) {
+        double[] b = new double[1 + xij[0].length];
+
+        b[0] = (double)regCoe.get("intercept");
+        System.arraycopy(regCoe.get("slope"), 0, b, 1, xij[0].length);
+        return calcAIC(b, xij);
+    }
+    private double[] initB(int xsie) {
+        double[] b = new double[1 + xsie];
+         
+        Arrays.fill(b, 0.0);
+        return b;
+    }
+    // q = b0 + b1 * x0
+    double regression(double[] b, double[] xi) {
+        double ret = 0.0;
+
+        for(int i = 0; i < xi.length; i++) {
+            ret += b[i] * xi[i];
+        }
+        return ret;
+    }
+    // p = 1 / (1 + exp( -q))
+    double linkFunc(double q) {
+        return 1.0 / (1.0 + Math.exp(-1.0 * q));
+    }
+    /*********************************/
+    /* 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;
+        }
+    }
+}
+

data/ext/num4regana/MultRegAna.java

@@ -12,33 +12,23 @@ public class MultRegAna {
         return regana;
     }
     public LineReg lineRegAna(double[] yi, double xij[][]) {
-        double[][] data = createData(yi, xij);
-        LineRegAna line = createLineRegAna(data);
+        LineRegAna line = createLineRegAna(yi, xij);
 
         return line.lineRegAna(yi, xij);
     }
     public double getR2(double[] yi, double xij[][]) {
-        double[][] data = createData(yi, xij);
-        LineRegAna line = createLineRegAna(data);
+        LineRegAna line = createLineRegAna(yi, xij);
 
         return line.getR2(yi, xij);
     }
     public double getAdjR2(double[] yi, double xij[][]) {
-        double[][] data = createData(yi, xij);
-        LineRegAna line = createLineRegAna(data);
+        LineRegAna line = createLineRegAna(yi, xij);
 
         return line.getAdjR2(yi, xij);
     }
-    private double[][] createData(double[] yi, double xij[][]) {
-        double[][] data = new double[yi.length][1 + xij[0].length];
+    private LineRegAna createLineRegAna(double[] yi, double xij[][]) {
+        double[][] data = createData(yi, xij);
 
-        for (int i = 0; i < yi.length; i++) {
-            data[i][0] = yi[i];
-            System.arraycopy(xij[i], 0, data[i], 1, xij[0].length);
-        }
-        return data;
-    }
-    private LineRegAna createLineRegAna(double data[][]) {
         // 等分散性の検定
         if (false == bartletTest(data)) {  // 等分散性
             return new OLSMultRegAna();
@@ -47,6 +37,15 @@ public class MultRegAna {
             return new GLSMultRegAna(data);
         }
     }
+    private double[][] createData(double[] yi, double xij[][]) {
+        double[][] data = new double[yi.length][1 + xij[0].length];
+
+        for (int i = 0; i < yi.length; i++) {
+            data[i][0] = yi[i];
+            System.arraycopy(xij[i], 0, data[i], 1, xij[0].length);
+        }
+        return data;
+    }
     private boolean bartletTest(double data[][]) {
         OneWayAnovaTest anova = new BartletTest();
         double statistic = anova.calcTestStatistic(data);

data/ext/num4regana/PoissonRegAna.java

@@ -0,0 +1,68 @@
+import java.util.Arrays;
+import java.util.Map;
+
+public class PoissonRegAna extends AbstratGLM {
+    private final int NUM = 1000;
+    private static PoissonRegAna regana = new PoissonRegAna();
+    public static PoissonRegAna getInstance() {
+        return regana;
+    }
+    public LineReg nonLineRegAna(double[] yi, double[][] xij) {
+        double[] b = initB(xij[0].length);
+
+        for  (int i = 0; i < NUM; i++) {
+            b = grand_metod(yi, b, xij);
+        }
+        return new LineReg(b);
+    }
+    public double getAIC(Map<String, Object> regCoe, double[][] xij) {
+        double[] b = new double[1 + xij[0].length];
+
+        b[0] = (double)regCoe.get("intercept");
+        System.arraycopy(regCoe.get("slope"), 0, b, 1, xij[0].length);
+        return calcAIC(b, xij);
+    }
+    private double[] initB(int xsie) {
+        double[] b = new double[1 + xsie];
+         
+        Arrays.fill(b, 0.0);
+        return b;
+    }
+    // q = b0 + b1 * x0
+    double regression(double[] b, double[] xi) {
+        double ret = 0.0;
+
+        for(int i = 0; i < xi.length; i++) {
+            ret += b[i] * xi[i];
+        }
+        return ret;
+    }
+    // p = exp(q)
+    double linkFunc(double q) {
+        return Math.exp(q);
+    }
+    /*********************************/
+    /* 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;
+        }
+    }
+}
+

data/ext/num4regana/ProBitRegAna.java

@@ -0,0 +1,71 @@
+import java.util.Arrays;
+import org.apache.commons.math3.distribution.NormalDistribution;
+import java.util.Map;
+
+public class ProBitRegAna extends AbstratGLM {
+    private final int NUM = 1000;
+    private static ProBitRegAna regana = new ProBitRegAna();
+    private NormalDistribution ndist = new NormalDistribution(0, 1);
+    
+    public static ProBitRegAna getInstance() {
+        return regana;
+    }
+    public LineReg nonLineRegAna(double[] yi, double[][] xij) {
+        double[] b = initB(xij[0].length);
+         
+        for  (int i = 0; i < NUM; i++) {
+            b = grand_metod(yi, b, xij);
+        }
+        return new LineReg(b);
+    }
+    public double getAIC(Map<String, Object> regCoe, double[][] xij) {
+        double[] b = new double[1 + xij[0].length];
+
+        b[0] = (double)regCoe.get("intercept");
+        System.arraycopy(regCoe.get("slope"), 0, b, 1, xij[0].length);
+        return calcAIC(b, xij);
+    }
+    private double[] initB(int xsie) {
+        double[] b = new double[1 + xsie];
+         
+        Arrays.fill(b, 0.0);
+        return b;
+    }
+    // q = b0 + b1 * x0
+    double regression(double[] b, double[] xi) {
+        double ret = 0.0;
+
+        for(int i = 0; i < xi.length; i++) {
+            ret += b[i] * xi[i];
+        }
+        return ret;
+    }
+    // 
+    double linkFunc(double q) {
+        return ndist.cumulativeProbability(q);
+    }
+    /*********************************/
+    /* 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;
+        }
+    }
+}
+

data/lib/num4glmregana.rb

@@ -0,0 +1,287 @@
+require 'java'
+require 'num4regana.jar'
+require 'commons-math3-3.6.1.jar'
+
+java_import 'LogitRegAna'
+java_import 'PoissonRegAna'
+java_import 'ProBitRegAna'
+java_import 'java.util.HashMap'
+
+# 一般化線形回帰分析
+#  (Apache commoms math3使用)
+module Num4GLMRegAnaLib
+    # (2項)ロジスティック回帰分析
+    class LogitRegAnaLib
+        def initialize
+            @multana = LogitRegAna.getInstance()
+        end
+        # (2項)ロジスティック回帰分析
+        #
+        # @overload non_line_reg_ana(yi, xij)
+        #   @param [Array] yi yの値(double[])
+        #   @param [Array] xij xの値(double[][])
+        #   @return [Hash] (intercept:定数項 slope:回帰係数)
+        # @example
+        #   glsyi = [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
+        #   glsxij = [
+        #        [1, 24],
+        #        [1, 18],
+        #        [0, 15],
+        #        [1, 16],
+        #        [0, 10],
+        #        [1, 26],
+        #        [1, 2],
+        #        [0, 24],
+        #        [1, 18],
+        #        [1, 22],
+        #        [1, 3],
+        #        [1, 6],
+        #        [0, 15],
+        #        [0, 12],
+        #        [1, 6],
+        #        [0, 6],
+        #        [1, 12],
+        #        [0, 12],
+        #        [1, 18],
+        #        [1, 3],
+        #        [1, 8],
+        #        [0, 9],
+        #        [0, 12],
+        #        [0, 6],
+        #        [0, 8],
+        #        [1, 12],
+        #   ]
+        #   regana = Num4RegAnaLib::LogitRegAnaLib.new
+        #   regana.non_line_reg_ana(glsyi, glsxij)
+        #   => 
+        #     {
+        #        "intercept":  -6.2313,              # 定数項
+        #        "slope":      [2.5995, 0.1652],     # 回帰係数
+        #     }
+        def non_line_reg_ana(yi, xij)
+            multRet = @multana.nonLineRegAna(yi.to_java(Java::double), xij.to_java(Java::double[]))
+            retRb = {
+                "intercept":  multRet.getIntercept(), # 定数項
+                "slope":      multRet.getSlope().to_a,     # 回帰係数
+            }
+            return retRb
+        end
+        # AIC
+        #
+        # @overload get_aic(regcoe, xij)
+        #   @param [Hash] regcoe 回帰係数(intercept:定数項 slope:回帰係数)
+        #   @param [Array] xij xの値(double[][])
+        #   @return double AIC値
+        # @example
+        #   reg = {
+        #        :intercept=>  -6.2313,    # 定数項
+        #        :slope=>      [2.5995, 0.1652],     # 回帰係数
+        #   }
+        #   xij = [
+        #        [1, 24],
+        #        [1, 18],
+        #        [0, 15],
+        #        [1, 16],
+        #        [0, 10],
+        #        [1, 26],
+        #        [1, 2],
+        #        [0, 24],
+        #        [1, 18],
+        #        [1, 22],
+        #        [1, 3],
+        #        [1, 6],
+        #        [0, 15],
+        #        [0, 12],
+        #        [1, 6],
+        #        [0, 6],
+        #        [1, 12],
+        #        [0, 12],
+        #        [1, 18],
+        #        [1, 3],
+        #        [1, 8],
+        #        [0, 9],
+        #        [0, 12],
+        #        [0, 6],
+        #        [0, 8],
+        #        [1, 12],
+        #   ]
+        #   regana = Num4RegAnaLib::LogitRegAnaLib.new
+        #   regana.get_aic(reg, xij)
+        #   => 155.612
+        def get_aic(regcoe, xij)
+            o = HashMap.new
+            o["intercept"] = regcoe[:intercept]
+            o["slope"]     = regcoe[:slope].to_java(Java::double)
+            @multana.getAIC(o, xij.to_java(Java::double[]))
+        end
+    end
+    # ポアソン回帰分析
+    class PoissonRegAnaLib
+        def initialize
+            @multana = PoissonRegAna.getInstance()
+        end
+        # ポアソン回帰分析
+        #
+        # @overload non_line_reg_ana(yi, xij)
+        #   @param [Array] yi yの値(double[])
+        #   @param [Array] xij xの値(double[][])
+        #   @return [Hash] (intercept:定数項 slope:回帰係数)
+        # @example
+        #   glsyi = [4, 10, 7, 14]
+        #   glsxij = [
+        #       [1],
+        #       [2],
+        #       [3],
+        #       [4],
+        #   ]
+        #   regana = Num4RegAnaLib::PoissonRegAnaLib.new
+        #   regana.non_line_reg_ana(glsyi, glsxij)
+        #   => 
+        #     {
+        #        "intercept":  1.3138,    # 定数項
+        #        "slope":      [0.3173],  # 回帰係数
+        #     }
+        def non_line_reg_ana(yi, xij)
+            multRet = @multana.nonLineRegAna(yi.to_java(Java::double), xij.to_java(Java::double[]))
+            retRb = {
+                "intercept":  multRet.getIntercept(), # 定数項
+                "slope":     multRet.getSlope().to_a,     # 回帰係数
+            }
+            return retRb
+        end
+        # AIC
+        #
+        # @overload get_aic(regcoe, xij)
+        #   @param [Hash] regcoe 回帰係数(intercept:定数項 slope:回帰係数)
+        #   @param [Array] xij xの値(double[][])
+        #   @return double AIC値
+        # @example
+        #   reg = {
+        #        :intercept => 1.3138,    # 定数項
+        #        :slope    =>  [0.3173],  # 回帰係数
+        #   }
+        #   xij = [
+        #       [1],
+        #       [2],
+        #       [3],
+        #       [4],
+        #   ]
+        #   regana = Num4RegAnaLib::PoissonRegAnaLib.new
+        #   regana.get_aic(reg, xij)
+        #   => -12.856
+        def get_aic(regcoe, xij)
+            o = HashMap.new
+            o["intercept"] = regcoe[:intercept]
+            o["slope"]     = regcoe[:slope].to_java(Java::double)
+            @multana.getAIC(o, xij.to_java(Java::double[]))
+        end
+    end
+    # プロビット回帰分析
+    class ProBitRegAnaLib
+        def initialize
+            @multana = ProBitRegAna.getInstance()
+        end
+        # プロビット回帰分析
+        #
+        # @overload non_line_reg_ana(yi, xij)
+        #   @param [Array] yi yの値(double[])
+        #   @param [Array] xij xの値(double[][])
+        #   @return [Hash] (intercept:定数項 slope:回帰係数)
+        # @example
+        #   glsyi = [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
+        #   glsxij = [
+        #        [1, 24],
+        #        [1, 18],
+        #        [0, 15],
+        #        [1, 16],
+        #        [0, 10],
+        #        [1, 26],
+        #        [1, 2],
+        #        [0, 24],
+        #        [1, 18],
+        #        [1, 22],
+        #        [1, 3],
+        #        [1, 6],
+        #        [0, 15],
+        #        [0, 12],
+        #        [1, 6],
+        #        [0, 6],
+        #        [1, 12],
+        #        [0, 12],
+        #        [1, 18],
+        #        [1, 3],
+        #        [1, 8],
+        #        [0, 9],
+        #        [0, 12],
+        #        [0, 6],
+        #        [0, 8],
+        #        [1, 12],
+        #   ]
+        #   regana = Num4RegAnaLib::ProBitRegAnaLib.new
+        #   regana.non_line_reg_ana(glsyi, glsxij)
+        #   => 
+        #     {
+        #        "intercept":  -5.0497,    # 定数項
+        #        "slope":      [2.2379, 0.2973],     # 回帰係数
+        #     }
+        def non_line_reg_ana(yi, xij)
+            multRet = @multana.nonLineRegAna(yi.to_java(Java::double), xij.to_java(Java::double[]))
+            retRb = {
+                "intercept":  multRet.getIntercept(), # 定数項
+                "slope":      multRet.getSlope().to_a,     # 回帰係数
+            }
+            return retRb
+        end
+        # AIC
+        #
+        # @overload get_aic(regcoe, xij)
+        #   @param [Hash] regcoe 回帰係数(intercept:定数項 slope:回帰係数)
+        #   @param [Array] xij xの値(double[][])
+        #   @return double AIC値
+        # @example
+        #   reg = {
+        #        :intercept=>  -5.0497,    # 定数項
+        #        :slope=>      [2.2379, 0.2973],     # 回帰係数
+        #   }
+        #   xij = [
+        #        [1, 24],
+        #        [1, 18],
+        #        [0, 15],
+        #        [1, 16],
+        #        [0, 10],
+        #        [1, 26],
+        #        [1, 2],
+        #        [0, 24],
+        #        [1, 18],
+        #        [1, 22],
+        #        [1, 3],
+        #        [1, 6],
+        #        [0, 15],
+        #        [0, 12],
+        #        [1, 6],
+        #        [0, 6],
+        #        [1, 12],
+        #        [0, 12],
+        #        [1, 18],
+        #        [1, 3],
+        #        [1, 8],
+        #        [0, 9],
+        #        [0, 12],
+        #        [0, 6],
+        #        [0, 8],
+        #        [1, 12],
+        #   ]
+        #   regana = Num4RegAnaLib::ProBitRegAnaLib.new
+        #   regana.get_aic(reg, xij)
+        #   => 119.674
+        def get_aic(regcoe, xij)
+            o = HashMap.new
+            o["intercept"] = regcoe[:intercept]
+            o["slope"]     = regcoe[:slope].to_java(Java::double)
+            @multana.getAIC(o, xij.to_java(Java::double[]))
+        end
+    end
+end
+
+
+

data/lib/num4lineregana.rb

@@ -0,0 +1,165 @@
+require 'java'
+require 'num4regana.jar'
+require 'commons-math3-3.6.1.jar'
+
+java_import 'SmplRegAna'
+java_import 'MultRegAna'
+# 線形回帰分析
+#  (Apache commoms math3使用)
+module Num4LineRegAnaLib
+    # 単回帰分析
+    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
+    # 重回帰分析(最小２乗法:等分散性checkあり)
+    class OLSMultRegAnaLib
+        def initialize
+            @multana = MultRegAna.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)
+            multRet = @multana.lineRegAna(yi.to_java(Java::double), xij.to_java(Java::double[]))
+
+            retRb = {
+                "intercept":  multRet.getIntercept(), # 定数項
+                "slope":      multRet.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 @multana.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 @multana.getAdjR2(yi.to_java(Java::double), xij.to_java(Java::double[]))
+        end
+    end
+end
+

data/lib/num4regana.rb

@@ -1,166 +1,10 @@
-require 'java'
-require 'num4regana.jar'
-require 'commons-math3-3.6.1.jar'
+require_relative('num4lineregana')
+require_relative('num4glmregana')
 
-java_import 'SmplRegAna'
-java_import 'MultRegAna'
 # 回帰分析
-#  (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
-    # 重回帰分析(最小２乗法:等分散性checkあり)
-    class OLSMultRegAnaLib
-        def initialize
-            @multana = MultRegAna.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)
-            multRet = @multana.lineRegAna(yi.to_java(Java::double), xij.to_java(Java::double[]))
-
-            retRb = {
-                "intercept":  multRet.getIntercept(), # 定数項
-                "slope":      multRet.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 @multana.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 @multana.getAdjR2(yi.to_java(Java::double), xij.to_java(Java::double[]))
-        end
-
-    end
+    include Num4LineRegAnaLib
+    include Num4GLMRegAnaLib
 end
 
+

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: num4regana
 version: !ruby/object:Gem::Version
-  version: 0.0.2
+  version: 0.0.4
 platform: java
 authors:
 - siranovel
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2024-08-08 00:00:00.000000000 Z
+date: 2024-09-13 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: rake
@@ -61,9 +61,16 @@ files:
 - Gemfile
 - LICENSE
 - Rakefile
+- ext/num4regana/AbstratGLM.java
+- ext/num4regana/AbstratGLMM.java
+- ext/num4regana/LogitRegAna.java
 - ext/num4regana/MultRegAna.java
+- ext/num4regana/PoissonRegAna.java
+- ext/num4regana/ProBitRegAna.java
 - ext/num4regana/SmplRegAna.java
 - lib/commons-math3-3.6.1.jar
+- lib/num4glmregana.rb
+- lib/num4lineregana.rb
 - lib/num4regana.rb
 homepage: http://github.com/siranovel/num4regana
 licenses: