num4tststatistic2 0.0.2 → 0.1.2

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 1cde00758d2f1f56ba425ce885da276d76272564848fd176dc78ab7e271ffdc5
-  data.tar.gz: c0631d73f042ab0bc7be1be74f408cc298cf0237a6ad56c7b9342edfbf8ef807
+  metadata.gz: cfd90c54dab3784be95fad71d373a338e24c17a41b0bb5ddde46e46d77a4e739
+  data.tar.gz: 6fa678137d1a9b7fa3336a8f61bdc3fe03d89cb67c6ece3ff5ddb27c5216cfb3
 SHA512:
-  metadata.gz: 80815f8945e899f8699bb3d06a27abc41c8b00967ab9f2b6b0badb2540c19a13f123f8a376c05358f277cc10c6947e6c4b028a5738d91b3f38002f6f614fd410
-  data.tar.gz: bd512533a00f76c4e6ed738c25bfb6cc61af2c634d5b28dac6236844e03e15088e5a846acead538b4dd8e1e4e38d99b258b5bc945fe68912f96999a519842374
+  metadata.gz: 20ec74c474a5d8524ab8c18eac95b33ac5c0bc96241a8aa8dea2aa3a2bb41d633c90aad4561611f25f421b68af754b5ebae32ac08a1ceb73b0ae672a1c9714ef
+  data.tar.gz: 2520663599be76415a0fb3a6106d096f5697be92fdbf1371affff34ac4f9aa1b6da0f24e36c5ef7150e6cd2eec05d94207fd38be17a97552ab78326d3e583e6a

data/CHANGELOG.md

@@ -2,6 +2,17 @@
 
 ## Unreleased
 
+## [0.1.2] - 2024-05-08
+
+### chg
+- change function name from diffPopulationMean2 to diffPopulationVarMean
+
+## [0.1.1] - 2024-05-06
+
+### add
+- add function of diffPopulationMean2.
+- add CorreFact.
+
 ## [0.0.2] - 2024-04-22
 
 ### add

data/lib/corrtest.rb

