Вращение граней в ggplot2

Думаю, у меня сложный случай. Я рисую эволюционные уровни заболеваемости растений во времени, используя geom_raster: x и y - произвольные координаты поля, а z - уровень болезни, измеренный в несколько временных точек, и я хочу, чтобы каждая дата была построена в другом аспекте. Пока нет проблем. Ниже приведен пример набора данных и кода:

library(ggplot2)
data <- data_frame(month=factor(rep(c("march","april","may","june"), each=100), levels=c("march","april","may","june")),
              x=rep(rep(1:10, each=10), 4),
              y=rep(rep(1:10, 10), 4),
              z=c(rnorm(100, 0.5, 1), rnorm(100, 3, 1.5), rnorm(100, 6, 2), rnorm(100, 9, 1)))
ggplot(data, aes(x=x, y=y, fill=z)) +
  geom_raster(color="white") +
  scale_fill_gradient2(low="white", mid=mean(range(dat$z)), high="red") +
  scale_x_discrete(limit=1:10, expand = c(0, 0)) +
  scale_y_discrete(limit=1:10, expand = c(0, 0)) +
  coord_equal() +
  facet_wrap(~month)

Но мне бы хотелось, чтобы каждая грань вращалась под определенным углом (например, 15 °), чтобы отразить тот факт, что мое поле отлично не ориентировано на север (т.е. вершина не является северной, а дно - не Юг). Есть ли возможность в ggplot2 или любых связанных с сетью инструментах сделать это автоматически? Даже автоматический способ сглаживания отдельных граней к изображениям, их вращение и печать поворотных изображений на новой странице будет достаточной для моих нужд. Вот пример изображения, который я хотел бы получить (графы повернули на 15 ° в редакторе изображений): http://imgur.com/RYJ3EaR повернутые грани

Ответ 1

Здесь можно вращать грани независимо. Мы создаем список, содержащий отдельный поворотный график для каждого уровня month, а затем используйте grid.arrange для компоновки четырех графиков. Я также удалил легенду с отдельных сюжетов и построил легенду отдельно. Следующий код содержит вспомогательную функцию для извлечения легенды.

Я извлекаю объект легенды в глобальную среду в функции lapply ниже (не говоря уже о повторении извлечения несколько раз). Вероятно, лучший способ, но этот способ был быстрым.

library(gridExtra)

# Helper function to extract the legend from a ggplot
# Source: http://stackoverflow.com/info/12539348/ggplot-separate-legend-and-plot
g_legend<-function(a.gplot){
  tmp <- ggplot_gtable(ggplot_build(a.gplot))
  leg <- which(sapply(tmp$grobs, function(x) x$name) == "guide-box")
  legend <- tmp$grobs[[leg]]
  legend
}

# Create a list containing a rotated plot for each level of month
pl = lapply(unique(data$month), function(m) {

  # Create a plot for the current level of month
  p1 = ggplot(data[data$month==m,], aes(x=x, y=y, fill=z)) +
    geom_raster(color="white") +
    scale_fill_gradient2(low="white", high="red", 
                         limits=c(floor(min(data$z)), ceiling(max(data$z)))) +
    scale_x_discrete(limit=1:10, expand = c(0, 0)) +
    scale_y_discrete(limit=1:10, expand = c(0, 0)) +
    coord_equal() +
    facet_wrap(~month) 

  # Extract legend into global environment
  leg <<- g_legend(p1)

  # Remove legend from plot
  p1 = p1 + guides(fill=FALSE)

  # Return rotated plot
  editGrob(ggplotGrob(p1), vp=viewport(angle=-20, width=unit(0.85,"npc"), 
                                       height=unit(0.85,"npc")))                    
})

# Lay out the rotated plots and the legend and save to a png file
png("rotated.png", 1100, 1000)
grid.arrange(do.call(arrangeGrob, c(pl, ncol=2)),
             leg, ncol=2, widths=c(0.9,0.1))
dev.off()

введите описание изображения здесь