生存分析R应用

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

# Load the library.
library("survival")
# Print first few rows.
print(head(pbc))
# Load the library.
library("survival")
# Create the survival object.
survfit(Surv(pbc$time,pbc$status == 2)~1)
# Give the chart file a name.
png(file = "survival.png")
# Plot the graph.
plot(survfit(Surv(pbc$time,pbc$status == 2)~1))
# Save the file.
dev.off()
为什么我的plot出不来图，麻烦大家帮我看一下

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

你呼叫library时，不应该用双引号，下面的code就能画出来图了
library(survival)
plot(survfit(Surv(pbc$time,pbc$status == 2)~ 1))

1. #请去掉这一行
   
2. png(file = "survival.png")

复制代码

路过学习一下

谢谢分享！

Did it work?

hifinecon 发表于 2018-6-21 07:21
Did it work?

yes,谢谢

hifinecon 发表于 2018-6-21 07:21
Did it work?

谢谢你的帮助

alright!

原先的代码中将绘制的图片以png格式保存在了工作目录中，你可以打开你的工作目录找文件名为survival.png的就是你想要的图片了,如果你不知道自己当前的工作目录，可以运行   > getwd()   命令就可以得到
例外，如果你想在R窗口中直接显示图片，直接运行plot(survfit(Surv(pbcstatus == 2)~ 1))这一行命令就可以了，不用该命令的上下两条命令。
最后，推荐使用RStudio，R集成开发环境（IDE），操作更加方便，附下载地址：https://www.rstudio.com/