num4anova 0.0.14-java → 0.0.16-java

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 844b392125f78a4966a67dd1162bef6676857f03ed9f49e48d4f00e9f294d6e9
-  data.tar.gz: da4999753cf5b2b73336b365c53e79ad5eb084fcb44151fdf6c05b5c691a40d2
+  metadata.gz: fb1dd347cd1156e1cc9e06ee7b4daff32442d6fc721e35553cbb406ae5d94cf1
+  data.tar.gz: 17d4c4084a40fe8e721829169e2551455d258f64dd4ef6c537b499ef9a5dab15
 SHA512:
-  metadata.gz: 6c8e9884146305d0db99633f36997de0dc3ec4593e2b497fa6f97c1d9cf957721ad233914454e2eb1f8334b86871fdb389f073d62492904432da76415f9a302f
-  data.tar.gz: 79b9cc886c7f2f8b794c204db714ca5663e9b9acb0388eee8afd60cef393b9f107e3f32f05855f848c53e39fc67f42212d352d87c0fdc2262718178f6d36d1d1
+  metadata.gz: 768b64803a235265e1dae41341428a0bda9fb12245e7f2fae54a3f75f8aeee607a4cd3709b578b1c7aeb6869fcc7f227b287d2142b99411ed8955248eeafedd4
+  data.tar.gz: da4c649d32a38c9f9a83adadc54d20605bfebe42f822fdd353df7ba8e00df8e3be8fbd8c4775474da5299957a477ec98c4912ec1ac5b8fb0d305c0f11ad91ff0

data/CHANGELOG.md

@@ -2,6 +2,16 @@
 
 ## Unreleased
 
+## [0.0.16] - 2024-04-16
+
+### add
+- add fuction of bonferrono_test in friedman_test.
+
+## [0.0.15] - 2024-04-13
+
+### add
+- add fuction of bonferrono_test in kruskalwallis_test.
+
 ## [0.0.14] - 2024-04-08
 
 ### add

data/ext/num4anova/OneWayLayout.java

@@ -30,6 +30,11 @@ import org.apache.commons.math3.stat.descriptive.DescriptiveStatistics;
 import org.apache.commons.math3.distribution.ChiSquaredDistribution;
 import org.apache.commons.math3.distribution.FDistribution;
 import java.util.Map;
+
+import org.apache.commons.math3.stat.ranking.NaNStrategy;
+import org.apache.commons.math3.stat.ranking.NaturalRanking;
+import org.apache.commons.math3.stat.ranking.TiesStrategy;
+import java.util.stream.IntStream;
 public class OneWayLayout {
     private static OneWayLayout oneWay = new OneWayLayout();
     public static OneWayLayout getInstance() {
@@ -72,6 +77,12 @@ public class OneWayLayout {
         double statistic = oneway.calcTestStatistic(xi);
         return oneway.execute_test(statistic, a);
     }
+    public boolean kruskalWallisTest(double[][] xi, double a) {
+        OneWayAnovaTest oneway = new KruskalWallisTest();
+
+        double statistic = oneway.calcTestStatistic(xi);
+        return oneway.execute_test(statistic, a);
+    }
     /*********************************/
     /* interface define              */
     /*********************************/
@@ -380,5 +391,66 @@ public class OneWayLayout {
             return (statistic >= f) ? true : false;
         }
     }
-}
+    // クラスカル・ウォリス検定
+    private class KruskalWallisTest implements OneWayAnovaTest {
+        private NaturalRanking naturalRanking;
+        private int[] ni = null;
+        public KruskalWallisTest() {
+            naturalRanking = new NaturalRanking(NaNStrategy.FIXED, 
+                                                 TiesStrategy.AVERAGE);
+        }
+        public double calcTestStatistic(double[][] xi) {
+            double[] z = concatSample(xi);            // 全てのデータをつなぐ
+            double[] ranks = naturalRanking.rank(z);  // rankに順位値に変換
+            double[] sumRankXi = calcSumRankXi(ranks);
+            double kw = 0.0;
+            int n = z.length;
+
+            for(int i = 0; i < sumRankXi.length; i++) {
+                kw += sumRankXi[i] * sumRankXi[i] / ni[i];
+            }
+            return 12.0 / (n * (n + 1.0)) * kw - 3.0 * (n + 1.0);
+        }
+        public boolean execute_test(double statistic, double a) {
+            ChiSquaredDistribution chi2Dist = new ChiSquaredDistribution(ni.length - 1);
+            double r_val = chi2Dist.inverseCumulativeProbability(1.0 - a);
+
+            return (r_val < statistic) ? true : false;
+        }
+        private int[] calcNi(double[][] xi) {
+            int[] ni = new int[xi.length];
+
+            for(int i = 0; i < ni.length; i++) {
+                ni[i] = xi[i].length;
+            }
+            return ni;
+        }
+        private double[] concatSample(double[][] xi) {
+            DescriptiveStatistics stat = new DescriptiveStatistics();
+            ni = calcNi(xi);
+            Arrays.stream(ni).forEach(stat::addValue);
+            int n = Double.valueOf(stat.getSum()).intValue();
+
+            double[] z = new double[n];
+            int idx = 0;
+            for(int cnt = 0; cnt < xi.length; cnt++) {
+                System.arraycopy(xi[cnt], 0, z, idx, ni[cnt]);
+                idx += ni[cnt];
+            }
+            return z;
+        }
+        private double[] calcSumRankXi(double[] ranks) {
+            double[] sumRi = new double[ni.length];
+            int idx = 0;
+
+            for(int cnt = 0; cnt < ni.length; cnt++) {
+                sumRi[cnt] = IntStream.range(idx, idx + ni[cnt])
+                                    .mapToDouble(i -> ranks[i])
+                                    .sum();
+                idx += ni[cnt];
+            }
+            return sumRi;
+        }
+    }
+ }
 

