num4anova 0.0.5-java → 0.0.6-java

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 7a089d52ead2c9726d9ca61db475f128925afddca7cbcf6f65f0da75e41f0b1b
-  data.tar.gz: 050e29e04bd0e8272da89d5aba231fe08294f50f090d261da7aa6822f5df5bc5
+  metadata.gz: 4dd60df363d899af7ca69bca4cb1da1689548f083abaacfd6441d5a589bf3501
+  data.tar.gz: 84469385e378cde1eddc3521582aefd2f5b2adce84e76d0c8a4c23bdd40d9374
 SHA512:
-  metadata.gz: f0cc84d4cebb416b04ed0360ef4cf72799aaebcf70104b6e2d2dc2bc471400d3ac313b71f5885fd0854ae8b3ae22dd66c753066233b1ed2848aabf586cee2c05
-  data.tar.gz: 88a74fb295e7ce0e0bdba405a14239a37fac4f3115bd29f358ce7b54fcc3b71a1533ff095e5aea7e2b141d867051eca703d179874f1ed349aef92575d8d0a186
+  metadata.gz: 7a036fe08363781f31b3583d7c7a454131680413e9c3aaeb0ebe135eec930912807e3f84be688d8bf70e1e56eafdefef0a3425f07eb18d6d042549dc2f7dd893
+  data.tar.gz: 7d7c136e1e457c08ac341fc542be1590888c233c161034eeff0f67d3d1466b4879ccd49c5f2a83e2f42f18f28a0279d66909d768d616afabf24b59be36a7e4e4

data/CHANGELOG.md

@@ -2,6 +2,11 @@
 
 ## Unreleased
 
+## [0.0.6] - 2024-02-04
+
+### add
+- add function of twoway2_anova.
+
 ## [0.0.5] - 2024-02-01
 
 ### add

data/ext/num4anova/TwoWayLayout.java

@@ -12,6 +12,12 @@ public class TwoWayLayout {
         double[] statistic = twoway.calcTestStatistic(xij);
         return twoway.execute_test(statistic, a);
     }
+    public boolean[] twoway2Anova(double[][] xij, double a) {
+        TwoWay2AnovaTest twoway = new TwoWay2Anova();
+
+        double[] statistic = twoway.calcTestStatistic(xij);
+        return twoway.execute_test(statistic, a);
+    }
     /*********************************/
     /* interface define              */
     /*********************************/
@@ -19,9 +25,14 @@ public class TwoWayLayout {
         double[] calcTestStatistic(double[][][] xij);
         boolean[] execute_test(double statistic[], double a);
     }
+    private interface TwoWay2AnovaTest {
+        double[] calcTestStatistic(double[][] xij);
+        boolean[] execute_test(double statistic[], double a);
+    }
     /*********************************/
     /* class define                  */
     /*********************************/
