Как показать иерархическую структуру на метках осей с помощью ggplot2?

Я хотел бы создать такой график с помощью ggplot2:

(Я допустил небольшую ошибку, из-за которой отображались не все ярлыки. Очевидно, я хочу, чтобы отображались все ярлыки :)) Однако мой поиск в Google ничего не дал. Самое близкое, что я получил, было с barNest от plotrix, однако мне не нравится качество получаемых графиков. Те, что ggplot2 выглядят намного лучше.

В идеальном случае я ожидал, что смогу сделать что-то вроде этого: scale_x_discrete ("Название оси", метки = Level1 + Level2 + Level3) И тогда Level1, Level2 и Level3 были бы столбцами фрейма данных, который использовался в функции qplot. Если я смогу заставить это работать, я хочу создать функцию, чтобы добавить этот тип функциональности к функции весов. Может, scale_x_hierarchical было бы лучше. Подобные данные могут быть сгенерированы с помощью следующего кода:

data <- data.frame(
  Subject = c(rep(1, 48), rep(2, 48), rep(3, 48), rep(4, 48)), 
  month = rep(c(rep(1, 12), rep(4, 12), rep(7, 12), rep(10, 12)), 4), 
  day = rep(c(rep(1, 4), rep(11, 4), rep(21, 4)), 16), 
  hour = rep(c(0, 6, 12, 18), 48), 
  data = rnorm(192))

Между прочим, изображение, используемое здесь, создано с помощью Spotfire.