num4normality 0.0.9-java → 0.0.10-java

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 28c1bcf4f9952a2f18587d5b23d301a43ea31cd109d6a2d99ce665f900877bbf
-  data.tar.gz: 29a902d3ef52739cb913bc523c5488ebcc14ac3288d8babc8d7e64fd28448cbf
+  metadata.gz: bc5dd36e009ab6ce0b3a044652184f723ce70b416bab8780d01a56952bcd7b8b
+  data.tar.gz: 46fa3aefe151fb4e6dd8a47a3c4a4e38670704dd74b21a9d6fb19f6a13349a82
 SHA512:
-  metadata.gz: 66fe7346cae038b6a399395d46863b2a7aaecff96687b451075e2affa08236c1f8eb5b18066eed22eafac9f1ac7c719989e43c07082fceb3d17b9b6d4310fe4b
-  data.tar.gz: 3b621221e49cfaabf960cb5f8d413aadb53ec1e280e7f3c65b1cb2d70dd8c960129e239aee017365d9910c609b84b35d875f7e12fcc5a3f0f84284db3a9d78b6
+  metadata.gz: f9933cc4c22ac645b17ae8e65b835d72fc80453e05c478903d1c748d3ff66615b9a0c631ced2ef4a2c8051297b2da16cdcd73b47d654032ed79599b0ae19f2bf
+  data.tar.gz: c74c186cc3caa4d511f8187ec9b7024d85ae8249d40692942800655c45fda60589be35745cb94659fd7d40038e65d4a0e3cb4ae91ef837d06302a6e180736e6d

data/CHANGELOG.md

@@ -2,6 +2,10 @@
 
 ## Unreleased
 
+## [0.0.10] - 2024-05-18
+### add
+- add function of jbtest.
+
 ## [0.0.9] - 2024-05-11
 ### add
 - add function of adtest.

data/ext/num4normality/Normality.java

@@ -94,6 +94,12 @@ public class Normality {
 
         return daigo.test(x, HypothesisTest.P);
     }
+    public static boolean jbTest(double[] xi) {
+        HypothesisTest daigo = new JBTest();
+        double x = daigo.calcTestStatistic(xi);
+
+        return daigo.test(x, HypothesisTest.P);
+    }
     /*********************************/
     /* interface define              */
     /*********************************/
@@ -644,5 +650,27 @@ public class Normality {
             return a2 * (1.0 + 0.75 / n + 2.25 / (n * n));
         }
     }
+    // Jarque-Bera検定
+    private static class JBTest implements HypothesisTest {
+        private HypothesisTest skewness = null;
+        private HypothesisTest kurtosis = null;
+        public JBTest() {
+            skewness = new SkewnessTest();
+            kurtosis = new KurtosisTest();
+        }
+        public double calcTestStatistic(double[] xi) {
+            int n = xi.length;
+            double x1 = skewness.calcTestStatistic(xi);
+            double x2 = kurtosis.calcTestStatistic(xi);
+
+            return n / 6.0 * (x1 *x1  + x2 * x2 / 4.0);
+        }
+        public boolean test(double statistic, double a) {
+            ChiSquaredDistribution chi2Dist = new ChiSquaredDistribution(2);
+            double p = chi2Dist.cumulativeProbability(statistic);
+
+            return (p < a) ? true  : false;
+        }
+    }
 }
 

data/lib/num4normality.rb

@@ -145,6 +145,18 @@ module Num4NormalityLib
         def adtest(xi)
              Normality.adTest(xi.to_java(Java::double))
         end
+        # Jarque-Bera検定
+        #
+        # @overload jbtest(xi)
+        #   @param  [Array]  xi データ(double[])
+        #   @return [boolean] 検定結果(true:棄却域内 false:棄却域外)
+        # @example
+        #   xi = [320, 240, 402, 325, 440, 286, 362, 281, 560, 212, 198, 209, 374]
+        #   Num4NormalityLib.jbtest(xi)
+        #   => false
+        def jbtest(xi)
+             Normality.jbTest(xi.to_java(Java::double))
+        end
     end
 end
 

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: num4normality
 version: !ruby/object:Gem::Version
-  version: 0.0.9
+  version: 0.0.10
 platform: java
 authors:
 - siranovel
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2024-05-11 00:00:00.000000000 Z
+date: 2024-05-18 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: rake