[英]Change Y-Axis to Percentage Points in plot()
我正在使用这个绘制图表:
### prepare survival object
jan.surv <- Surv(time=data1$start_year,time2=data1$time_since_conflict,data1$instrastate_war_recurrence==1)
### Binning the Governance.Score variable
data1$gov_score_binned <- ifelse(data1$Governance.Score<3,0,
ifelse(data1$Governance.Score>7,2,1))
### Run and plot Kaplan Meier
surv.kaplan1 <- survfit(jan.surv ~ gov_score_binned,data=data1)
plot(surv.kaplan1, col="blue", bg="lightblue", xlab = "Years since conflict", ylab = "Probability of Peace")
我得到这个 plot:或: https://docs.google.com/document/d/1kQDOeErpaS-5oA-hUtD0s3Xj7ig6entkRcGF8pfnGp8/edit?usp=sharing
如何将 Y 轴的标签从小数更改为百分比?
提前致谢: :)
数据: https://docs.google.com/spreadsheets/d/1S2288ASa0qzSsrR6BLhXDkBVAGBMH_qVY7a-IFh3Jz4/edit?usp=sharing
使用yaxt = 'n'
绘制 plot 时需要抑制 y 轴,然后在适当的点绘制另一个带有标签的 y 轴:
plot(surv.kaplan1, col = "blue", bg = "lightblue",
xlab = "Years since conflict",
ylab = "Probability of Peace", yaxt = 'n')
axis(2, at = 0:5 / 5, labels = paste0(0:5 * 20, '%'), las = 1)
使用的数据
显然,我们没有数据来重新创建这个 plot,所以我创建了一个具有相似值和相同名称的可重现 model 用于本示例的目的:
set.seed(1)
data1 <- data.frame(EP = c(rep(1, 10), rep(0, 4)),
time = c(rexp(10, 0.1), rep(60, 4)))
data1$Surv <- Surv(data1$time, event = data1$EP)
surv.kaplan1 <- survfit(Surv ~ 1, data = data1)
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.