dxy logo
首页丁香园病例库全部版块
搜索
登录

R语言绘制森林图报错,求指导!

发布于 2020-05-09 · 浏览 1346 · 来自 Android · IP 湖南湖南
这个帖子发布于 4 年零 359 天前,其中的信息可能已发生改变或有所发展。

rt=read.table("uniCox.xls",header=T,sep="\t",row.names=1,check.names=F)

data=as.matrix(rt)

HR=data[,1:3]

hr=sprintf("%.3f",HR[,"HR"])

hrLow=sprintf("%.3f",HR[,"HR.95L"])

hrHigh=sprintf("%.3f",HR[,"HR.95H"])

pVal=data[,"pvalue"]

pVal=ifelse(pVal<0.001, "<0.001", sprintf("%.3f", pVal))

tabletext <- 

 list(c(NA, rownames(HR)),

    append("pvalue", pVal),

    append("Hazard ratio",paste0(hr,"(",hrLow,"-",hrHigh,")")) )      

pdf(file="forest.pdf",

    width = 60,       

    height = 40,       

    )

forestplot(tabletext, 

      rbind(rep(NA, 3), HR),

      col=clrs,

      graphwidth=unit(200, "mm"),

      xlog=T,

      lwd.ci=2,

      boxsize=1.6,

      xlab="Hazard ratio"

      )

dev.off()

然后报错显示:

Error in valid.viewport(x, y, width, height, just, gp, clip, xscale, yscale, : 

 invalid 'xscale' in viewport

请问哪里出了错,需要怎么调整?谢谢!

最后编辑于 2022-10-09 · 浏览 1346

2 收藏2

全部讨论0

默认最新
avatar
2
分享帖子
share-weibo分享到微博
share-weibo分享到微信
认证
返回顶部