@@ -0,0 +1,151 @@
+require 'num4tststatistic'
+require 'hypothTest3'
+
+# 相関検定
+module CorrTestLib
+    # 無相関の検定
+    class DecorrTestLib
+        def initialize
+            @corr = CorrStatisticLib.new
+            @hypothTest = Num4HypothTestLib::DecorrTestLib.new
+        end
+        # ピアソン相関係数
+        #
+        # @overload pearsoCorrelation(x, y, a)
+        #   @param [Array] x xのデータ(double[])
+        #   @param [Array] y yのデータ(double[])
+        #   @param  [double] a  有意水準
+        #   @return [boolean] 検定結果(true:棄却域内 false:棄却域外)
+        # @example
+        #   x = [113, 64, 16, 45, 28, 19, 30, 82, 76]
+        #   y = [31, 5, 2, 17, 18, 2, 9, 25, 13]
+        #   corrTest = CorrTestLib::DecorrTestLib.new
+        #   corrTest.pearsoCorrelation(x, y, 0.05)
+        #   => true
+        def pearsoCorrelation(x, y, a)
+            df = x.size - 2
+            statistic = @corr.pearsoCorrelation(x, y)
+            return @hypothTest.twoSideTest(statistic, df, a)
+        end
+        # スピアマンの順位相関係数
+        #
+        # @overload spearmanscorr(x, y, a)
+        #   @param [Array] x xのデータ(double[])
+        #   @param [Array] y yのデータ(double[])
+        #   @param  [double] a  有意水準
+        #   @return [boolean] 検定結果(true:棄却域内 false:棄却域外)
+        # @example
+        #   x = [113, 64, 16, 45, 28, 19, 30, 82, 76]
+        #   y = [31, 5, 2, 17, 18, 2, 9, 25, 13]
+        #   corrTest = CorrTestLib::DecorrTestLib.new
+        #   corrTest.spearmanscorr(x, y, 0.05)
+        #   => true
+        def spearmanscorr(x, y, a)
+            df = x.size - 2
+            statistic = @corr.spearmanscorr(x, y)
+            return @hypothTest.twoSideTest(statistic, df, a)
+        end
+        # ケンドールの順位相関係数
+        #
+        # @overload kendallscorr(x, y, a)
+        #   @param [Array] x xのデータ(double[])
+        #   @param [Array] y yのデータ(double[])
+        #   @param  [double] a  有意水準
+        #   @return [boolean] 検定結果(true:棄却域内 false:棄却域外)
+        # @example
+        #   x = [113, 64, 16, 45, 28, 19, 30, 82, 76]
+        #   y = [31, 5, 2, 17, 18, 2, 9, 25, 13]
+        #   corrTest = CorrTestLib::DecorrTestLib.new
+        #   corrTest.kendallscorr(x, y, 0.05)
+        #   => false
+        def kendallscorr(x, y, a)
+            df = x.size - 2
+            statistic = @corr.kendallscorr(x, y)
+            return @hypothTest.twoSideTest(statistic, df, a)
+        end
+    end
+    # 母相関係数の検定
+    class CorreFactLib
+        def initialize(hypothTest3)
+            @hypothTest3 = hypothTest3
+            @corr = CorrStatisticLib.new
+        end
+        # ピアソン相関係数
+        #
+        # @overload pearsoCorrelation(x, y, rth0, a)
+        #   @param [Array] x xのデータ(double[])
+        #   @param [Array] y yのデータ(double[])
+        #   @param [double] rth0 母相関係数
+        #   @param  [double] a  有意水準
+        #   @return [boolean] 検定結果(true:棄却域内 false:棄却域外)
+        # @example
+        #   x = [113, 64, 16, 45, 28, 19, 30, 82, 76]
+        #   y = [31, 5, 2, 17, 18, 2, 9, 25, 13]
+        #   hypothTest = Num4HypothTestLib::TwoSideTestLib.new
+        #   corrTest = CorrTestLib::CorreFactLib.new(hypothTest)
+        #   corrTest.pearsoCorrelation(x, y,  -0.3, 0.05)
+        #   => true
+        def pearsoCorrelation(x, y, rth0, a)
+            raise TypeError unless @hypothTest3.kind_of?(HypothTest3IF)
+            statistic = @corr.pearsoCorrelation(x, y)
+            return @hypothTest3.populationCorre(statistic, x.size, rth0, a)
+        end
+        # スピアマンの順位相関係数
+        #
+        # @overload spearmanscorr(x, y, rth0, a)
+        #   @param [Array] x xのデータ(double[])
+        #   @param [Array] y yのデータ(double[])
+        #   @param [double] rth0 母相関係数
+        #   @param  [double] a  有意水準
+        #   @return [boolean] 検定結果(true:棄却域内 false:棄却域外)
+        # @example
+        #   x = [113, 64, 16, 45, 28, 19, 30, 82, 76]
+        #   y = [31, 5, 2, 17, 18, 2, 9, 25, 13]
+        #   hypothTest = Num4HypothTestLib::TwoSideTestLib.new
+        #   corrTest = CorrTestLib::CorreFactLib.new(hypothTest)
+        #   corrTest.spearmanscorr(x, y, -0.3, 0.05)
+        #   => true
+        def spearmanscorr(x, y, rth0, a)
+            raise TypeError unless @hypothTest3.kind_of?(HypothTest3IF)
+            statistic = @corr.spearmanscorr(x, y)
+            return @hypothTest3.populationCorre(statistic, x.size, rth0, a)
+        end
+        # ケンドールの順位相関係数
+        #
+        # @overload kendallscorr(x, y, rth0, a)
+        #   @param [Array] x xのデータ(double[])
+        #   @param [Array] y yのデータ(double[])
+        #   @param [double] rth0 母相関係数
+        #   @param  [double] a  有意水準
+        #   @return [boolean] 検定結果(true:棄却域内 false:棄却域外)
+        # @example
+        #   x = [113, 64, 16, 45, 28, 19, 30, 82, 76]
+        #   y = [31, 5, 2, 17, 18, 2, 9, 25, 13]
+        #   hypothTest = Num4HypothTestLib::TwoSideTestLib.new
+        #   corrTest = CorrTestLib::CorreFactLib.new(hypothTest)
+        #   corrTest.kendallscorr(x, y, -0.3, 0.05)
+        #   => true
+        def kendallscorr(x, y, rth0, a)
+            raise TypeError unless @hypothTest3.kind_of?(HypothTest3IF)
+            statistic = @corr.kendallscorr(x, y)
+            return @hypothTest3.populationCorre(statistic, x.size, rth0, a)
+        end
+    end
+
+    class CorrStatisticLib
+        def initialize
+            @paraTest = Num4TstStatisticLib::ParametrixTestLib.new
+            @nonParaTest = Num4TstStatisticLib::NonParametrixTestLib.new
+        end
+        def pearsoCorrelation(x, y)
+            return @paraTest.pearsoCorrelation(x, y)
+        end
+        def spearmanscorr(x, y)
+            return @nonParaTest.spearmanscorr(x, y)
+        end
+        def kendallscorr(x, y)
+            return  @nonParaTest.kendallscorr(x, y)
+        end
+    end
+    private_constant :CorrStatisticLib
+end

