Étiquettes d'axe sur deux lignes avec variables x imbriquées (année sous mois)

Question

Je voudrais afficher les mois (sous forme abrégée) le long de l'axe horizontal, avec l'année correspondante imprimée une fois. Je sais afficher mois-année:

La répétition inutile de l'année encombre les étiquettes. Au lieu de cela, je voudrais quelque chose comme ça:

sauf que l'année serait imprimée sous les mois.

J'ai imprimé l'année au-dessus des étiquettes des axes, car c'est le mieux que je pouvais faire. Cela suit une limitation de la fonction annotate(), qui est tronquée si elle se trouve en dehors de la zone de traçage. Je suis conscient des solutions de contournement possibles basées sur annotate_custom(), mais je n'ai pas pu les faire fonctionner avec des objets de date (je n'ai pas essayé de convertir les dates en nombres et de nouveau en dates, car cela semblait plus compliqué que si tout va bien nécessaire)

Je me demande si la nouvelle dup_axis() pourrait être détournée à cet effet. Si au lieu d'envoyer l'axe dupliqué du côté opposé du panneau, il pourrait lui envoyer quelques lignes en dessous de l'axe dupliqué, alors il s'agirait peut-être simplement de configurer un axe avec panel.grid.major Effacé et les étiquettes définies sur %b, tandis que l'autre axe aurait panel.grid.minor effacé et les étiquettes définies sur %Y. (un défi supplémentaire est que les étiquettes de l'année seraient décalées en octobre au lieu de janvier)

Ces questions sont liées. Cependant, pour autant que je sache, les fonctions annotate_custom() et textGrob() ne fonctionnent pas bien avec les dates.

comment-puis-je-ajouter-des annotations-en-dessous-de-l'axe-x-dans-ggplot2

affichage du texte sous le tracé généré par ggplot2

