num4anova 0.0.2-java → 0.0.3-java

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 3987057ab3576c0da8a66256901cf7548ad7bf4bef7aaafccf128b083655333d
-  data.tar.gz: 4b35674f595de9f05b4056ec6c20d671549341aad88083e5dfd928a881de5c7a
+  metadata.gz: 89f91d9c210bb998f14a726f1c7147470cba19fba721724625b7ec87e4c88ad9
+  data.tar.gz: a79f6e3ca4adc650193065ae62037cfed0580e5f7b9c5e1d06109d74f17a05d6
 SHA512:
-  metadata.gz: 75fd29b81f0c79fa4f72520ced191bccc543d1cf86875425f9408a7c25aeeb8dc0b3417a3af3838bbc2dadee641da79f1fd43b46dc3a02f7e920ffafc1f6e1e7
-  data.tar.gz: 4a9226cc5d20f86fda993eecc91048edc77a7157da8b2d44c305fdcbf9849aeb0c2c6d4c270b623d62a061daee67aff8864e039280282fc8c2a415453f8e8334
+  metadata.gz: 97a1409dfb791766292e8d6c0a49502555911529632d45b72aa6c68bc6fda7e82ed813fff4c005c9c76745a2ad05fe79d1131e23879da18afcc88d0652617c29
+  data.tar.gz: 4cd28879971daf44d20ff8122120327e8a74519f658a22e052724af316387c057cfb315a6d1cb23324bbd08634c066c1045594d6c8804f1f87e12d31b6acd0fb

data/CHANGELOG.md

@@ -2,6 +2,11 @@
 
 ## Unreleased
 
+## [0.0.3] - 2024-01-26
+
+### add
+- add function of dunnet_test.
+
 ## [0.0.2] - 2024-01-23
 
 ### add

data/ext/num4anova/MultiComp.java

@@ -23,6 +23,28 @@ public class MultiComp {
             double[][] statistic = hypoth.calcTestStatistic(xi);
             return hypoth.executeTest(statistic, a * 0.5);
         }
+        public boolean scheffe_test(double[][] xi, double a) {
+            return false;
+        }
+
+        public boolean[][] twosideTest(double[][] xi, double a) {
+            HypothesisTest hypoth = new TwoSideTest();
+            double[][] statistic = hypoth.calcTestStatistic(xi);
+
+            return hypoth.executeTest(statistic, a / 2.0);
+        }
+        public boolean[][] rightsideTest(double[][] xi, double a) {
+            HypothesisTest hypoth = new RightSideTest();
+            double[][] statistic = hypoth.calcTestStatistic(xi);
+
+            return hypoth.executeTest(statistic, a);
+        }
+        public boolean[][] leftsideTest(double[][] xi, double a) {
+            HypothesisTest hypoth = new LeftSideTest();
+            double[][] statistic = hypoth.calcTestStatistic(xi);
+
+            return hypoth.executeTest(statistic, a);
+        }
         /*********************************/
         /* interface define              */
         /*********************************/
@@ -145,7 +167,6 @@ public class MultiComp {
                 TDistribution tDist = new TDistribution(v);
                 double t = 
                     tDist.inverseCumulativeProbability(p);
-
                 return Math.sqrt(2) * t;
             }
         }
@@ -220,6 +241,127 @@ public class MultiComp {
                 return sumSq / na;
             } 
         }
