请教两个有关于rcspline.plot的问题，敬请指教！

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

file.xls (966 KB)

2011-8-14 18:57:14 上传

我希望利用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轴的名称，还请各位大侠指教！
万分感谢！！

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：Spline plot line Lin CSP library

我把你画图所需要的资料

通通抓出来

剩下的就由你自由发挥了.

######

library(Hmisc)
A<-read.csv("file.csv",header=TRUE)
attach(A)
out=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)

x=factor
xrange=c(0,3)
xe <- seq(xrange[1],xrange[2],length=600)
ylim=c(-1,2)
main <- "Estimated Spline Transformation"
xlim <- range(pretty(xe))
ylabl <- "log Relative Hazard"
xlab <- label(x)
lty=1
plot(out$x,out$xbeta,type="n",main=main,xlab=xlab,ylab=ylabl,xlim=xlim,ylim=ylim)
lines(xe,out$xbeta,lty=lty)
#######

  

万分感谢高手的指点

jsntxt 发表于 2011-8-15 19:34
万分感谢高手的指点

你好，能否把你的file.csv文件发到我的邮箱，我想看一看你的数据录入形式，因为也要做这个图形，但是我是R语言初学者，不太懂，所以想参考下你的，我的邮箱是sweettian2014@hotmail.com.谢谢高手了

jsntxt 发表于 2011-8-15 19:34
万分感谢高手的指点

你好，能否把你的file.csv文件发到我的邮箱，我想看一看你的数据录入形式，因为也要做这个图形，但是我是R语言初学者，不太懂，所以想参考下你的，我的邮箱是sweettian2014@hotmail.com.谢谢高手了

jsntxt 发表于 2011-8-15 19:34
万分感谢高手的指点

你好，能否把你的file.csv文件发到我的邮箱，我想看一看你的数据录入形式，因为也要做这个图形，但是我是R语言初学者，不太懂，所以想参考下你的，我的邮箱是sweettian2014@hotmail.com.谢谢高手了

你好，能否把你的file.csv文件发到我的邮箱，我想看一看你的数据录入形式，因为也要做这个图形，但是我是R语言初学者，不太懂，所以想参考下你的，我的邮箱是sweettian2014@hotmail.com.谢谢高手了

谢谢 已经成功下载了

您好~我是新人，能不能给我发一下您的数据的格式呀，我看不了。。。544426114@qq.com   谢谢您

你好，我想请教一下rcs如何取nkots，以及如何检验所取的nkots是否合适？我是用R语言做的。