在R中的x轴上显示较小的对数刻度

Question

我有一个正态分布图和一个直方图,x轴以对数刻度显示0,10 ^ 0,10 ^ 1 ...我想在主要的之间包括小刻度.实际上,我能够使用前一个问题中给出的解决方案将主要的刻度格式从1,2,3等改为10 ^ 0,10 ^ 1,10 ^ 2,10 ^ 3 .我使用以下代码进行主要滴答:

major.ticks <- axTicks(1)
labels <- sapply(major.ticks,function(i)
            as.expression(bquote(10^ .(i)))
          )
axis(1,at=major.ticks,labels=labels)

这可以编辑为只标记次要刻度而不标记它们吗？

Answer 1

minor.tick包装中有一个功能Hmisc,但是这个功能与对数刻度很差.我使用以下函数来获得遵循对数标度的次要标记.ax是您使用它的轴(与函数相同axis),n是次要刻度的数量(默认为9),t.ratio是主刻度和次刻度之间的比率,并且...您可以将额外的参数传递给axis

编辑:评论中的好主意,所以我编辑了我的功能.有两个额外的参数,mn和mx用于最小和最大对数刻度(mn=0从而表示最小为10 ^ 0或1!)

功能:

minor.ticks.axis <- function(ax,n,t.ratio=0.5,mn,mx,...){

  lims <- par("usr")
  if(ax %in%c(1,3)) lims <- lims[1:2] else lims[3:4]

  major.ticks <- pretty(lims,n=5)
  if(missing(mn)) mn <- min(major.ticks)
  if(missing(mx)) mx <- max(major.ticks)

  major.ticks <- major.ticks[major.ticks >= mn & major.ticks <= mx]

  labels <- sapply(major.ticks,function(i)
            as.expression(bquote(10^ .(i)))
          )
  axis(ax,at=major.ticks,labels=labels,...)

  n <- n+2
  minors <- log10(pretty(10^major.ticks[1:2],n))-major.ticks[1]
  minors <- minors[-c(1,n)]

  minor.ticks = c(outer(minors,major.ticks,`+`))
  minor.ticks <- minor.ticks[minor.ticks > mn & minor.ticks < mx]


  axis(ax,at=minor.ticks,tcl=par("tcl")*t.ratio,labels=FALSE)
}

这可以如下应用:

x <- 10^(0:8)
y <- 1:9
plot(log10(x),y,xaxt="n",xlab="x",xlim=c(0,9))
minor.ticks.axis(1,9,mn=0,mx=8)

给:

Answer 2

这是一个简单的功能：

log10.axis <- function(side, at, ...) {
    at.minor <- log10(outer(1:9, 10^(min(at):max(at))))
    lab <- sapply(at, function(i) as.expression(bquote(10^ .(i))))
    axis(side=side, at=at.minor, labels=NA, tcl=par("tcl")*0.5, ...)
    axis(side=side, at=at, labels=lab, ...)
}

这是一个例子：

x <- exp(rnorm(200, 5))
hist(log(x), 20, xaxt="n", xlim=c(0, 8))
log10.axis(1, at=seq(0, 8, 2))

给出：