+        // ダネット法
+        private class DunnetTest{
+            private int k = 0;
+            private int v = 0;
+            private double[] mean = null;
+            private double[] n = null;
+            protected int getK() { return k;}
+            protected int getV() { return v;}
+            public double[][] calcTestStatistic(double[][] xi) {
+                k = xi.length;
+                mean = new double[k];
+                n = new double[k];
+                double[][] statistic = new double[k][k];
+                double ve = calcVe(xi);
+
+                for(int i = 0; i < k; i++) {
+                    for(int j = 0; j < k; j++) {
+                        statistic[i][j] = (mean[j] - mean[i]) 
+                                        / Math.sqrt(ve * (1.0 / n[j] + 1.0 / n[i]));
+                    }
+                }                
+                return statistic;
+            }
+            private double calcVe(double[][] xi) {
+                double sumSq = 0.0;
+                int sumN = 0;
+                for(int i = 0; i < k; i++) {
+                    DescriptiveStatistics stat = new DescriptiveStatistics();
+                    Arrays.stream(xi[i]).forEach(stat::addValue);
+                    mean[i] = stat.getMean();
+                    n[i] = stat.getN();
+                    sumSq += (n[i] - 1) * stat.getVariance();
+                    sumN += n[i];
+                    stat.clear();
+                }
+                v = sumN - k;
+                return sumSq / v;
+            }
+        }        
+        private class TwoSideTest extends DunnetTest 
+                                  implements HypothesisTest {
+            public boolean[][] executeTest(double[][] statistic, double a) {
+                int v = super.getV();
+                int k = super.getK();
+                double den = k - 1;
+                double p = 1.0 - a / den;
+                TDistribution tDist = new TDistribution(v);
+                double l_val = tDist.inverseCumulativeProbability(a  / den);
+                double r_val = tDist.inverseCumulativeProbability(1.0 - a / den);
+                boolean[][] ret = new boolean[k][k];
+
+                for(int i = 0; i < k; i++) {
+                    for(int j = 0; j < k; j++) {
+                       ret[i][j] = evaluation(statistic[i][j], l_val, r_val ); 
+                    }
+                }
+                return ret;
+            }
+            private boolean evaluation(double statistic, double l_val, double r_val) {
+                boolean ret = true;
+
+                if ((l_val < statistic) && (statistic < r_val)) {
+                    ret = false;
+                }
+                return ret;
+            }
+        }
+        private class RightSideTest extends DunnetTest 
+                                  implements HypothesisTest {
+            public boolean[][] executeTest(double[][] statistic, double a) {
+                int v = super.getV();
+                int k = super.getK();
+                double den = k - 1;
+                double p = 1.0 - a / den;
+                TDistribution tDist = new TDistribution(v);
+                double r_val = tDist.inverseCumulativeProbability(1.0 - a);
+                boolean[][] ret = new boolean[k][k];
+
+                for(int i = 0; i < k; i++) {
+                    for(int j = 0; j < k; j++) {
+                       ret[i][j] = evaluation(statistic[i][j], r_val ); 
+                    }
+                }
+                return ret;
+            }
+            private boolean evaluation(double statistic, double r_val) {
+                boolean ret = true;
+
+                if (statistic < r_val) {
+                    ret = false;
+                }
+                return ret;
+            }
+        }
+        private class LeftSideTest extends DunnetTest 
+                                  implements HypothesisTest {
+            public boolean[][] executeTest(double[][] statistic, double a) {
+                int v = super.getV();
+                int k = super.getK();
+                double den = k - 1;
+                double p = a / den;
+                TDistribution tDist = new TDistribution(v);
+                double l_val = tDist.inverseCumulativeProbability(a);
+                boolean[][] ret = new boolean[k][k];
+
+                for(int i = 0; i < k; i++) {
+                    for(int j = 0; j < k; j++) {
+                       ret[i][j] = evaluation(statistic[i][j], l_val ); 
+                    }
+                }
+                return ret;
+            }
+            private boolean evaluation(double statistic, double l_val) {
+                boolean ret = true;
+
+                if (l_val < statistic) {
+                    ret = false;
+                }
+                return ret;
+            }
+        }
     }
 }
 

data/lib/dunnet.rb

