plot errbar {Hmisc} x ось, когда x является фактором R

Question

Я написал этот код:

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

Он строит сюжет так:

Но моя цель - нарисовать Chr в xlab, а не в 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)

Надеюсь, это кому-нибудь поможет.