+    // 二元配置の分散分析(繰り返し数が等しい時)
     private class TwoWayAnova implements TwoWayAnovaTest {
         private int a = 0;
         private int b = 0;
@@ -32,6 +43,7 @@ public class TwoWayLayout {
         private int en = 0;
         public double[] calcTestStatistic(double[][][] xij) {
             double statistic[] = new double[3];
+
             a = xij.length;
             b = xij[0].length;
             n = xij[0][0].length;
@@ -45,7 +57,6 @@ public class TwoWayLayout {
             double[] meanBn = calcMeanBn(meanXij);
             double meanABn = calcMeanABn(meanAn);
 
-            double allDrift = calcAllDrift(xij, meanABn); // 全変動
             double anDrift  = calcAnDrift(meanAn, meanABn); // 水準Ai間変動
             double bnDrift  = calcBnDrift(meanBn, meanABn); // 水準Bj間変動
                                                             // 交互作用の変動
@@ -76,20 +87,17 @@ public class TwoWayLayout {
         }
         private double[] calcMeanAn(double[][] meanXij) {
             double[] an = new double[a];
-            DescriptiveStatistics stat = new DescriptiveStatistics();
 
             for(int i = 0; i < a; i++) {
                 double sumSa = 0.0;
                 for(int j = 0; j < b; j++) {
-                    sumSa += meanXij[i][j];
+                    an[i] += meanXij[i][j] / b;
                 }
-                an[i] = sumSa / b;
             }
             return an;
         }
         private double[] calcMeanBn(double[][] meanXij) {
             double[] bn = new double[b];
-            double[] sumA = new double[b];
 
             for(int i = 0; i < a; i++) {
                 for(int j = 0; j < b; j++) {
@@ -104,20 +112,6 @@ public class TwoWayLayout {
             Arrays.stream(meanAn).forEach(stat::addValue);
             return stat.getMean();
         }
-        // 全変動
-        private double calcAllDrift(double[][][] xij, double meanABn) {
-            double sumDrift = 0.0;
-
-            for(int i = 0; i < a; i++) {
-                for(int j = 0; j < b; j++) {
-                    for(int k = 0; k < xij[i][j].length; k++) {
-                        double diffXijk = xij[i][j][k] - meanABn;
-                        sumDrift += diffXijk * diffXijk; 
-                    }
-                }
-            }
-            return sumDrift;
-        }
         // 水準Ai間変動
         private double calcAnDrift(double[] meanAn, double meanABn) {
             double sumDrift = 0.0;
@@ -179,7 +173,113 @@ public class TwoWayLayout {
         private boolean evaluation(FDistribution fDist, double statistic, double a) {
             double r_val = fDist.inverseCumulativeProbability(1.0 - a);
 
-            return (statistic < r_val) ? false : true;
+            return (statistic >= r_val) ? true : false;
+        }
+    }
+    // 二元配置の分散分析(繰り返しのない時)
+    private class TwoWay2Anova implements TwoWay2AnovaTest {
+        private int a = 0;
+        private int b = 0;
+        private int an = 0;
+        private int bn = 0;
+        private int en = 0;
+        public double[] calcTestStatistic(double[][] xij) {
+            double statistic[] = new double[2];
+
+            a = xij.length;
+            b = xij[0].length;
+            an = a- 1;
+            bn = b - 1;
+            en = (a- 1) * (b - 1);
+
+            double[] meanAn = calcMeanAn(xij);
+            double[] meanBn = calcMeanBn(xij);
+            double meanAB = calcMeanAB(meanAn);
+
+            double anDrift  = calcAnDrift(meanAn, meanAB); // 水準Ai間変動
+            double bnDrift  = calcBnDrift(meanBn, meanAB); // 水準Bj間変動
+            double benchDrift  = calcBenchDrift(xij, meanAn, meanBn, meanAB);   // 水準内変動
+            double va = anDrift / an;
+            double vb = bnDrift / bn;
+            double ve = benchDrift / en;
+
+            statistic[0] = va / ve;
+            statistic[1] = vb / ve;
+            return statistic;
+        }
+        private double[] calcMeanAn(double[][] xij) {
+            double[] an = new double[a];
+
+            for(int i = 0; i < a; i++) {
+                for(int j = 0; j < b; j++) {
+                    an[i] += xij[i][j] / b;
+                }
+            }
+            return an;
+        }
+        private double[] calcMeanBn(double[][] xij) {
+            double[] bn = new double[b];
+
+            for(int i = 0; i < a; i++) {
+                for(int j = 0; j < b; j++) {
+                    bn[j] += xij[i][j] / a;
+                }
+            }
+            return bn;
+        }
+        private double calcMeanAB(double[] meanAn) {
+            DescriptiveStatistics stat = new DescriptiveStatistics();
+
+            Arrays.stream(meanAn).forEach(stat::addValue);
+            return stat.getMean();
+        }
+        // 水準Ai間変動
+        private double calcAnDrift(double[] meanAn, double meanAB) {
+            double sumDrift = 0.0;
+
+            for(int i =0; i < meanAn.length; i++) {
+                double diffXi = meanAn[i] - meanAB;
+
+                sumDrift += diffXi * diffXi;
+            }
+            return b * sumDrift;      
+        }
+        // 水準Bj間変動
+        private double calcBnDrift(double[] meanBn, double meanAB) {
+            double sumDrift = 0.0;
+
+            for(int j = 0; j < meanBn.length; j++) {
+                double diffXj = meanBn[j] - meanAB;
+
+                sumDrift += diffXj * diffXj;
+            }
+            return a * sumDrift;
+        }
+        // 水準内変動
+        private double calcBenchDrift(double[][] xij, double[] meanAn, double[] meanBn, double meanAB) {
+            double sumDrift = 0.0;
+
+            for(int i = 0; i < a; i++) {
+                for(int j = 0; j < b; j++) {
+                    double diffXj = xij[i][j] - meanAn[i] - meanBn[j] + meanAB;
+
+                    sumDrift += diffXj * diffXj;
+                }
+            }
+            return sumDrift;
+        }
+
+        public boolean[] execute_test(double statistic[], double a) {
+            boolean[] ret = new boolean[2];
+
+            ret[0] = evaluation(new FDistribution(an, en), statistic[0], a);
+            ret[1] = evaluation(new FDistribution(bn, en), statistic[1], a);
+            return ret;
+        }
+        private boolean evaluation(FDistribution fDist, double statistic, double a) {
+            double r_val = fDist.inverseCumulativeProbability(1.0 - a);
+
+            return (statistic >= r_val) ? true : false;
         }
     }
 }

data/lib/num4anova.rb

@@ -158,6 +158,7 @@ module Num4AnovaLib
             @twoWay = TwoWayLayout.getInstance()
         end
         # 二元配置の分散分析
+        #  (繰り返し数が等しい時)
         #
         # @overload twoway_anova(xij, a)
         #   @param [array]  xij データ(double[][][])
@@ -194,6 +195,28 @@ module Num4AnovaLib
             ret = @twoWay.twowayAnova(xij.to_java(Java::double[][]), a)
             return ret.to_a
         end
+        # 二元配置の分散分析
+        #  (繰り返しのない時)
+        #
+        # @overload twoway2_anova(xij, a)
+        #   @param [array]  xij データ(double[][])
+        #   @param [double] a         有意水準
+        #   @return [Array] 検定結果(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.twoway2_anova(xij, 0.05)
+        #   =>
+        #     [true, true]
+        def twoway2_anova(xij, a)
+            ret = @twoWay.twoway2Anova(xij.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.5
+  version: 0.0.6
 platform: java
 authors:
 - siranovel
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2024-02-01 00:00:00.000000000 Z
+date: 2024-02-04 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: rake