num4tststatistic 0.0.2-java → 0.0.4-java

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 7cf4040a8f70a0ee55f1bcbb70397220d88aaab31d3edaf96dc2b8d4e9589bea
-  data.tar.gz: b1eb211d1b876bbd7a801b4a9eb8398eca8ba1d3dd1b4eee72783d73b0b691b2
+  metadata.gz: 8b8a5b476ae860590efbcc631f0851e19d12b37f7f43e127d0637d7754912d4c
+  data.tar.gz: 2cd3007bd6872e547a7bc594d4d9d9af1b4b1ebcd8a45c12c69c5568752ff923
 SHA512:
-  metadata.gz: b2802b5d55555c08234609c5d1cf29a5e6ce223aab56cf2b3a66b34ee3cd1078fd94ecc4ef05a069709d56c03ac51abcad2e05aa9a8fe0c88fa34e57a1428c5b
-  data.tar.gz: fa23a4ebe88f4dd924651a6976956948bf3ad9494dbdd7ea77a26123b5d45bc86dd40e7b88bffe4c31e231cfe2a02ebd6726f78b140277d061a517361c8f54b5
+  metadata.gz: 1cf442842fe6bd7fd8447c12910c214aabdbe63aa0ab3811e984e289e5fd14276d99391672898c14c4a61368d403165f7a5e9d3c3cfccc9b640f095cb1f2f209
+  data.tar.gz: 7af0e381b6833ae5215d1048e7bdcb4fc8080d3fd0013eb3c7aabe868b5e44ebbb9ba5d89e5030a2428dd39d62e3f7d44f35772a162c4a6f3543c95aa9a6d25f

data/CHANGELOG.md

@@ -2,6 +2,16 @@
 
 ## Unreleased
 
+## [0.0.4] - 2023-12-18
+
+### add
+- add function of wilcoxon.
+
+## [0.0.3] - 2023-11-16
+
+### add
+- add function of grubbs.
+
 ## [0.0.2] - 2023-11-15
 
 ### add

data/ext/num4tststatistic/TstStatistic.java

@@ -2,7 +2,7 @@ import java.util.Arrays;
 
 import org.apache.commons.math3.stat.descriptive.SummaryStatistics;
 import org.apache.commons.math3.stat.correlation.PearsonsCorrelation;
-
+import org.apache.commons.math3.stat.inference.WilcoxonSignedRankTest;
 public class TstStatistic {
     public static double populationMean(double[] xi, double m0) {
         int n = xi.length;
@@ -172,6 +172,25 @@ public class TstStatistic {
 
         return t;
     }
+    public static double grubbs(double[] xi, double xk) {
+        SummaryStatistics stat = new SummaryStatistics();
+
+        Arrays.stream(xi).forEach(stat::addValue);
+        double m = stat.getMean();     // 平均
+        double sd = stat.getStandardDeviation();// 標準偏差
+        double min = stat.getMin();
+        double max = stat.getMax();
+        double t = 0.0;
+
+        if (xk == min) { t = (m - xk) / sd; }
+        if (xk == max) { t = (xk - m) / sd; }
+        return t;
+    }
+    public static double wilcoxon(double[] x, double[] y) {
+        WilcoxonSignedRankTest wilcox = new WilcoxonSignedRankTest();
+
+        return wilcox.wilcoxonSignedRank(x, y); 
+    }
 }
 
 

data/lib/num4tststatistic.rb

@@ -214,7 +214,7 @@ module Num4TstStatisticLib
         # 独立性の検定量
         #
         # @overload independency(fij)
-        #   @param [Array] fij 実測度数(double[])
+        #   @param [Array] fij 実測度数(double[][])
         #   @return [double] 検定統計量
         # @example
         #   fij = [
@@ -228,6 +228,37 @@ module Num4TstStatisticLib
         def independency(fij)
             return TstStatistic.independency(fij.to_java(Java::double[]))
         end
+        # グラプス・スミルノフの外れ値の検定量
+        #
+        # @overload grubbs(xi, xk)
+        #   @param [Array] xi xiのデータ(double[])
+        #   @param [double] xk 外れ値
+        #   @return [double] 検定統計量
+        # @example
+        #   xi = [3.4, 3.5, 3.3, 2.2, 3.3, 3.4, 3.6, 3.2]
+        #   Num4TstStatisticLib.grubbs(xi, 2.2)
+        #   => 2.3724
+        # @note
+        #   グラプス・スミルノフの数表に従う
+        def grubbs(xi, xk)
+            return TstStatistic.grubbs(xi.to_java(Java::double), xk)
+        end
+        # ウィルコクス符号付き順位検定
+        #
+        # @overload wilcoxon(x, y)
+        #   @param [Array] x xのデータ(double[])
+        #   @param [Array] y yのデータ(double[])
+        #   @return [double] 検定統計量
+        # @example
+        #   x = [28, 25, 29, 28, 30, 20, 31, 27, 24, 26, 35, 23, 27, 32]
+        #   y = [32, 30, 31, 27, 35, 25, 40, 30, 45, 28, 32, 30, 30, 38]
+        #   Num4TstStatisticLib.wilcoxon(x, y)
+        #   => 99
+        # @note
+        #   標準正規分布 N(0,1*1)に従う
+        def wilcoxon(x, y)
+            return TstStatistic.wilcoxon(x.to_java(Java::double), y.to_java(Java::double))
+        end
     end
 end
 

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: num4tststatistic
 version: !ruby/object:Gem::Version
-  version: 0.0.2
+  version: 0.0.4
 platform: java
 authors:
 - siranovel
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2023-11-15 00:00:00.000000000 Z
+date: 2023-12-18 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: rake