data/lib/num4tststatistic2.rb

@@ -3,6 +3,7 @@ require 'hypothTest3'
 
 # 統計的仮説検定
 module Num4TstStatistic2Lib
+    # パラメトリック検定
     class ParametrixTestLib
         def initialize(hypothTest3)
             @hypothTest3 = hypothTest3
@@ -68,6 +69,29 @@ module Num4TstStatistic2Lib
             statistic = @paraTest.populationRatio(m, n, p0)
             return @hypothTest3.normDistTest(statistic, a)
         end
+        # 2つの母平均の差の検定
+        #
+        # @overload diffPopulationVarMean(xi1, xi2, a)
+        #   @param [Array] xi1 x1のデータ(double[])
+        #   @param [Array] xi2 x2のデータ(double[])
+        #   @param  [double] a  有意水準
+        #   @return [boolean] 検定結果(true:棄却域内 false:棄却域外)
+        # @example
+        #   xi1 = [165, 130, 182, 178, 194, 206, 160, 122, 212, 165, 247, 195]
+        #   xi2 = [180, 180, 235, 270, 240, 285, 164, 152]
+        #   hypothTest = Num4HypothTestLib::TwoSideTestLib.new
+        #   paraTest = Num4TstStatistic2Lib::ParametrixTestLib.new(hypothTest)
+        #   paraTest.diffPopulationVarMean(xi1, xi2, 0.05)
+        #   => false
+        def diffPopulationVarMean(xi1, xi2, a)
+            bRet = diffPopulationVar(xi1, xi2, a)
+
+            if bRet == true # 等分散ではない
+                return diffPopulationMean2UnEquVar(xi1, xi2, a)
+            else            # 等分散性
+                return diffPopulationMean2EquVar(xi1, xi2, a)
+            end
+        end
         # 2つの母平均の差の検定(等分散性を仮定)
         #
         # @overload diffPopulationMean2EquVar(xi1, xi2, a)
@@ -222,6 +246,7 @@ module Num4TstStatistic2Lib
             return @hypothTest3.chi2DistTest(statistic, df, a)
         end
     end
+    # ノンパラメトリック検定
     class NonParametrixTestLib
         def initialize(hypothTest3)
             @hypothTest3 = hypothTest3
@@ -287,6 +312,7 @@ module Num4TstStatistic2Lib
             return @nonParaTest.ks2test(xi1, xi2, a)
         end
     end
+    # 外れ値検定
     class OutlierLib
         def initialize
             @outlier = Num4TstStatisticLib::OutlierLib.new
@@ -317,73 +343,11 @@ module Num4TstStatistic2Lib
         # @example
         #   xi = [3.4, 3.5, 3.3, 2.2, 3.3, 3.4, 3.6, 3.2]
         #   outlier = Num4TstStatistic2Lib::OutlierLib.new
