svykmで失敗/累積発生率/イベント生存率プロットを作成する

Question

twangパッケージを使用して私のデータの傾向スコアを測定し、調査パッケージ機能を使用しています。私がオンラインに見てきた

km <- svykm(Surv(rec_time,rec)~procedure, design=design.mnps) 
plot(km) 

、それがsvykmlistオブジェクトでなかった場合は、次のように動作します：

私は、次のコードを使用して問題ない標準生存曲線を作成することができています

km <- svykm(Surv(rec_time,rec)~procedure, design=design.mnps, function(x) 1-x) 
plot(km) 

私は次のエラーを取得する：

1: In plot.window(...) : "FUN" is not a graphical parameter 

など

また、私はFUN = "イベント"を試しました。これは従来のサーベイオブジェクトではうまく動作しますが、svykmオブジェクトではうまく動作しません。私はまた、失敗した調査対象を操作しようとしました。

まだ、svykmlistオブジェクトで動作することはありませんでしたか？

おかげで、 マット

Answer 1

あなたがSOに投稿し、最小限の再現性の例を使用してください。オブジェクトdesign.mnpsは役に立ちません。あなたの質問を編集して、なぜこのような公式であなたの関数を反転させていないのかを明確にしてください。おかげ

例から？svykm

library(survey) 
data(pbc, package="survival") 
pbc$randomized <- with(pbc, !is.na(trt) & trt>0) 
biasmodel<-glm(randomized~age*edema,data=pbc) 
pbc$randprob<-fitted(biasmodel) 
dpbc<-svydesign(id=~1, prob=~randprob, strata=~edema, data=subset(pbc,randomized)) 

s2<-svykm(Surv(time,status>0)~I(bili>6), design=dpbc) 
plot(s2) 

sz<-svykm(Surv(time,status>0)~I(1 - I(bili>6)), design=dpbc) 
plot(sz)