「統計学関連なんでもあり」の過去ログ--- 042

No.10678　Kaplan-Meier推定量　　【初心者子】　2009/08/22(Sat) 10:02

統計を始めたばかりで初歩的な質問ですが
よろしくお願いいたします。

現在，コマンドがまだわからず
Rコマンダーで解析しております。
メインの暴露（2群）と環境（2群），性別による生存解析を行いたいです。
そこで
1，Kaplan-Meier推定量はallで求めるのか
　　この時点で比較したい群（男性環境A暴露と男性環境A非暴露）を抽出して求め，
　　順にやっていく（男性環境B暴露と男性環境B非暴露）のでしょうか？
　　（上記の内容ですと8群になりますが，層別因子選択のところは2群の設定のようなので）

2，もし抽出するなら条件をかえて何度も抽出して解析していいのか？

3，8群一緒に図に載せたいが可能か

ネットをいろいろ見てみましたが
因子の数が多いのが見当たりませんでした。
ご教授お願いいたします。

No.10679　Re: Kaplan-Meier推定量　　【青木繁伸】　2009/08/22(Sat) 13:45

比較する必要がある群分けで何枚も図を描かなければならないこともあるでしょう。

Rコマンダーはわかりませんが，以下のようにすれば，複数の群分けの生存曲線を一枚の図に描くことはできます。見やすいかどうかは保証できませんが。
set.seed(123)
n <- 100
d <- data.frame(time=round(rexp(n, 0.05)),
exposure=sample(c("-", "+"), n, replace=TRUE),
environment=sample(LETTERS[1:2], n, replace=TRUE),
sex=sample(c("M", "F"), n, replace=TRUE),
status=sample(0:1, n, replace=TRUE, prob=c(7, 3)))
head(d)
d2 <- split(d, list(d$exposure, d$environment, d$sex))
plot(c(0,130), c(0,1), type="n")
mapply(function(x, no) lines(survfit(Surv(time, status)~1, data=x), col=no), d2, 1:8)

No.10682　Re: Kaplan-Meier推定量　　【初心者子】　2009/08/23(Sun) 09:47

ありがとうございます。
やってみます。

N数が少ないわけではないですが
Kaplan-Meierと生命表はそんなに違いがないということですし
Kaplan-Meierでやってみます。