num4anova 0.0.12-java → 0.0.14-java

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 87febfd8031fd924430223c93fad6b0606aa212ab109e1b758db2c1d23b6dace
-  data.tar.gz: 48b610b93d2c5341f05c25ef230acea09b02fa75c924129f5dafba56fcb6051f
+  metadata.gz: 844b392125f78a4966a67dd1162bef6676857f03ed9f49e48d4f00e9f294d6e9
+  data.tar.gz: da4999753cf5b2b73336b365c53e79ad5eb084fcb44151fdf6c05b5c691a40d2
 SHA512:
-  metadata.gz: 979fe122d0cf17086e80a45ff6c44b67153822a8d89d83bceb84a25d3d68190f087c3ed8bf4b192e9afd48a38b3ec0d1456fa0821c84256113154e04d5a6e461
-  data.tar.gz: 0bcad9f7a234a96759697440f540f0caddf70c1a25cf46409dead11478faa0017cfe9e98d2a23feaad2648bea641d5bca80cb73adab5001ec46d93ef8899ce70
+  metadata.gz: 6c8e9884146305d0db99633f36997de0dc3ec4593e2b497fa6f97c1d9cf957721ad233914454e2eb1f8334b86871fdb389f073d62492904432da76415f9a302f
+  data.tar.gz: 79b9cc886c7f2f8b794c204db714ca5663e9b9acb0388eee8afd60cef393b9f107e3f32f05855f848c53e39fc67f42212d352d87c0fdc2262718178f6d36d1d1

data/CHANGELOG.md

@@ -2,6 +2,16 @@
 
 ## Unreleased
 
+## [0.0.14] - 2024-04-08
+
+### add
+- add fuction of bonferrono_test in NonParametrixTestLib.
+
+## [0.0.13] - 2024-04-04
+
+### chg
+- chg fuction of interval_estim.
+
 ## [0.0.12] - 2024-03-27
 
 ### add

data/ext/num4anova/Ancova.java

@@ -24,7 +24,7 @@ public class Ancova {
         double statistic = hypoth.calcTestStatistic(xi);
         return hypoth.executeTest(statistic, a);
     }
