从上图中可以看出,第二条和第三条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))仅绘制第二和第三个绘图,第一个绘制得很好(这是图中绘图中的第二个),第二个绘图再次超过了框。
有谁知道为什么会发生这种情况以及如何解决这个问题?
由于 维基