data/ext/num4anova/TwoWayLayout.java

@@ -1,6 +1,11 @@
 import org.apache.commons.math3.stat.descriptive.DescriptiveStatistics;
 import org.apache.commons.math3.distribution.FDistribution;
+import org.apache.commons.math3.distribution.ChiSquaredDistribution;
 import java.util.Arrays;
+
+import org.apache.commons.math3.stat.ranking.NaNStrategy;
+import org.apache.commons.math3.stat.ranking.NaturalRanking;
+import org.apache.commons.math3.stat.ranking.TiesStrategy;
 public class TwoWayLayout {
     private static TwoWayLayout twoWay = new TwoWayLayout();
     public static TwoWayLayout getInstance() {
@@ -18,16 +23,26 @@ public class TwoWayLayout {
         double[] statistic = twoway.calcTestStatistic(xij);
         return twoway.execute_test(statistic, a);
     }
+    public boolean friedmanTest(double[][] xij, double a) {
+        OneWayAnovaTest twoway = new FriedmanTest();
+
+        double statistic = twoway.calcTestStatistic(xij);
+        return twoway.execute_test(statistic, a);
+    }
     /*********************************/
     /* interface define              */
     /*********************************/
     private interface TwoWayAnovaTest {
         double[] calcTestStatistic(double[][][] xij);
-        boolean[] execute_test(double statistic[], double a);
+        boolean[] execute_test(double[] statistic, double a);
     }
     private interface TwoWay2AnovaTest {
         double[] calcTestStatistic(double[][] xij);
-        boolean[] execute_test(double statistic[], double a);
+        boolean[] execute_test(double[] statistic, double a);
+    }
+    private interface OneWayAnovaTest {
+        double calcTestStatistic(double[][] xi);
+        boolean execute_test(double statistic, double a);
     }
     /*********************************/
     /* class define                  */
@@ -162,7 +177,7 @@ public class TwoWayLayout {
             }
             return sumDrift;
         }
-        public boolean[] execute_test(double statistic[], double a) {
+        public boolean[] execute_test(double[] statistic, double a) {
             boolean[] ret = new boolean[3];
 
             ret[0] = evaluation(new FDistribution(an, en), statistic[0], a);
@@ -269,7 +284,7 @@ public class TwoWayLayout {
             return sumDrift;
         }
 
-        public boolean[] execute_test(double statistic[], double a) {
+        public boolean[] execute_test(double[] statistic, double a) {
             boolean[] ret = new boolean[2];
 
             ret[0] = evaluation(new FDistribution(an, en), statistic[0], a);
@@ -282,5 +297,64 @@ public class TwoWayLayout {
             return (statistic >= r_val) ? true : false;
         }
     }
+    // フリードマンの検定
+    private class FriedmanTest implements  OneWayAnovaTest{
+        private NaturalRanking naturalRanking;
+        private int n = 0;
+        private int k = 0;
+        public FriedmanTest()  {
+            naturalRanking = new NaturalRanking(NaNStrategy.FIXED, 
+                                                 TiesStrategy.AVERAGE);
+        }
+        public double calcTestStatistic(double[][] xij) {
+            n = xij[0].length;
+            k = xij.length;
+            double[][] z = concatSample(xij);
+            double[][] ranks = calcRankij(z);
+            double[] sumRankXi = calcSumRankXi(ranks);
+            double kw = 0.0;
+
+            for(int i = 0; i < sumRankXi.length; i++) {
+                kw += sumRankXi[i] * sumRankXi[i];
+            }
+            return 12.0 / (n * k * (k + 1.0)) * kw - 3.0 * n * (k + 1.0);
+        }
+        public boolean execute_test(double statistic, double a) {
+            ChiSquaredDistribution chi2Dist = new ChiSquaredDistribution(k - 1);
+
+            double r_val = chi2Dist.inverseCumulativeProbability(1.0 - a);
+
+            return (r_val < statistic) ? true : false;
+        }
+        private double[][] concatSample(double[][] xi) {
+            double[][] z = new double[n][k];
+
+            for(int i = 0; i < k; i++) {
+                for(int j = 0; j < n; j++) {
+                    z[j][i] = xi[i][j];
+                }
+            }
+            return z;
+        }
+        private double[][] calcRankij(double[][] z) {
+            double[][] ranks = new double[n][k];
+
+            for(int i = 0; i < n; i++) {
+                ranks[i] = naturalRanking.rank(z[i]);
+            }
+            return ranks;
+        }
+        private double[] calcSumRankXi(double[][] ranks) {
+            double[] sumRi = new double[k];
+
+            for(int i = 0; i < k; i++) {
+                sumRi[i] = 0.0;
+                for(int j = 0; j < n; j++) {
+                    sumRi[i] += ranks[j][i];
+                }
+            }
+            return sumRi;
+        }
+    }
 }
 

data/lib/num4anova.rb

@@ -151,6 +151,24 @@ module Num4AnovaLib
         def replicate_test(xi, a)
             return @oneWay.replicateTest(xi.to_java(Java::double[]), a)
         end
+        # クラスカル・ウォリスの検定
+        #
+        # @overload kruskalwallis_test(xi, a)
+        #   @param [array]  xi データ(double[][])
+        #   @param [double] a         有意水準
+        #   @return [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],
+        #   ]
+        #   oneWay = Num4AnovaLib::OneWayLayoutLib.new 
+        #   oneWay.kruskalwallis_test(xi, 0.05)
+        #   => true
+        def kruskalwallis_test(xi, a)
+            return @oneWay.kruskalWallisTest(xi.to_java(Java::double[]), a)
+        end
     end
 
     # 二元配置の分散分析
@@ -217,6 +235,27 @@ module Num4AnovaLib
         def twoway2_anova(xij, a)
             ret = @twoWay.twoway2Anova(xij.to_java(Java::double[]), a)
             return ret.to_a
+        end
+        # フリードマン検定
+        #
+        # @overload friedman_test(xij, a)
+        #   @param [array]  xij データ(double[][])
+        #   @param [double] a         有意水準
+        #   @return [boolean] 検定結果(boolean true:棄却域内 false:棄却域外)
+        # @example
+        #   xij = [
+        #       [13.6, 15.6, 9.2],
+        #       [22.3, 23.3, 13.3],
+        #       [26.7, 28.8, 15.0],
+        #       [28.0, 31.2, 15.8],
+        #   ]
+        #   twoWay = Num4AnovaLib::TwoWayLayoutLib.new 
+        #   twoWay.friedman_test(xij, 0.05)
+        #   =>
+        #     true
+        def friedman_test(xij, a)
+            ret = @twoWay.friedmanTest(xij.to_java(Java::double[]), a)
+            return ret
         end        
     end
     # 共分散分析

metadata

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