Заполнение строк при использовании facet_grid и geom_raster в ggplot2

Question

Допустим, я провел четыре измерения / эксперимента (M1..M4), получая значения интенсивности, которые варьируются в 5 местах (loc_1...5). Я наблюдал различные классы элементов (n=7). Теперь я хотел бы суммировать результаты в одну диаграмму, используя facet_grid и geom_raster из ggplot2. Я заканчиваю диаграммой ниже (см. Изображение).

Теперь вопрос заключается в том, существует ли простое решение, чтобы заставить строки заполнить пространство в каждом фасете (т.е. удалить неиспользуемые строки в каждой панели).

Спасибо

rm(list=ls())
library(ggplot2)
library(reshape2)

set.seed(123)

# let's create a fake dataset
nb.mesure <- 4
n.row <- 200
n.col <- 5
nb.class <- 7


d <- matrix(round(runif(n.row * n.col),2), nc=n.col) 
colnames(d) <- sprintf("Loc_%02d", 1:5)

# These strings will be the row names of each heatmap
# in the subsequent facet plot
elements <- sample(replicate(n.row/2, 1:100))

# let's create a data.frame d
d <- data.frame(d, 
                mesure = sort(rep(c("M1","M2","M3", "M4"), n.row/4)),
                elements= elements,
                class=sample(nb.class, 
                             length(elements), 
                             rep=T,
                             prob = seq(0.01, 0.25, length.out=7))
                             )



# Data are melt
dm <- melt(d, id.var=c( "mesure", "elements", "class"))
colnames(dm) <- c("mesure","elements", "class", "pos", "intensity")

# Plotting
p <- ggplot(dm, aes(x = pos, y = elements, fill = intensity))
p <- p + geom_raster()                                        
p <- p + facet_grid(mesure~class , scales = "free", space="free_y")
p <- p + theme_bw()
p <- p + theme(text = element_text(size=8))
p <- p + theme(text = element_text(family = "mono", face = "bold"))
p <- p + theme(axis.text.y = element_blank(), 
               axis.ticks.y=element_blank(),
               axis.text.x = element_text(colour="grey20",
                                           size=6,angle=45,
                                        vjust = 0.3))
print(p)

Answer 1

Если ваша цель представить NA значения (или что-то еще) для всех отсутствующих комбинаций mesure, elements, class а также pos Вы можете использовать complete функция от tidyr пакет вроде так:

library(tidyr)

dm <- complete(dm, mesure, elements, class, pos, fill = list(intensity = NA))

ОБНОВИТЬ

Если вы хотите иметь неNA Значения расширяются, чтобы заполнить каждый аспект, от которого вы должны отойти facet_grid и переключиться на facet_wrap,

p <- p + facet_wrap(mesure~class , scales = "free_y", nrow = 4)