1

我正在使用生存包为两个群体生成 KM 曲线，但是当我尝试定义 xmax 参数时，我无法将 x 轴范围从 0-15 扩展到 0-50。

这是我的代码：

 library(sm)
library(survival)
comb = read.csv("S:/Novocure/ChemoSKTime2Resp.csv", head = TRUE)
comb
survobj = with(comb, Surv(SK.response.time,Censored))
survobj

# Compare the survival distributions of nonresponders(0) and responders(1)

fitresp = survfit(survobj ~ response, data=comb)
fitresp
summary(fitresp)

# plot the survival distributions by response

plot(fitresp, xlab="Tumor Progression Time in Months", 
  ylab="% Progressed", yscale=100, col=c("red","blue"),
  main="Chemo Time to Tumor Progression Distribution", xmax = 50) 
  legend("topright", title="Chi-sq = 5.2, P = 0.0222", 
  c("Non-responders", "Responders"), fill=c("red", "blue"))
4

1 回答 1

3

to的xmax参数plot.survfit似乎只缩短了 x 轴的范围。要扩展它，请改用参数xlim：

plot(fitresp, xlab="Tumor Progression Time in Months", 
  ylab="% Progressed", yscale=100, col=c("red","blue"),
  main="Chemo Time to Tumor Progression Distribution", xlim = c(0,50))
于 2012-11-01T17:12:46.607 回答