+
+
+
+
+
+
+
+
+
+
+Autour du SARS-CoV-2 (Covid-19)
+Elia Pérennès
+ +
+
+
+
+
+
+
+Le but est ici de reproduire des graphes semblables à ceux du South China Morning Post (SCMP), sur la page The Coronavirus Pandemic et qui montrent pour différents pays le nombre cumulé (c’est-à-dire le nombre total de cas depuis le début de l’épidémie) de personnes atteintes de la maladie à coronavirus 2019. Les données que nous utiliserons sont compilées par le Johns Hopkins University Center for Systems Science and Engineering (JHU CSSE).
+Tout d’abord, on nettoie l’espace de travail sur R et on charge les librairies nécessaires à l’exécution de notre code :
+ + + + + + + +
+
+Etude du nombre de cas
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+Préparation des données
+Nous téléchargeons les données recensant le nombre de cas de Covid-19 disponibles sur le Github du JHU CSSE.
+ + + +data_url = "https://raw.githubusercontent.com/CSSEGISandData/COVID-19/master/csse_covid_19_data/csse_covid_19_time_series/time_series_covid19_confirmed_global.csv"
+data = read.csv(data_url, stringsAsFactors=FALSE)Regardons nos données :
+ + + + + + +
+
+
+
+
+
+On retraite dans un premier temps les noms des colonnes de dates en enlevant le X devant chaque nom de colonne :
+ + + + + + + [1] "Province.State" "Country.Region" "Lat" "Long" "1.22.20"
+ [6] "1.23.20" "1.24.20" "1.25.20" "1.26.20" "1.27.20" On sépare les données sur Hong-Kong de celles du reste de la Chine :
+ + + + + + + +On filtre la base de données pour ne garder que les pays suivants : la Belgique (Belgium), la Chine - toutes les provinces sauf Hong-Kong (China), Hong Kong (China, Hong-Kong), la France métropolitaine (France), l’Allemagne (Germany), l’Iran (Iran), l’Italie (Italy), le Japon (Japan), la Corée du Sud (Korea, South), la Hollande sans les colonies (Netherlands), le Portugal (Portugal), l’Espagne (Spain), le Royaume-Unis sans les colonies (United Kingdom) ainsi que les États-Unis (US).
+ + + +data = data %>% filter(Country.Region %in% c("Belgium","China",
+ "Hong Kong","France","Germany","Iran","Italy",
+ "Japan","Korea, South","Netherlands","Portugal","Spain",
+ "United Kingdom", "US") )
+ On enlève les DOM-TOM, ainsi que les colonies du Royaume-Uni et des Pays-Bas :
+ + + +data = data %>% filter(Country.Region == "China" | Country.Region == "Hong Kong" | (Country.Region != "China" & Province.State == ""))On supprime les colonnes inutiles pour la suite de l’analyse :
+ + + + + + + +On regarde s’il y a des valeurs manquantes dans chaque colonne :
+ + + + + + +[1] 0Aucune valeur manquante à signaler.
+On regroupe les différentes provinces de Chine au sein d’une seule observation :
+ + + + + + + +On transpose la base de données pour avoir les dates en lignes et les noms de pays en colonnes, puis on retraite les noms des lignes pour qu’ils soient au format date :
+ + + +data = as.data.frame(data) # pour éviter que le df soit au format tibble
+rownames(data) = data$Country.Region
+data$Country.Region = NULL
+tdata = as.data.frame(t(data))
+tdata$Date = as.Date(rownames(tdata), format = "%m.%d.%y")On visualise le rendu final de nos données :
+ + + +rownames(tdata) = c()
+tdata = tdata %>%
+ select(Date, everything()) # On met la colonne de Date en première position
+tdata
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+Visualisation des données
+Avant de réaliser nos graphiques, on change le format de la base de données pour le passer de “Wide” à “Long”:
+ + + + + + +
+
+
+
+
+
+Puis on trace un graphique animé représentant la série temporelle du nombre cumulé de cas de Covid-19 pour les différents pays de notre base de données :
+ + + +library(ggplot2)
+library(gganimate)
+ylab = c(0,1,2,3,4,5,6)
+p = ggplot(meltdf,aes(x=Date,y=value,colour=variable,group=variable)) + geom_line(size = 1) + scale_x_date(date_breaks = "1 month", date_labels = "%b-%y") + labs(y="Nombre de cas en cumulé", x = "Mois") + scale_y_continuous(labels = paste0(ylab, "M"), breaks = 10^6 * ylab) +
+ labs(color='Pays') + transition_reveal(Date) + ggtitle("Evolution du nombre de cas de Covid-19 par pays")
+animate(p, duration = 8, fps = 20, width = 800, height = 500, renderer = gifski_renderer())
+anim_save('output.gif')Afin de pouvoir visualiser plus clairement les évolutions des différents pays, on réalise maintenant le même graphique mais avec une échelle transformée via la fonction Sinus hyperbolique inverse (la transformation logarithmique n’étant pas appropriée à cause des valeurs nulles) :
+ + + +asinh_trans = scales::trans_new(
+ "inverse_hyperbolic_sine",
+ transform = function(x) {asinh(x)},
+ inverse = function(x) {sinh(x)}
+ )
+p = ggplot(meltdf,aes(x=Date,y=value,colour=variable,group=variable)) +
+ geom_line(size = 1) + scale_x_date(date_breaks = "1 month", date_labels = "%b-%y") + labs(y="Nombre de cas en cumulé (échelle transformée)", x = "Mois") +
+ scale_y_continuous(trans = asinh_trans) + labs(color='Pays') +
+ ggtitle("Evolution du nombre de cas de Covid-19 par pays") + theme(axis.text=element_text(size=6.5),
+ axis.title=element_text(size=8),
+ plot.title = element_text(size=10),
+ legend.title = element_text(size = 8),
+ legend.text = element_text(size = 6.5))
+
+
+Etude du nombre de décès
+Nous réalisons maintenant la même étude que précédemment sur le nombre de décès dus à la Covid-19. Les données sont toujours téléchargées sur le Github du JHU CSSE.
+
+
+
+
+
+
+
+
+
+
+Préparation des données
+ + + +data_url = "https://raw.githubusercontent.com/CSSEGISandData/COVID-19/master/csse_covid_19_data/csse_covid_19_time_series/time_series_covid19_deaths_global.csv"
+data = read.csv(data_url, stringsAsFactors=FALSE)Regardons nos données :
+ + + + + + + +Nous effectuons les mêmes retraitements que précédemment :
+ + + +names(data)[5:ncol(data)] = substring(names(data)[5:ncol(data)],2)
+data[data$Province.State == "Hong Kong","Country.Region"] = "Hong Kong"
+data = data %>% filter(Country.Region %in% c("Belgium","China",
+ "Hong Kong","France","Germany","Iran","Italy",
+ "Japan","Korea, South","Netherlands","Portugal","Spain",
+ "United Kingdom", "US") )
+data = data %>% filter(Country.Region == "China" | Country.Region == "Hong Kong" | (Country.Region != "China" & Province.State == ""))
+data = data %>% select(-c('Province.State', 'Lat', 'Long'))
+data = data %>% group_by(Country.Region) %>% summarise_all(sum)
+data = as.data.frame(data) # pour éviter que le df soit au format tibble
+rownames(data) = data$Country.Region
+data$Country.Region = NULL
+tdata = as.data.frame(t(data))
+tdata$Date = as.Date(rownames(tdata), format = "%m.%d.%y")On visualise le rendu final de nos données :
+ + + + + + + +
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+Visualisation des données
+On change le format de la base de données pour le passer de “Wide” à “Long”:
+ + + + + + + +Puis on trace un graphique animé représentant la série temporelle du nombre cumulé de décès dus à la Covid-19 pour les différents pays de notre base de données :
+ + + + +
+Afin de pouvoir visualiser plus clairement les évolutions des différents pays, on réalise maintenant le même graphique mais avec une échelle transformée via la fonction Sinus hyperbolique inverse (la transformation logarithmique n’étant pas appropriée à cause des valeurs nulles) :
+ + + +p = ggplot(meltdf,aes(x=Date,y=value,colour=variable,group=variable)) +
+ geom_line(size = 1) + scale_x_date(date_breaks = "1 month", date_labels = "%b-%y") + labs(y="Nombre de décès en cumulé (échelle transformée)", x = "Mois") +
+ scale_y_continuous(trans = asinh_trans) + labs(color='Pays') +
+ ggtitle("Evolution du nombre de décès dus à la Covid-19 par pays") + theme(axis.text=element_text(size=6.5),
+ axis.title=element_text(size=8),
+ plot.title = element_text(size=10),
+ legend.title = element_text(size = 8),
+ legend.text = element_text(size = 6.5))
LS0tCnRpdGxlOiAiQXV0b3VyIGR1IFNBUlMtQ29WLTIgKENvdmlkLTE5KSIKYXV0aG9yOiAiRWxpYSBQw6lyZW5uw6hzIgpvdXRwdXQ6CiAgaHRtbF9ub3RlYm9vayA6CiAgICB0b2M6IHRydWUKICAgIHRvY19kZXB0aDogMgogICAgbnVtYmVyX3NlY3Rpb25zOiBmYWxzZQogICAgdGhlbWU6IGpvdXJuYWwKICAgIGhpZ2hsaWdodDogdGFuZ28KZG9jdW1lbnRjbGFzczogYXJ0aWNsZQpjbGFzc29wdGlvbjogYTRwYXBlcgpoZWFkZXItaW5jbHVkZXM6Ci0gXHVzZXBhY2thZ2VbZnJlbmNoXXtiYWJlbH0KLSBcdXNlcGFja2FnZVt1cHJpZ2h0XXtmb3VyaWVyfQotIFxoeXBlcnNldHVwe2NvbG9ybGlua3M9dHJ1ZSxwYWdlYmFja3JlZj10cnVlfQotLS0KCmBgYHtyIHNldHVwLCBpbmNsdWRlPUZBTFNFfQprbml0cjo6b3B0c19jaHVuayRzZXQoZWNobyA9IFRSVUUsIGZpZy5hbGlnbj0iY2VudGVyIikgCmBgYAoKTGUgYnV0IGVzdCBpY2kgZGUgcmVwcm9kdWlyZSBkZXMgZ3JhcGhlcyBzZW1ibGFibGVzIMOgIGNldXggZHUgW1NvdXRoIENoaW5hIE1vcm5pbmcgUG9zdCAoU0NNUCldKGh0dHBzOi8vd3d3LnNjbXAuY29tLyksIHN1ciBsYSBwYWdlIFtUaGUgQ29yb25hdmlydXMgUGFuZGVtaWNdKGh0dHBzOi8vd3d3LnNjbXAuY29tL2Nvcm9uYXZpcnVzP3NyYz1ob21lcGFnZV9jb3ZpZF93aWRnZXQpIGV0IHF1aSBtb250cmVudCBwb3VyIGRpZmbDqXJlbnRzIHBheXMgbGUgbm9tYnJlIGN1bXVsw6kgKGMnZXN0LcOgLWRpcmUgbGUgbm9tYnJlIHRvdGFsIGRlIGNhcyBkZXB1aXMgbGUgZMOpYnV0IGRlIGwnw6lwaWTDqW1pZSkgZGUgcGVyc29ubmVzIGF0dGVpbnRlcyBkZSBsYSBtYWxhZGllIMOgIGNvcm9uYXZpcnVzIDIwMTkuIExlcyBkb25uw6llcyBxdWUgbm91cyB1dGlsaXNlcm9ucyBzb250IGNvbXBpbMOpZXMgcGFyIGxlIFsqSm9obnMgSG9wa2lucyBVbml2ZXJzaXR5IENlbnRlciBmb3IgU3lzdGVtcyBTY2llbmNlIGFuZCBFbmdpbmVlcmluZyAoSkhVIENTU0UpKl0oaHR0cHM6Ly9zeXN0ZW1zLmpodS5lZHUvKS4KClRvdXQgZCdhYm9yZCwgb24gbmV0dG9pZSBsJ2VzcGFjZSBkZSB0cmF2YWlsIHN1ciBSIGV0IG9uIGNoYXJnZSBsZXMgbGlicmFpcmllcyBuw6ljZXNzYWlyZXMgw6AgbCdleMOpY3V0aW9uIGRlIG5vdHJlIGNvZGUgOgpgYGB7cn0Kcm0obGlzdCA9IGxzKCkpCmxpYnJhcnkoZHBseXIpCmxpYnJhcnkoZ2dwbG90MikKbGlicmFyeShnZ2FuaW1hdGUpCmxpYnJhcnkocmVzaGFwZTIpCmBgYAoKCiMgRXR1ZGUgZHUgbm9tYnJlIGRlIGNhcyAKCiMjIFByw6lwYXJhdGlvbiBkZXMgZG9ubsOpZXMKCk5vdXMgdMOpbMOpY2hhcmdlb25zIGxlcyBkb25uw6llcyByZWNlbnNhbnQgbGUgbm9tYnJlIGRlIGNhcyBkZSBDb3ZpZC0xOSBkaXNwb25pYmxlcyBzdXIgbGUgW0dpdGh1YiBkdSBKSFUgQ1NTRV0oaHR0cHM6Ly9naXRodWIuY29tL0NTU0VHSVNhbmREYXRhL0NPVklELTE5KS4KYGBge3J9CmRhdGFfdXJsID0gImh0dHBzOi8vcmF3LmdpdGh1YnVzZXJjb250ZW50LmNvbS9DU1NFR0lTYW5kRGF0YS9DT1ZJRC0xOS9tYXN0ZXIvY3NzZV9jb3ZpZF8xOV9kYXRhL2Nzc2VfY292aWRfMTlfdGltZV9zZXJpZXMvdGltZV9zZXJpZXNfY292aWQxOV9jb25maXJtZWRfZ2xvYmFsLmNzdiIKZGF0YSA9IHJlYWQuY3N2KGRhdGFfdXJsLCBzdHJpbmdzQXNGYWN0b3JzPUZBTFNFKQpgYGAKCgpSZWdhcmRvbnMgbm9zIGRvbm7DqWVzIDoKYGBge3J9CmRhdGEKYGBgCgpPbiByZXRyYWl0ZSBkYW5zIHVuIHByZW1pZXIgdGVtcHMgbGVzIG5vbXMgZGVzIGNvbG9ubmVzIGRlIGRhdGVzIGVuIGVubGV2YW50IGxlICpYKiBkZXZhbnQgY2hhcXVlIG5vbSBkZSBjb2xvbm5lIDoKYGBge3J9Cm5hbWVzKGRhdGEpWzU6bmNvbChkYXRhKV0gPSBzdWJzdHJpbmcobmFtZXMoZGF0YSlbNTpuY29sKGRhdGEpXSwyKQpuYW1lcyhkYXRhKVsxOjEwXQpgYGAKCk9uIHPDqXBhcmUgbGVzIGRvbm7DqWVzIHN1ciBIb25nLUtvbmcgZGUgY2VsbGVzIGR1IHJlc3RlIGRlIGxhIENoaW5lIDoKCmBgYHtyfQpkYXRhW2RhdGEkUHJvdmluY2UuU3RhdGUgPT0gIkhvbmcgS29uZyIsIkNvdW50cnkuUmVnaW9uIl0gPSAiSG9uZyBLb25nIgpgYGAKCgpPbiBmaWx0cmUgbGEgYmFzZSBkZSBkb25uw6llcyBwb3VyIG5lIGdhcmRlciBxdWUgbGVzIHBheXMgc3VpdmFudHMgOiBsYSBCZWxnaXF1ZSAoQmVsZ2l1bSksIGxhIENoaW5lIC0gdG91dGVzIGxlcyBwcm92aW5jZXMgc2F1ZiBIb25nLUtvbmcgKENoaW5hKSwgSG9uZyBLb25nIChDaGluYSwgSG9uZy1Lb25nKSwgbGEgRnJhbmNlIG3DqXRyb3BvbGl0YWluZSAoRnJhbmNlKSwgbOKAmUFsbGVtYWduZSAoR2VybWFueSksIGzigJlJcmFuIChJcmFuKSwgbOKAmUl0YWxpZSAoSXRhbHkpLCBsZSBKYXBvbiAoSmFwYW4pLCBsYSBDb3LDqWUgZHUgU3VkIChLb3JlYSwgU291dGgpLCBsYSBIb2xsYW5kZSBzYW5zIGxlcyBjb2xvbmllcyAoTmV0aGVybGFuZHMpLCBsZSBQb3J0dWdhbCAoUG9ydHVnYWwpLCBs4oCZRXNwYWduZSAoU3BhaW4pLCBsZSBSb3lhdW1lLVVuaXMgc2FucyBsZXMgY29sb25pZXMgKFVuaXRlZCBLaW5nZG9tKSBhaW5zaSBxdWUgbGVzIMOJdGF0cy1VbmlzIChVUykuCgpgYGB7cn0KZGF0YSA9IGRhdGEgJT4lIGZpbHRlcihDb3VudHJ5LlJlZ2lvbiAlaW4lIGMoIkJlbGdpdW0iLCJDaGluYSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJIb25nIEtvbmciLCJGcmFuY2UiLCJHZXJtYW55IiwiSXJhbiIsIkl0YWx5IiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIkphcGFuIiwiS29yZWEsIFNvdXRoIiwiTmV0aGVybGFuZHMiLCJQb3J0dWdhbCIsIlNwYWluIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlVuaXRlZCBLaW5nZG9tIiwgIlVTIikgKQogIApgYGAKCk9uIGVubMOodmUgbGVzIERPTS1UT00sIGFpbnNpIHF1ZSBsZXMgY29sb25pZXMgZHUgUm95YXVtZS1VbmkgZXQgZGVzIFBheXMtQmFzIDogCmBgYHtyfQpkYXRhID0gZGF0YSAlPiUgZmlsdGVyKENvdW50cnkuUmVnaW9uID09ICJDaGluYSIgfCBDb3VudHJ5LlJlZ2lvbiA9PSAiSG9uZyBLb25nIiB8IChDb3VudHJ5LlJlZ2lvbiAhPSAiQ2hpbmEiICYgUHJvdmluY2UuU3RhdGUgPT0gIiIpKQpgYGAKCk9uIHN1cHByaW1lIGxlcyBjb2xvbm5lcyBpbnV0aWxlcyBwb3VyIGxhIHN1aXRlIGRlIGwnYW5hbHlzZSA6CmBgYHtyfQpkYXRhID0gZGF0YSAlPiUgc2VsZWN0KC1jKCdQcm92aW5jZS5TdGF0ZScsICdMYXQnLCAnTG9uZycpKQpgYGAKCk9uIHJlZ2FyZGUgcydpbCB5IGEgZGVzIHZhbGV1cnMgbWFucXVhbnRlcyBkYW5zIGNoYXF1ZSBjb2xvbm5lIDoKYGBge3J9CnN1bShjb2xTdW1zKGlzLm5hKGRhdGEpKSkKYGBgCkF1Y3VuZSB2YWxldXIgbWFucXVhbnRlIMOgIHNpZ25hbGVyLgoKT24gcmVncm91cGUgbGVzIGRpZmbDqXJlbnRlcyBwcm92aW5jZXMgZGUgQ2hpbmUgYXUgc2VpbiBkJ3VuZSBzZXVsZSBvYnNlcnZhdGlvbiA6CmBgYHtyfQpkYXRhID0gZGF0YSAlPiUgZ3JvdXBfYnkoQ291bnRyeS5SZWdpb24pICU+JSBzdW1tYXJpc2VfYWxsKHN1bSkKYGBgCgpPbiB0cmFuc3Bvc2UgbGEgYmFzZSBkZSBkb25uw6llcyBwb3VyIGF2b2lyIGxlcyBkYXRlcyBlbiBsaWduZXMgZXQgbGVzIG5vbXMgZGUgcGF5cyBlbiBjb2xvbm5lcywgcHVpcyBvbiByZXRyYWl0ZSBsZXMgbm9tcyBkZXMgbGlnbmVzIHBvdXIgcXUnaWxzIHNvaWVudCBhdSBmb3JtYXQgZGF0ZSA6CmBgYHtyfQpkYXRhID0gYXMuZGF0YS5mcmFtZShkYXRhKSAjIHBvdXIgw6l2aXRlciBxdWUgbGUgZGYgc29pdCBhdSBmb3JtYXQgdGliYmxlCnJvd25hbWVzKGRhdGEpID0gZGF0YSRDb3VudHJ5LlJlZ2lvbgpkYXRhJENvdW50cnkuUmVnaW9uID0gTlVMTAp0ZGF0YSA9IGFzLmRhdGEuZnJhbWUodChkYXRhKSkKdGRhdGEkRGF0ZSA9IGFzLkRhdGUocm93bmFtZXModGRhdGEpLCBmb3JtYXQgPSAiJW0uJWQuJXkiKQpgYGAKCk9uIHZpc3VhbGlzZSBsZSAqKnJlbmR1IGZpbmFsIGRlIG5vcyBkb25uw6llcyoqIDoKYGBge3J9CnJvd25hbWVzKHRkYXRhKSA9IGMoKQp0ZGF0YSA9IHRkYXRhICU+JQogIHNlbGVjdChEYXRlLCBldmVyeXRoaW5nKCkpICMgT24gbWV0IGxhIGNvbG9ubmUgZGUgRGF0ZSBlbiBwcmVtacOocmUgcG9zaXRpb24KdGRhdGEKYGBgCgoKIyMgVmlzdWFsaXNhdGlvbiBkZXMgZG9ubsOpZXMKCkF2YW50IGRlIHLDqWFsaXNlciBub3MgZ3JhcGhpcXVlcywgb24gY2hhbmdlIGxlIGZvcm1hdCBkZSBsYSBiYXNlIGRlIGRvbm7DqWVzIHBvdXIgbGUgcGFzc2VyIGRlICJXaWRlIiDDoCAiTG9uZyI6CmBgYHtyfQptZWx0ZGYgPSBtZWx0KHRkYXRhLCBpZCA9ICJEYXRlIikKbWVsdGRmCmBgYAoKUHVpcyBvbiB0cmFjZSB1biBncmFwaGlxdWUgYW5pbcOpIHJlcHLDqXNlbnRhbnQgbGEgc8OpcmllIHRlbXBvcmVsbGUgZHUgbm9tYnJlIGN1bXVsw6kgZGUgY2FzIGRlIENvdmlkLTE5IHBvdXIgbGVzIGRpZmbDqXJlbnRzIHBheXMgZGUgbm90cmUgYmFzZSBkZSBkb25uw6llcyA6CmBgYHtyLCByZXN1bHRzID0gImhpZGUifQp5bGFiID0gYygwLDEsMiwzLDQsNSw2KQpwID0gZ2dwbG90KG1lbHRkZixhZXMoeD1EYXRlLHk9dmFsdWUsY29sb3VyPXZhcmlhYmxlLGdyb3VwPXZhcmlhYmxlKSkgKyBnZW9tX2xpbmUoc2l6ZSA9IDEpICsgc2NhbGVfeF9kYXRlKGRhdGVfYnJlYWtzID0gIjEgbW9udGgiLCBkYXRlX2xhYmVscyA9ICIlYi0leSIpICsgbGFicyh5PSJOb21icmUgZGUgY2FzIGVuIGN1bXVsw6kiLCB4ID0gIk1vaXMiKSArIHNjYWxlX3lfY29udGludW91cyhsYWJlbHMgPSBwYXN0ZTAoeWxhYiwgIk0iKSwgYnJlYWtzID0gMTBeNiAqIHlsYWIpICsKICBsYWJzKGNvbG9yPSdQYXlzJykgKyB0cmFuc2l0aW9uX3JldmVhbChEYXRlKSArIGdndGl0bGUoIkV2b2x1dGlvbiBkdSBub21icmUgZGUgY2FzIGRlIENvdmlkLTE5IHBhciBwYXlzIikKYW5pbWF0ZShwLCBkdXJhdGlvbiA9IDgsIGZwcyA9IDIwLCB3aWR0aCA9IDgwMCwgaGVpZ2h0ID0gNTAwLCByZW5kZXJlciA9IGdpZnNraV9yZW5kZXJlcigpKQphbmltX3NhdmUoJ291dHB1dC5naWYnKQpgYGAKCjxjZW50ZXI+CiFbXShvdXRwdXQuZ2lmKXtoZWlnaHQ9IjUwMCIgd2lkdGg9IjgwMCJ9CjwvY2VudGVyPgoKCkFmaW4gZGUgcG91dm9pciB2aXN1YWxpc2VyIHBsdXMgY2xhaXJlbWVudCBsZXMgw6l2b2x1dGlvbnMgZGVzIGRpZmbDqXJlbnRzIHBheXMsIG9uIHLDqWFsaXNlIG1haW50ZW5hbnQgbGUgbcOqbWUgZ3JhcGhpcXVlIG1haXMgYXZlYyB1bmUgw6ljaGVsbGUgdHJhbnNmb3Jtw6llIHZpYSBsYSBmb25jdGlvbiBTaW51cyBoeXBlcmJvbGlxdWUgaW52ZXJzZSAobGEgdHJhbnNmb3JtYXRpb24gbG9nYXJpdGhtaXF1ZSBuJ8OpdGFudCBwYXMgYXBwcm9wcmnDqWUgw6AgY2F1c2UgZGVzIHZhbGV1cnMgbnVsbGVzKSAgOgoKYGBge3IsIGZpZy5hbGlnbj0nY2VudGVyJ30KYXNpbmhfdHJhbnMgPSBzY2FsZXM6OnRyYW5zX25ldygKICAiaW52ZXJzZV9oeXBlcmJvbGljX3NpbmUiLAogIHRyYW5zZm9ybSA9IGZ1bmN0aW9uKHgpIHthc2luaCh4KX0sCiAgaW52ZXJzZSA9IGZ1bmN0aW9uKHgpIHtzaW5oKHgpfQogICkKcCA9IGdncGxvdChtZWx0ZGYsYWVzKHg9RGF0ZSx5PXZhbHVlLGNvbG91cj12YXJpYWJsZSxncm91cD12YXJpYWJsZSkpICsgCiAgZ2VvbV9saW5lKHNpemUgPSAxKSArIHNjYWxlX3hfZGF0ZShkYXRlX2JyZWFrcyA9ICIxIG1vbnRoIiwgZGF0ZV9sYWJlbHMgPSAiJWItJXkiKSArIGxhYnMoeT0iTm9tYnJlIGRlIGNhcyBlbiBjdW11bMOpICjDqWNoZWxsZSB0cmFuc2Zvcm3DqWUpIiwgeCA9ICJNb2lzIikgKyAKICBzY2FsZV95X2NvbnRpbnVvdXModHJhbnMgPSBhc2luaF90cmFucykgKyBsYWJzKGNvbG9yPSdQYXlzJykgKyAKICBnZ3RpdGxlKCJFdm9sdXRpb24gZHUgbm9tYnJlIGRlIGNhcyBkZSBDb3ZpZC0xOSBwYXIgcGF5cyIpICsgdGhlbWUoYXhpcy50ZXh0PWVsZW1lbnRfdGV4dChzaXplPTYuNSksCiAgICAgICAgYXhpcy50aXRsZT1lbGVtZW50X3RleHQoc2l6ZT04KSwKICAgICAgICBwbG90LnRpdGxlID0gZWxlbWVudF90ZXh0KHNpemU9MTApLAogICAgICAgIGxlZ2VuZC50aXRsZSA9IGVsZW1lbnRfdGV4dChzaXplID0gOCksCiAgICAgICAgbGVnZW5kLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDYuNSkpCmBgYAoKPGNlbnRlcj4KYGBge3IgZmlnLndpZHRoPTguNSwgZmlnLmhlaWdodD01LCBlY2hvID0gRkFMU0V9CnAKYGBgCjwvY2VudGVyPgoKIyBFdHVkZSBkdSBub21icmUgZGUgZMOpY8OocwoKTm91cyByw6lhbGlzb25zIG1haW50ZW5hbnQgbGEgbcOqbWUgw6l0dWRlIHF1ZSBwcsOpY8OpZGVtbWVudCBzdXIgbGUgbm9tYnJlIGRlIGTDqWPDqHMgZHVzIMOgIGxhIENvdmlkLTE5LiBMZXMgZG9ubsOpZXMgc29udCB0b3Vqb3VycyB0w6lsw6ljaGFyZ8OpZXMgc3VyIGxlIFtHaXRodWIgZHUgSkhVIENTU0VdKGh0dHBzOi8vZ2l0aHViLmNvbS9DU1NFR0lTYW5kRGF0YS9DT1ZJRC0xOSkuCgojIyBQcsOpcGFyYXRpb24gZGVzIGRvbm7DqWVzCgpgYGB7cn0KZGF0YV91cmwgPSAiaHR0cHM6Ly9yYXcuZ2l0aHVidXNlcmNvbnRlbnQuY29tL0NTU0VHSVNhbmREYXRhL0NPVklELTE5L21hc3Rlci9jc3NlX2NvdmlkXzE5X2RhdGEvY3NzZV9jb3ZpZF8xOV90aW1lX3Nlcmllcy90aW1lX3Nlcmllc19jb3ZpZDE5X2RlYXRoc19nbG9iYWwuY3N2IgpkYXRhID0gcmVhZC5jc3YoZGF0YV91cmwsIHN0cmluZ3NBc0ZhY3RvcnM9RkFMU0UpCmBgYAoKClJlZ2FyZG9ucyBub3MgZG9ubsOpZXMgOgpgYGB7cn0KZGF0YQpgYGAKCk5vdXMgZWZmZWN0dW9ucyBsZXMgbcOqbWVzIHJldHJhaXRlbWVudHMgcXVlIHByw6ljw6lkZW1tZW50IDoKYGBge3J9Cm5hbWVzKGRhdGEpWzU6bmNvbChkYXRhKV0gPSBzdWJzdHJpbmcobmFtZXMoZGF0YSlbNTpuY29sKGRhdGEpXSwyKQpkYXRhW2RhdGEkUHJvdmluY2UuU3RhdGUgPT0gIkhvbmcgS29uZyIsIkNvdW50cnkuUmVnaW9uIl0gPSAiSG9uZyBLb25nIgpkYXRhID0gZGF0YSAlPiUgZmlsdGVyKENvdW50cnkuUmVnaW9uICVpbiUgYygiQmVsZ2l1bSIsIkNoaW5hIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIkhvbmcgS29uZyIsIkZyYW5jZSIsIkdlcm1hbnkiLCJJcmFuIiwiSXRhbHkiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiSmFwYW4iLCJLb3JlYSwgU291dGgiLCJOZXRoZXJsYW5kcyIsIlBvcnR1Z2FsIiwiU3BhaW4iLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiVW5pdGVkIEtpbmdkb20iLCAiVVMiKSApCmRhdGEgPSBkYXRhICU+JSBmaWx0ZXIoQ291bnRyeS5SZWdpb24gPT0gIkNoaW5hIiB8IENvdW50cnkuUmVnaW9uID09ICJIb25nIEtvbmciIHwgKENvdW50cnkuUmVnaW9uICE9ICJDaGluYSIgJiBQcm92aW5jZS5TdGF0ZSA9PSAiIikpCmRhdGEgPSBkYXRhICU+JSBzZWxlY3QoLWMoJ1Byb3ZpbmNlLlN0YXRlJywgJ0xhdCcsICdMb25nJykpCmRhdGEgPSBkYXRhICU+JSBncm91cF9ieShDb3VudHJ5LlJlZ2lvbikgJT4lIHN1bW1hcmlzZV9hbGwoc3VtKQpkYXRhID0gYXMuZGF0YS5mcmFtZShkYXRhKSAjIHBvdXIgw6l2aXRlciBxdWUgbGUgZGYgc29pdCBhdSBmb3JtYXQgdGliYmxlCnJvd25hbWVzKGRhdGEpID0gZGF0YSRDb3VudHJ5LlJlZ2lvbgpkYXRhJENvdW50cnkuUmVnaW9uID0gTlVMTAp0ZGF0YSA9IGFzLmRhdGEuZnJhbWUodChkYXRhKSkKdGRhdGEkRGF0ZSA9IGFzLkRhdGUocm93bmFtZXModGRhdGEpLCBmb3JtYXQgPSAiJW0uJWQuJXkiKQpgYGAKCk9uIHZpc3VhbGlzZSBsZSAqKnJlbmR1IGZpbmFsIGRlIG5vcyBkb25uw6llcyoqIDoKYGBge3J9CnJvd25hbWVzKHRkYXRhKSA8LSBjKCkKdGRhdGEgPSB0ZGF0YSAlPiUKICBzZWxlY3QoRGF0ZSwgZXZlcnl0aGluZygpKSAjIE9uIG1ldCBsYSBjb2xvbm5lIGRlIERhdGUgZW4gcHJlbWnDqHJlIHBvc2l0aW9uCnRkYXRhCmBgYAoKCiMjIFZpc3VhbGlzYXRpb24gZGVzIGRvbm7DqWVzCgpPbiBjaGFuZ2UgbGUgZm9ybWF0IGRlIGxhIGJhc2UgZGUgZG9ubsOpZXMgcG91ciBsZSBwYXNzZXIgZGUgIldpZGUiIMOgICJMb25nIjoKYGBge3J9Cm1lbHRkZiA9IG1lbHQodGRhdGEsIGlkID0gIkRhdGUiKQptZWx0ZGYKYGBgCgpQdWlzIG9uIHRyYWNlIHVuIGdyYXBoaXF1ZSBhbmltw6kgcmVwcsOpc2VudGFudCBsYSBzw6lyaWUgdGVtcG9yZWxsZSBkdSBub21icmUgY3VtdWzDqSBkZSBkw6ljw6hzIGR1cyDDoCBsYSBDb3ZpZC0xOSBwb3VyIGxlcyBkaWZmw6lyZW50cyBwYXlzIGRlIG5vdHJlIGJhc2UgZGUgZG9ubsOpZXMgOgpgYGB7ciwgcmVzdWx0cyA9ICJoaWRlIn0KcCA9IGdncGxvdChtZWx0ZGYsYWVzKHg9RGF0ZSx5PXZhbHVlLGNvbG91cj12YXJpYWJsZSxncm91cD12YXJpYWJsZSkpICsgZ2VvbV9saW5lKHNpemUgPSAxKSArIHNjYWxlX3hfZGF0ZShkYXRlX2JyZWFrcyA9ICIxIG1vbnRoIiwgZGF0ZV9sYWJlbHMgPSAiJWItJXkiKSArIGxhYnMoeT0iTm9tYnJlIGRlIGTDqWPDqHMgZW4gY3VtdWzDqSIsIHggPSAiTW9pcyIpICsgbGFicyhjb2xvcj0nUGF5cycpICsgdHJhbnNpdGlvbl9yZXZlYWwoRGF0ZSkgKyBnZ3RpdGxlKCJFdm9sdXRpb24gZHUgbm9tYnJlIGRlIGTDqWPDqHMgZHVzIMOgIGxhIENvdmlkLTE5IHBhciBwYXlzIikKYW5pbWF0ZShwLCBkdXJhdGlvbiA9IDgsIGZwcyA9IDIwLCB3aWR0aCA9IDgwMCwgaGVpZ2h0ID0gNTAwLCByZW5kZXJlciA9IGdpZnNraV9yZW5kZXJlcigpKQphbmltX3NhdmUoJ291dHB1dGRlYXRocy5naWYnKQpgYGAKCjxjZW50ZXI+CiFbXShvdXRwdXRkZWF0aHMuZ2lmKXtoZWlnaHQ9IjUwMCIgd2lkdGg9IjgwMCJ9CjwvY2VudGVyPgoKCkFmaW4gZGUgcG91dm9pciB2aXN1YWxpc2VyIHBsdXMgY2xhaXJlbWVudCBsZXMgw6l2b2x1dGlvbnMgZGVzIGRpZmbDqXJlbnRzIHBheXMsIG9uIHLDqWFsaXNlIG1haW50ZW5hbnQgbGUgbcOqbWUgZ3JhcGhpcXVlIG1haXMgYXZlYyB1bmUgw6ljaGVsbGUgdHJhbnNmb3Jtw6llIHZpYSBsYSBmb25jdGlvbiBTaW51cyBoeXBlcmJvbGlxdWUgaW52ZXJzZSAobGEgdHJhbnNmb3JtYXRpb24gbG9nYXJpdGhtaXF1ZSBuJ8OpdGFudCBwYXMgYXBwcm9wcmnDqWUgw6AgY2F1c2UgZGVzIHZhbGV1cnMgbnVsbGVzKSAgOgoKYGBge3J9CnAgPSBnZ3Bsb3QobWVsdGRmLGFlcyh4PURhdGUseT12YWx1ZSxjb2xvdXI9dmFyaWFibGUsZ3JvdXA9dmFyaWFibGUpKSArIAogIGdlb21fbGluZShzaXplID0gMSkgKyBzY2FsZV94X2RhdGUoZGF0ZV9icmVha3MgPSAiMSBtb250aCIsIGRhdGVfbGFiZWxzID0gIiViLSV5IikgKyBsYWJzKHk9Ik5vbWJyZSBkZSBkw6ljw6hzIGVuIGN1bXVsw6kgKMOpY2hlbGxlIHRyYW5zZm9ybcOpZSkiLCB4ID0gIk1vaXMiKSArICAKICBzY2FsZV95X2NvbnRpbnVvdXModHJhbnMgPSBhc2luaF90cmFucykgKyBsYWJzKGNvbG9yPSdQYXlzJykgKyAKICBnZ3RpdGxlKCJFdm9sdXRpb24gZHUgbm9tYnJlIGRlIGTDqWPDqHMgZHVzIMOgIGxhIENvdmlkLTE5IHBhciBwYXlzIikgICsgdGhlbWUoYXhpcy50ZXh0PWVsZW1lbnRfdGV4dChzaXplPTYuNSksCiAgICAgICAgYXhpcy50aXRsZT1lbGVtZW50X3RleHQoc2l6ZT04KSwKICAgICAgICBwbG90LnRpdGxlID0gZWxlbWVudF90ZXh0KHNpemU9MTApLAogICAgICAgIGxlZ2VuZC50aXRsZSA9IGVsZW1lbnRfdGV4dChzaXplID0gOCksCiAgICAgICAgbGVnZW5kLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDYuNSkpCmBgYAoKPGNlbnRlcj4KYGBge3IgZmlnLndpZHRoPTguNSwgZmlnLmhlaWdodD01LCBlY2hvID0gRkFMU0V9CnAKYGBgCjwvY2VudGVyPgoKCg==
+
+
+
+