2011-01-18 4 views
1

ワンサンプルテストにsurvdiff関数を使用したいと思います。survdiff one-sampleテスト - 式でオフセット項を使用するにはどうすればよいですか?

私は、オフセット用語をどこに含めるべきか理解していませんでした。

surdiffは次のようになります。EXについて

survdiff(formula, data, subset, na.action, rho=0) 

を、K-試料試験において、式は次のようになります。Surv(futime, fustat)と予測因子である:rx。そして、我々は次のようにsurdiffを使用します。

survdiff(Surv(futime, fustat) ~ rx,data=ovarian) 

(この例は、Rの例から取られている。)

私は理解していないですが、私は予測因子でオフセット項を含めるべきか1つのサンプルのテストを実行したい場合は、 Rの説明では、spが確率ベクトルであるところでoff​​set(sp)を使うべきだと言う。

誰にでも正しい構文が何であるか教えてもらえますか?それは「インクルード」ですか?

ありがとうございます! 路地

答えて

1

あなたは、モデル式のRHS上のみものとしてoffset(var)が含まれ、var生存確率、後者はsurvexp機能を使用して計算することが可能でなければなりません。あなたは?survdiffでさらに上のビットを読んだ場合

、例まで、あなたはこの見ていると思います:あなたを示し

## Expected survival for heart transplant patients based on 
## US mortality tables 
expect <- survexp(futime ~ ratetable(age=(accept.dt - birth.dt), 
        sex = 1, year = accept.dt, race = "white"), jasa, 
        cohort = FALSE, ratetable = survexp.usr) 

## actual survival is much worse (no surprise) 
survdiff(Surv(jasa$futime, jasa$fustat) ~ offset(expect)) 

正確で単一offset用語を含める方法式:survdiff

+0

本当にありがとう、私はその元を逃した。私は1サンプルのテストを見たいと思っていました。 については、私はそれを使用しようとしましたが、なぜ年齢の値を取得する必要があるのか​​分かりませんでした。このパラメータはrmapに属しています。これはオプションのパラメータです。さらに、私の実験では、私はこの「年齢」を全く使用しません。 「年齢」パラメータを持たずにsurvexpベクトルを使用する方法をアドバイスできれば、非常に良いでしょう。あなたの助けをもう一度ありがとう! Alley – user573347

+0

ニースの返信、Gavin。ここでは、ワンサンプル法(TherneauがMayoのウェブサイトで予想される生存予想のすべてのリソースに加えて)の良い引用があります。http://jnci.oxfordjournals.org/content/95/19/1434完全 –

+0

これは素晴らしいです、私は間違いなくDWinが必要です、ありがとう! – user573347

関連する問題