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

No.15545　Rによるforest plot作成の件で　　【さとる】　2011/10/25(Tue) 00:35

いつもお世話になっております。
Rでのforest plot作成についてご質問させていただきます。
libraryからmetaforをインストールして，以下のプログラムでエクセルのCSVファイルにある下記のデータからforest plotを作成できました。

dat <- read.csv("forest plot.csv", header=TRUE)
library(metafor)
dat.escalc <- escalc(measure="OR", ai=a, n1i=n1, ci=c, n2i=n0, data=dat)
res.reml <- rma.uni(yi, vi, method="REML", data=dat.escalc)
forest(res.reml, atransf=exp, at=log(c(.05,.25,1,4,20)))

　a 　n1 　c 　n0 　d
　6 　16 　22 　36 　1
　20 　38 　8 　14 　2
　17 　31 　11 　21 　3
　13 　23 　4 　8 　4
　5 　9 　12 　22 　5

この場合は2×2の表におけるデータを用いたplotであり，各項目（5項目）におけるオッズ比の95％信頼区間が表示されます。

それで今回ご質問したいのは，下記のような2×2の表ではない時間的要素を含むデータを用いてハザード比のforest plotを作ろうとした場合，どのようなプログラムを用いればよいのでしょうか？簡単にですが，以下のデータであればCDDP（シスプラチン）， TXL（タキサン）の有る無しでPFS（無増悪生存期間）のハザード比のforest plotを作成しようとした場合，どのようにしたらよいのか具体的に教えて頂けたら本当に助かります。いろいろ調べましたが，どこにも記載がなく，こちらでご質問させていただいた次第です。よろしくお願いいたします。

PFS event 　CDDP　　　TXL　　
46 　　1 　　1　　　　　0
47 　　1 　　1 　　 1
118 　0 　　0 　　 0
5 　　1 　　1 　　0
61 　　1 　　0 　　 1
32 　　0 　　1 　　 1
207 　1 　　0 　　 0