num4tststatistic2 0.0.2 → 0.1.1

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 1cde00758d2f1f56ba425ce885da276d76272564848fd176dc78ab7e271ffdc5
-  data.tar.gz: c0631d73f042ab0bc7be1be74f408cc298cf0237a6ad56c7b9342edfbf8ef807
+  metadata.gz: 60785711542acd092cb1688f03e95bb31b372dab71e52048ec6adec7fcc2bcf5
+  data.tar.gz: fc3dd84abfa51c591061cde8e4ebad3e56985dbb0ae30928f8456a2488aa9c42
 SHA512:
-  metadata.gz: 80815f8945e899f8699bb3d06a27abc41c8b00967ab9f2b6b0badb2540c19a13f123f8a376c05358f277cc10c6947e6c4b028a5738d91b3f38002f6f614fd410
-  data.tar.gz: bd512533a00f76c4e6ed738c25bfb6cc61af2c634d5b28dac6236844e03e15088e5a846acead538b4dd8e1e4e38d99b258b5bc945fe68912f96999a519842374
+  metadata.gz: bccb6876cf7adb09a192dd34faf6fbd4f29652974e1cced0ae49bde5e30ce19c9c9a171db5744e2cd4128a3bbaf2a9d56c43720825693348c056bfe84b81e68b
+  data.tar.gz: 1f8c945bc94f924990b7451d809ed883c0e7dee6ab4b25ab1c08c822314d7a4c9e1a4ba5d8f0c749369efd319e579aedc83d3cd7c1abaa37dc54ea2acf69c9b3

data/CHANGELOG.md

@@ -2,6 +2,12 @@
 
 ## Unreleased
 
+## [0.1.1] - 2024-05-06
+
+### add
+- add function of diffPopulationMean2.
+- add CorreFact.
+
 ## [0.0.2] - 2024-04-22
 
 ### add

data/lib/decorrtest.rb

@@ -0,0 +1,22 @@
+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

data/lib/num4tststatistic2.rb

@@ -1,8 +1,10 @@
 require 'num4tststatistic'
 require 'hypothTest3'
+require_relative('decorrtest')
 
 # 統計的仮説検定
 module Num4TstStatistic2Lib
+    # パラメトリック検定
     class ParametrixTestLib
         def initialize(hypothTest3)
             @hypothTest3 = hypothTest3
@@ -68,6 +70,28 @@ module Num4TstStatistic2Lib
             statistic = @paraTest.populationRatio(m, n, p0)
             return @hypothTest3.normDistTest(statistic, a)
         end
+        # 2つの母平均の差の検定
+        #
+        # @overload diffPopulationMean2(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.diffPopulationMean2(xi1, xi2, 0.05)
+        #   => false
+        def diffPopulationMean2(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
@@ -323,8 +349,11 @@ module Num4TstStatistic2Lib
             @outlier.errbar(dname, xi)
         end
     end
+end
+# 相関検定
+module DecorrTestLib
     # 無相関の検定
-    class DecorrTestLib
+    class UnDecorrTestLib < DecorrTestIF
         def initialize
             @paraTest = Num4TstStatisticLib::ParametrixTestLib.new
             @nonParaTest = Num4TstStatisticLib::NonParametrixTestLib.new
@@ -340,8 +369,8 @@ module Num4TstStatistic2Lib
         # @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)
+        #   decorrTest = DecorrTestLib::UnDecorrTestLib.new
+        #   decorrTest.pearsoCorrelation(x, y, 0.05)
         #   => true
         def pearsoCorrelation(x, y, a)
             df = x.size - 2
@@ -358,8 +387,8 @@ module Num4TstStatistic2Lib
         # @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)
+        #   decorrTest = DecorrTestLib::UnDecorrTestLib.new
+        #   decorrTest.spearmanscorr(x, y, 0.05)
         #   => true
         def spearmanscorr(x, y, a)
             df = x.size - 2
@@ -376,8 +405,8 @@ module Num4TstStatistic2Lib
         # @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)
+        #   decorrTest = DecorrTestLib::UnDecorrTestLib.new
+        #   decorrTest.kendallscorr(x, y, 0.05)
         #   => false
         def kendallscorr(x, y, a)
             df = x.size - 2
@@ -385,5 +414,72 @@ module Num4TstStatistic2Lib
             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.0.2
+  version: 0.1.1
 platform: ruby
 authors:
 - siranovel
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2024-04-22 00:00:00.000000000 Z
+date: 2024-05-06 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: num4tststatistic
@@ -59,8 +59,9 @@ files:
 - CHANGELOG.md
 - Gemfile
 - LICENSE
+- lib/decorrtest.rb
 - lib/num4tststatistic2.rb
-homepage: http://github.com/siranovel/num4tststatistic2
+homepage: https://github.com/siranovel/num4tststatistic2
 licenses:
 - MIT
 metadata: {}