Вопрос по r, plot – Как инвертировать ось Y на графике

12

Я хотел бы знать, как сделать график в R, где ось Y инвертирована так, что данные на графике появляются в том, что будет четвертый квадрант (IV) декартовой плоскости, в отличие от первого (I) квадранта.

Для справки, сюжет, который я пытаюсь сделать, выглядит очень похоже на следующее (источник):

a plot with inverted y-axis

Я нашел несколько вопросов в Интернете, касающихся измененияnumbering на оси Y, но все они по-прежнему отображают данные в первом квадранте. Кто-нибудь может подсказать, как мне подготовить сюжет, аналогичный показанному выше?

Я перехожу ко второму комментарию Тима: это не похоже на то, что вы действительно "переворачиваете" что-нибудь. Вы просто строите ось х сверху, а не снизу. joran
в своем сюжете звоните, укажитеaxes = FALSEи после использованияaxis() пару раз, указавside аргумент, чтобы получить ваши оси, где вы хотите их. tim riffe

Ваш Ответ

3   ответа
4

Отвечая на вопрос в заголовке, лучший / самый простой способ инвертировать ось - это перевернутьlimit переменные вокруг:

> plot(1:10, xlim=c(1,10));

standard plot

> plot(1:10, xlim=c(10,1));

plot with inverted X axis

> plot(1:10, ylim=c(10,1));

plot with inverted Y axis

Выполнение этого означает, что вам не нужно возиться с осями, которые отличаются от координат изображения.

Это может быть объединено с "xaxt =" n "; параметр и дополнительныйaxis Команда для размещения оси на другой стороне:

> plot(1:10, ylim=c(10,1), xaxt="n"); axis(3);

Inverted graph with X axis on top

Не работает с scatterplot3d.
10

Просто чтобы дать выработанный ответ, следуя комментариям @timriffe и @joran ...

Используйте функцию для небольших тиков журналаэтот ответ:

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 <- 1:8
y <- 10^(sort(runif(8, 1, 10), decreasing = TRUE))

Участок без осей:

plot(x, log10(y), # function to plot
xlab="",          # suppress x labels
type = 'l',       # specify line graph
xlim = c(min(x), (max(x)*1.3)),  # extend axis limits to give space for text annotation
ylim = c(0, max(log10(y))),      # ditto
axes = FALSE)    # suppress both axes

Добавьте причудливую ось журнала и поверните метки вправо вверх (спасибо @joran!):

minor.ticks.axis(2, 9, mn=0, mx=10, las=1)

Добавьте ось X вверх:

axis(3)

Добавьте метку оси X (спасибо за подсказку, @WojciechSobala)

mtext("x", side = 3, line = 2)

И добавить аннотацию в конец строки

text(max(x), min(log10(y)), "Example", pos = 1)

Вот результат:

enter image description here

С mtext вы можете добавить метку оси X, чтобы получить желаемый результат.
@WojciechSobala, спасибо за подсказку, я добавил это, и аннотацию в конце строки, чтобы более близко походить на пример OP.
+1 Отлично сделано. Вам не нужно загружать другое изображение, чтобы доказать это, но, поскольку ОП обязательно спросит, проходит лиlas = 1 вminor.ticks.axis повернуть метки оси Y, чтобы они были горизонтальными?
@joran, хорошая идея, готово!
0

Теперь довольно легко повернуть ось Y с помощьюscale_y_reverse и указатьposition = "top" для оси х вggplot2

Example

library(ggplot2)
library(scales)

set.seed(99)
Date <- seq(from = as.Date("2017-12-01"), to = as.Date("2017-12-15"),
             by = "days")
Flux <- runif(length(Date), 1, 10000)
Flux_df <- data.frame(Date, Flux)

p1 <- ggplot(Flux_df, aes(Date, Flux)) +
  geom_col() +
  xlab("") +
  scale_x_date(position = "top", breaks = pretty_breaks(), expand = c(0, 0)) +
  scale_y_reverse(expand = expand_scale(mult = c(0.2, 0))) +
  theme_bw(base_size = 16) +
  theme(panel.border = element_blank(),
        panel.grid.major.x = element_blank(),
        panel.grid.minor = element_blank(),
        axis.line = element_line()) +
  theme(legend.position = "none")
p1

enter image description here

Если нам нужна логарифмическая и обратная оси, нам нужен обходной путьВот какggplot2 не имеет этой опции банкомат

reverselog_trans <- function(base = exp(1)) { 
  trans <- function(x) -log(x, base) 
  inv <- function(x) base^(-x) 
  scales::trans_new(paste0("reverselog-", format(base)), trans, inv, 
                    scales::log_breaks(base = base), domain = c(1e-100, Inf)) 
}


p1 + scale_y_continuous(trans = reverselog_trans(10),
                        breaks = scales::trans_breaks("log10", function(x) 10^x),
                        labels = scales::trans_format("log10", scales::math_format(10^.x)),
                        expand = expand_scale(mult = c(0.2, 0))) +
  annotation_logticks()

enter image description here

Похожие вопросы