num4tststatistic2 0.1.1 → 0.2.1

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 60785711542acd092cb1688f03e95bb31b372dab71e52048ec6adec7fcc2bcf5
-  data.tar.gz: fc3dd84abfa51c591061cde8e4ebad3e56985dbb0ae30928f8456a2488aa9c42
+  metadata.gz: 7959285c4da926ed687eb9d7292df6e49fdbf6e6f478f747060d4313b6b75a86
+  data.tar.gz: 26f6275766687e48dbd58ef413efc57edc2d29fea9bf200b9d8d845987f87fba
 SHA512:
-  metadata.gz: bccb6876cf7adb09a192dd34faf6fbd4f29652974e1cced0ae49bde5e30ce19c9c9a171db5744e2cd4128a3bbaf2a9d56c43720825693348c056bfe84b81e68b
-  data.tar.gz: 1f8c945bc94f924990b7451d809ed883c0e7dee6ab4b25ab1c08c822314d7a4c9e1a4ba5d8f0c749369efd319e579aedc83d3cd7c1abaa37dc54ea2acf69c9b3
+  metadata.gz: 97b33ebc741f49d7ac1e4b5ca9b9811d67458a5a58d1cb06353eca28422ea329d9656e402703f41a6a7f3a44ed93a5c31790d1e1cc4fff4c1add6a127af69220
+  data.tar.gz: d08688f8e75d54eafdda166e7b7055a008719a823bb2a62a0dc4f4c00fbc12aa166ef76d68fb1f07397c65eae762fead655613f4101ecbb5623b2857ab7028e0

data/CHANGELOG.md

@@ -2,6 +2,15 @@
 
 ## Unreleased
 
+## [0.2.1] - 2024-05-15
+### change
+- change module name to Num4CorrTestLib from CorrTestLib.
+
+## [0.1.2] - 2024-05-08
+
+### chg
+- change function name from diffPopulationMean2 to diffPopulationVarMean
+
 ## [0.1.1] - 2024-05-06
 
 ### add

data/lib/num4corrtest.rb

@@ -0,0 +1,151 @@
+require 'num4tststatistic'
+require 'hypothTest3'
+
+# 相関検定
+module Num4CorrTestLib
+    # 無相関の検定
+    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

@@ -1,6 +1,5 @@
 require 'num4tststatistic'
 require 'hypothTest3'
-require_relative('decorrtest')
 
 # 統計的仮説検定
 module Num4TstStatistic2Lib
@@ -70,9 +69,9 @@ module Num4TstStatistic2Lib
             statistic = @paraTest.populationRatio(m, n, p0)
             return @hypothTest3.normDistTest(statistic, a)
         end
-        # 2つの母平均の差の検定
+        # 2つの母平均の差の検定(等分散性check有り)
         #
-        # @overload diffPopulationMean2(xi1, xi2, a)
+        # @overload diffPopulationVarMean(xi1, xi2, a)
         #   @param [Array] xi1 x1のデータ(double[])
         #   @param [Array] xi2 x2のデータ(double[])
         #   @param  [double] a  有意水準
@@ -82,10 +81,11 @@ module Num4TstStatistic2Lib
         #   xi2 = [180, 180, 235, 270, 240, 285, 164, 152]
         #   hypothTest = Num4HypothTestLib::TwoSideTestLib.new
         #   paraTest = Num4TstStatistic2Lib::ParametrixTestLib.new(hypothTest)
-        #   paraTest.diffPopulationMean2(xi1, xi2, 0.05)
+        #   paraTest.diffPopulationVarMean(xi1, xi2, 0.05)
         #   => false
-        def diffPopulationMean2(xi1, xi2, a)
+        def diffPopulationVarMean(xi1, xi2, a)
             bRet = diffPopulationVar(xi1, xi2, a)
+
             if bRet == true # 等分散ではない
                 return diffPopulationMean2UnEquVar(xi1, xi2, a)
             else            # 等分散性