-        #   outlier.grubbs("LDH", xi)
+        #   outlier.errbar("LDH", xi)
         #   => errbar.jpeg
         def errbar(dname, xi)
             @outlier.errbar(dname, xi)
         end
     end
-    # 無相関の検定
-    class DecorrTestLib
-        def initialize
-            @paraTest = Num4TstStatisticLib::ParametrixTestLib.new
-            @nonParaTest = Num4TstStatisticLib::NonParametrixTestLib.new
-            @hypothTest = Num4HypothTestLib::DecorrTestLib.new
-        end
-        # ピアソン相関係数
-        #
-        # @overload pearsoCorrelation(x, y, a)
-        #   @param [Array] x xのデータ(double[])
-        #   @param [Array] y yのデータ(double[])
-        #   @param  [double] a  有意水準
-        #   @return [boolean] 検定結果(true:棄却域内 false:棄却域外)
-        # @example
-        #   x = [113, 64, 16, 45, 28, 19, 30, 82, 76]
-        #   y = [31, 5, 2, 17, 18, 2, 9, 25, 13]
-        #   paraTest = Num4TstStatistic2Lib::DecorrTestLib.new
-        #   paraTest.pearsoCorrelation(x, y, 0.05)
-        #   => true
-        def pearsoCorrelation(x, y, a)
-            df = x.size - 2
-            statistic = @paraTest.pearsoCorrelation(x, y)
-            return @hypothTest.twoSideTest(statistic, df, a)
-        end
-        # スピアマンの順位相関係数
-        #
-        # @overload spearmanscorr(x, y, a)
-        #   @param [Array] x xのデータ(double[])
-        #   @param [Array] y yのデータ(double[])
-        #   @param  [double] a  有意水準
-        #   @return [boolean] 検定結果(true:棄却域内 false:棄却域外)
-        # @example
-        #   x = [113, 64, 16, 45, 28, 19, 30, 82, 76]
-        #   y = [31, 5, 2, 17, 18, 2, 9, 25, 13]
-        #   paraTest = Num4TstStatistic2Lib::DecorrTestLib.new
-        #   paraTest.spearmanscorr(x, y, 0.05)
-        #   => true
-        def spearmanscorr(x, y, a)
-            df = x.size - 2
-            statistic = @nonParaTest.spearmanscorr(x, y)
-            return @hypothTest.twoSideTest(statistic, df, a)
-        end
-        # ケンドールの順位相関係数
-        #
-        # @overload kendallscorr(x, y, a)
-        #   @param [Array] x xのデータ(double[])
-        #   @param [Array] y yのデータ(double[])
-        #   @param  [double] a  有意水準
-        #   @return [boolean] 検定結果(true:棄却域内 false:棄却域外)
-        # @example
-        #   x = [113, 64, 16, 45, 28, 19, 30, 82, 76]
-        #   y = [31, 5, 2, 17, 18, 2, 9, 25, 13]
-        #   paraTest = Num4TstStatistic2Lib::DecorrTestLib.new
-        #   paraTest.kendallscorr(x, y, 0.05)
-        #   => false
-        def kendallscorr(x, y, a)
-            df = x.size - 2
-            statistic = @nonParaTest.kendallscorr(x, y)
-            return @hypothTest.twoSideTest(statistic, df, a)
-        end
-    end
 end
 

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: num4tststatistic2
 version: !ruby/object:Gem::Version
-  version: 0.0.2
+  version: 0.1.2
 platform: ruby
 authors:
 - siranovel
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2024-04-22 00:00:00.000000000 Z
+date: 2024-05-08 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: num4tststatistic
@@ -59,8 +59,9 @@ files:
 - CHANGELOG.md
 - Gemfile
 - LICENSE
+- lib/corrtest.rb
 - lib/num4tststatistic2.rb
-homepage: http://github.com/siranovel/num4tststatistic2
+homepage: https://github.com/siranovel/num4tststatistic2
 licenses:
 - MIT
 metadata: {}