@@ -0,0 +1,89 @@
+# Dunnet検定
+#  (Apache commoms math3使用)
+module DunnetTestLib
+    # Dunnet検定の両側検定
+    #
+    # @overload twoside_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],
+    #    ]
+    #    paraTest = MultiCompLib::ParametrixTestLib.new
+    #    paraTest.twoside_test(xi, 0.05)
+    #    => 
+    #      res = [
+    #          [false, false, false, true, false],
+    #          [false, false, false, true, false],
+    #          [false, false, false, false, false],
+    #          [true,  true,  false, false, false],
+    #          [false, false, false, false, false],
+    #      ]
+    def twoside_test(xi, a)
+        ret = @paramTest.twosideTest(xi.to_java(Java::double[]), a)
+        return ret.to_a
+    end
+    # Dunnet検定の右側検定
+    #
+    # @overload rightside_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],
+    #    ]
+    #    paraTest = MultiCompLib::ParametrixTestLib.new
+    #    paraTest.rightside_test(xi, 0.05)
+    #    => 
+    #      res = [
+    #          [false, false, true,  true,  true],
+    #          [false, false, false, true,  false],
+    #          [false, false, false, true,  false],
+    #          [false, false, false, false, false],
+    #          [false, false, false, true, false],
+    #      ]
+    def rightside_test(xi, a)
+        ret = @paramTest.rightsideTest(xi.to_java(Java::double[]), a)
+        return ret.to_a
+    end
+    # Dunnet検定の左側検定
+    #
+    # @overload leftside_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],
+    #    ]
+    #    paraTest = MultiCompLib::ParametrixTestLib.new
+    #    paraTest.leftside_test(xi, 0.05)
+    #    => 
+    #      res = [
+    #          [false, false, false, false, false],
+    #          [false, false, false, false, false],
+    #          [true,  false, false, false, false],
+    #          [true,  true,  true,  false, true],
+    #          [true,  false, false, false, false],
+    #      ]
+    def leftside_test(xi, a)
+        ret = @paramTest.leftsideTest(xi.to_java(Java::double[]), a)
+        return ret.to_a
+    end
+end
+

data/lib/multicomp.rb

@@ -3,12 +3,15 @@ require 'num4anova.jar'
 require 'jfreechart-1.5.4.jar'
 require 'commons-math3-3.6.1.jar'
 
+require_relative('dunnet')
+
 java_import 'MultiComp'
 # 多重比較を行う
 #  (Apache commoms math3使用)
 module MultiCompLib
     # パラメトリック検定
     class ParametrixTestLib
+        include DunnetTestLib
         def initialize
             @paramTest = MultiComp::ParametrixTest.getInstance()
         end
@@ -26,7 +29,8 @@ module MultiCompLib
         #        [26.4, 32.5, 31.3],
         #        [24.5, 21.2, 22.4],
         #    ]
-        #    paraTest.turkey_test(xi, a)
+        #    paraTest = MultiCompLib::ParametrixTestLib.new
+        #    paraTest.turkey_test(xi, 0.05)
         #    => 
         #      [
         #        [false, false, false, true, false],
@@ -53,7 +57,8 @@ module MultiCompLib
         #        [26.4, 32.5, 31.3],
         #        [24.5, 21.2, 22.4],
         #    ]
-        #    paraTest.bonferrono_test(xi, a)
+        #    paraTest = MultiCompLib::ParametrixTestLib.new
+        #    paraTest.bonferrono_test(xi, 0.05)
         #    => 
         #      [
         #        [false, false, false, true, false],

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: num4anova
 version: !ruby/object:Gem::Version
-  version: 0.0.2
+  version: 0.0.3
 platform: java
 authors:
 - siranovel
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2024-01-23 00:00:00.000000000 Z
+date: 2024-01-26 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: rake
@@ -64,6 +64,7 @@ files:
 - ext/num4anova/MultiComp.java
 - ext/num4anova/OneWayLayout.java
 - lib/commons-math3-3.6.1.jar
+- lib/dunnet.rb
 - lib/jcommon-1.0.23.jar
 - lib/jfreechart-1.5.4.jar
 - lib/multicomp.rb