num4tststatistic 0.0.1-java

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA256:
+  metadata.gz: a8b0c9145ecc7379e5619ad7f8ffb33b61d5560bf92c594e42f53e5c0e288a1b
+  data.tar.gz: e8b836c2974645183f78e3e9c6533b105b3d364f2405c10c68ab18372f06e685
+SHA512:
+  metadata.gz: 20fa87418cb04465c0c62d33dd64016e67a57ad9f0c0dccddf2572b8f6aa6618ddbc76ff6111c12baa414103e24cc8dc74abcec87dfac38768e28a3ae80233f5
+  data.tar.gz: a01a432d8a7530b3be8f568e047faed5a42d9d79dde68f60426bcda8a751c166212e772495bd74e7b4277d3e7045e9443a5993c19a5af36b451306486dce528f

data/CHANGELOG.md

@@ -0,0 +1,11 @@
+# Change Log
+
+## Unreleased
+
+## [0.0.1] - 2023-11-11
+
+### Fixed
+- fix first fixed.
+
+
+

data/Gemfile

@@ -0,0 +1,3 @@
+source "https://rubygems.org"
+
+gemspec

data/LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2023 siranovel
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.

data/Rakefile

@@ -0,0 +1,7 @@
+require  'rake/javaextensiontask'
+
+Rake::JavaExtensionTask.new(name='num4tststatistic')  do | ext |
+  ext.release          = '11'
+  ext.classpath        = 'lib/commons-math3-3.6.1.jar'
+end
+task :default => [:compile]

data/ext/num4tststatistic/TstStatistic.java

