悬赏 1 个论坛币 未解决
我希望利用Hmisc中rcspline.plot画出一条rcs曲线,数据置于附件中(请自行转换为CSV文件)
library(Hmisc)
A<-read.csv("file.csv",header=TRUE)
attach(A)
rcspline.plot(factor,Time,model="cox",xrange=c(0,3),ylim=c(-1,2),event=event,nk=4,knots=c(0.8, 1.0,1.5,
2.0),showknots=TRUE,plotcl=FALSE,statloc="none",subset=SEX=="2",lty=2)
问题1:显示的图形的纵轴为log relative Hazard的值,有什么办法可以使其变为relative Hazard的值,也就是希望取其以e为底的指数值。这个问题我在一个国外网站上得到这样的解答:You will have to learn some R programming to know how to change the code inside rcspline.plot. Better: get into the R rms。但究竟应该如何修改,还请各位大侠指教!
问题2 在Hmisc包的PDF文档对于rcspline.plot的说明中,是可以修改xlab和ylab的,但我发现怎么只能修改x轴名称而无法修改y轴的名称,还请各位大侠指教!
万分感谢!!
有人给出答案,但是最后两行命令没看懂,貌似不能运行,求大神指点!