R: назначить цвета значениям / палитре цветового градиента

У меня есть образец фрейма данных, который выглядит так:

reg1 <- structure(list(REGION = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L,  1L, 1L, 2L, 2L, 2L, 2L, 2L), .Label = c("REG1", "REG2"), class = "factor"),STARTYEAR = c(1959L, 1960L, 1961L, 1962L, 1963L, 1964L, 1965L, 1966L, 1967L, 1945L, 1946L, 1947L, 1948L, 1949L), ENDYEAR = c(1960L, 1961L, 1962L, 1963L, 1964L, 1965L, 1966L, 1967L, 1968L, 1946L, 1947L, 1948L, 1949L, 1950L), Y_START = c(0L, 0L, 0L, 0L,  0L, 0L, 0L, 0L, 0L, 2L, 2L, 2L, 2L, 2L), Y_END = c(1L, 1L,  1L, 1L, 1L, 1L, 1L, 1L, 1L, 3L, 3L, 3L, 3L, 3L), COLOR_VALUE = c(-969L, -712L, -574L, -312L, -12L, 1L, 0L, -782L, -999L, -100L, 23L, 45L, NA, 999L)), .Names = c("REGION", "STARTYEAR", "ENDYEAR", "Y_START", "Y_END", "COLOR_VALUE"), class = "data.frame", row.names = c(NA, -14L))

   REGION STARTYEAR ENDYEAR Y_START Y_END COLOR_VALUE
1    REG1      1959    1960       0     1        -969
2    REG1      1960    1961       0     1        -712
3    REG1      1961    1962       0     1        -574
4    REG1      1962    1963       0     1        -312
5    REG1      1963    1964       0     1         -12
6    REG1      1964    1965       0     1           1
7    REG1      1965    1966       0     1           0
8    REG1      1966    1967       0     1        -782
9    REG1      1967    1968       0     1        -999
10   REG2      1945    1946       2     3        -100
11   REG2      1946    1947       2     3          23
12   REG2      1947    1948       2     3          45
13   REG2      1948    1949       2     3          NA
14   REG2      1949    1950       2     3         999

Я создаю график с функцией rect(), которая отлично работает.

xx = unlist(reg1[, c(2, 3)])
yy = unlist(reg1[, c(4, 5)])
png(width=1679, height=1165, res=150)
if(any(xx < 1946)) {my_x_lim <- c(min(xx), 2014)} else {my_x_lim <- c(1946, 2014)} 
plot(xx, yy, type='n', xlim = my_x_lim)
apply(reg1, 1, function(y) 
  rect(y[2], y[4], y[3], y[5]))
dev.off()

В моих данных reg1 у меня есть 6-й столбец, который содержит значения от +1000 до -1000. Мне было интересно, есть ли метод, с помощью которого я мог бы раскрасить прямоугольники на моем графике в соответствии с моими значениями цвета. Низкие значения должны быть синими, значения около 0 должны приводить к белому, а высокие значения - к красному (если значение отсутствует или NA, тогда должен быть нанесен серый цвет).

Мой вопрос: Как я могу создать цветовую палитру со значениями от 1000 до -1000 (от красного над белым до синего) и применить ее к моему графику, чтобы каждый прямоугольник окрашивался в соответствии со значением цвета ?