@@ -0,0 +1,135 @@
+import java.util.Arrays;
+
+import org.apache.commons.math3.stat.descriptive.SummaryStatistics;
+import org.apache.commons.math3.stat.correlation.PearsonsCorrelation;
+
+public class TstStatistic {
+    public static double populationMean(double[] xi, double m0) {
+        int n = xi.length;
+        SummaryStatistics stat = new SummaryStatistics();
+
+        Arrays.stream(xi).forEach(stat::addValue);
+        double m = stat.getMean();     // 平均
+        double s2 = stat.getVariance();// 分散
+        return (m - m0) / Math.sqrt(s2 / n);
+    }
+    public static double populationVar(double[] xi, double sig0) {
+        int n = xi.length;
+        SummaryStatistics stat = new SummaryStatistics();
+        Arrays.stream(xi).forEach(stat::addValue);
+        double s2 = stat.getVariance();// 分散
+        
+        return (n - 1) * s2 / sig0;
+    }
+    public static double populationRatio(int m, int n, double p0) {
+        double p = (double)m / (double)n;
+
+        return (p - p0) / Math.sqrt(p0 * (1-p0) / n);
+    }
+    public static double diffPopulationMean2EquVar(double[] xi1, double[] xi2) {
+        int n1 = xi1.length;
+        int n2 = xi2.length;
+        SummaryStatistics stat1 = new SummaryStatistics();
+        SummaryStatistics stat2 = new SummaryStatistics();
+        Arrays.stream(xi1).forEach(stat1::addValue);
+        Arrays.stream(xi2).forEach(stat2::addValue);
+
+        double m1 = stat1.getMean();
+        double m2 = stat2.getMean();
+        double s12 = stat1.getVariance();// 分散
+        double s22 = stat2.getVariance();// 分散
+        double s2 = ((n1 - 1) * s12 + (n2 - 1) * s22) / (n1 + n2 - 2);
+
+        return (m1 - m2) / Math.sqrt((1.0 / n1 + 1.0 / n2) * s2);
+    }
+    public static double diffPopulationMean2UnEquVar(double[] xi1, double[] xi2) {
+        int n1 = xi1.length;
+        int n2 = xi2.length;
+        SummaryStatistics stat1 = new SummaryStatistics();
+        SummaryStatistics stat2 = new SummaryStatistics();
+        Arrays.stream(xi1).forEach(stat1::addValue);
+        Arrays.stream(xi2).forEach(stat2::addValue);
+
+        double m1 = stat1.getMean();
+        double m2 = stat2.getMean();
+        double s12 = stat1.getVariance();// 分散
+        double s22 = stat2.getVariance();// 分散
+
+        return (m1 - m2) / Math.sqrt(s12 / n1 + s22 / n2);
+    }
+    public static int df4welch(double[] xi1, double[] xi2) {
+        int n1 = xi1.length;
+        int n2 = xi2.length;
+        SummaryStatistics stat1 = new SummaryStatistics();
+        SummaryStatistics stat2 = new SummaryStatistics();
+        Arrays.stream(xi1).forEach(stat1::addValue);
+        Arrays.stream(xi2).forEach(stat2::addValue);
+
+        double s12 = stat1.getVariance();// 分散
+        double s22 = stat2.getVariance();// 分散
+        double s14 = s12 * s12;
+        double s24 = s22 * s22;
+        int n12 = n1 * n1;
+        int n22 = n2 * n2;
+        double ns = (s12 / n1) + (s22 / n2);
+
+        return (int)
+        (
+           (ns * ns) / 
+           (
+             s14 / (n12 * (n1 - 1)) + s24 / (n22 * (n2 - 1))
+           )
+         );
+    }
+    public static double diffPopulationMean(double[] xi1, double[] xi2) {
+        int n = xi1.length;
+        SummaryStatistics stat = new SummaryStatistics();
+
+        for(int i = 0; i < n; i++) {
+            stat.addValue(xi1[i] - xi2[i]);
+        }
+        double m = stat.getMean();
+        double s2 = stat.getVariance();// 分散
+        
+        return (m - 0) / Math.sqrt(s2/n);
+    }
+    public static double diffPopulationVar(double[] xi1, double[] xi2) {
+        SummaryStatistics stat1 = new SummaryStatistics();
+        SummaryStatistics stat2 = new SummaryStatistics();
+        Arrays.stream(xi1).forEach(stat1::addValue);
+        Arrays.stream(xi2).forEach(stat2::addValue);
+
+        double s12 = stat1.getVariance();// 分散
+        double s22 = stat2.getVariance();// 分散
+        return s12 / s22;
+    }
+    public static double diffPopulationRatio(int m1, int n1, int m2, int n2) {
+        double p1 = (double)m1 / (double)n1;
+        double p2 = (double)m2 / (double)n2;
+        double p = (double)(m1 + m2) / (double)(n1 + n2);
+
+        return (p1 - p2) / Math.sqrt(p * (1 - p) * (1.0 / n1 + 1.0 / n2));
+    }
+    public static double unCorrelation(double[] x, double[] y) {
+        int n = x.length;
+        PearsonsCorrelation corel = new PearsonsCorrelation();
+
+        double r = corel.correlation(x, y);
+
+        return r * Math.sqrt(n - 2.0) / Math.sqrt(1.0 - r * r);
+    }
+    public static double populationCorre(double[] x, double[] y, double rth0) {
+        int n = x.length;
+        PearsonsCorrelation corel = new PearsonsCorrelation();
+
+        double r = corel.correlation(x, y);
+
+        return Math.sqrt(n-3.0) * 
+               (
+               0.5 * Math.log((1.0 + r) / (1.0 - r)) 
+             - 0.5 * Math.log((1.0 + rth0) / (1.0 - rth0))
+               );
+    }
+}
+
+

data/lib/num4tststatistic.rb

