Rを使用してカプランマイヤー曲線を作成したいです。

Question

前提・実現したいこと

R　を使用しています。
カプランマイヤー曲線が表示されません。
何がエラーなのかご教示いただきたいです。

発生している問題・エラーメッセージ

#####(Logrank検定)#####

library(survival, pos=18)

km <- NULL

km.summary.table <- NULL

km <- survfit(Surv((PFS/1),PFS.イベント１.打ち切り０==1)~GnP.FFX.0.S.1.GEM.1, data=Dataset,

• na.action = na.omit, conf.type="log-log")

summary(km)

Length Class Mode
0 NULL NULL

len <- nchar("GnP.FFX.0.S.1.GEM.1")

legend <- substring(names(km$strata), len+2)

quartz(width=7, height=7); par(lwd=1, las=1, family="sans", cex=1, mgp=c(2.5,1,0))

mar <- par("mar")

mar[1] <- mar[1] + length(km$strata) + 0.5

mar[2] <- mar[2] + 2

par(mar=mar)

opar <- par(mar = mar)

on.exit(par(opar))

plot(km, bty="l", col=1:32, lty=1, lwd=1, conf.int=FALSE, mark.time=TRUE, xlab="PFS",

• ylab="Probability")

xticks <- axTicks(1)

n.atrisk <- nrisk(km, xticks)

for (i in 1:length(km$strata)){axis(1, at = xticks, labels = n.atrisk[i,], line=3+i, tick

• = FALSE)}

for (i in 1:length(km$strata)){mtext(legend[i], at=-(xticks[2]-xticks[1])/2, side=1,

• line=4+i, cex=1)}

title(xlab = "Number at risk", line = 3.5, adj = 0)

legend ("topright", legend, col=1:32, lty=1, lwd=1, box.lty=0,

• title="GnP.FFX.0.S.1.GEM.1")

res <- NULL

(res <- survdiff(Surv(PFS,PFS.イベント１.打ち切り０==1)~GnP.FFX.0.S.1.GEM.1, data=Dataset, rho=0,

• na.action = na.omit))

km.summary.table <- summary.km(survfit=km, survdiff=res)

km.summary.table

Comments

[KojiDoi]

グラフ画像出力用の作図デバイスの設定はどうなっているのでしょうか。