如何从R图导出X和Y图值？

Question

使用R中的eha和生存包，我使用以下代码创建以下生存图：

with(dataset, plot(Surv(enter, exit, event), ylim= c(0.87, 1.0), fn = "surv", strata = Gender))

除了导出图的图像，我如何导出该图上每条线的X和Y值？ 我试图使用以下代码，但是它只会创建一个空白的.csv文件：

write.table((with(dataset, plot(Surv(enter, exit, event), ylim= c(0.87, 1.0), fn = "surv", strata = Gender))), file="testfit.csv", sep=',')

感谢您的任何帮助。

在第二点上，有没有人碰巧知道如何使该图在ggplot2中看起来更好？

Answer 1

您可以使用gridGraphics包将绘图转换为grob，然后从grob对象中获取点。 我喜欢此解决方案，因为它适用于各种问题。

使用此处创建的图形作为示例图形：

# plot a graph
library(survival)
library(grid)
library(gridGraphics)
data(lung)
lung.surv <- survfit(Surv(time,status) ~ 1, data = lung)
plot(lung.surv)

# capture the plotted output as a grob
grid.echo()
grid.grab() -> k

# pull out the data from the grob..
k$children$`graphics-plot-1-points-1`$x -> x
k$children$`graphics-plot-1-points-1`$y -> y

我注意到您的情节按性别分层...请确保为每种性别提取正确的点子集。