Données et code de base ci-dessous:

 library("ggplot2") library("scales") ggplot(data = df, aes(x = Date, y = value)) + geom_line() + scale_x_date(date_breaks = "2 month", date_minor_breaks = "1 month", labels = date_format("%b %Y")) + xlab(NULL) ggplot(data = df, aes(x = Date, y = value)) + geom_line() + scale_x_date(date_minor_breaks = "2 month", labels = date_format("%b")) + annotate(geom = "text", x = as.Date("1719-10-01"), y = 0, label = "1719") + annotate(geom = "text", x = as.Date("1720-10-01"), y = 0, label = "1720") + xlab(NULL) # data df <- structure(list(Date = structure(c(-91455, -91454, -91453, -91452, -91451, -91450, -91448, -91447, -91446, -91445, -91444, -91443, -91441, -91440, -91439, -91438, -91437, -91436, -91434, -91433, -91431, -91430, -91429, -91427, -91426, -91425, -91424, -91423, -91422, -91420, -91419, -91418, -91417, -91416, -91415, -91413, -91412, -91411, -91410, -91409, -91408, -91406, -91405, -91404, -91403, -91402, -91401, -91399, -91398, -91397, -91396, -91395, -91394, -91392, -91391, -91390, -91389, -91388, -91387, -91385, -91384, -91382, -91381, -91380, -91379, -91377, -91376, -91375, -91374, -91373, -91372, -91371, -91370, -91369, -91368, -91367, -91366, -91364, -91363, -91362, -91361, -91360, -91359, -91357, -91356, -91355, -91354, -91353, -91352, -91350, -91349, -91348, -91347, -91346, -91345, -91343, -91342, -91341, -91340, -91339, -91338, -91336, -91335, -91334, -91333, -91332, -91331, -91329, -91328, -91327, -91326, -91325, -91324, -91322, -91321, -91320, -91319, -91315, -91314, -91313, -91312, -91311, -91310, -91308, -91307, -91306, -91305, -91304, -91303, -91301, -91300, -91299, -91298, -91297, -91296, -91294, -91293, -91292, -91291, -91290, -91289, -91287, -91286, -91285, -91284, -91283, -91282, -91280, -91279, -91278, -91277, -91276, -91275, -91273, -91272, -91271, -91270, -91269, -91268, -91266, -91265, -91264, -91263, -91262, -91261, -91259, -91258, -91257, -91256, -91255, -91254, -91252, -91251, -91250, -91249, -91248, -91247, -91245, -91244, -91243, -91242, -91241, -91240, -91238, -91237, -91236, -91235, -91234, -91233, -91231, -91230, -91229, -91228, -91227, -91226, -91224, -91223, -91222, -91221, -91220, -91219, -91217, -91216, -91215, -91214, -91213, -91212, -91210, -91209, -91208, -91207, -91205, -91201, -91200, -91199, -91198, -91196, -91195, -91194, -91193, -91192, -91191, -91189, -91188, -91187, -91186, -91185, -91184, -91182, -91181, -91180, -91179, -91178, -91177, -91175, -91174, -91173, -91172, -91171, -91170, -91168, -91167, -91166, -91165, -91164, -91163, -91161, -91160, -91159, -91158, -91157, -91156, -91154, -91153, -91152, -91151, -91150, -91149, -91147, -91146, -91145, -91144, -91143, -91142, -91140, -91139, -91138, -91131, -91130, -91129, -91128, -91126, -91125, -91124, -91123, -91122, -91121, -91119, -91118, -91117, -91116, -91115, -91114, -91112, -91111, -91110, -91109, -91108, -91107, -91104, -91103, -91102, -91101, -91100, -91099, -91097, -91096, -91095, -91094, -91093, -91091, -91090, -91089, -91088, -91087, -91086, -91084, -91083, -91082, -91081, -91080, -91079, -91077, -91076, -91075, -91074, -91073, -91072, -91070, -91069, -91068, -91065, -91063, -91062, -91061, -91060, -91059, -91058, -91056, -91055, -91054, -91053, -91052, -91051, -91049, -91048, -91047, -91046, -91045, -91044, -91042, -91041, -91040, -91039, -91038, -91037, -91035, -91034, -91033, -91032, -91031, -91030, -91028, -91027, -91026, -91025, -91024, -91023, -91021, -91020, -91019, -91018, -91017, -91016, -91014, -91013, -91012, -91011, -91010, -91009, -91007, -91006, -91005, -91004, -91003, -91002, -91000, -90999, -90998, -90997, -90996, -90995, -90993, -90992, -90991, -90990, -90989, -90988, -90986, -90985, -90984, -90983, -90982), class = "Date"), value = c(113, 113, 113, 113, 114, 114, 114, 115, 115, 115, 116, 116, 116, 116, 117, 117, 117, 117, 116, 117, 116, 116, 116, 117, 117, 117, 117, 117, 117, 117, 116, 117, 116, 116, 116, 117, 117, 117, 117, 117, 117, 117, 116, 116, 117, 117, 117, 117, 117, 117, 117, 117, 117, 117, 117, 118, 118, 118, 118, 117, 118, 117, 117, 117, 117, 117, 117, 118, 116, 116, 116, 116, 116, 116, 116, 117, 117, 118, 118, 118, 118, 118, 119, 120, 120, 119, 119, 120, 120, 121, 121, 122, 124, 124, 122, 123, 124, 123, 123, 123, 123, 123, 124, 124, 126, 126, 126, 126, 126, 125, 125, 126, 127, 126, 126, 125, 126, 126, 126, 128, 128, 128, 130, 133, 131, 133, 134, 134, 134, 136, 136, 136, 135, 135, 135, 136, 136, 136, 136, 135, 135, 135, 135, 130, 129, 129, 130, 131, 136, 138, 155, 157, 161, 170, 174, 168, 165, 169, 171, 181, 184, 182, 179, 181, 179, 175, 177, 177, 174, 170, 174, 173, 178, 173, 178, 179, 182, 184, 184, 180, 181, 182, 182, 184, 184, 188, 195, 198, 220, 255, 275, 350, 310, 315, 320, 320, 316, 300, 310, 310, 320, 317, 313, 312, 310, 297, 285, 285, 286, 288, 315, 328, 338, 344, 345, 352, 352, 342, 335, 343, 340, 342, 339, 337, 336, 336, 342, 347, 352, 352, 351, 352, 352, 351, 352, 352, 355, 375, 400, 452, 487, 476, 475, 473, 485, 500, 530, 595, 720, 720, 770, 750, 770, 750, 735, 740, 745, 735, 700, 700, 750, 760, 755, 755, 760, 760, 765, 950, 950, 950, 875, 875, 875, 880, 880, 880, 900, 900, 900, 880, 880, 890, 895, 890, 880, 870, 870, 870, 870, 870, 860, 860, 860, 860, 850, 840, 810, 820, 810, 810, 805, 810, 805, 820, 815, 820, 805, 790, 800, 780, 760, 765, 750, 740, 820, 810, 800, 800, 775, 750, 810, 750, 740, 700, 705, 660, 630, 640, 595, 590, 570, 565, 535, 440, 400, 410, 400, 405, 390, 370, 300, 300, 180, 200, 310, 290, 260, 260, 275, 260, 270, 265, 255, 250, 210, 210, 200, 195, 210, 215, 240, 240, 220, 220, 220, 220, 210, 212, 208, 220, 210, 212, 208, 220, 215, 220, 214, 214, 213, 212, 210, 210, 195, 195, 160, 160, 175, 205, 210, 208, 197, 181, 185)), .Names = c("Date", "value"), row.names = c(NA, 393L ), class = "data.frame")

