第二和第三个Kaplan Meier图超过设定的x限制

时间:2015-07-23 12:35:05

标签: r plot survival-analysis

enter image description here

从上图中可以看出,第二条和第三条Kaplan Meier曲线超出了我设置的x限制和框。

我对所有三个图使用了完全相同的代码,每次只更改实际的生存对象,以及y轴标签和范围。每个图的x限制和标签都相同。代码示例如下:

par(mfrow=c(2,2))

surv.obj <- survfit(Surv(Start, Stop, Event) ~ factor(Characteristic),
    data=dataset)

par( mai=c(2,3,1.5,0.42) )

colours <- c("gray0", "gray75")


# x-axis
max( unlist(unname(summary(surv.obj)["time"] )) ) # 1274 


plot(surv.obj, mark.time=F, col=colours, lty=1, lwd=2, yscale=100,
xscale=1,
fun="event",

ylab="",
xlab="Follow-Up (days)",
xlim=c(0, 1280),
xaxt="n",

axes=F,
cex.lab=1.8,
cex.main=1.8
)

box("plot", lty=1)

axis(2, cex.axis=1.5, at=c(0, 0.04, 0.08, 0.12, 0.16) )
axis(1, cex.axis=1.5, at=c(0, 300, 600, 900, 1200) )

对于接下来的两条曲线,再重复两次。

我知道这不是特定的情节&#39;问题的代码,因为当我单独绘制每个(并忽略par(mfrow = c(2,2)))时,它们按照我希望的方式生成(生存曲线截止约1280天) )。此外,我尝试使用par(mfrow = c(1,2))仅绘制第二和第三个绘图,第一个绘制得很好(这是图中绘图中的第二个),第二个绘图再次超过了框。

有谁知道为什么会发生这种情况以及如何解决这个问题?

由于 维基

0 个答案:

没有答案