Построить непрерывные переменные в качестве цветового градиента на фоне ggplot2 - PullRequest
Новая
       40

Построить непрерывные переменные в качестве цветового градиента на фоне ggplot2

0 голосов
/ 24 марта 2020

У меня есть набор данных о температуре окружающей среды в море, который выглядит следующим образом:

Дата и время SW
01.11.2017 02:00:00 AM 17.8804
01.11.2017 04 : 00:00:00 17.8106
01.11.2017 06:00:00 AM 17.752
01.11.2017 08:00:00 17.7315
01.11.2017 10:00:00 AM 17.9066
01.11.2017 12:00:00 18.1229
11.11.2017 14:00:00 18.5551
11.11.2017 16:00:00 19.2719
11 / 1/2017 18:00:00 18.6102
01.11.2017 08:00:00 18.0809
01.11.2017 22:00:00 17.9975
02.11.2017 02 : 00:00 - 17.6566
02.11.2017, 04:00:00, 17.4234
, 02.11.2017, 06:00:00, 17.6084
, 02.11.2017, 08:00:00 17.5954 ...

Я бы хотел, чтобы температура была градиентом цвета на фоне другого графика ggplot2, показывающего скорости роста за тот же период времени, чтобы найти корреляции между температурой и скоростью роста двух переменных. Подобная проблема была поднята здесь: Отображение цвета фона в градиенте , но я не выяснил, как использовать непрерывные переменные для кодирования градиента.

Большое спасибо

1 Ответ

0 голосов
/ 24 марта 2020

То, что вы искали, было на самом деле ниже в ответах на пост, который вы связали. Это функция grad_by_val() из ответа jkd . 

library(ggplot2)

# Load your data
dat <- "Date, SW
11/1/2017 02:00:00 AM, 17.8804
11/1/2017 04:00:00 AM, 17.8106
11/1/2017 06:00:00 AM, 17.752
11/1/2017 08:00:00 AM, 17.7315
11/1/2017 10:00:00 AM, 17.9066
11/1/2017 12:00:00 PM, 18.1229
11/1/2017 02:00:00 PM, 18.5551
11/1/2017 04:00:00 PM, 19.2719
11/1/2017 06:00:00 PM, 18.6102
11/1/2017 08:00:00 PM, 18.0809
11/1/2017 10:00:00 PM, 17.9975
11/2/2017 02:00:00 AM, 17.6566
11/2/2017 04:00:00 AM, 17.4234
11/2/2017 06:00:00 AM, 17.6084
11/2/2017 08:00:00 AM, 17.5954"
dat <- read.table(text = dat, header = T, sep = ",", as.is = T)
# convert to date format
dat$Date <- as.POSIXct(dat$Date, format = "%m/%d/%Y %I:%M:%S %p")

# function from jkd's post
grad_by_val <- function(x, y, cols = blues9) {
  require(grid)
  y <- y[order(x)]
  ys <- (y - min(y)) / diff(range(y))
  cols <- colorRamp(cols)(ys) / 256
  colnames(cols) <- c("red", "green", "blue")
  cols <- apply(cols, 1, function(z) do.call(rgb, as.list(z)))
  mat <- matrix(cols, ncol = length(x))
  rasterGrob(
    image = mat,
    width = unit(1, "npc"),
    height = unit(1, "npc"),
    interpolate = TRUE
  )
}

По вашим данным это дает: 

ggplot(dat, aes(x = Date, y = SW)) +      
  annotation_custom(
    grob = grad_by_val(dat$Date, dat$SW, cols = heat.colors(10)),
    xmin = -Inf,
    xmax = Inf,
    ymin = -Inf,
    ymax = Inf) +
  geom_point() # make sure you plot stuff AFTER the gradient (or it'll be below it)

Вывод:

enter image description here

...