eipi10 · Accepted Answer

Le code ci-dessous fournit deux options possibles pour ajouter des étiquettes d'année.

Option 1a: facettage

Vous pouvez utiliser des facettes pour marquer les années. Par exemple:

library(ggplot2) library(lubridate) ggplot(df, aes(Date, value)) + geom_line() + scale_x_date(date_labels="%b", date_breaks="month", expand=c(0,0)) + facet_grid(~ year(Date), space="free_x", scales="free_x", switch="x") + theme_bw() + theme(strip.placement = "outside", strip.background = element_rect(fill=NA,colour="grey50"), panel.spacing=unit(0,"cm"))

Notez qu'avec cette approche, s'il y a des dates manquantes au début ou à la fin d'une année (par "manquant", je veux dire que les lignes pour ces dates ne sont même pas présentes dans les données) alors l'axe x commencera/se terminera au première/dernière date dans les données de cette année, plutôt que de passer du 1er janvier au 31 décembre. Dans ce cas, vous devez ajouter des lignes pour les dates manquantes et NA pour value ou interpoler value. De plus, avec cette méthode, il n'y a ni espace ni ligne entre le 31 décembre d'une année et le 1er janvier de l'année suivante, il y a donc une discontinuité à chaque année.

Option 1b: facettes + étiquettes de mois centrées

Pour répondre au commentaire de @ AF7. Vous pouvez centrer les étiquettes du mois en ajoutant des espaces avant chaque étiquette. Mais vous devez choisir le nombre d'espaces manuellement, en fonction de la taille physique du tracé lorsque vous l'imprimez sur un appareil. (Il y a probablement un moyen de centrer les étiquettes par programme sur la base des mesures internes de grob, mais je ne sais pas comment le faire.) J'ai également supprimé les lignes de grille verticales mineures et allégé la ligne entre les années.

ggplot(df, aes(Date, value)) + geom_line() + scale_x_date(date_labels=paste(c(rep(" ",11), "%b"), collapse=""), date_breaks="month", expand=c(0,0)) + facet_grid(~ year(Date), space="free_x", scales="free_x", switch="x") + theme_bw() + theme(strip.placement = "outside", strip.background = element_blank(), panel.grid.minor.x = element_blank(), panel.border = element_rect(colour="grey70"), panel.spacing=unit(0,"cm"))

Option 2a: modifier le groupe d'étiquettes de l'axe X

Voici une méthode plus complexe et capricieuse (bien qu'elle puisse probablement être automatisée par quelqu'un qui comprend mieux que moi la structure et les espacements unitaires des graphiques en grille) qui évite les pièges de la méthode de facettage décrite ci-dessus:

library(grid) # Fake data with an extra year added for illustration set.seed(2) df = data.frame(Date=seq(as.Date("1718-03-01"),as.Date("1721-09-20"), by="1 day")) df$value = cumsum(rnorm(nrow(df))) # The plot we'll start with p = ggplot(df, aes(Date, value)) + geom_vline(xintercept=as.numeric(df$Date[yday(df$Date)==1]), colour="grey60") + geom_line() + scale_x_date(date_labels="%b", date_breaks="month", expand=c(0,0)) + theme_bw() + theme(panel.grid.minor.x = element_blank()) + labs(x="")

Maintenant, nous voulons ajouter les valeurs de l'année ci-dessous et entre juin et juillet de chaque année. Le code ci-dessous fait cela en modifiant le groupe d'étiquettes de l'axe x et est adapté de this SO answer by @SandyMuspratt.

# Get the grob g <- ggplotGrob(p) # Get the y axis index <- which(g$layout$name == "axis-b") # Which grob xaxis <- g$grobs[[index]] # Get the ticks (labels and marks) ticks <- xaxis$children[[2]] # Get the labels ticksB <- ticks$grobs[[2]] # Edit x-axis label grob # Find every index of Jun in the x-axis labels and add a newline and # then a year label junes = which(ticksB$children[[1]]$label == "Jun") ticksB$children[[1]]$label[junes] = paste0(ticksB$children[[1]]$label[junes], "
 ", unique(year(df$Date))) # Put the edited labels back into the plot ticks$grobs[[2]] <- ticksB xaxis$children[[2]] <- ticks g$grobs[[index]] <- xaxis # Draw the plot grid.newpage() grid.draw(g)

