num4tststatistic 0.0.1-java → 0.0.2-java

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: a8b0c9145ecc7379e5619ad7f8ffb33b61d5560bf92c594e42f53e5c0e288a1b
-  data.tar.gz: e8b836c2974645183f78e3e9c6533b105b3d364f2405c10c68ab18372f06e685
+  metadata.gz: 7cf4040a8f70a0ee55f1bcbb70397220d88aaab31d3edaf96dc2b8d4e9589bea
+  data.tar.gz: b1eb211d1b876bbd7a801b4a9eb8398eca8ba1d3dd1b4eee72783d73b0b691b2
 SHA512:
-  metadata.gz: 20fa87418cb04465c0c62d33dd64016e67a57ad9f0c0dccddf2572b8f6aa6618ddbc76ff6111c12baa414103e24cc8dc74abcec87dfac38768e28a3ae80233f5
-  data.tar.gz: a01a432d8a7530b3be8f568e047faed5a42d9d79dde68f60426bcda8a751c166212e772495bd74e7b4277d3e7045e9443a5993c19a5af36b451306486dce528f
+  metadata.gz: b2802b5d55555c08234609c5d1cf29a5e6ce223aab56cf2b3a66b34ee3cd1078fd94ecc4ef05a069709d56c03ac51abcad2e05aa9a8fe0c88fa34e57a1428c5b
+  data.tar.gz: fa23a4ebe88f4dd924651a6976956948bf3ad9494dbdd7ea77a26123b5d45bc86dd40e7b88bffe4c31e231cfe2a02ebd6726f78b140277d061a517361c8f54b5

data/CHANGELOG.md

@@ -2,6 +2,11 @@
 
 ## Unreleased
 
+## [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,48 @@ 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;
+    }
 }
 
 

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,39 @@ 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
     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.2
 platform: java
 authors:
 - siranovel
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2023-11-11 00:00:00.000000000 Z
+date: 2023-11-15 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: rake