num4tststatistic 0.0.1-java → 0.0.3-java

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: a8b0c9145ecc7379e5619ad7f8ffb33b61d5560bf92c594e42f53e5c0e288a1b
-  data.tar.gz: e8b836c2974645183f78e3e9c6533b105b3d364f2405c10c68ab18372f06e685
+  metadata.gz: 8c015fe7102176f7c7c2314cc1dfb6622c07bc8fe1e5b36c3417c0c1df3e255a
+  data.tar.gz: f348504d4e1245853d4cc93a19e4dbe16fa1c4ac9ef6be75f53d8f4b7d1df767
 SHA512:
-  metadata.gz: 20fa87418cb04465c0c62d33dd64016e67a57ad9f0c0dccddf2572b8f6aa6618ddbc76ff6111c12baa414103e24cc8dc74abcec87dfac38768e28a3ae80233f5
-  data.tar.gz: a01a432d8a7530b3be8f568e047faed5a42d9d79dde68f60426bcda8a751c166212e772495bd74e7b4277d3e7045e9443a5993c19a5af36b451306486dce528f
+  metadata.gz: a93a631e82744d7d8130e32951caf07a82d6b0530d74f2fb54fda082dfb8269f8ccdcad69a99e3a9b81e14399077b10ec42b8b6f6b2e535704da87bf993c4a5b
+  data.tar.gz: 19cc9ecf0e5bcc36ec32bfed81ffa8b3301588617267edac4e1be18c5e4ddae8f8f8d693a7b67503c9c8c79c227c7104f324ef56a45486ee6d1fc7f812084a57

data/CHANGELOG.md

@@ -2,6 +2,16 @@
 
 ## Unreleased
 
+## [0.0.3] - 2023-11-16
+
+### add
+- add function of grubbs.
+
+## [0.0.2] - 2023-11-15
+
+### add
+- add function of fidelity and independency.
+
 ## [0.0.1] - 2023-11-11
 
 ### Fixed

data/ext/num4tststatistic/TstStatistic.java

@@ -130,6 +130,62 @@ public class TstStatistic {
              - 0.5 * Math.log((1.0 + rth0) / (1.0 - rth0))
                );
     }
+    public static double fidelity(double[] fi, double[] pi) {
+        double[] e = new double[fi.length];
+        double t = 0.0;
+        SummaryStatistics stat = new SummaryStatistics();
+
+        Arrays.stream(fi).forEach(stat::addValue);
+        double s = stat.getSum();
+
+        for (int i = 0; i < fi.length; i++) {
+            e[i] = s * pi[i];
+            double f_e = fi[i] - e[i];
+            t += f_e * f_e / e[i];
+        }
+        return t;
+    }
+    public static double independency(double[][] fij) {
+        double t = 0.0;
+        long n = 0;
+        double[] fa = new double[fij.length];
+        double[] fb = new double[fij[0].length];
+
+        // 各セルの計算
+        for (int i = 0; i < fij.length; i++) {
+            fa[i] = 0.0;
+            for (int j = 0; j < fij[i].length; j++)  {
+                fa[i] += fij[i][j];
+                fb[j] += fij[i][j];
+                n += fij[i][j];
+            }
+        }
+
+        // 検定統計量計算
+        for (int i = 0; i < fij.length; i++) {
+            for (int j = 0; j < fij[i].length; j++)  {
+                double f_e = n * fij[i][j] - fa[i] * fb[j];
+
+                t += f_e * f_e / (n * fa[i] * fb[j]);
+            }
+        }
+
+        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;
+    }
 }
 
 

data/lib/num4tststatistic.rb

@@ -58,7 +58,7 @@ module Num4TstStatisticLib
         #
         # @overload diffPopulationMean2EquVar(xi1, xi2)
         #   @param [Array] xi1 x1のデータ(double[])
-        #   @param [Array] xi2 x2のデータ)double[])
+        #   @param [Array] xi2 x2のデータ(double[])
         #   @return [double] 検定統計量
         # @example
         #   xi1 = [165, 130, 182, 178, 194, 206, 160, 122, 212, 165, 247, 195]
