2007年12月10日月曜日

データにエラー?

Multi-stateモデル解析用のRプログラムでは、個々の状態について

  1. その状態に推移した時間
  2. その状態から別の状態に移った時間
を引数に取り、1.<2. となることを事前要件としている。
一方、(自分で導出したので怪しい限りだが)ID

  • 3905
  • 4999
  • 6490
については1.=2.となっている。

追跡中止・治療中止の管理が甘いのがこのデータの欠点で、Time-to-event解析にとっては結構痛いエラーなのだが・・・

1 件のコメント:

Isao Tsumiyama さんのコメント...

事実を確認した所、「中止日」=「イベント日」というデータであった。

欲を言えばモニタリングで工夫してもらいたかった所だが、もう10年以上も前のことだし、現実的には十分起こりうるデータである。

イベント時間の同順位(調整法は数種類存在)ではなく、時間の起点と終点が同じという、生存時間分析ではまず想定しない状況である。そりゃsurvivalパッケージ作成者でも気が付かないだろう。