当x为因子R时绘制errbar {Hmisc} x轴

Question

我写了下面的代码：

with (data = puresa, 
      expr = errbar(Chr, Purity, Purity+sd, Purity-sd, add=F, pch=1, cap=.015))

它会像这样绘制一个图：

​

​

但我的目标是在xlab中绘制Chr，而不是像上面的图那样在ylab中绘制。问题是"Chr“是一个因素。因为当我把它打印成数字时，我得到了我的目标：

with (data = puresa,
      expr =  errbar(as.numeric(Chr), Purity, Purity+sd, Purity-sd, cap=.015))

​

​

最后一个图的问题是我不能将xlab值替换为Chr名称。

我正在尝试这个代码，但它不工作...

with (data = puresa, 
      expr =  errbar(as.numeric(Chr), Purity, Purity+sd, Purity-sd,
                     xlab=if(is.character(Chr) || is.character(Chr)) "" else (substitute(Chr))))

​

​

非常感谢您的帮助！

Answer 1

好了，我成功了！

par(las=2)
plot(as.numeric(puresa$Chr), puresa$Purity, type="n", xaxt = "n", ylim=c(0,1), main="Tumor purity", xlab="Chromosome", ylab="Purity")
errbar(as.numeric(puresa$Chr), puresa$Purity, puresa$Purity+puresa$sd, puresa$Purity-puresa$sd, add=T)
axis(1, at=1:25, labels=levels(puresa$Chr), cex.axis=0.8)

​

​

希望对某些人有帮助。