@@ -152,6 +152,7 @@ module Num4TstStatistic2Lib
         #   => true
         def diffPopulationMean(xi1, xi2, a)
             raise TypeError unless @hypothTest3.kind_of?(HypothTest3IF)
+            raise RangeError unless xi1.size == xi2.size
 
             n = xi1.size
             df = n - 1
@@ -343,143 +344,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
 end
-# 相関検定
-module DecorrTestLib
-    # 無相関の検定
-    class UnDecorrTestLib < DecorrTestIF
-        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]
-        #   decorrTest = DecorrTestLib::UnDecorrTestLib.new
-        #   decorrTest.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]
-        #   decorrTest = DecorrTestLib::UnDecorrTestLib.new
-        #   decorrTest.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]
-        #   decorrTest = DecorrTestLib::UnDecorrTestLib.new
-        #   decorrTest.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
-    # 相関係数の検定
-    class CorreFactLib < CorreFactIF
-        def initialize(hypothTest3)
-            @hypothTest3 = hypothTest3
-            @paraTest = Num4TstStatisticLib::ParametrixTestLib.new
-            @nonParaTest = Num4TstStatisticLib::NonParametrixTestLib.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
-        #   decorrTest = DecorrTestLib::CorreFactLib.new(hypothTest)
-        #   decorrTest.pearsoCorrelation(x, y,  -0.3, 0.05)
-        #   => true
-        def pearsoCorrelation(x, y, rth0, a)
-            raise TypeError unless @hypothTest3.kind_of?(HypothTest3IF)
-            statistic = @paraTest.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
-        #   decorrTest = DecorrTestLib::CorreFactLib.new(hypothTest)
-        #   decorrTest.spearmanscorr(x, y, -0.3, 0.05)
-        #   => true
-        def spearmanscorr(x, y, rth0, a)
-            raise TypeError unless @hypothTest3.kind_of?(HypothTest3IF)
-            statistic = @nonParaTest.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
-        #   decorrTest = DecorrTestLib::CorreFactLib.new(hypothTest)
-        #   decorrTest.kendallscorr(x, y, -0.3, 0.05)
-        #   => true
-        def kendallscorr(x, y, rth0, a)
-            raise TypeError unless @hypothTest3.kind_of?(HypothTest3IF)
-            statistic = @nonParaTest.kendallscorr(x, y)
-            return @hypothTest3.populationCorre(statistic, x.size, rth0, a)
-        end
-    end
-end
+

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: num4tststatistic2
 version: !ruby/object:Gem::Version
-  version: 0.1.1
+  version: 0.2.1
 platform: ruby
 authors:
 - siranovel
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2024-05-06 00:00:00.000000000 Z
+date: 2024-05-15 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: num4tststatistic
@@ -59,7 +59,7 @@ files:
 - CHANGELOG.md
 - Gemfile
 - LICENSE
-- lib/decorrtest.rb
+- lib/num4corrtest.rb
 - lib/num4tststatistic2.rb
 homepage: https://github.com/siranovel/num4tststatistic2
 licenses:

data/lib/decorrtest.rb

@@ -1,22 +0,0 @@
-class DecorrTestIF
-    def pearsoCorrelation(x, y, a)
-        raise NotImplementedError.new("#{self.class}##{__method__} が実装されていません")
-    end
-    def spearmanscorr(x, y, a)
-        raise NotImplementedError.new("#{self.class}##{__method__} が実装されていません")
-    end
-    def kendallscorr(x, y, a)
-        raise NotImplementedError.new("#{self.class}##{__method__} が実装されていません")
-    end
-end
-class CorreFactIF
-    def pearsoCorrelation(x, y, rth0, a)
-        raise NotImplementedError.new("#{self.class}##{__method__} が実装されていません")
-    end
-    def spearmanscorr(x, y, rth0, a)
-        raise NotImplementedError.new("#{self.class}##{__method__} が実装されていません")
-    end
-    def kendallscorr(x, y, rth0, a)
-        raise NotImplementedError.new("#{self.class}##{__method__} が実装されていません")
-    end
-end