构造hist和boxplot的对数标度

Question

我需要你的帮助，因为我希望我的直方图和箱形图有一个对数刻度。 问题是我有很多观察集中在一个小间隔，我认为使用对数标度将解决我的问题。 例如，

x <- c(1:5000)
y <- c(rnorm(x),0)
z <- c(exp(y^2 +5*y),0)

par(mfrow=c(1,2))
hist(z, main="Histogramme des durées",
     ylab="Fréquences")
boxplot(z, main="Boîtes des durées")
par(mfrow=c(1,1))

分析我的数据非常困难，所以我尝试添加log =＆＃34; x＆＃34;但它不能做我想做的事！

提前谢谢。

Answer 1

您可以绘制数据的对数，而不是缩放绘图轴。

par(mfrow=c(1,2))
hist(log(z), main="Histogramme des durées",ylab="Fréquences")
boxplot(log(z), main="Boîtes des durées")

如果希望轴刻度为原始单位（不是对数单位），则可以更改轴刻度标签。

par(mfrow=c(1,2))
hist(log10(z), xaxt='n',xlab='',main="Histogramme des durées",ylab="Fréquences")
atx <- axTicks(1)
labels <- sapply(atx,function(i)
        as.expression(bquote(10^.(i) ))
       )
axis(1,at=atx,labels=labels)

boxplot(log10(z),yaxt='n', ylab='z',main="Boîtes des durées")
aty <- axTicks(2)
labels <- sapply(aty,function(i)
        as.expression(bquote(10^.(i) ))
      )
axis(2,at=aty,labels=labels)