Как получить 100 номеров для каждой ячейки в ggplot2 с geom_bin2d

Question

Прежде чем спросить, я прочитал это post, но мой более конкретный.

library(ggplot2) 
library(scales) 

set.seed(1) 
dat <- data.frame(x = rnorm(1000), y = rnorm(1000)) 

я заменяю свои реальные данные с dat, область х и у [-4,4] в этой случайной семени, и я разделить площадь на 256 (16 * 16) клеток, интервал который равен 0,5. Для каждой ячейки я хочу получить номера счетчиков.

Да, это довольно легко, geom_bin2d может решить.

# plot 
p <- ggplot(dat, aes(x = x, y = y)) + geom_bin2d() 

# Get data - this includes counts and x,y coordinates 
newdat <- ggplot_build(p)$data[[1]] 

# add in text labels 
p + geom_text(data=newdat, aes((xmin + xmax)/2, (ymin + ymax)/2, 
        label=count), col="white") 

До сих пор так хорошо, но я только хочу, чтобы получить топ 100 подсчета числа и сюжет в ПОС, как рис ниже.

После прочтения ?geom_bin2d, drop = TRUE удаляет только все ячейки с 0 подсчетов, и мое беспокойство топ-100 на счету. Что мне делать, это вопрос 1.

Просьба посмотреть еще на legend второго рисунка, количество отсчетов мало и близко, а если оно 10 000, 20 000, 30 000.

Метод используется trans в scale_fill_gradient, встроенная функция exp, log, sqrt и т. Д., Но я хочу разделить 1000. Затем я нашел trans_new() в пакете scales и попробовал, но отрицательный.

sci_trans <- function(){ trans_new('sci', function(x) x/1000, function(x) x*1000)} 
p + scale_fill_gradient(trans='sci') 

И это вопрос 2. Я гугл много, но не может найти способ решить эту проблему, спасибо много для тех, кто делает мне одолжение, спасибо!

Answer 1

Видимо Вы не можете получить выпускные лоты или количество от stat_bin2d или stat_summary_2d; в соответствии с соответствующим вопросом: How to use stat_bin2d() to compute counts labels in ggplot2?, где комментарий @MrFlick цитирует Хэдли с 2010 года: «он в основном говорит, что вы не можете использовать stat_bin2d, вам нужно будет сделать сводку самостоятельно».

Итак, обходной путь: создайте координатные ячейки вручную самостоятельно, получите значения 2D, затем возьмите верхний n. Например, с помощью dplyr:

dat %>% mutate(x_binned=some_fn(x), y_binned=some_fn(y)) %>% 
     group_by(x_binned,y_binned) %>% # maybe can skip this line 
     summarize(count = count()) %>% # NOTE: no need to sort() or order() 
     top_n(..., 100) 

Вы, возможно, придется тыкать в stat_bin2d, чтобы скопировать (или вызова) их точное координатно-биннинговые код. UPDATE: Вот source for stat-bin2d.r

StatBin2d <- ggproto("StatBin2d", Stat, 
    default_aes = aes(fill = ..count..), 
    required_aes = c("x", "y"), 

    compute_group = function(data, scales, binwidth = NULL, bins = 30, 
          breaks = NULL, origin = NULL, drop = TRUE) { 

    origin <- dual_param(origin, list(NULL, NULL)) 
    binwidth <- dual_param(binwidth, list(NULL, NULL)) 
    breaks <- dual_param(breaks, list(NULL, NULL)) 
    bins <- dual_param(bins, list(x = 30, y = 30)) 

    xbreaks <- bin2d_breaks(scales$x, breaks$x, origin$x, binwidth$x, bins$x) 
    ybreaks <- bin2d_breaks(scales$y, breaks$y, origin$y, binwidth$y, bins$y) 

    xbin <- cut(data$x, xbreaks, include.lowest = TRUE, labels = FALSE) 
    ybin <- cut(data$y, ybreaks, include.lowest = TRUE, labels = FALSE) 

    ... 

    } 

bin2d_breaks <- function(scale, breaks = NULL, origin = NULL, binwidth = NULL, 
         bins = 30, right = TRUE) { 
    ... 

(Но это кажется достойным повышения запрос на ggplot2, если она еще не была подана.)