num4anova 0.0.13-java → 0.0.14-java

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 84f7f49327e10a3da88a363daa4d6d5dec1dd4da00827d6f8e5c60ee6db4bc08
-  data.tar.gz: 495495a96d218034014971770b654431bb0df61aec47f8f3171fe351285022c4
+  metadata.gz: 844b392125f78a4966a67dd1162bef6676857f03ed9f49e48d4f00e9f294d6e9
+  data.tar.gz: da4999753cf5b2b73336b365c53e79ad5eb084fcb44151fdf6c05b5c691a40d2
 SHA512:
-  metadata.gz: 1b46cc7670940eb88343959c80bad55209a981388acdd823c3d9eeafa72a629d2319cfb49707f3f998d72f0d080480785129c948c8ebedb388693476be59add1
-  data.tar.gz: 83585b2c13c1f964b8bf92e507a005deb4d358a612e23b939cb9d9911175006a1378e16719c1b73073834aa510f245b4432f0aba9f5dfcdceb78ee2cf5638ba3
+  metadata.gz: 6c8e9884146305d0db99633f36997de0dc3ec4593e2b497fa6f97c1d9cf957721ad233914454e2eb1f8334b86871fdb389f073d62492904432da76415f9a302f
+  data.tar.gz: 79b9cc886c7f2f8b794c204db714ca5663e9b9acb0388eee8afd60cef393b9f107e3f32f05855f848c53e39fc67f42212d352d87c0fdc2262718178f6d36d1d1

data/CHANGELOG.md

@@ -2,6 +2,11 @@
 
 ## Unreleased
 
+## [0.0.14] - 2024-04-08
+
+### add
+- add fuction of bonferrono_test in NonParametrixTestLib.
+
 ## [0.0.13] - 2024-04-04
 
 ### chg

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

metadata

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