使用R在直方图中进行y轴对数

Question

嗨,我正在使用R进行直方图,但Y轴的数量是如此之大,我需要将其转换为对数.请参阅我的脚本下方:

hplot<-read.table("libl")
hplot
pdf("first_end")
hist(hplot$V1, breaks=24, xlim=c(0,250000000), ylim=c(0,2000000),main="first end mapping", xlab="Coordinates")
dev.off()

那么我应该如何更改脚本？谢谢

Answer 1

在绘图之前,您可以保存直方图数据以进行调整:

set.seed(12345)
x = rnorm(1000)

hist.data = hist(x, plot=F)
hist.data$counts = log10(hist.data$counts)

dev.new(width=4, height=4)
hist(x)

dev.new(width=4, height=4)
plot(hist.data, ylab='log10(Frequency)')