R中标记对数刻度显示

Question

在 R 中绘制直方图、散点图和其他轴缩放为对数刻度的图时，如何使用 10^-1 10^0 10^1 10^2 10^3 等标签而不是显示的轴只是 -1, 0, 1, 2, 3 等。 应该向 hist()、plot() 等命令添加哪些参数？

Answer 1

除了ggplot2的解决方案（参见gsk3的评论）之外，我想补充一点，当使用正确的参数时，这也会在plot()中自动发生，例如：

x <- 1:10
y <- exp(1:10)
plot(x,y,log="y")

您可以使用参数log="x" 表示 X 轴，或 log="xy" 表示两者。

如果您想格式化数字，或者您有日志格式的数据，您可以使用 axis() 进行解决。一些有趣的函数：

• axTicks(x) 为您提供 X 轴 (x=1) 或 Y 轴 (x=2) 上刻度的位置
• bquote() 将表达式转换为语言，但可以用变量的值替换变量。有关问题中的 bquote() 的更多信息 R 中绘图标签中的乳胶和变量？ .
• as.expression() 使来自 bquote() 的语言对象成为一个表达式。这允许 axis() 进行格式化，如 ?plotmath 中所述。它无法对语言对象执行此操作。

良好格式设置的示例：

x <- y <- 1:10
plot(x,y,yaxt="n")
aty <- axTicks(2)
labels <- sapply(aty,function(i)
            as.expression(bquote(10^ .(i)))
          )
axis(2,at=aty,labels=labels)

"https://i.sstatic.net/chb6E.png" alt="在此处输入图像描述">

Answer 2

以下是绘制此类轴的不同方法：

plot(NA, xlim=c(0,10), ylim=c(1, 10^4), xlab="x", ylab="y", log="y", yaxt="n")
at.y <- outer(1:9, 10^(0:4))
lab.y <- ifelse(log10(at.y) %% 1 == 0, at.y, NA)
axis(2, at=at.y, labels=lab.y, las=1)

编辑：这也在 latticeExtra 中通过 scale.components

Answer 3

在 ggplot2 中，您只需添加一个

... + 
scale_x_log10() + 
scale_y_log10(limits = c(1e-4,1), breaks=c(1e-4,1e-3,1e-2,0.1,1)) + ...

来缩放轴，为它们添加标签并添加自定义中断。