生存分析の検閲を考慮したプロットを取得

Question

（私は英語のため申し訳ありません）私は現在、行政の病気の葉の長さを研究しようとしていますので、長さ、日付などを入れたデータベースがあります。何千もの葉についての情報。私はその種の分析の初心者です。結局、私は、病気の休暇のx日目を過ぎて起こる病気の日の割合を知ることができるはずです。私は私の仕事を簡素化するために生存分析パッケージ（生存）を使用しようとしています。これまでのところ、私は私の質問に答えることができますが、私は打ち切りを考慮に入れません。私は、検閲されたデータについては0、そうでないデータに対しては1の列を持っています。何かの前に、私は打ち切りが訂正されたデータの生存プロットを取得しようとしています。これは、原理が正しく理解されていれば、プロットはすべての病欠が終わった予測されたものであることを意味します。ここで私が使用しているコードがありますが、私はRを正しい方法で使用していないと確信しています。私はcoxや他の何かを使うべきかどうかもわかりません。

私は葉の長さを持っていますが、event2013には検閲済みデータが0、そうでないデータが1です。すべてのデータを含むデータフレームはmydataと呼ばれます。

このコードを使って私が探しているプロットは何ですか？

s <- Surv(time2013,event2013) 
fKM <- survfit(s ~ 1,data=mydata) 
plot(fKM,conf.int=FALSE) 

Answer 1

まあ、それは....これに

1 ***おそらく私がしようとするだろう最初のものだが、そこ側面に非常に多くのより多くの事）あなたは何をしようとしていますか？それは純粋に探検ですか、あるいはいくつかの仮定/理論/仮説をテストしていますか？

2）検診フィールドに病気のない人が含まれているのか、それとも返されていない人だけが含まれていますか？データ期間の最終日の前に死亡した場合、病気休暇中にまだ生存している人が同じ方法で検閲されたと仮定すると、それはどのように影響を及ぼしますか？

3）その他の情報はありますか？また、公式の右側に共変量を入れることができますか。その人の性別と年齢は価値があるように見えるかもしれません。勤続年数、年功序列、以前の病気休暇期間の数はありますか？

私は私があなたの質問に答えるために、この記事を編集したかったが、私が行われた後、私は問題を実現:-)