問題タブ [hazard]
For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.
r - muhaz パッケージ ハザード関数の信頼区間
このパッケージは、カーネル平滑化法を使用して、右打ち切りデータからハザード関数muhazを推定します。私の質問は、計算するハザード関数の信頼区間を取得する方法はありますか?muhaz
上記の例では、muhaz.object fitにいくつかのエントリfit1$msemin,がありますがfit1$var.min、fit1$haz.estそれらの長さは の半分ですfit1$haz.est。
ハザード関数の信頼区間を抽出することが可能であれば、何か考えはありますか?
編集：@ user20650が提案したことに基づいて、次のことを試しました
max.time の設定は機能しているようで、すべてのブートストラップ サンプルには同じ推定グリッド ポイントがあります。ただし、得られたCIはほとんど意味がありません。通常、間隔は t=0 では狭く、時間とともに広くなる (情報が少なくなり、不確実性が増す) と予想されますが、得られた間隔は、時間とともに多かれ少なかれ一定であるように見えます。
r - 標準正規分布のハザード関数の計算
Mathematica UUPDE データベースにある式に基づいて 、R の標準正規分布のハザード関数をプロットしました。
特定の範囲では正しいようです。数値の問題は、より大きな値で発生します。添付の図を参照してください。以下は、完全な R コードです。
コメントをお待ちしております。
r - ハザード比のCoxph方向
次のデータセットがあるとします。
私がしたことは、参照 A にレベルを変更して実行したことです。
結果を見て、提供されたサンプル プロトコルと比較したところ、私の exp(-coeff) はケース B と C の両方で exp(coeff) に等しいことがわかりました。
したがって、参照 B を使用してコードを再度実行し、参照 C を使用してコードを 1 回実行したところ、下限と上限の両方がサンプルと一致することがわかりました。ただし、さらにログランクのp値が必要であり、これはこの方法では提供されないため、調べたいと思います:
逆ハザード比と conf. 制限と p 値が表示されます (つまり、逆モデルの場合)/ または、何か他の問題が発生したと思いますか?