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

    Function mcreg

    • y (検査法) の x (基準法) への測定法比較の回帰です: Deming、重みつき Deming、または Passing-Bablok の回帰で、解析的、ジャックナイフ、またはブートストラップの区間を伴います。R の mcr::mcreg() と同等です。

      • Deming ("Deming"): 傾きは、誤差分散による尺度での二乗距離を最小にするものです。λ = errorRatio = σ²(x)/σ²(y) です (mcr の規約。λ → 0 のとき、y の x への最小二乗の傾きになります)。解析的な標準誤差は Strike (1991) によるもので、t 区間は自由度 n - 2 を用います。
      • 重みつき Deming ("WDeming"): 誤差が水準に比例する場合 (変動係数が一定) の Linnet (1990) の反復です。mcr と同じく、両方の係数の変化が threshold (絶対値) より小さくなると止めます。iterMax 回の反復に達すると警告を出します。解析的な区間はありません (既定はジャックナイフ)。値は正でなければなりません。threshold は絶対値のため、非常に大きな大きさのデータ (約 1e150 程度) では収束しないことがあり、非常に小さな大きさのデータ (約 1e-150 程度) では最初の反復の後で止まります。threshold はデータの単位で選んでください。
      • Passing-Bablok ("PaBa"): n(n - 1)/2 個の対ごとの傾きのシフトされた中央値と、Passing と Bablok (1983) の順位区間です。切片は y - b x の中央値で、その限界は傾きの限界と同じです。相対で 1e-12 より小さい差は同順位として扱われます (mcr の規則。x だけが同順位の対は垂直です)。そのため、大きなオフセット (1e12 + O(1)) を持つデータでは、差の大部分が 0 になります。中央の 2 つの傾きの平均は、平均の角度の tan (slopeMeasure: "radian"、mcr の既定。y を再スケールすると変わります。角度は尺度変換されないためです) か、平均の傾き ("tangent") です。区間が定義されない場合 (点が少なすぎる場合、n <= 4 では常に)、傾きと切片の限界は -Infinity と Infinity になり、警告を出します。値は負であってはいけません。結果には線形性の CUSUM 統計量 (cusum。検定の判定は行いません) が含まれます。
      • ジャックナイフ: 1 つ抜きの推定値、標準誤差 sd(擬似値)/sqrt(n)、および区間の推定値 ± t se (自由度 n - 2、Linnet (1993)) です。推定値は全データによるものです。
      • ブートストラップ: 対の nsamples 個の再標本を復元抽出します。"quantile" はパーセンタイル区間、"Student" は推定値 ± t sd(B)、"BCa" は偏りを補正し加速を加えた区間 (ジャックナイフによる加速)、"tBoot" は解析的な標準誤差を用いたブートストラップ t 区間 (Deming の場合)、または "nestedbootstrap" の場合は、各再標本の内側の nnested 個の再標本の標準偏差です。x または y が定数の再標本は、mcr の警告とともに除かれます (5% を超える場合はエラーです)。再標本は、現在の乱数ストリームから R の sample(1:n, n, replace = TRUE) として引かれるか、mcr の rng.seed として set.seed(rngSeed) から引かれます (生成器の状態は、あった場合は復元されます)。そのため、同じ種に対して値は mcr のものと一致します。点が非常に少ない場合 (n = 5) には、相異なる 2 点か 3 点の再標本の標準誤差が 0 になるため、ブートストラップ t の限界が無限大や非常に大きな値になることがあります。

      mcr と異なり、既定の区間は、ランダムなブートストラップではなく "analytical" (重みつき Deming では "jackknife") です。confLevel は mcr の alpha (= 1 - confLevel) に代わるものです。2 点以下、x または y が定数、無限大の値、Deming の共分散が 0 のデータは、エラーになります。Passing-Bablok の奇数位置の傾きは dy/dx です (mcr では角度の tan で、数 ulp ずれます)。PaBaLarge、LinReg、WLinReg、PBequi、TS は提供していません。

      Parameters

      • x: Sample

        基準法による測定値です。

      • y: Sample

        検査法による測定値で、x と対になっています。

      • options: McregOptions = {}

        McregOptions を参照してください。

      戻り値 McregResult

      McregResult を参照してください。

      x と y の長さが異なる場合は TypeError です。naRm を指定しない欠測値、無限大の値、3 対未満、x または y が定数、負の値 (Passing-Bablok)、正でない値 (重みつき Deming)、5793 点を超える場合 (Passing-Bablok)、重みつき Deming の解析的な区間、Passing-Bablok または重みつき Deming で入れ子の再標本のないブートストラップ t 区間、退化した再標本が多すぎる場合、共分散が 0 の場合 (Deming)、または無効なオプションの場合は RangeError です。

      const x = [1, 2, 3, 4, 5, 6, 7];
      const y = [1.2, 1.9, 3.2, 3.8, 5.1, 6.4, 6.6];
      const d = mcreg(x, y, { methodReg: "Deming" });
      d.slope.estimate; // 0.97547950725034482
      d.slope.confInt; // [0.83267195599246846, 1.1182870585082212]
      const pb = mcreg(x, y); // Passing-Bablok
      pb.slope.estimate; // 0.95
      pb.intercept.confInt; // [-0.8524794459077718, 0.80365513709218339]