num4anova 0.0.5-java → 0.0.7-java

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 7a089d52ead2c9726d9ca61db475f128925afddca7cbcf6f65f0da75e41f0b1b
-  data.tar.gz: 050e29e04bd0e8272da89d5aba231fe08294f50f090d261da7aa6822f5df5bc5
+  metadata.gz: 72ca9e237cf35e5d7c29cef5603044824f821b53974275cf2025db717cc74947
+  data.tar.gz: 8898ddd64a38258ac4208906fed49b9b15f5f0328150b570b91a7c48af8ff96b
 SHA512:
-  metadata.gz: f0cc84d4cebb416b04ed0360ef4cf72799aaebcf70104b6e2d2dc2bc471400d3ac313b71f5885fd0854ae8b3ae22dd66c753066233b1ed2848aabf586cee2c05
-  data.tar.gz: 88a74fb295e7ce0e0bdba405a14239a37fac4f3115bd29f358ce7b54fcc3b71a1533ff095e5aea7e2b141d867051eca703d179874f1ed349aef92575d8d0a186
+  metadata.gz: f0d19ec6277e2551d90268319d6866ea15fa467ed1da1ad8ad776db3e4c57d4b2ab4ae3f57abaf3fe976ee5adc68a0b1a2e615a1b7c88b5c841252bdd1bae292
+  data.tar.gz: 759dab8f82a35ca0a65a1d81a45b670fc5476799df44517ff52f379fb83f29a661fd75be82dc79fa0835b3f9edb8bb06a4746869a12eee66e06d77c68c504517

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.7
 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