rStats.js v1.3.0
    検索インデックスを準備しています...

    Function fisherTest

    • 度数データに対する Fisher の正確検定です。R の fisher.test() と同等です。

      周辺和が固定された r x c 分割表の独立性を検定します。2 x 2 表では、オッズ比の条件付き最尤推定値とその信頼区間も返します。R の数値を再現します (R の求根の許容誤差も含みます)。両側 p 値は、観測された表の確からしさ以下の表の確率を合計したものです (相対許容誤差 1e-7、R と同じ)。非整数の度数は、最も近い整数に丸められます (端数が 0.5 のときは偶数へ)。その際、R と同じ警告が出ます ("'x' has been rounded to integer: Mean relative difference: 0.0625")。これは度数と丸めた度数を R の all.equal で比較したときの文面です。差が約 1.5e-8 未満のときは出ません。

      simulatePValue: true のとき、2 x 2 より大きい表の p 値は、観測された周辺和をもつ B 個の乱数表から計算されます。乱数は setSeed の乱数ストリームから、R とまったく同じ方法で抽出されます。

      midp: true (2 x 2 のみ) のとき、exact2x2 の fisher.exact() の mid-p 版になります。p 値と区間は mid-p の裾から求められます (両側: central の方法で、小さい方の裾の 2 倍。最大 1)。限界は、mid-p の裾の根を倍精度の精度で求めたものです (exact2x2 の探索は、そこから最大 1% ずれたところで止まります)。method は "Central Fisher's Exact Test (mid-p version)" または "One-sided Fisher's Exact Test (mid-p version)" になります。

      未実装: R の hybrid オプション。

      Parameters

      戻り値 FisherTestResult

      検定の結果です。詳細は FisherTestResult を参照してください。

      const r = fisherTest([[3, 1], [1, 3]]);
      r.pValue; // 0.4857142857142857
      r.estimate; // 6.408319658199675 (conditional MLE of the odds ratio; R: 6.408309)
      r.confInt; // [0.2117355954465785, 626.2435305888127] (R: 0.2117329, 621.93)
      const m = fisherTest([[3, 1], [1, 3]], { midp: true });
      m.pValue; // 0.2571428571428571
      m.confInt; // [0.31005488171907974, 308.5567715811547] (exact2x2: 0.3100451, 311.4363)