Отображение незначительных логарифмических тиков по оси X в R

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

В пакете Hmisc есть функция minor.tick, но она плохо работает с логарифмическими шкалами.Я использую следующую функцию для получения небольших тиков, которые следуют за логарифмической шкалой.ax - это ось, на которой вы ее используете (так же, как и для функции axis), n - количество младших тиков (по умолчанию 9), t.ratio - соотношение между старшими и младшими тиками,а с помощью ... вы можете передать дополнительные параметры в axis

edit: Хорошая идея в комментариях, поэтому я отредактировал свою функцию.Есть два дополнительных параметра, 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))

Дает:

Answer 3

Попробуйте magaxis в упаковке magicaxis .

Answer 4

Метки оси x на первом графике на этой странице ошибочны. Незначительные тики распределены неправильно.

Answer 5

В ggplot2 мы можем использовать annotation_logticks вместе с scales::trans_breaks и scales::trans_format. Ниже приведен пример, взятый по ссылке выше.

    library(ggplot2)

    a <- ggplot(msleep, aes(bodywt, brainwt)) +
      geom_point(na.rm = TRUE) +
      scale_x_log10(
        breaks = scales::trans_breaks("log10", function(x) 10^x),
        labels = scales::trans_format("log10", scales::math_format(10^.x))
      ) +
      scale_y_log10(
        breaks = scales::trans_breaks("log10", function(x) 10^x),
        labels = scales::trans_format("log10", scales::math_format(10^.x))
      ) +
      theme_bw()

    a + annotation_logticks() # Default: log ticks on bottom and left

Answer 6

Произошла небольшая ошибка,

lims<-lims[3:4] отсутствует

   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 <- 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)
}

Answer 7

Используйте "" для надписей второстепенных тиков.