-    public Interval intervalEstim(double[][][] xi, double a) {
+    public Interval[] intervalEstim(double[][][] xi, double a) {
         Estim estim = new IntervalEstim();
 
         return estim.calcInterval(xi, a);
@@ -37,7 +37,7 @@ public class Ancova {
         boolean executeTest(double statistic, double a);
     }
     private interface Estim {
-        Interval calcInterval(double[][][] xi, double a);
+        Interval[] calcInterval(double[][][] xi, double a);
     }
     /*********************************/
     /* class define                  */
@@ -346,7 +346,8 @@ public class Ancova {
         private int n = 0;
         private int[] ni = null;
         private double sumex = 0.0;
-        public Interval calcInterval(double[][][] xi, double a) {
+        public Interval[] calcInterval(double[][][] xi, double a) {
+            Interval[] ret = new Interval[xi.length];
             ni = calcNi(xi);
             int sumn = calcSumn(xi);
             n = sumn - xi.length - 1;
@@ -360,12 +361,16 @@ public class Ancova {
 
             TDistribution tDist = new TDistribution(n);
             double t = tDist.inverseCumulativeProbability(1.0 - a / 2.0);
-            double wk = (meanxi[0] - meanx);
-            double wk2 = t * Math.sqrt((1/ni[0] + wk * wk / sumex) * ve);
-            double min = meanyi[0] - b * wk - wk2;
-            double max = meanyi[0] - b * wk + wk2;
+            for(int i = 0; i < ret.length; i++) {
+                double wk = (meanxi[i] - meanx);
+                double wk2 = t * Math.sqrt((1/ni[i] + wk * wk / sumex) * ve);
+                double min = meanyi[i] - b * wk - wk2;
+                double max = meanyi[i] - b * wk + wk2;
 
-            return new Interval(min, max);
+                ret[i] = new Interval(min, max);
+            }
+            
+            return ret;
         }
         private int[] calcNi(double[][][] xi) {
             int[] ni = new int[xi.length];

data/ext/num4anova/MultiComp.java

@@ -3,9 +3,17 @@ import java.util.Arrays;
 import java.util.List;
 import java.util.ArrayList;
 import org.apache.commons.math3.distribution.TDistribution;
+import org.apache.commons.math3.stat.inference.MannWhitneyUTest;
 
 import org.apache.commons.math3.util.Combinations;
 public class MultiComp {
+    /*********************************/
+    /* interface define              */
+    /*********************************/
+    private interface HypothesisTest {
+        double[][] calcTestStatistic(double[][] xi);
+        boolean[][] executeTest(double[][] statistic, double a);
+    }
     public static class ParametrixTest {
         private static ParametrixTest paramTest = new ParametrixTest();
         public static ParametrixTest getInstance() {
@@ -46,13 +54,6 @@ public class MultiComp {
             return hypoth.executeTest(statistic, a);
         }
         /*********************************/
-        /* interface define              */
-        /*********************************/
-        private interface HypothesisTest {
-            double[][] calcTestStatistic(double[][] xi);
-            boolean[][] executeTest(double[][] statistic, double a);
-        }
-        /*********************************/
         /* Class define                  */
         /*********************************/
         // turkey法
@@ -366,5 +367,55 @@ public class MultiComp {
             }
         }
     }
+    public static class NonParametrixTest {
+        private static NonParametrixTest nonParamTest = new NonParametrixTest();
+        public static NonParametrixTest getInstance() {
+            return nonParamTest;
+        }
+        public boolean[][] bonferronoTest(double[][] xi, double a) {
+            HypothesisTest hypoth = new BonferroniTest();
+
+            double[][] statistic = hypoth.calcTestStatistic(xi);
+            return hypoth.executeTest(statistic, a * 0.5);
+        }
+        /*********************************/
+        /* Class define                  */
+        /*********************************/
+        // ボンフェロー法
+        private class BonferroniTest implements HypothesisTest {
+            private int n = 0;
+            private int k = 0;
+            public double[][] calcTestStatistic(double[][] xi) {
+                n = xi.length;
+                double[][] statistic = new double[n][n];
+                //
+                Combinations c = new Combinations(n, 2);
+                List<int[]> al = new ArrayList<>();
+                for(int[] iterate : c) {
+                    al.add(iterate);
+                }
+                k = al.size();
+                //
+                MannWhitneyUTest utest = new MannWhitneyUTest();
+                for(int[] array : al) {
+                    int i = array[0];
+                    int j = array[1];
+
+                    statistic[i][j] = utest.mannWhitneyUTest(xi[i], xi[j]);
+                }
+                return statistic;
+            }
+            public boolean[][] executeTest(double[][] statistic, double a) {
+                boolean[][] ret = new boolean[n][n];
+
+                for(int i = 0; i < n; i++) {
+                    for(int j = i+1; j < n; j++) {
+                        ret[i][j] = (statistic[i][j] < a * k) ? true : false;
+                    }
+                }
+                return ret;
+            }
+        }        
+    }
 }
 

data/lib/multicomp.rb

@@ -43,7 +43,7 @@ module MultiCompLib
             ret = @paramTest.turkeyTest(xi.to_java(Java::double[]), a)
             return ret.to_a
         end
-        # ボンフェロー二の不等式による多重比較
+        # ボンフェロー二の不等式による多重比較(T検定)
         #
         # @overload bonferrono_test(xi, a)
         #   @param [array]  xi データ(double[][])
@@ -74,5 +74,37 @@ module MultiCompLib
     end
     # ノンパラメトリック検定
     class NonParametrixTestLib
+        def initialize
+            @nonParamTest = MultiComp::NonParametrixTest.getInstance()
+        end
+        # ボンフェロー二の不等式による多重比較(マン・ホイットニーU検定)
+        #
+        # @overload bonferrono_test(xi, a)
+        #   @param [array]  xi データ(double[][])
+        #   @param [double] a         有意水準
+        #   @return [Array] 検定結果(boolean[][] true:棄却域内 false:棄却域外)
+        # @example
+        #    xi = [
+        #        [12.2, 18.8, 18.2],
+        #        [22.2, 20.5, 14.6],
+        #        [20.8, 19.5, 26.3],
+        #        [26.4, 32.5, 31.3],
+        #        [24.5, 21.2, 22.4],
+        #    ]
+        #    nonParaTest = MultiCompLib::NonParametrixTestLib.new
+        #    nonParaTest.bonferrono_test(xi, 0.05)
+        #    => 
+        #      [
+        #        [false, false, true,  true,  true],
+        #        [false, false, false, true,  true],
+        #        [false, false, false, true,  false],
+        #        [false, false, false, false, true],
+        #        [false, false, false, false, false],
+        #     ]
+        def bonferrono_test(xi, a)
+            ret = @nonParamTest.bonferronoTest(xi.to_java(Java::double[]), a)
+            return ret.to_a
+        end
+
     end
 end

data/lib/num4anova.rb

@@ -341,14 +341,24 @@ module Num4AnovaLib
         #    ancova = Num4AnovaLib::Num4AncovaLib.new
         #    ancova.interval_estim(xi, 0.05)
         #    =>
-        #      {:min=>4.466605469341916, :max=>7.1909253948556096}
+        #       {:min=>4.466605469341916,  :max=>7.1909253948556096}
+        #       {:min=>5.05699825110459,   :max=>6.386335082228742}
+        #       {:min=>2.510804295684195,  :max=>4.250430272217034}
+        #       {:min=>2.8089257316042135, :max=>2.9566298239513418}
+        #       {:min=>-6.303283147572267, :max=>-0.6577045067487104}
         def interval_estim(xi, a)
+            retRb = []
             retJava = @ancova.intervalEstim(xi.to_java(Java::double[][]), a)
-            retHash = {
-                "min": retJava.getMin(),
-                "max": retJava.getMax()
-            }
-            return retHash
+            sz = retJava.size
+            sz.times do |i|
+                retRb.push(
+                  {
+                    "min": retJava[i].getMin(),
+                    "max": retJava[i].getMax()
+                  }
+                )
+            end
+            return retRb
         end
     end
 end

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: num4anova
 version: !ruby/object:Gem::Version
-  version: 0.0.12
+  version: 0.0.14
 platform: java
 authors:
 - siranovel
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2024-03-27 00:00:00.000000000 Z
+date: 2024-04-08 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: rake