Legenda R dla wykresu rozrzutu gęstości koloru utworzonego za pomocą smoothScatter

Question

Generuję wykres rozproszenia gęstości kolorów w R za pomocą funkcji smoothScatter().

Przykład:

## A largish data set 
n <- 10000 
x1 <- matrix(rnorm(n), ncol = 2) 
x2 <- matrix(rnorm(n, mean = 3, sd = 1.5), ncol = 2) 
x <- rbind(x1, x2) 
oldpar <- par(mfrow = c(2, 2)) 
smoothScatter(x, nrpoints = 0) 

wyjściowa:

Problem mam jest to, że jestem pewien jak dodać skalę legenda/kolorów, który opisuje względną różnicę w kategoriach liczbowych między różnymi odcieniami. Na przykład, nie ma sposobu, aby stwierdzić, czy najciemniejszy niebieski na powyższym rysunku jest 2 razy, 10 razy lub 100 razy gęściej niż najjaśniejszy niebieski bez jakiejś legendy lub skali kolorów. Czy istnieje jakaś metoda w R, aby uzyskać wymagane informacje, aby stworzyć taką skalę, lub wszystko, co jest wbudowane, może automatycznie wygenerować skalę kolorów?

Answer 1

Oto odpowiedź, która opiera się na fields::imageplot a niektóre błahy z par(mar) uzyskać marże poprawne

fudgeit <- function(){ 
    xm <- get('xm', envir = parent.frame(1)) 
    ym <- get('ym', envir = parent.frame(1)) 
    z <- get('dens', envir = parent.frame(1)) 
    colramp <- get('colramp', parent.frame(1)) 
    fields::image.plot(xm,ym,z, col = colramp(256), legend.only = T, add =F) 
} 

par(mar = c(5,4,4,5) + .1) 
smoothScatter(x, nrpoints = 0, postPlotHook = fudgeit) 

Można pobawić się wokół z image.plot aby dostać to, czego chcesz, i spojrzeć na ?bkde2D i argument transformation do smoothScatter, aby uzyskać pojęcie o tym, co reprezentują kolory.