@@ -0,0 +1,200 @@
+require 'java'
+require 'num4tststatistic.jar'
+require 'commons-math3-3.6.1.jar'
+
+java_import 'TstStatistic'
+
+# 検定統計量を計算
+#  (Apache commoms math3使用)
+module Num4TstStatisticLib
+    class << self
+        # 正規母集団の母平均の検定量
+        #
+        # @overload populationMean(xi, m0)
+        #   @param  [Array]  xi データ(double[])
+        #   @param  [double] m0 母平均
+        #   @return [double] 検定統計量
+        # @example
+        #   xi = [15.5, 15.7, 15.4, 15.4, 15.6, 15.4, 15.6, 15.5, 15.4]
+        #   Num4TstStatisticLib.populationMean(xi, 15.4)
+        #   => 2.683
+        # @note
+        #   自由度(N-1)のt分布に従う
+        def populationMean(xi, m0)
+            return TstStatistic.populationMean(xi.to_java(Java::double), m0)
+        end
+        # 正規母集団の母分散の検定量
+        #
+        # @overload populationVar(xi, sig0)
+        #   @param  [Array]  xi データ(double[])
+        #   @param  [double] sig0 母分散
+        #   @return [double] 検定統計量
+        # @example
+        #   xi = xi = [35.2, 34.5, 34.9, 35.2, 34.8, 35.1, 34.9, 35.2, 34.9, 34.8]
+        #   sd = 0.4
+        #   Num4TstStatisticLib.populationVar(xi, sd*sd)
+        #   => 2.906
+        # @note
+        #   自由度(N-1)の階2乗分布に従う
+        def populationVar(xi, sig0)
+            return TstStatistic.populationVar(xi.to_java(Java::double), sig0)
+        end
+        # 母比率の検定量
+        # @overload populationRatio(m, n, p0)
+        #   @param [int]     m m値
+        #   @param [int]     n N値
+        #   @param [double]  p0 母比率
+        #   @return [double] 検定統計量
+        # @example
+        #   Num4TstStatisticLib.populationRatio(29, 346, 0.12)
+        #   => -2.071
+        # @note
+        #   標準正規分布 N(0,1*1)に従う
+        def populationRatio(m, n, p0)
+            return TstStatistic.populationRatio(m, n, p0)
+        end
+        # 2つの母平均の差の検定量
+        # (等分散性を仮定)
+        #
+        # @overload diffPopulationMean2EquVar(xi1, xi2)
+        #   @param [Array] xi1 x1のデータ(double[])
+        #   @param [Array] xi2 x2のデータ)double[])
+        #   @return [double] 検定統計量
+        # @example
+        #   xi1 = [165, 130, 182, 178, 194, 206, 160, 122, 212, 165, 247, 195]
+        #   xi2 = [180, 180, 235, 270, 240, 285, 164, 152]
+        #   Num4TstStatisticLib.diffPopulationMean2EquVar(xi1, xi2)
+        #   => -1.765
+        # @note
+        #   N1+N2-2のt分布に従う
+        def diffPopulationMean2EquVar(xi1, xi2)
+            return TstStatistic.diffPopulationMean2EquVar(
+                xi1.to_java(Java::double), xi2.to_java(Java::double)
+            )
+        end
+        # 2つの母平均の差の検定量
+        # (不等分散性を仮定)
+        #
+        # @overload diffPopulationMean2UnEquVar(xi1, xi2)
+        #   @param [Array] xi1 x1のデータ(double[])
+        #   @param [Array] xi2 x2のデータ)double[])
+        #   @return [double] 検定統計量
+        # @example
+        #   xi1 = [165, 130, 182, 178, 194, 206, 160, 122, 212, 165, 247, 195]
+        #   xi2 = [180, 180, 235, 270, 240, 285, 164, 152]
+        #   Num4TstStatisticLib.diffPopulationMean2UnEquVar(xi1, xi2)
+        #   => -1.636
+        # @note
+        #   df4welch関数で求めた自由度のt分布に従う
+        def diffPopulationMean2UnEquVar(xi1, xi2)
+            return TstStatistic.diffPopulationMean2UnEquVar(
+                xi1.to_java(Java::double), xi2.to_java(Java::double)
+            )
+        end
+        # ウェルチ検定の為の自由度
+        # @overload df4welch(xi1, xi2)
+        #   @param [Array] xi1 x1のデータ(double[])
+        #   @param [Array] xi2 x2のデータ)double[])
+        #   @return [int] 自由度
+        # @example
+        #   xi1 = [165, 130, 182, 178, 194, 206, 160, 122, 212, 165, 247, 195]
+        #   xi2 = [180, 180, 235, 270, 240, 285, 164, 152]
+        #   Num4TstStatisticLib.df4welch(xi1, xi2)
+        #   => 11
+        def df4welch(xi1, xi2)
+            return TstStatistic.df4welch(
+                xi1.to_java(Java::double), xi2.to_java(Java::double)
+            )
+        end
+        # 対応のある2つの母平均の差の検定量
+        #
+        # @overload diffPopulationMean(xi1, xi2)
+        #   @param [Array] xi1 x1のデータ(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]
+        #   xi2 = [36.8, 36.6, 36.5, 37.0, 36.0, 36.5, 36.6, 37.1, 36.4, 36.7]
+        #   Num4TstStatisticLib.diffPopulationMean(xi1, xi2)
+        #   => 2.283
+        # @note
+        #   自由度(N-1)のt分布に従う
+         def diffPopulationMean(xi1, xi2)
+            return TstStatistic.diffPopulationMean(
+                xi1.to_java(Java::double), xi2.to_java(Java::double)
+            )
+        end
+        # 2つの母分散の差の検定量
+        #
+        # @oerload diffPopulationVar(xi1, xi2)
+        #   @param [Array] xi1 x1のデータ(double[])
+        #   @param [Array] xi2 x2のデータ)double[])
+        #   @return [double] 検定統計量
+        # @example
+        #   xi1 = [165, 130, 182, 178, 194, 206, 160, 122, 212, 165, 247, 195]
+        #   xi2 = [180, 180, 235, 270, 240, 285, 164, 152]
+        #   Num4TstStatisticLib.diffPopulationVar(xi1, xi2)
+        #   => 0.4727
+        # @note
+        #   自由度(N1-1,N2-1)のF分布に従う
+        def diffPopulationVar(xi1, xi2)
+            return TstStatistic.diffPopulationVar(
+                xi1.to_java(Java::double), xi2.to_java(Java::double)
+            )
+        end
+        # 2つの母比率の差の検定量
+        #
+        # @overload diffPopulationRatio(m1, n1, m2, n2)
+        #   @param [int] m1 m1値
+        #   @param [int] n1 N1値
+        #   @param [int] m2 m2値
+        #   @param [int] n2 N2値
+        #   @return [double] 検定統計量
+        # @example
+        #   Num4TstStatisticLib.diffPopulationRatio(469, 1200, 308, 900)
+        #   => 2.283
+        # @note
+        #   標準正規分布 N(0,1*1)に従う
+        def diffPopulationRatio(m1, n1, m2, n2)
+            return TstStatistic.diffPopulationRatio(m1, n1, m2, n2)
+        end
+        # 無相関の検定量
+        #
+        # @overload unCorrelation(x, y)
+        #   @param [Array] x xのデータ(double[])
+        #   @param [Array] y yのデータ)double[])
+        #   @return [double] 検定統計量
+        # @example
+        #   x = [113, 64, 16, 45, 28, 19, 30, 82, 76]
+        #   y = [31, 5, 2, 17, 18, 2, 9, 25, 13]
+        #   Num4TstStatisticLib.unCorrelation(x, y)
+        #   => 3.1035
+        # @note
+        #   自由度(N-2)t分布に従う
+        def unCorrelation(x, y)
+            return TstStatistic.unCorrelation(
+                x.to_java(Java::double), y.to_java(Java::double)
+            )
+        end
+        # 母相関係数の検定量
+        #
+        # @overload populationCorre(x, y, rth0)
+        #   @param [Array] x xのデータ(double[])
+        #   @param [Array] y yのデータ)double[])
+        #   @param [double] rth0 母相関係数
+        #   @return [double] 検定統計量
+        # @example
+        #   x = [2750, 2956, 2675, 3198, 1816, 2233, 2375, 2288, 1932, 2036, 2183, 2882]
+        #   y = [249, 713, 1136, 575, 5654, 2107, 915, 4193, 7225, 3730, 472, 291]
+        #   Num4TstStatisticLib.populationCorre(x, y, -0.3)
+        #   => -2.107168
+        # @note
+        #   標準正規分布 N(0,1*1)に従う
+        def populationCorre(x, y, rth0)
+            return TstStatistic.populationCorre(
+                x.to_java(Java::double), y.to_java(Java::double), rth0
+            )
+        end
+    end
+end
+

metadata

@@ -0,0 +1,90 @@
+--- !ruby/object:Gem::Specification
+name: num4tststatistic
+version: !ruby/object:Gem::Version
+  version: 0.0.1
+platform: java
+authors:
+- siranovel
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2023-11-11 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: rake
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '12.3'
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 12.3.3
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '12.3'
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 12.3.3
+- !ruby/object:Gem::Dependency
+  name: rake-compiler
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.2'
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 1.2.5
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.2'
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 1.2.5
+description: numerical solution for test statistic
+email: siranovel@gmail.com
+executables: []
+extensions:
+- Rakefile
+extra_rdoc_files: []
+files:
+- CHANGELOG.md
+- Gemfile
+- LICENSE
+- Rakefile
+- ext/num4tststatistic/TstStatistic.java
+- lib/commons-math3-3.6.1.jar
+- lib/num4tststatistic.rb
+homepage: http://github.com/siranovel/num4tststatistic
+licenses:
+- MIT
+metadata: {}
+post_install_message: 
+rdoc_options: []
+require_paths:
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+required_rubygems_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+requirements: []
+rubygems_version: 3.3.7
+signing_key: 
+specification_version: 4
+summary: num for test statistic!
+test_files: []