全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
6255 16
2011-08-14
file.xls
大小:(966 KB)

 马上下载


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

扫码加我 拉你入群

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

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

全部回复
2011-8-15 18:59:42

我把你画图所需要的资料

通通抓出来

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

######

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)
#######

   rcspline.plot.jpeg

二维码

扫码加我 拉你入群

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

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

2011-8-15 19:34:23
万分感谢高手的指点
二维码

扫码加我 拉你入群

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

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

2015-12-11 11:14:29
jsntxt 发表于 2011-8-15 19:34
万分感谢高手的指点
你好,能否把你的file.csv文件发到我的邮箱,我想看一看你的数据录入形式,因为也要做这个图形,但是我是R语言初学者,不太懂,所以想参考下你的,我的邮箱是sweettian2014@hotmail.com.谢谢高手了
二维码

扫码加我 拉你入群

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

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

2015-12-11 11:15:28
jsntxt 发表于 2011-8-15 19:34
万分感谢高手的指点
你好,能否把你的file.csv文件发到我的邮箱,我想看一看你的数据录入形式,因为也要做这个图形,但是我是R语言初学者,不太懂,所以想参考下你的,我的邮箱是sweettian2014@hotmail.com.谢谢高手了
二维码

扫码加我 拉你入群

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

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

2015-12-11 11:15:31
jsntxt 发表于 2011-8-15 19:34
万分感谢高手的指点
你好,能否把你的file.csv文件发到我的邮箱,我想看一看你的数据录入形式,因为也要做这个图形,但是我是R语言初学者,不太懂,所以想参考下你的,我的邮箱是sweettian2014@hotmail.com.谢谢高手了
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群