Option 2b: modifier le groupe d'étiquettes de l'axe des x et centrer les étiquettes du mois

Vous trouverez ci-dessous la seule modification à apporter à l'option 2a pour centrer les étiquettes de mois, mais, encore une fois, le nombre d'espaces doit être modifié manuellement.

# Make the edit # Center the month labels between ticks ticksB$children[[1]]$label = paste0(paste(rep(" ",7),collapse=""), ticksB$children[[1]]$label) # Find every index of Jun in the x-axis labels and a year label junes = grep("Jun", ticksB$children[[1]]$label) ticksB$children[[1]]$label[junes] = paste0(ticksB$children[[1]]$label[junes], "
 ", unique(year(df$Date)))

Mike H. · Answer

Si vous voulez essayer de pirater ensemble un sous-label, vous pouvez le convertir en grob. J'ai édité cela à partir du message d'origine pour créer une fonction qui ajoute les sous-étiquettes et renvoie un objet gtable. Notez que l'entrée sublabs doit être de la même longueur que vos ruptures d'axe x:

library(grid) library(gtable) library(gridExtra) add_sublabs <- function(plot, sublabs){ gg <- ggplotGrob(plot) axis_num <- which(gg$layout[,"name"] == "axis-b") xbreaks <- gg[["grobs"]][[axis_num]][["children"]][[2]][["grobs"]][[2]][["children"]][[1]]$x if(length(xbreaks) != length(sublabs)) stop("Sub-labels must be the same length as the x-axis breaks") to_breaks <- c(as.numeric(xbreaks),1)[which(!duplicated(sublabs, fromLast = TRUE))+1] sublabs_x <- diff(c(0,to_breaks)) sublabs_labels <- sublabs[!duplicated(sublabs, fromLast = TRUE)] tg <- tableGrob(matrix(sublabs_labels, nrow = 1)) tg$widths = unit(sublabs_x, attr(xbreaks,"unit")) pos <- gg$layout[axis_num,c("t","l")] gg2 <- gtable_add_rows(gg, heights = sum(tg$heights)+unit(4,"mm"), pos = pos$t) gg3 <- gtable_add_grob(gg2, tg, t = pos$t+1, l = pos$l) return(gg3) } #Plot and sublabels p <- ggplot(data = df, aes(x = Date, y = value)) + geom_line() + scale_x_date(date_breaks = "2 month", date_minor_breaks = "1 month", labels = date_format("%b")) + xlab(NULL) sublabs <- c(rep("1719",2),rep("1720",6)) #Draw grid.draw(add_sublabs(p, sublabs))

pasipasi · Answer

Je suis tombé sur cette question et j'ai pensé que je pourrais peut-être ajouter une solution. Nous pouvons afficher le mois et l'année dans le premier mois affiché de chaque année en utilisant une condition simple. Vous pouvez jouer avec le date_breaks Pour supprimer janvier des étiquettes, et cela fonctionnera toujours. J'utilise month() et year() de lubridate.

library(tidyverse) library(lubridate) df %>% ggplot(aes(Date, value)) + geom_line() + scale_x_date(date_breaks = "2 months", labels = function(x) if_else(is.na(lag(x)) | !year(lag(x)) == year(x), paste(month(x, label = TRUE), "
", year(x)), paste(month(x, label = TRUE))))

G. Grothendieck · Answer

Une façon d'éviter les complexités serait de modifier la production requise afin que janvier soit remplacé par l'année.

La fonction lab renvoie les étiquettes en fonction des ruptures. De manière inattendue, ggplot lui transmettra des NA, donc dans la première ligne du corps de la fonction, nous remplaçons ceux par une date - peu importe la date puisque ces valeurs ne sont pas utilisées par la suite par ggplot. Enfin, nous formaterons la date en année ou en mois abrégé selon que le mois est janvier (ce qui correspond au composant POSIXlt mon égal à 0) ou non.

library(ggplot2) library(scales) lab <- function(b) { b[is.na(b)] <- Sys.Date() format(b, ifelse(as.POSIXlt(b)$mon == 0, "%Y", "%b")) } ggplot(df, aes(Date, value)) + geom_line() + scale_x_date(date_breaks = "month", labels = lab)

Remarque: J'ai ajouté problème 2182 à la liste des problèmes de github ggplot2 concernant les NA qui sont passées à la fonction d'étiquette. Si les versions suivantes de ggplot2 ne passent plus les NA, la première ligne du corps de lab pourrait être omise.

Mise à jour: corrigé.