@@ -77,7 +77,7 @@ module Num4TstStatisticLib
         #
         # @overload diffPopulationMean2UnEquVar(xi1, xi2)
         #   @param [Array] xi1 x1のデータ(double[])
-        #   @param [Array] xi2 x2のデータ)double[])
+        #   @param [Array] xi2 x2のデータ(double[])
         #   @return [double] 検定統計量
         # @example
         #   xi1 = [165, 130, 182, 178, 194, 206, 160, 122, 212, 165, 247, 195]
@@ -94,7 +94,7 @@ module Num4TstStatisticLib
         # ウェルチ検定の為の自由度
         # @overload df4welch(xi1, xi2)
         #   @param [Array] xi1 x1のデータ(double[])
-        #   @param [Array] xi2 x2のデータ)double[])
+        #   @param [Array] xi2 x2のデータ(double[])
         #   @return [int] 自由度
         # @example
         #   xi1 = [165, 130, 182, 178, 194, 206, 160, 122, 212, 165, 247, 195]
@@ -110,7 +110,7 @@ module Num4TstStatisticLib
         #
         # @overload diffPopulationMean(xi1, xi2)
         #   @param [Array] xi1 x1のデータ(double[])
-        #   @param [Array] xi2 x2のデータ)double[])
+        #   @param [Array] xi2 x2のデータ(double[])
         #   @return [double] 検定統計量
         # @example
         #   xi1 = [37.1, 36.2, 36.6, 37.4, 36.8, 36.7, 36.9, 37.4, 36.6, 36.7]
@@ -126,9 +126,9 @@ module Num4TstStatisticLib
         end
         # 2つの母分散の差の検定量
         #
-        # @oerload diffPopulationVar(xi1, xi2)
+        # @overload diffPopulationVar(xi1, xi2)
         #   @param [Array] xi1 x1のデータ(double[])
-        #   @param [Array] xi2 x2のデータ)double[])
+        #   @param [Array] xi2 x2のデータ(double[])
         #   @return [double] 検定統計量
         # @example
         #   xi1 = [165, 130, 182, 178, 194, 206, 160, 122, 212, 165, 247, 195]
@@ -162,7 +162,7 @@ module Num4TstStatisticLib
         #
         # @overload unCorrelation(x, y)
         #   @param [Array] x xのデータ(double[])
-        #   @param [Array] y yのデータ)double[])
+        #   @param [Array] y yのデータ(double[])
         #   @return [double] 検定統計量
         # @example
         #   x = [113, 64, 16, 45, 28, 19, 30, 82, 76]
@@ -180,7 +180,7 @@ module Num4TstStatisticLib
         #
         # @overload populationCorre(x, y, rth0)
         #   @param [Array] x xのデータ(double[])
-        #   @param [Array] y yのデータ)double[])
+        #   @param [Array] y yのデータ(double[])
         #   @param [double] rth0 母相関係数
         #   @return [double] 検定統計量
         # @example
@@ -195,6 +195,54 @@ module Num4TstStatisticLib
                 x.to_java(Java::double), y.to_java(Java::double), rth0
             )
         end
+        # 適合度の検定量
+        #
+        # @overload fidelity(fi, pi)
+        #   @param [Array] fi 実測度数(double[])
+        #   @param [Array] pi 比率(double[])
+        #   @return [double] 検定統計量
+        # @example
+        #   fi = [57, 33, 46, 14]
+        #   pi = [0.4, 0.2, 0.3, 0.1]
+        #   Num4TstStatisticLib.fidelity(fi, pi)
+        #   => 0.5389
+        # @note
+        #   自由度(n-1)の階２乗分布に従う
+        def fidelity(fi, pi)
+            return TstStatistic.fidelity(fi.to_java(Java::double), pi.to_java(Java::double))
+        end
+        # 独立性の検定量
+        #
+        # @overload independency(fij)
+        #   @param [Array] fij 実測度数(double[][])
+        #   @return [double] 検定統計量
+        # @example
+        #   fij = [
+        #     [57, 33, 46, 14],
+        #     [89, 24, 75, 12],
+        #   ]
+        #   Num4TstStatisticLib.independency(fij)
+        #   => 8.5711
+        # @note
+        #   自由度(m-1)(n-1)の階２乗分布に従う
+        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
     end
 end
 

metadata

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