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

    Function coxph

    • Cox 比例ハザード回帰です。R の survival::coxph() と同等です。右側打ち切りのデータについて、coxph(Surv(time, status) ~ x [+ strata(strata)], weights, ties, init, robust, cluster, control) の形で呼び出した場合と同じく、summary(fit) の値を得ます。

      robust を指定したとき (既定では、cluster を指定したとき、または観測値の重みに整数でないものがあるとき) は、分散はロバストな (Lin-Wei のサンドイッチ) 分散になります。これは、クラスターごとに合計したものです (cluster がなければ、観測値ごとに独自のクラスター)。R と同じく、var はロバスト分散、naiveVar はモデルに基づく分散です。se (と stdError) はモデルに基づく標準誤差のままです (R の se(coef))。一方、z、p、区間 lower/upper、wald、waldtest はロバストなもの (robustSe) を使います。robscore は init におけるロバストスコア検定です。尤度比検定、スコア検定、R^2 は、クラスター化を無視します。

      時刻、status、共変量、重み、層のいずれかが欠測している行は除外されます。互いの差が sqrt(eps) 以内の時刻は、同順位として扱われます。nocenter に値がすべて含まれる共変量は中心化されません。そのほかの共変量は、重み付き平均で中心化されます。これは線形予測子をずらすだけです。R の警告 (非収束、係数が無限大になる可能性) は、warnings で返されます。

      Parameters

      • time: readonly number[]

        追跡時間。

      • status: readonly number[]

        イベントなら 1、打ち切りなら 0。

      • x: readonly (readonly number[])[]

        共変量の行列。観測値ごとに 1 行 (列は少なくとも 1 つ)。

      • options: CoxphOptions = {}

        CoxphOptions を参照してください。

      戻り値 CoxphResult

      robust が真偽値でない場合は TypeError です。不正な引数に対しては RangeError です (無限大の重み、2 つ未満のクラスターでのロバスト分散を含みます)。

      const fit = coxph([5, 3, 8, 1, 9, 4], [1, 1, 0, 1, 0, 1], [[1.2], [0.3], [-0.5], [2], [0.7], [0.1]]);
      fit.coefficients; // [1.0324]
      fit.logtest.pValue; // 0.18817