問題タブ [survival-analysis]

このサバイバル モデルを作成しました (複数の曲線を使用)。

..そして結果のフィットオブジェクトは結果を に保存しますfit$surv。ただし、私の場合、生存推定の長さは、さまざまな曲線/グループで等しくありません。さまざまな曲線の長さは にありfit$strataます。

基本的に、データは次のようになります。

data.tableこのデータを同じグループ長に抽出する必要があります..

これを行う簡単な方法はありますか?または、明らかなことを完全に見逃していますか?

Rで生存関数の導関数を取り、それを新しい関数に保存しようとしています。

これまでの私のコードは次のとおりです。

生存パッケージのSurvおよび関数を使用して生存分析を実行しようとしています。survfitほとんどのデータが切り捨てられたままになっていて、Surv正しく機能します。私の応答変数は、橋が欠陥があると分類されてから崩壊するまでの時間 (年単位で測定) です。2012 年から 1992 年までさかのぼって各橋の欠陥状況を追跡できますが、それ以上は追跡できません。検閲が行われるのは、多くの橋が崩壊時から 1992 年にさかのぼって欠陥があると分類されたため、橋が欠陥になった正確な時期がわからないため、橋の真の「寿命」(橋梁からの年数) がわからないためです。崩壊する不完全な分類）。たとえば、橋が 1995 年に崩壊し、1995 年、1994 年、1993 年、1992 年に欠陥があると分類されたとします。最初に欠陥があると分類されたのは 1992 年である可能性もあり、欠陥があると分類された可能性もあります。 1984年以来。

いくつかのサンプルデータ:

8 つの橋のうち 4 つ (橋 1、2、4、および 8) が切り捨てられたままになっていることに注意してください。それらが最初に欠陥として分類された時期を正確に観察することはできません。ブリッジ 3、5、6、および 7 については、1992 年以降に欠陥が生じて以来、正確な寿命を知っているため、これらの観察は打ち切られていません。

次に、以下のモデルに適合します。

このモデルが正しいかどうかはわかりません。一つには、ドキュメントでは、「時間」は右打ち切りデータまたは間隔打ち切りデータの開始時間であると書かれています。別の理由として、このモデルが打ち切られていない観察をどのように説明しているかわかりません。これが正しいかどうか、そうでない場合は、何を変更する必要があるのか​​、その理由を教えてください。どんな助けでも大歓迎です。本当にありがとう！

cv.glm が glm オブジェクトでどのように機能するかはわかりますが、適合生存モデルはどうでしょうか?

たくさんのモデルがあります (Weibull、Gompertz、lognormal など)。交差検証を使用して予測誤差を評価したいと考えています。Rでこれを行うことができるパッケージ/機能はどれですか?

生存分析オブジェクト (suvfit/coxph) の AIC 値を計算したいと考えています。私がそれをしようとすると、次のように表示されます。

私が理解していることは、ソフトウェアの制限です。手でAIC値を計算せずにこの問題を解決するのを手伝ってくれる人はいますか?

トピックの基本を扱うのに問題があるようです。文献/サイト/例の推奨事項はありますか?

「survfit」の構文とその引数、および KM 曲線のプロットに特に問題がありました。

ありがとうリラン

R のサバイバル パッケージを使用して、生存率の中央値でログ ランク テストを実行しようとしていますが、それは可能ですか? 私はコードを使用しています：

ここで、年齢 = 診断時の年齢。Survival = 彼らが生き残った年数です。治療 = 治療を受けたかどうか (1 または 2 としてコード化)

これが正しい場合、なぜエラーが発生するのですか:

治療は、1 または 2 ではなく 0 または 1 としてコード化する必要がありますか?