Périphérisation et périphérisation infrastructurelle en Bulgarie
Raphaëlle Segond
2026-01-07
0.Préalables à la réalisation de l’analyse
0.A.Chargement des packages et options
options("digits"=2,"scipen" =100)
ggrepel.max.overlaps=Inf
library(shiny)
library(rsconnect)##
## Attachement du package : 'rsconnect'## L'objet suivant est masqué depuis 'package:shiny':
##
## serverInfolibrary(rmarkdown)
library(dplyr)##
## Attachement du package : 'dplyr'## Les objets suivants sont masqués depuis 'package:stats':
##
## filter, lag## Les objets suivants sont masqués depuis 'package:base':
##
## intersect, setdiff, setequal, unionlibrary(FactoMineR)
library(factoextra)## Le chargement a nécessité le package : ggplot2## Welcome! Want to learn more? See two factoextra-related books at https://goo.gl/ve3WBalibrary(psych)##
## Attachement du package : 'psych'## Les objets suivants sont masqués depuis 'package:ggplot2':
##
## %+%, alphalibrary(Rcpp)
library(data.table)##
## Attachement du package : 'data.table'## Les objets suivants sont masqués depuis 'package:dplyr':
##
## between, first, lastlibrary(questionr)##
## Attachement du package : 'questionr'## L'objet suivant est masqué depuis 'package:psych':
##
## describelibrary(explor)
library(ade4)##
## Attachement du package : 'ade4'## L'objet suivant est masqué depuis 'package:FactoMineR':
##
## reconstlibrary(cluster)
library(corrplot)## corrplot 0.95 loadedlibrary(missMDA)
library(stats)
library(reshape2)##
## Attachement du package : 'reshape2'## Les objets suivants sont masqués depuis 'package:data.table':
##
## dcast, meltlibrary(ggplot2)
library(readxl)
library(tibble)
library(scales)##
## Attachement du package : 'scales'## Les objets suivants sont masqués depuis 'package:psych':
##
## alpha, rescalelibrary(tidyverse)## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ forcats 1.0.1 ✔ readr 2.1.6
## ✔ lubridate 1.9.4 ✔ stringr 1.6.0
## ✔ purrr 1.2.0 ✔ tidyr 1.3.2## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ psych::%+%() masks ggplot2::%+%()
## ✖ scales::alpha() masks psych::alpha(), ggplot2::alpha()
## ✖ data.table::between() masks dplyr::between()
## ✖ readr::col_factor() masks scales::col_factor()
## ✖ purrr::discard() masks scales::discard()
## ✖ dplyr::filter() masks stats::filter()
## ✖ data.table::first() masks dplyr::first()
## ✖ lubridate::hour() masks data.table::hour()
## ✖ lubridate::isoweek() masks data.table::isoweek()
## ✖ lubridate::isoyear() masks data.table::isoyear()
## ✖ dplyr::lag() masks stats::lag()
## ✖ data.table::last() masks dplyr::last()
## ✖ lubridate::mday() masks data.table::mday()
## ✖ lubridate::minute() masks data.table::minute()
## ✖ lubridate::month() masks data.table::month()
## ✖ lubridate::quarter() masks data.table::quarter()
## ✖ lubridate::second() masks data.table::second()
## ✖ purrr::transpose() masks data.table::transpose()
## ✖ lubridate::wday() masks data.table::wday()
## ✖ lubridate::week() masks data.table::week()
## ✖ lubridate::yday() masks data.table::yday()
## ✖ lubridate::year() masks data.table::year()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errorslibrary(ggpubr)
library(rstatix)##
## Attachement du package : 'rstatix'
##
## L'objet suivant est masqué depuis 'package:stats':
##
## filter0.B. Importation du tableur de données et du tableur de métadonnées
0.B.a. Tableur de données
Le tableur de données initial (“ViK”) contient 28 individus (les 28
régions et ViKs régionaux) et 36 variables. Une variable, la part du
prix de l’eau potable dans le budget des ViKs sera calculée par la
suite, pour l’analyse. Ces variables sont toutes décrites dans le
tableur des métadonnées (“META”).
Les ViKs municipaux n’ont pas été pris en compte dans cette analyse. Les
données sont essentiellement issues de la statistique nationale (INS) et
des rapports du régulateur, le KEVR pour les années 2013 (données
anciennes) et 2021 à 2024 (données actuelles). L’année 2021 est une
année particulière à prendre en compte puisque c’est celle où la
pandémie de COVID et la guerre en Ukraine entraînent une forte
inflation, notamment des prix de l’énergie, et où beaucoup de travaux de
maintenance sont reportés. Parmi les indicateurs de performance, seuls
trois ont été retenus parmi l’ensemble des KPIs du KEVER car ils
concernent directement la qualité du service d’approvisionnement en eau.
L’indicateur de la qualité de l’eau aurait pu être retenu mais il est
très homogène et globalement élevé.
ViK<-read.csv2("C:/Users/Raph/Documents/thesebulgarie/STATISTIQUES/analysesR/VIK_INEG/STAT_VIK_R.csv")row.names(ViK) <- ViK$REG_FRViK## REG_FR CODE clusCAH RevHab_24 Pop65_24 TxMig_14.24
## Blagoevgrad Blagoevgrad BGR.1_1 CLUS 4 25997 23 -20.0
## Montana Montana BGR.10_1 CLUS 1 26217 28 -28.2
## Pazardzhik Pazardzhik BGR.11_1 CLUS 1 24736 26 -17.9
## Pernik Pernik BGR.12_1 CLUS 1 30570 30 37.5
## Pleven Pleven BGR.13_1 CLUS 1 26815 31 -36.1
## Plovdiv Plovdiv BGR.14_1 CLUS 2 24742 26 49.8
## Razgrad Razgrad BGR.15_1 CLUS 4 21316 27 -23.3
## Ruse Ruse BGR.16_1 CLUS 1 22413 27 -5.2
## Shumen Shumen BGR.17_1 CLUS 4 27406 26 38.3
## Silistra Silistra BGR.18_1 CLUS 4 25220 29 -14.6
## Sliven Sliven BGR.19_1 CLUS 4 23765 23 -44.4
## Burgas Burgas BGR.2_1 CLUS 2 26769 29 69.0
## Smoljan Smoljan BGR.20_1 CLUS 1 29927 26 -87.8
## Sofia region Sofia region BGR.21_1 CLUS 3 18820 24 59.1
## Stara Zagora Stara Zagora BGR.22_1 CLUS 3 24704 26 6.7
## Targovishte Targovishte BGR.23_1 CLUS 4 23476 25 -3.5
## Varna Varna BGR.24_1 CLUS 2 23076 28 54.1
## Veliko Tarnovo Veliko Tarnovo BGR.25_1 CLUS 1 22549 25 -20.4
## Vidin Vidin BGR.26_1 CLUS 4 21539 23 -46.7
## Vraca Vraca BGR.27_1 CLUS 3 23540 23 -48.2
## Jambol Jambol BGR.28_1 CLUS 1 17821 27 -32.2
## Dobrich Dobrich BGR.3_1 CLUS 1 20742 28 -11.0
## Gabrovo Gabrovo BGR.4_1 CLUS 1 28608 31 -10.0
## Sofia capitale Sofia capitale BGR.5_1 CLUS 5 31829 19 41.2
## Haskovo Haskovo BGR.6_1 CLUS 1 20292 27 8.5
## Kardzhali Kardzhali BGR.7_1 CLUS 6 21024 25 188.4
## Kjustendil Kjustendil BGR.8_1 CLUS 1 23491 28 -18.2
## Lovech Lovech BGR.9_1 CLUS 4 17555 22 -19.6
## DensPop_24 PEVO_POP_01.21 NonBG_21 AgriVA_23 IndusVA_23
## Blagoevgrad 44 -11.4 18.9 5.534 25
## Montana 31 -31.6 17.0 11.619 24
## Pazardzhik 50 -17.8 19.6 5.953 39
## Pernik 46 -23.0 5.9 3.486 31
## Pleven 46 -30.5 11.5 5.763 28
## Plovdiv 107 -11.2 19.1 2.427 35
## Razgrad 41 -29.7 62.3 11.602 29
## Ruse 65 -27.5 23.1 5.187 31
## Shumen 44 -25.5 45.9 10.169 32
## Silistra 33 -30.6 47.2 16.487 16
## Sliven 48 -20.5 33.1 6.727 24
## Burgas 50 -10.0 26.3 2.801 23
## Smoljan 28 -30.6 9.8 5.659 36
## Sofia region 32 -14.2 11.8 3.293 54
## Stara Zagora 56 -19.6 19.1 2.083 66
## Targovishte 35 -30.3 52.7 10.452 33
## Varna 115 -6.3 18.4 1.904 26
## Veliko Tarnovo 43 -28.8 13.9 5.545 28
## Vidin 23 -41.1 9.6 11.928 13
## Vraca 40 -31.6 10.0 4.561 60
## Jambol 31 -28.9 17.3 8.961 29
## Dobrich 31 -29.6 27.4 12.360 25
## Gabrovo 46 -31.1 9.1 3.843 44
## Sofia capitale 965 8.1 16.9 0.092 13
## Haskovo 37 -23.1 27.2 7.257 22
## Kardzhali 47 -13.6 73.5 9.823 42
## Kjustendil 35 -33.2 9.0 6.365 27
## Lovech 27 -31.7 11.1 6.698 32
## PIBhab_23 PopRural_24 TxRisq_24 FUSION_ViK
## Blagoevgrad 16449 39.8 26 Petric
## Montana 15127 36.3 24
## Pazardzhik 16057 38.3 32 Panagjurishte, Strelcha
## Pernik 15608 21.7 37
## Pleven 17703 32.5 32 Knezha
## Plovdiv 22209 25.2 35
## Razgrad 17081 52.6 48 Kurbat
## Ruse 21678 22.0 23
## Shumen 17129 38.1 19
## Silistra 13114 55.5 27
## Sliven 13493 35.0 40
## Burgas 22169 23.4 30
## Smoljan 17016 42.5 16
## Sofia region 24251 39.8 23 Botevgrad
## Stara Zagora 36470 28.6 37
## Targovishte 18802 43.9 33
## Varna 26946 15.5 29
## Veliko Tarnovo 20187 29.9 34 Svishtov
## Vidin 15071 36.8 37
## Vraca 27113 40.7 30
## Jambol 17865 30.1 27
## Dobrich 17075 30.2 27
## Gabrovo 24297 20.6 25 Sevlievo
## Sofia capitale 61833 4.9 25
## Haskovo 13012 28.5 25 Dimitrovgrad + Stambolovo
## Kardzhali 15593 59.3 30
## Kjustendil 15030 31.3 25 Dupnica
## Lovech 18024 38.3 24
## REV_VIK_21 EVO_REVcon_VIK_pourc_13.21 EVO_REVcon_VIK_31.21
## Blagoevgrad 24782754 75.4 9118111
## Montana 12392993 33.5 2664439
## Pazardzhik 10761841 37.3 2501905
## Pernik 12066065 24.6 2036535
## Pleven 27415887 25.9 4835506
## Plovdiv 64806330 31.8 13395872
## Razgrad 9907945 10.7 820830
## Ruse 27851183 10.5 2275191
## Shumen 20580338 -5.5 -1016848
## Silistra 12022800 22.0 1855493
## Sliven 15631200 26.4 2794877
## Burgas 92973600 71.6 33216685
## Smoljan 8781623 35.9 1987513
## Sofia region 15990648 18.3 2117623
## Stara Zagora 37340883 8.0 2362934
## Targovishte 9997000 -6.4 -589925
## Varna 69827463 41.9 17642787
## Veliko Tarnovo 25490400 46.0 6871973
## Vidin 7185600 9.3 521055
## Vraca 19227600 19.2 2654986
## Jambol 11711968 17.7 1511370
## Dobrich 23638800 20.2 3400699
## Gabrovo 11879122 60.0 3814481
## Sofia capitale 196243690 36.8 45179858
## Haskovo 19172873 25.0 3278131
## Kardzhali 11143000 46.2 3016240
## Kjustendil 7603100 44.8 2015503
## Lovech 11235860 2.9 269743
## Km_reseau_eau_21 dens_pop_reseau21 Rev.VIK_milliers.KM_21
## Blagoevgrad 1770 150 14.0
## Montana 1824 61 6.8
## Pazardzhik 1150 155 9.4
## Pernik 1481 82 8.1
## Pleven 3255 74 8.4
## Plovdiv 4344 154 14.9
## Razgrad 1428 60 6.9
## Ruse 2629 82 10.6
## Shumen 2367 73 8.7
## Silistra 1861 59 6.5
## Sliven 1857 100 8.4
## Burgas 4304 95 21.6
## Smoljan 1512 70 5.8
## Sofia region 3632 62 4.4
## Stara Zagora 3475 90 10.7
## Targovishte 2212 50 4.5
## Varna 4334 109 16.1
## Veliko Tarnovo 2299 103 11.1
## Vidin 1684 50 4.3
## Vraca 2865 56 6.7
## Jambol 1785 67 6.6
## Dobrich 3404 51 6.9
## Gabrovo 1958 55 6.1
## Sofia capitale 3699 359 53.1
## Haskovo 3193 73 6.0
## Kardzhali 1628 89 6.8
## Kjustendil 1413 91 5.4
## Lovech 1803 52 6.2
## POP_VIK_21 EVO_POP_VIK_13.21 PEVO_POP_VIK_13.21 REV_POP_VIK_21
## Blagoevgrad 265892 55768 26.5 93
## Montana 111151 -6552 -5.6 111
## Pazardzhik 178564 32474 22.2 60
## Pernik 120880 -5771 -4.6 100
## Pleven 240380 -9165 -3.7 114
## Plovdiv 668334 -14693 -2.2 97
## Razgrad 85825 -21508 -20.0 115
## Ruse 215865 -11650 -5.1 129
## Shumen 172367 -27459 -13.7 119
## Silistra 109271 -7925 -6.8 110
## Sliven 185191 -10503 -5.4 84
## Burgas 410331 -4732 -1.1 227
## Smoljan 105421 -3384 -3.1 83
## Sofia region 226385 13767 6.5 71
## Stara Zagora 313396 -50287 -13.8 119
## Targovishte 111597 26083 30.5 90
## Varna 471252 -38558 -7.6 148
## Veliko Tarnovo 237628 26029 12.3 107
## Vidin 84865 -68201 -44.6 85
## Vraca 160829 -32143 -16.7 120
## Jambol 118897 -9535 -7.4 99
## Dobrich 173657 -12522 -6.7 136
## Gabrovo 108481 22157 25.7 110
## Sofia capitale 1328790 19156 1.5 148
## Haskovo 234354 67725 40.6 82
## Kardzhali 145218 9640 7.1 77
## Kjustendil 128673 45072 53.9 59
## Lovech 94650 -10230 -9.8 119
## EVO_REVcon_POP_13.21 PRIX_ViK_23 PEVO_PRIX_cou_13.23
## Blagoevgrad 22.2 2.2 38.5
## Montana 28.0 3.4 42.8
## Pazardzhik 5.7 1.4 5.1
## Pernik 20.0 2.0 25.9
## Pleven 23.0 3.7 53.5
## Plovdiv 21.4 2.5 46.4
## Razgrad 27.5 4.5 50.1
## Ruse 15.7 3.7 49.7
## Shumen 9.0 2.9 22.9
## Silistra 22.2 4.4 54.7
## Sliven 18.2 2.0 14.2
## Burgas 82.2 3.5 48.3
## Smoljan 20.5 3.9 51.0
## Sofia region 6.0 2.1 12.1
## Stara Zagora 20.6 3.8 40.9
## Targovishte -30.3 3.8 29.5
## Varna 44.2 4.2 52.3
## Veliko Tarnovo 21.2 3.1 52.7
## Vidin 35.7 2.9 39.8
## Vraca 30.8 3.3 41.5
## Jambol 18.0 3.7 55.7
## Dobrich 26.1 2.5 5.5
## Gabrovo 20.1 2.9 37.2
## Sofia capitale 32.6 2.4 39.9
## Haskovo -8.8 3.6 31.4
## Kardzhali 17.6 3.2 57.0
## Kjustendil -3.2 1.7 0.0
## Lovech 12.5 3.4 40.8
## EVO_PRIX_cou_13.23 EAUprel_21 EAUprelHab_21
## Blagoevgrad 0.85 27368811 103
## Montana 1.48 19510132 176
## Pazardzhik 0.07 19844284 111
## Pernik 0.52 24727882 205
## Pleven 1.99 27109307 113
## Plovdiv 1.18 80112632 120
## Razgrad 2.25 11059338 129
## Ruse 1.84 18441714 85
## Shumen 0.66 36739201 213
## Silistra 2.42 10280450 94
## Sliven 0.28 39440815 213
## Burgas 1.69 60083418 146
## Smoljan 1.96 7924761 75
## Sofia region 0.26 23889253 106
## Stara Zagora 1.57 31325423 100
## Targovishte 1.11 12009848 108
## Varna 2.18 65314429 139
## Veliko Tarnovo 1.64 27424738 115
## Vidin 1.16 6102371 72
## Vraca 1.36 18876730 117
## Jambol 2.05 20653219 174
## Dobrich 0.14 40835687 235
## Gabrovo 1.08 12564286 116
## Sofia capitale 0.98 142388515 107
## Haskovo 1.14 17299143 74
## Kardzhali 1.83 10237822 70
## Kjustendil 0.00 13299419 103
## Lovech 1.37 9049870 96
## Perteseau_Pourc_VIK21 Energie.kwh.m3_21ViK Continuite_21_VIK
## Blagoevgrad 58 0.06 1.85
## Montana 76 0.50 15.28
## Pazardzhik 61 0.38 0.06
## Pernik 80 0.29 5.43
## Pleven 62 0.81 1.27
## Plovdiv 58 0.60 2.41
## Razgrad 69 1.53 0.16
## Ruse 45 1.08 0.18
## Shumen 81 1.03 0.73
## Silistra 58 1.48 0.20
## Sliven 85 0.52 8.28
## Burgas 54 0.40 0.18
## Smoljan 45 0.59 0.59
## Sofia region 56 0.56 0.29
## Stara Zagora 57 1.08 3.53
## Targovishte 69 0.95 0.24
## Varna 64 0.46 0.78
## Veliko Tarnovo 58 0.56 1.38
## Vidin 47 0.78 4.79
## Vraca 60 0.58 0.12
## Jambol 76 0.67 2.75
## Dobrich 74 1.33 1.01
## Gabrovo 57 0.31 9.21
## Sofia capitale 40 0.04 0.83
## Haskovo 57 2.34 2.33
## Kardzhali 49 0.53 0.46
## Kjustendil 66 0.38 0.91
## Lovech 55 0.82 5.40
## Avaries100Km_VIK_21 Conso_eau_hab_23
## Blagoevgrad 86 125
## Montana 76 90
## Pazardzhik 79 115
## Pernik 143 111
## Pleven 71 100
## Plovdiv 81 102
## Razgrad 65 88
## Ruse 63 109
## Shumen 68 88
## Silistra 60 96
## Sliven 84 75
## Burgas 55 115
## Smoljan 72 84
## Sofia region 93 98
## Stara Zagora 87 95
## Targovishte 55 74
## Varna 79 110
## Veliko Tarnovo 69 105
## Vidin 127 92
## Vraca 67 103
## Jambol 147 91
## Dobrich 78 96
## Gabrovo 63 105
## Sofia capitale 95 126
## Haskovo 114 83
## Kardzhali 109 81
## Kjustendil 72 120
## Lovech 96 990.B.b.Création d’une base de données sans Sofia.
Pour de nombreux indicateurs Sofia ne représente pas une valeur extrême, mais lors de la réalisation de l’ACP la base de données (BD) sans Sofia permet de tester la pertinence des axes avec ou sans la ville de loin la plus peuplée, la plus riche et où le ViK a de très loin le plus d’abonné·es du pays
ViKnoSOF <- ViK[-24, ]
ViKnoSOF## REG_FR CODE clusCAH RevHab_24 Pop65_24 TxMig_14.24
## Blagoevgrad Blagoevgrad BGR.1_1 CLUS 4 25997 23 -20.0
## Montana Montana BGR.10_1 CLUS 1 26217 28 -28.2
## Pazardzhik Pazardzhik BGR.11_1 CLUS 1 24736 26 -17.9
## Pernik Pernik BGR.12_1 CLUS 1 30570 30 37.5
## Pleven Pleven BGR.13_1 CLUS 1 26815 31 -36.1
## Plovdiv Plovdiv BGR.14_1 CLUS 2 24742 26 49.8
## Razgrad Razgrad BGR.15_1 CLUS 4 21316 27 -23.3
## Ruse Ruse BGR.16_1 CLUS 1 22413 27 -5.2
## Shumen Shumen BGR.17_1 CLUS 4 27406 26 38.3
## Silistra Silistra BGR.18_1 CLUS 4 25220 29 -14.6
## Sliven Sliven BGR.19_1 CLUS 4 23765 23 -44.4
## Burgas Burgas BGR.2_1 CLUS 2 26769 29 69.0
## Smoljan Smoljan BGR.20_1 CLUS 1 29927 26 -87.8
## Sofia region Sofia region BGR.21_1 CLUS 3 18820 24 59.1
## Stara Zagora Stara Zagora BGR.22_1 CLUS 3 24704 26 6.7
## Targovishte Targovishte BGR.23_1 CLUS 4 23476 25 -3.5
## Varna Varna BGR.24_1 CLUS 2 23076 28 54.1
## Veliko Tarnovo Veliko Tarnovo BGR.25_1 CLUS 1 22549 25 -20.4
## Vidin Vidin BGR.26_1 CLUS 4 21539 23 -46.7
## Vraca Vraca BGR.27_1 CLUS 3 23540 23 -48.2
## Jambol Jambol BGR.28_1 CLUS 1 17821 27 -32.2
## Dobrich Dobrich BGR.3_1 CLUS 1 20742 28 -11.0
## Gabrovo Gabrovo BGR.4_1 CLUS 1 28608 31 -10.0
## Haskovo Haskovo BGR.6_1 CLUS 1 20292 27 8.5
## Kardzhali Kardzhali BGR.7_1 CLUS 6 21024 25 188.4
## Kjustendil Kjustendil BGR.8_1 CLUS 1 23491 28 -18.2
## Lovech Lovech BGR.9_1 CLUS 4 17555 22 -19.6
## DensPop_24 PEVO_POP_01.21 NonBG_21 AgriVA_23 IndusVA_23
## Blagoevgrad 44 -11.4 18.9 5.5 25
## Montana 31 -31.6 17.0 11.6 24
## Pazardzhik 50 -17.8 19.6 6.0 39
## Pernik 46 -23.0 5.9 3.5 31
## Pleven 46 -30.5 11.5 5.8 28
## Plovdiv 107 -11.2 19.1 2.4 35
## Razgrad 41 -29.7 62.3 11.6 29
## Ruse 65 -27.5 23.1 5.2 31
## Shumen 44 -25.5 45.9 10.2 32
## Silistra 33 -30.6 47.2 16.5 16
## Sliven 48 -20.5 33.1 6.7 24
## Burgas 50 -10.0 26.3 2.8 23
## Smoljan 28 -30.6 9.8 5.7 36
## Sofia region 32 -14.2 11.8 3.3 54
## Stara Zagora 56 -19.6 19.1 2.1 66
## Targovishte 35 -30.3 52.7 10.5 33
## Varna 115 -6.3 18.4 1.9 26
## Veliko Tarnovo 43 -28.8 13.9 5.5 28
## Vidin 23 -41.1 9.6 11.9 13
## Vraca 40 -31.6 10.0 4.6 60
## Jambol 31 -28.9 17.3 9.0 29
## Dobrich 31 -29.6 27.4 12.4 25
## Gabrovo 46 -31.1 9.1 3.8 44
## Haskovo 37 -23.1 27.2 7.3 22
## Kardzhali 47 -13.6 73.5 9.8 42
## Kjustendil 35 -33.2 9.0 6.4 27
## Lovech 27 -31.7 11.1 6.7 32
## PIBhab_23 PopRural_24 TxRisq_24 FUSION_ViK
## Blagoevgrad 16449 40 26 Petric
## Montana 15127 36 24
## Pazardzhik 16057 38 32 Panagjurishte, Strelcha
## Pernik 15608 22 37
## Pleven 17703 33 32 Knezha
## Plovdiv 22209 25 35
## Razgrad 17081 53 48 Kurbat
## Ruse 21678 22 23
## Shumen 17129 38 19
## Silistra 13114 55 27
## Sliven 13493 35 40
## Burgas 22169 23 30
## Smoljan 17016 42 16
## Sofia region 24251 40 23 Botevgrad
## Stara Zagora 36470 29 37
## Targovishte 18802 44 33
## Varna 26946 15 29
## Veliko Tarnovo 20187 30 34 Svishtov
## Vidin 15071 37 37
## Vraca 27113 41 30
## Jambol 17865 30 27
## Dobrich 17075 30 27
## Gabrovo 24297 21 25 Sevlievo
## Haskovo 13012 28 25 Dimitrovgrad + Stambolovo
## Kardzhali 15593 59 30
## Kjustendil 15030 31 25 Dupnica
## Lovech 18024 38 24
## REV_VIK_21 EVO_REVcon_VIK_pourc_13.21 EVO_REVcon_VIK_31.21
## Blagoevgrad 24782754 75.4 9118111
## Montana 12392993 33.5 2664439
## Pazardzhik 10761841 37.3 2501905
## Pernik 12066065 24.6 2036535
## Pleven 27415887 25.9 4835506
## Plovdiv 64806330 31.8 13395872
## Razgrad 9907945 10.7 820830
## Ruse 27851183 10.5 2275191
## Shumen 20580338 -5.5 -1016848
## Silistra 12022800 22.0 1855493
## Sliven 15631200 26.4 2794877
## Burgas 92973600 71.6 33216685
## Smoljan 8781623 35.9 1987513
## Sofia region 15990648 18.3 2117623
## Stara Zagora 37340883 8.0 2362934
## Targovishte 9997000 -6.4 -589925
## Varna 69827463 41.9 17642787
## Veliko Tarnovo 25490400 46.0 6871973
## Vidin 7185600 9.3 521055
## Vraca 19227600 19.2 2654986
## Jambol 11711968 17.7 1511370
## Dobrich 23638800 20.2 3400699
## Gabrovo 11879122 60.0 3814481
## Haskovo 19172873 25.0 3278131
## Kardzhali 11143000 46.2 3016240
## Kjustendil 7603100 44.8 2015503
## Lovech 11235860 2.9 269743
## Km_reseau_eau_21 dens_pop_reseau21 Rev.VIK_milliers.KM_21
## Blagoevgrad 1770 150 14.0
## Montana 1824 61 6.8
## Pazardzhik 1150 155 9.4
## Pernik 1481 82 8.1
## Pleven 3255 74 8.4
## Plovdiv 4344 154 14.9
## Razgrad 1428 60 6.9
## Ruse 2629 82 10.6
## Shumen 2367 73 8.7
## Silistra 1861 59 6.5
## Sliven 1857 100 8.4
## Burgas 4304 95 21.6
## Smoljan 1512 70 5.8
## Sofia region 3632 62 4.4
## Stara Zagora 3475 90 10.7
## Targovishte 2212 50 4.5
## Varna 4334 109 16.1
## Veliko Tarnovo 2299 103 11.1
## Vidin 1684 50 4.3
## Vraca 2865 56 6.7
## Jambol 1785 67 6.6
## Dobrich 3404 51 6.9
## Gabrovo 1958 55 6.1
## Haskovo 3193 73 6.0
## Kardzhali 1628 89 6.8
## Kjustendil 1413 91 5.4
## Lovech 1803 52 6.2
## POP_VIK_21 EVO_POP_VIK_13.21 PEVO_POP_VIK_13.21 REV_POP_VIK_21
## Blagoevgrad 265892 55768 26.5 93
## Montana 111151 -6552 -5.6 111
## Pazardzhik 178564 32474 22.2 60
## Pernik 120880 -5771 -4.6 100
## Pleven 240380 -9165 -3.7 114
## Plovdiv 668334 -14693 -2.2 97
## Razgrad 85825 -21508 -20.0 115
## Ruse 215865 -11650 -5.1 129
## Shumen 172367 -27459 -13.7 119
## Silistra 109271 -7925 -6.8 110
## Sliven 185191 -10503 -5.4 84
## Burgas 410331 -4732 -1.1 227
## Smoljan 105421 -3384 -3.1 83
## Sofia region 226385 13767 6.5 71
## Stara Zagora 313396 -50287 -13.8 119
## Targovishte 111597 26083 30.5 90
## Varna 471252 -38558 -7.6 148
## Veliko Tarnovo 237628 26029 12.3 107
## Vidin 84865 -68201 -44.6 85
## Vraca 160829 -32143 -16.7 120
## Jambol 118897 -9535 -7.4 99
## Dobrich 173657 -12522 -6.7 136
## Gabrovo 108481 22157 25.7 110
## Haskovo 234354 67725 40.6 82
## Kardzhali 145218 9640 7.1 77
## Kjustendil 128673 45072 53.9 59
## Lovech 94650 -10230 -9.8 119
## EVO_REVcon_POP_13.21 PRIX_ViK_23 PEVO_PRIX_cou_13.23
## Blagoevgrad 22.2 2.2 38.5
## Montana 28.0 3.4 42.8
## Pazardzhik 5.7 1.4 5.1
## Pernik 20.0 2.0 25.9
## Pleven 23.0 3.7 53.5
## Plovdiv 21.4 2.5 46.4
## Razgrad 27.5 4.5 50.1
## Ruse 15.7 3.7 49.7
## Shumen 9.0 2.9 22.9
## Silistra 22.2 4.4 54.7
## Sliven 18.2 2.0 14.2
## Burgas 82.2 3.5 48.3
## Smoljan 20.5 3.9 51.0
## Sofia region 6.0 2.1 12.1
## Stara Zagora 20.6 3.8 40.9
## Targovishte -30.3 3.8 29.5
## Varna 44.2 4.2 52.3
## Veliko Tarnovo 21.2 3.1 52.7
## Vidin 35.7 2.9 39.8
## Vraca 30.8 3.3 41.5
## Jambol 18.0 3.7 55.7
## Dobrich 26.1 2.5 5.5
## Gabrovo 20.1 2.9 37.2
## Haskovo -8.8 3.6 31.4
## Kardzhali 17.6 3.2 57.0
## Kjustendil -3.2 1.7 0.0
## Lovech 12.5 3.4 40.8
## EVO_PRIX_cou_13.23 EAUprel_21 EAUprelHab_21
## Blagoevgrad 0.85 27368811 103
## Montana 1.48 19510132 176
## Pazardzhik 0.07 19844284 111
## Pernik 0.52 24727882 205
## Pleven 1.99 27109307 113
## Plovdiv 1.18 80112632 120
## Razgrad 2.25 11059338 129
## Ruse 1.84 18441714 85
## Shumen 0.66 36739201 213
## Silistra 2.42 10280450 94
## Sliven 0.28 39440815 213
## Burgas 1.69 60083418 146
## Smoljan 1.96 7924761 75
## Sofia region 0.26 23889253 106
## Stara Zagora 1.57 31325423 100
## Targovishte 1.11 12009848 108
## Varna 2.18 65314429 139
## Veliko Tarnovo 1.64 27424738 115
## Vidin 1.16 6102371 72
## Vraca 1.36 18876730 117
## Jambol 2.05 20653219 174
## Dobrich 0.14 40835687 235
## Gabrovo 1.08 12564286 116
## Haskovo 1.14 17299143 74
## Kardzhali 1.83 10237822 70
## Kjustendil 0.00 13299419 103
## Lovech 1.37 9049870 96
## Perteseau_Pourc_VIK21 Energie.kwh.m3_21ViK Continuite_21_VIK
## Blagoevgrad 58 0.06 1.85
## Montana 76 0.50 15.28
## Pazardzhik 61 0.38 0.06
## Pernik 80 0.29 5.43
## Pleven 62 0.81 1.27
## Plovdiv 58 0.60 2.41
## Razgrad 69 1.53 0.16
## Ruse 45 1.08 0.18
## Shumen 81 1.03 0.73
## Silistra 58 1.48 0.20
## Sliven 85 0.52 8.28
## Burgas 54 0.40 0.18
## Smoljan 45 0.59 0.59
## Sofia region 56 0.56 0.29
## Stara Zagora 57 1.08 3.53
## Targovishte 69 0.95 0.24
## Varna 64 0.46 0.78
## Veliko Tarnovo 58 0.56 1.38
## Vidin 47 0.78 4.79
## Vraca 60 0.58 0.12
## Jambol 76 0.67 2.75
## Dobrich 74 1.33 1.01
## Gabrovo 57 0.31 9.21
## Haskovo 57 2.34 2.33
## Kardzhali 49 0.53 0.46
## Kjustendil 66 0.38 0.91
## Lovech 55 0.82 5.40
## Avaries100Km_VIK_21 Conso_eau_hab_23
## Blagoevgrad 86 125
## Montana 76 90
## Pazardzhik 79 115
## Pernik 143 111
## Pleven 71 100
## Plovdiv 81 102
## Razgrad 65 88
## Ruse 63 109
## Shumen 68 88
## Silistra 60 96
## Sliven 84 75
## Burgas 55 115
## Smoljan 72 84
## Sofia region 93 98
## Stara Zagora 87 95
## Targovishte 55 74
## Varna 79 110
## Veliko Tarnovo 69 105
## Vidin 127 92
## Vraca 67 103
## Jambol 147 91
## Dobrich 78 96
## Gabrovo 63 105
## Haskovo 114 83
## Kardzhali 109 81
## Kjustendil 72 120
## Lovech 96 990.B.b. métadonnées
META<-read.csv2("C:/Users/Raph/Documents/thesebulgarie/STATISTIQUES/analysesR/VIK_INEG/Metadonnees.csv",stringsAsFactors=FALSE, fileEncoding="latin1")
META## Titre.variable
## 1 REG_FR
## 2 CODE
## 3 clusCAH
## 4 RevHab_24
## 5 Pop65_24
## 6 TxMig_14-24
## 7 DensPop_24
## 8 PEVO_POP_01-21
## 9 NonBG_21
## 10 AgriVA_23
## 11 IndusVA_23
## 12 PIBhab_23
## 13 PopRural_24
## 14 TxRisq_24
## 15 FUSION_ViK
## 16 REV_VIK_21
## 17 EVO_REVcon_VIK_pourc_13-21
## 18 EVO_REVcon_VIK_31-21
## 19 Km_reseau_eau_21
## 20 dens_pop_reseau21
## 21 Rev.VIK_milliers.KM_21
## 22 POP_VIK_21
## 23 EVO_POP_VIK_13-21
## 24 PEVO_POP_VIK_13-21
## 25 REV_POP_VIK_21
## 26 EVO_REVcon_POP_13-21
## 27 PRIX_ViK_23
## 28 PEVO_PRIX_cou_13-23
## 29 EVO_PRIX_cou_13-23
## 30 EAUprel_21
## 31 EAUprelHab_21
## 32 Perteseau_Pourc_VIK21
## 33 Energie-kwh.m3_21ViK
## 34 Continuite_21_VIK
## 35 Avaries100Km_VIK_21
## 36 Conso_eau_hab_23
## Description
## 1 Nom de la région en français
## 2 Code de la région par rapport aux couches GADM
## 3 Groupe de région auquel appartient le ViK étudié selon la classification ascendante hiérarchique réalisée dans la présentation des terrains de l'enquête
## 4 Revenu moyen de la population par habitant dans la région
## 5 Population de plus de 65 ans
## 6 Taux d'accroissement migratoire entre 2014 et 2024\nCorrespond à l'évolution de la population sur une période donnée imputable aux mouvements migratoires. Il est exprimé pour 1000 habitants. (Solde migratoire pour 1000 hab)/(Population moyenne de la région)\nDepuis 2007 l\u0092accroissement migratoire inclue non seulement les migrations intérieures mais aussi internationales.
## 7 Densité de la population dans les différents oblasts
## 8 Évolution de la population entre 2001 et 2021
## 9 % de population se définissant comme non bulgare dans le recensement de 2021, soit Turque, Rom et Autre
## 10 Part de l'agriculture dans la valeur ajouté
## 11 Part de l'industrie dans la valeur ajouté
## 12 PIB par habitant de la région
## 13 Pourcentage de population vivant en zone rurale
## 14 Taux de risque de pauvreté et d'exclusion sociale. \nPersonnes qui sont menacées par la pauvreté et/ou vivent dans le dénuement matériel et/ou vivent dans des ménages à très faible intensité de travail. \n\u0095 Les personnes qui sont menacées par la pauvreté sont les personnes dont le revenu disponible équivalent se situe en dessous du seuil de pauvreté fixé à 60 % du revenu disponible équivalent médian national (après transferts sociaux). \n\u0095 Les personnes en situation de privation matérielle grave ont des conditions de vie fortement affectées par le manque de ressources; quatre des neuf indicateurs de privation s\u0092appliquent à elles: impossibilité i) de régler le loyer ou les factures pour les services d\u0092utilité publique, ii) de chauffer convenablement le domicile, iii) de faire face à des dépenses imprévues, iv) de consommer de la viande, du poisson ou un équivalent protéiné tous les deux jours, v) de partir en vacances hors du domicile une semaine par an, vi) d\u0092acheter une voiture, vii) d\u0092acheter une machine à laver le linge, viii) d\u0092acheter une télévision couleur ou ix) de payer une connexion téléphonique. \n\u0095 Les personnes vivant dans des ménages à très faible intensité de travail sont les personnes âgées de 0 à 59 ans vivant dans des ménages dans lesquels les adultes ont travaillé 20 % ou moins de leur potentiel de travail total au cours de l\u0092année écoulée. \n\n
## 15 Fusions d'opérateurs de réseau d'eau potable entre 2013 et 2021 : permet d'expliquer l'augmeentation de la population de certains ViKs
## 16 Revenu gagné par les ViKs. Attention : ce n'est pas le bénéfice net.
## 17 Évolution du revenu des ViKs en leva constants
## 18 Évolution du revenu des ViKs
## 19 Kilomètres de réseau d'eau de chaque ViK
## 20 Densité de population sur le réseau
## 21 Revenu des ViKs pour 1000 km de réseau
## 22 Population desservie par les opérateurs Viks
## 23 Évolution de la population desservie par les opérateurs ViKs
## 24 Évolution de la population desservie par les opérateurs ViKs
## 25 Revenu généré par les ViKs par rapport à la population desservie en 2021
## 26 Évolution du revenu par habitant des ViKs entre 2013 et 2021
## 27 Prix de l'eau des ViKs
## 28 Évolution des prix de l'eau entre 2013 et 2023
## 29 Évolution du prix de l'eau entre 2013 et 2023
## 30 Eau prélevée à la source par les opérateurs ViKs
## 31 Eau prélevée à la source par rapport à la population desservie
## 32 Pourcentage de pertes en eau (pertes physiques + commerciales) par rapport à l'eau prélevée
## 33 Consommation énergétique des ViK en KwH par mètre cube d'eau produit
## 34 Ratio de continuité du service de distribution de l'eau. Il est déterminé par le rapport entre 1) le produit du nombre total de personnes touchées par des interruptions de l'approvisionnement en eau pendant la durée correspondante (en heures) et 2) le produit du nombre total de personnes approvisionnées en eau, multiplié par 24 et par le nombre de jours de la période considérée.
## 35 Nombre de pannes pour 100km de réseau
## 36 Consommation en eau du réseau en litres par jour et par habitant en 2021. Indique la consommation domestique moyenne de la population. Ne prend pas en compte les consommations industrielles. Est mesurée au compteur d'eau.
## Unité.de.mesure Années
## 1 Qualitative
## 2 Qualitative
## 3 qualitative 2025
## 4 leva courant/hab. 2024 (dernière année disponible)
## 5 % 2024 (dernière année disponible)
## 6 \u0089 2014-2024 (à partir de la dernière année disponible)
## 7 hab./km2 2024 (dernière année disponible)
## 8 % 2001-2021 (entre les deux recensements sur 20 ans)
## 9 % 2021 (dernier recensement)
## 10 % 2023 (dernière année disponible)
## 11 % 2023 (dernière année disponible)
## 12 leva/hab. 2023 (dernière année disponible)
## 13 % 2024 (dernière année disponible)
## 14 % 2024 (dernière année disponible)
## 15 Qualitative 2013-2021
## 16 leva courant 2021
## 17 % 2013-2021 (laps de temps le plus important disponible)
## 18 leva constants 2013-2021 (laps de temps le plus important disponible)
## 19 km 2021
## 20 hab/km 2021
## 21 leva courant/1000km 2021
## 22 hab. 2021
## 23 hab. 2013-2021
## 24 % 2013-2021
## 25 leva courant/hab. 2021
## 26 leva constant/hab 2013-2021
## 27 leva courant 2023
## 28 % (leva courant) 2013-2023
## 29 leva courant 2013-2023
## 30 m3 2021
## 31 m3/hab 2021
## 32 % 2021
## 33 KwH/m3 2021
## 34 Ratio 2021
## 35 pannes/100km 2021
## 36 L/hab./J 2023
## Source
## 1
## 2 GADM
## 3 R.Segond
## 4 INS
## 5 INS
## 6 INS
## 7 INS
## 8 INS
## 9 INS
## 10 INS
## 11 INS
## 12 INS
## 13 INS
## 14 INS
## 15 KEVR
## 16 KEVR
## 17 KEVR
## 18 KEVR
## 19 KEVR
## 20 KEVR
## 21 KEVR
## 22 KEVR
## 23 KEVR
## 24 KEVR
## 25 KEVR
## 26 KEVR
## 27 KEVR
## 28 KEVR
## 29 KEVR
## 30 KEVR
## 31 KEVR
## 32 KEVR
## 33 KEVR
## 34 KEVR
## 35 KEVR
## 36 INS
## Détail.source
## 1 <NA>
## 2 <NA>
## 3 R.Segond à partir d'un ensemble de données statistiques. Réalisation d'une CAH
## 4 Enquête budget des ménages
## 5 Enquête population
## 6 Enquête population
## 7 Enquête population
## 8 Recensements (calculs de l'autrice)
## 9 Recensements (calculs de l'autrice)
## 10 PIB et VAB par secteur économique et région (calculs de l'autrice)
## 11 PIB et VAB par secteur économique et région (calculs de l'autrice)
## 12 PIB et VAB par secteur économique et région + Enquête population (calculs de l'autrice)
## 13 Enquête
## 14 Enquête SILS
## 15 Registres du KEVR (compilation de l'autrice)
## 16 Rapports annuels du KEVR
## 17 Rapports annuels du KEVR
## 18 Rapports annuels du KEVR
## 19 Business Plans des ViKs (compilation de l'autricele)
## 20 Rapports annuels du KEVR et business plans (calcul de l'autrice)
## 21 Calcul de l'autrice à partir des business plans des ViKs et des rapports annuels du KEVR
## 22 Rapports annuels du KEVR
## 23 Rapports annuels du KEVR
## 24 Calcul de l'autrice à partir des rapports annuels du KEVR
## 25 Calcul de l'autrice à partir des rapports annuels du KEVR
## 26 Calcul de l'autrice à partir des rapports annuels du KEVR
## 27 Registre des prix 2023
## 28 Registre des prix 2013 et 2023 (calcul de l'autrice)
## 29 Registre des prix 2013 et 2023 (calcul de l'autrice)
## 30 Rapports annuels du KEVR
## 31 Rapports annuels du KEVR (calcul de l'autrice)
## 32 Rapport d'évaluation en ligne du KEVR (compilation de l'autrice)
## 33 Rapport d'évaluation en ligne du KEVR (compilation de l'autrice). Selon le KEVR données peu fiables pour : Stara Zagora, Sliven et Dobric \n
## 34 Rapport d'évaluation en ligne du KEVR (compilation de l'autrice). Selon le KEVR données peu fiables pour : Vraca, Plovdiv, Montana, Sliven, Burgas, Gabrovo
## 35 Rapport d'évaluation en ligne du KEVR (compilation de l'autrice). Selon le KEVR données peu fiables pour : Vraca
## 36 Enquête Environnementc. matrices de corrélations
J’enlève les variables qualitatives et j’utilise le coefficient de Spearman car au regard des différents nuages de points, la majorité des variables sont non linéaires
ViKmatrice<-ViK[,c(4:14,16:36)]
Matrice_globale_cor<- cor(ViKmatrice, method=c("spearman"))
Matrice_globale_cor## RevHab_24 Pop65_24 TxMig_14.24 DensPop_24
## RevHab_24 1.000 0.241 0.0213 0.379
## Pop65_24 0.241 1.000 0.1177 0.056
## TxMig_14.24 0.021 0.118 1.0000 0.517
## DensPop_24 0.379 0.056 0.5167 1.000
## PEVO_POP_01.21 0.135 -0.145 0.6844 0.737
## NonBG_21 -0.229 -0.032 0.3191 0.224
## AgriVA_23 -0.339 0.077 -0.4007 -0.667
## IndusVA_23 -0.026 -0.101 0.1232 0.110
## PIBhab_23 0.090 -0.118 0.3098 0.441
## PopRural_24 -0.238 -0.326 -0.3372 -0.529
## TxRisq_24 -0.045 -0.024 -0.1338 0.248
## REV_VIK_21 0.243 0.063 0.4313 0.643
## EVO_REVcon_VIK_pourc_13.21 0.367 0.158 0.1297 0.363
## EVO_REVcon_VIK_31.21 0.273 0.094 0.2983 0.608
## Km_reseau_eau_21 0.036 0.052 0.4811 0.410
## dens_pop_reseau21 0.319 -0.114 0.3273 0.755
## Rev.VIK_milliers.KM_21 0.362 0.032 0.3311 0.790
## POP_VIK_21 0.171 -0.072 0.4899 0.719
## EVO_POP_VIK_13.21 0.129 -0.030 0.1319 -0.034
## PEVO_POP_VIK_13.21 0.144 0.012 0.2923 0.177
## REV_POP_VIK_21 0.181 0.224 0.1544 0.296
## EVO_REVcon_POP_13.21 0.283 0.118 -0.1494 0.147
## PRIX_ViK_23 -0.121 0.273 -0.1029 -0.081
## PEVO_PRIX_cou_13.23 -0.045 0.160 -0.0991 0.102
## EVO_PRIX_cou_13.23 -0.097 0.244 -0.1363 0.033
## EAUprel_21 0.286 0.068 0.3842 0.604
## EAUprelHab_21 0.219 0.368 0.0044 0.137
## Perteseau_Pourc_VIK21 0.039 0.350 -0.2157 -0.136
## Energie.kwh.m3_21ViK -0.475 0.091 -0.1210 -0.343
## Continuite_21_VIK 0.078 0.014 -0.1762 -0.149
## Avaries100Km_VIK_21 -0.297 -0.369 0.0635 -0.105
## Conso_eau_hab_23 0.327 0.131 0.1922 0.418
## PEVO_POP_01.21 NonBG_21 AgriVA_23 IndusVA_23
## RevHab_24 0.1352 -0.229 -0.3394 -0.0257
## Pop65_24 -0.1448 -0.032 0.0772 -0.1007
## TxMig_14.24 0.6844 0.319 -0.4007 0.1232
## DensPop_24 0.7369 0.224 -0.6667 0.1100
## PEVO_POP_01.21 1.0000 0.370 -0.5686 0.0025
## NonBG_21 0.3696 1.000 0.3897 -0.1062
## AgriVA_23 -0.5686 0.390 1.0000 -0.3005
## IndusVA_23 0.0025 -0.106 -0.3005 1.0000
## PIBhab_23 0.3129 -0.173 -0.7126 0.4346
## PopRural_24 -0.3498 0.306 0.5928 0.2326
## TxRisq_24 0.1246 0.146 0.0079 -0.0567
## REV_VIK_21 0.6296 0.116 -0.6152 -0.1418
## EVO_REVcon_VIK_pourc_13.21 0.3430 -0.167 -0.3454 -0.1790
## EVO_REVcon_VIK_31.21 0.5782 -0.033 -0.5605 -0.2271
## Km_reseau_eau_21 0.4440 0.130 -0.4915 -0.0553
## dens_pop_reseau21 0.7479 0.063 -0.5950 -0.0843
## Rev.VIK_milliers.KM_21 0.7109 0.234 -0.5260 -0.1303
## POP_VIK_21 0.7889 0.138 -0.6661 -0.0892
## EVO_POP_VIK_13.21 0.1087 -0.077 -0.0405 -0.0640
## PEVO_POP_VIK_13.21 0.2439 -0.042 -0.1949 0.0186
## REV_POP_VIK_21 0.1177 0.102 -0.2419 -0.1653
## EVO_REVcon_POP_13.21 0.0808 -0.120 -0.1757 -0.3996
## PRIX_ViK_23 -0.1911 0.278 0.1401 0.0060
## PEVO_PRIX_cou_13.23 0.0183 0.124 -0.0192 -0.0383
## EVO_PRIX_cou_13.23 -0.0819 0.185 0.0525 -0.0575
## EAUprel_21 0.6953 0.110 -0.5287 -0.1593
## EAUprelHab_21 0.1385 0.071 -0.0203 -0.0575
## Perteseau_Pourc_VIK21 -0.1361 0.177 0.3826 -0.0761
## Energie.kwh.m3_21ViK -0.3160 0.466 0.5071 0.0120
## Continuite_21_VIK -0.1708 -0.398 0.0027 -0.1959
## Avaries100Km_VIK_21 0.2357 -0.236 -0.0783 -0.0755
## Conso_eau_hab_23 0.2680 -0.462 -0.6219 -0.0953
## PIBhab_23 PopRural_24 TxRisq_24 REV_VIK_21
## RevHab_24 0.090 -0.2381 -0.0446 0.2430
## Pop65_24 -0.118 -0.3262 -0.0241 0.0635
## TxMig_14.24 0.310 -0.3372 -0.1338 0.4313
## DensPop_24 0.441 -0.5287 0.2483 0.6431
## PEVO_POP_01.21 0.313 -0.3498 0.1246 0.6296
## NonBG_21 -0.173 0.3060 0.1464 0.1160
## AgriVA_23 -0.713 0.5928 0.0079 -0.6152
## IndusVA_23 0.435 0.2326 -0.0567 -0.1418
## PIBhab_23 1.000 -0.3815 -0.0400 0.5654
## PopRural_24 -0.381 1.0000 -0.0071 -0.5988
## TxRisq_24 -0.040 -0.0071 1.0000 -0.0208
## REV_VIK_21 0.565 -0.5988 -0.0208 1.0000
## EVO_REVcon_VIK_pourc_13.21 -0.066 -0.2485 -0.0495 0.1888
## EVO_REVcon_VIK_31.21 0.310 -0.5463 0.0862 0.7121
## Km_reseau_eau_21 0.638 -0.4778 -0.0736 0.8385
## dens_pop_reseau21 0.134 -0.4154 0.1330 0.5183
## Rev.VIK_milliers.KM_21 0.397 -0.4910 0.2140 0.7947
## POP_VIK_21 0.435 -0.5260 0.0802 0.8725
## EVO_POP_VIK_13.21 -0.252 0.0411 -0.2833 -0.1828
## PEVO_POP_VIK_13.21 -0.128 -0.1144 -0.2077 -0.0252
## REV_POP_VIK_21 0.561 -0.3941 -0.0490 0.6377
## EVO_REVcon_POP_13.21 0.262 -0.2031 0.2808 0.4067
## PRIX_ViK_23 0.164 0.1511 -0.0161 0.0301
## PEVO_PRIX_cou_13.23 0.165 0.0750 0.0777 0.1648
## EVO_PRIX_cou_13.23 0.146 0.1040 0.0679 0.0963
## EAUprel_21 0.425 -0.5939 0.1645 0.8303
## EAUprelHab_21 0.152 -0.2906 0.2223 0.2731
## Perteseau_Pourc_VIK21 -0.322 0.0854 0.2751 -0.1379
## Energie.kwh.m3_21ViK -0.123 0.2995 0.0141 -0.0895
## Continuite_21_VIK -0.166 -0.3881 0.0742 -0.0022
## Avaries100Km_VIK_21 -0.188 -0.1401 0.0884 -0.0794
## Conso_eau_hab_23 0.357 -0.5047 -0.0712 0.4140
## EVO_REVcon_VIK_pourc_13.21 EVO_REVcon_VIK_31.21
## RevHab_24 0.3673 0.273
## Pop65_24 0.1582 0.094
## TxMig_14.24 0.1297 0.298
## DensPop_24 0.3634 0.608
## PEVO_POP_01.21 0.3430 0.578
## NonBG_21 -0.1675 -0.033
## AgriVA_23 -0.3454 -0.560
## IndusVA_23 -0.1790 -0.227
## PIBhab_23 -0.0657 0.310
## PopRural_24 -0.2485 -0.546
## TxRisq_24 -0.0495 0.086
## REV_VIK_21 0.1888 0.712
## EVO_REVcon_VIK_pourc_13.21 1.0000 0.739
## EVO_REVcon_VIK_31.21 0.7395 1.000
## Km_reseau_eau_21 -0.0454 0.537
## dens_pop_reseau21 0.5878 0.615
## Rev.VIK_milliers.KM_21 0.3558 0.652
## POP_VIK_21 0.3492 0.755
## EVO_POP_VIK_13.21 0.5309 0.205
## PEVO_POP_VIK_13.21 0.5824 0.363
## REV_POP_VIK_21 -0.1554 0.257
## EVO_REVcon_POP_13.21 0.2529 0.470
## PRIX_ViK_23 -0.2841 -0.173
## PEVO_PRIX_cou_13.23 0.1144 0.147
## EVO_PRIX_cou_13.23 -0.0613 -0.011
## EAUprel_21 0.2386 0.641
## EAUprelHab_21 0.0099 0.154
## Perteseau_Pourc_VIK21 -0.1560 -0.214
## Energie.kwh.m3_21ViK -0.7203 -0.468
## Continuite_21_VIK 0.0443 0.095
## Avaries100Km_VIK_21 -0.0635 0.000
## Conso_eau_hab_23 0.4830 0.487
## Km_reseau_eau_21 dens_pop_reseau21
## RevHab_24 0.036 0.319
## Pop65_24 0.052 -0.114
## TxMig_14.24 0.481 0.327
## DensPop_24 0.410 0.755
## PEVO_POP_01.21 0.444 0.748
## NonBG_21 0.130 0.063
## AgriVA_23 -0.492 -0.595
## IndusVA_23 -0.055 -0.084
## PIBhab_23 0.638 0.134
## PopRural_24 -0.478 -0.415
## TxRisq_24 -0.074 0.133
## REV_VIK_21 0.839 0.518
## EVO_REVcon_VIK_pourc_13.21 -0.045 0.588
## EVO_REVcon_VIK_31.21 0.537 0.615
## Km_reseau_eau_21 1.000 0.167
## dens_pop_reseau21 0.167 1.000
## Rev.VIK_milliers.KM_21 0.426 0.758
## POP_VIK_21 0.708 0.765
## EVO_POP_VIK_13.21 -0.262 0.247
## PEVO_POP_VIK_13.21 -0.060 0.372
## REV_POP_VIK_21 0.554 -0.043
## EVO_REVcon_POP_13.21 0.314 0.049
## PRIX_ViK_23 0.168 -0.321
## PEVO_PRIX_cou_13.23 0.105 0.038
## EVO_PRIX_cou_13.23 0.100 -0.106
## EAUprel_21 0.652 0.620
## EAUprelHab_21 0.186 0.048
## Perteseau_Pourc_VIK21 -0.212 -0.099
## Energie.kwh.m3_21ViK 0.148 -0.516
## Continuite_21_VIK -0.043 -0.065
## Avaries100Km_VIK_21 -0.137 0.157
## Conso_eau_hab_23 0.139 0.506
## Rev.VIK_milliers.KM_21 POP_VIK_21 EVO_POP_VIK_13.21
## RevHab_24 0.362 0.1708 0.129
## Pop65_24 0.032 -0.0722 -0.030
## TxMig_14.24 0.331 0.4899 0.132
## DensPop_24 0.790 0.7192 -0.034
## PEVO_POP_01.21 0.711 0.7889 0.109
## NonBG_21 0.234 0.1385 -0.077
## AgriVA_23 -0.526 -0.6661 -0.041
## IndusVA_23 -0.130 -0.0892 -0.064
## PIBhab_23 0.397 0.4351 -0.252
## PopRural_24 -0.491 -0.5260 0.041
## TxRisq_24 0.214 0.0802 -0.283
## REV_VIK_21 0.795 0.8725 -0.183
## EVO_REVcon_VIK_pourc_13.21 0.356 0.3492 0.531
## EVO_REVcon_VIK_31.21 0.652 0.7553 0.205
## Km_reseau_eau_21 0.426 0.7077 -0.262
## dens_pop_reseau21 0.758 0.7646 0.247
## Rev.VIK_milliers.KM_21 1.000 0.7559 -0.155
## POP_VIK_21 0.756 1.0000 0.074
## EVO_POP_VIK_13.21 -0.155 0.0744 1.000
## PEVO_POP_VIK_13.21 -0.042 0.2748 0.922
## REV_POP_VIK_21 0.545 0.2600 -0.532
## EVO_REVcon_POP_13.21 0.423 0.1795 -0.453
## PRIX_ViK_23 -0.063 -0.1784 -0.313
## PEVO_PRIX_cou_13.23 0.211 0.0033 -0.246
## EVO_PRIX_cou_13.23 0.117 -0.1045 -0.345
## EAUprel_21 0.786 0.8342 -0.138
## EAUprelHab_21 0.346 0.1264 -0.272
## Perteseau_Pourc_VIK21 -0.011 -0.1516 -0.123
## Energie.kwh.m3_21ViK -0.278 -0.2385 -0.417
## Continuite_21_VIK -0.087 -0.0892 -0.049
## Avaries100Km_VIK_21 -0.076 0.0733 -0.020
## Conso_eau_hab_23 0.495 0.4419 0.213
## PEVO_POP_VIK_13.21 REV_POP_VIK_21
## RevHab_24 0.144 0.181
## Pop65_24 0.012 0.224
## TxMig_14.24 0.292 0.154
## DensPop_24 0.177 0.296
## PEVO_POP_01.21 0.244 0.118
## NonBG_21 -0.042 0.102
## AgriVA_23 -0.195 -0.242
## IndusVA_23 0.019 -0.165
## PIBhab_23 -0.128 0.561
## PopRural_24 -0.114 -0.394
## TxRisq_24 -0.208 -0.049
## REV_VIK_21 -0.025 0.638
## EVO_REVcon_VIK_pourc_13.21 0.582 -0.155
## EVO_REVcon_VIK_31.21 0.363 0.257
## Km_reseau_eau_21 -0.060 0.554
## dens_pop_reseau21 0.372 -0.043
## Rev.VIK_milliers.KM_21 -0.042 0.545
## POP_VIK_21 0.275 0.260
## EVO_POP_VIK_13.21 0.922 -0.532
## PEVO_POP_VIK_13.21 1.000 -0.524
## REV_POP_VIK_21 -0.524 1.000
## EVO_REVcon_POP_13.21 -0.459 0.608
## PRIX_ViK_23 -0.366 0.379
## PEVO_PRIX_cou_13.23 -0.278 0.288
## EVO_PRIX_cou_13.23 -0.400 0.359
## EAUprel_21 0.022 0.432
## EAUprelHab_21 -0.242 0.386
## Perteseau_Pourc_VIK21 -0.166 -0.022
## Energie.kwh.m3_21ViK -0.427 0.147
## Continuite_21_VIK -0.043 -0.102
## Avaries100Km_VIK_21 -0.068 -0.332
## Conso_eau_hab_23 0.246 0.265
## EVO_REVcon_POP_13.21 PRIX_ViK_23 PEVO_PRIX_cou_13.23
## RevHab_24 0.283 -0.121 -0.0454
## Pop65_24 0.118 0.273 0.1604
## TxMig_14.24 -0.149 -0.103 -0.0991
## DensPop_24 0.147 -0.081 0.1024
## PEVO_POP_01.21 0.081 -0.191 0.0183
## NonBG_21 -0.120 0.278 0.1242
## AgriVA_23 -0.176 0.140 -0.0192
## IndusVA_23 -0.400 0.006 -0.0383
## PIBhab_23 0.262 0.164 0.1653
## PopRural_24 -0.203 0.151 0.0750
## TxRisq_24 0.281 -0.016 0.0777
## REV_VIK_21 0.407 0.030 0.1648
## EVO_REVcon_VIK_pourc_13.21 0.253 -0.284 0.1144
## EVO_REVcon_VIK_31.21 0.470 -0.173 0.1472
## Km_reseau_eau_21 0.314 0.168 0.1051
## dens_pop_reseau21 0.049 -0.321 0.0383
## Rev.VIK_milliers.KM_21 0.423 -0.063 0.2113
## POP_VIK_21 0.180 -0.178 0.0033
## EVO_POP_VIK_13.21 -0.453 -0.313 -0.2458
## PEVO_POP_VIK_13.21 -0.459 -0.366 -0.2775
## REV_POP_VIK_21 0.608 0.379 0.2879
## EVO_REVcon_POP_13.21 1.000 0.233 0.4078
## PRIX_ViK_23 0.233 1.000 0.7099
## PEVO_PRIX_cou_13.23 0.408 0.710 1.0000
## EVO_PRIX_cou_13.23 0.388 0.880 0.9425
## EAUprel_21 0.315 -0.296 -0.1319
## EAUprelHab_21 0.247 -0.204 -0.1888
## Perteseau_Pourc_VIK21 -0.106 -0.134 -0.2611
## Energie.kwh.m3_21ViK -0.161 0.580 0.2015
## Continuite_21_VIK 0.027 -0.226 -0.1270
## Avaries100Km_VIK_21 -0.122 -0.339 -0.1067
## Conso_eau_hab_23 0.276 -0.377 -0.1235
## EVO_PRIX_cou_13.23 EAUprel_21 EAUprelHab_21
## RevHab_24 -0.097 0.286 0.2189
## Pop65_24 0.244 0.068 0.3678
## TxMig_14.24 -0.136 0.384 0.0044
## DensPop_24 0.033 0.604 0.1368
## PEVO_POP_01.21 -0.082 0.695 0.1385
## NonBG_21 0.185 0.110 0.0706
## AgriVA_23 0.053 -0.529 -0.0203
## IndusVA_23 -0.057 -0.159 -0.0575
## PIBhab_23 0.146 0.425 0.1516
## PopRural_24 0.104 -0.594 -0.2906
## TxRisq_24 0.068 0.164 0.2223
## REV_VIK_21 0.096 0.830 0.2731
## EVO_REVcon_VIK_pourc_13.21 -0.061 0.239 0.0099
## EVO_REVcon_VIK_31.21 -0.011 0.641 0.1538
## Km_reseau_eau_21 0.100 0.652 0.1861
## dens_pop_reseau21 -0.106 0.620 0.0476
## Rev.VIK_milliers.KM_21 0.117 0.786 0.3465
## POP_VIK_21 -0.105 0.834 0.1264
## EVO_POP_VIK_13.21 -0.345 -0.138 -0.2720
## PEVO_POP_VIK_13.21 -0.400 0.022 -0.2425
## REV_POP_VIK_21 0.359 0.432 0.3864
## EVO_REVcon_POP_13.21 0.388 0.315 0.2469
## PRIX_ViK_23 0.880 -0.296 -0.2036
## PEVO_PRIX_cou_13.23 0.943 -0.132 -0.1888
## EVO_PRIX_cou_13.23 1.000 -0.223 -0.1866
## EAUprel_21 -0.223 1.000 0.6065
## EAUprelHab_21 -0.187 0.606 1.0000
## Perteseau_Pourc_VIK21 -0.207 0.230 0.7482
## Energie.kwh.m3_21ViK 0.377 -0.289 -0.1725
## Continuite_21_VIK -0.184 0.099 0.1631
## Avaries100Km_VIK_21 -0.212 0.049 -0.1730
## Conso_eau_hab_23 -0.210 0.389 0.0564
## Perteseau_Pourc_VIK21 Energie.kwh.m3_21ViK
## RevHab_24 0.039 -0.475
## Pop65_24 0.350 0.091
## TxMig_14.24 -0.216 -0.121
## DensPop_24 -0.136 -0.343
## PEVO_POP_01.21 -0.136 -0.316
## NonBG_21 0.177 0.466
## AgriVA_23 0.383 0.507
## IndusVA_23 -0.076 0.012
## PIBhab_23 -0.322 -0.123
## PopRural_24 0.085 0.300
## TxRisq_24 0.275 0.014
## REV_VIK_21 -0.138 -0.090
## EVO_REVcon_VIK_pourc_13.21 -0.156 -0.720
## EVO_REVcon_VIK_31.21 -0.214 -0.468
## Km_reseau_eau_21 -0.212 0.148
## dens_pop_reseau21 -0.099 -0.516
## Rev.VIK_milliers.KM_21 -0.011 -0.278
## POP_VIK_21 -0.152 -0.238
## EVO_POP_VIK_13.21 -0.123 -0.417
## PEVO_POP_VIK_13.21 -0.166 -0.427
## REV_POP_VIK_21 -0.022 0.147
## EVO_REVcon_POP_13.21 -0.106 -0.161
## PRIX_ViK_23 -0.134 0.580
## PEVO_PRIX_cou_13.23 -0.261 0.202
## EVO_PRIX_cou_13.23 -0.207 0.377
## EAUprel_21 0.230 -0.289
## EAUprelHab_21 0.748 -0.172
## Perteseau_Pourc_VIK21 1.000 0.035
## Energie.kwh.m3_21ViK 0.035 1.000
## Continuite_21_VIK 0.178 -0.176
## Avaries100Km_VIK_21 -0.080 -0.151
## Conso_eau_hab_23 -0.216 -0.609
## Continuite_21_VIK Avaries100Km_VIK_21
## RevHab_24 0.0777 -0.297
## Pop65_24 0.0137 -0.369
## TxMig_14.24 -0.1762 0.063
## DensPop_24 -0.1489 -0.105
## PEVO_POP_01.21 -0.1708 0.236
## NonBG_21 -0.3979 -0.236
## AgriVA_23 0.0027 -0.078
## IndusVA_23 -0.1959 -0.076
## PIBhab_23 -0.1664 -0.188
## PopRural_24 -0.3881 -0.140
## TxRisq_24 0.0742 0.088
## REV_VIK_21 -0.0022 -0.079
## EVO_REVcon_VIK_pourc_13.21 0.0443 -0.063
## EVO_REVcon_VIK_31.21 0.0947 0.000
## Km_reseau_eau_21 -0.0427 -0.137
## dens_pop_reseau21 -0.0651 0.157
## Rev.VIK_milliers.KM_21 -0.0870 -0.076
## POP_VIK_21 -0.0892 0.073
## EVO_POP_VIK_13.21 -0.0493 -0.020
## PEVO_POP_VIK_13.21 -0.0432 -0.068
## REV_POP_VIK_21 -0.1024 -0.332
## EVO_REVcon_POP_13.21 0.0274 -0.122
## PRIX_ViK_23 -0.2261 -0.339
## PEVO_PRIX_cou_13.23 -0.1270 -0.107
## EVO_PRIX_cou_13.23 -0.1839 -0.212
## EAUprel_21 0.0985 0.049
## EAUprelHab_21 0.1631 -0.173
## Perteseau_Pourc_VIK21 0.1779 -0.080
## Energie.kwh.m3_21ViK -0.1758 -0.151
## Continuite_21_VIK 1.0000 0.500
## Avaries100Km_VIK_21 0.5003 1.000
## Conso_eau_hab_23 -0.1043 -0.068
## Conso_eau_hab_23
## RevHab_24 0.327
## Pop65_24 0.131
## TxMig_14.24 0.192
## DensPop_24 0.418
## PEVO_POP_01.21 0.268
## NonBG_21 -0.462
## AgriVA_23 -0.622
## IndusVA_23 -0.095
## PIBhab_23 0.357
## PopRural_24 -0.505
## TxRisq_24 -0.071
## REV_VIK_21 0.414
## EVO_REVcon_VIK_pourc_13.21 0.483
## EVO_REVcon_VIK_31.21 0.487
## Km_reseau_eau_21 0.139
## dens_pop_reseau21 0.506
## Rev.VIK_milliers.KM_21 0.495
## POP_VIK_21 0.442
## EVO_POP_VIK_13.21 0.213
## PEVO_POP_VIK_13.21 0.246
## REV_POP_VIK_21 0.265
## EVO_REVcon_POP_13.21 0.276
## PRIX_ViK_23 -0.377
## PEVO_PRIX_cou_13.23 -0.123
## EVO_PRIX_cou_13.23 -0.210
## EAUprel_21 0.389
## EAUprelHab_21 0.056
## Perteseau_Pourc_VIK21 -0.216
## Energie.kwh.m3_21ViK -0.609
## Continuite_21_VIK -0.104
## Avaries100Km_VIK_21 -0.068
## Conso_eau_hab_23 1.000Sans Sofia
ViKnoSOFMatrice<-ViKnoSOF[,c(4:14,16:36)]
Matrice_global_nosof<-cor(ViKnoSOFMatrice,method=c("spearman"))
Matrice_global_nosof## RevHab_24 Pop65_24 TxMig_14.24 DensPop_24
## RevHab_24 1.0000 0.384 -0.043 0.307
## Pop65_24 0.3840 1.000 0.204 0.178
## TxMig_14.24 -0.0427 0.204 1.000 0.484
## DensPop_24 0.3071 0.178 0.484 1.000
## PEVO_POP_01.21 0.0354 -0.046 0.658 0.707
## NonBG_21 -0.2198 -0.070 0.345 0.256
## AgriVA_23 -0.2631 -0.029 -0.352 -0.628
## IndusVA_23 0.0745 -0.225 0.191 0.222
## PIBhab_23 -0.0153 -0.016 0.256 0.376
## PopRural_24 -0.1502 -0.479 -0.297 -0.474
## TxRisq_24 -0.0067 -0.070 -0.120 0.303
## REV_VIK_21 0.1557 0.186 0.390 0.602
## EVO_REVcon_VIK_pourc_13.21 0.3553 0.219 0.082 0.335
## EVO_REVcon_VIK_31.21 0.1886 0.220 0.237 0.563
## Km_reseau_eau_21 -0.0556 0.139 0.437 0.347
## dens_pop_reseau21 0.2399 -0.012 0.275 0.726
## Rev.VIK_milliers.KM_21 0.2888 0.151 0.283 0.766
## POP_VIK_21 0.0751 0.035 0.450 0.687
## EVO_POP_VIK_13.21 0.0861 0.026 0.115 -0.096
## PEVO_POP_VIK_13.21 0.1319 0.065 0.280 0.158
## REV_POP_VIK_21 0.1001 0.338 0.095 0.218
## EVO_REVcon_POP_13.21 0.2265 0.217 -0.239 0.069
## PRIX_ViK_23 -0.0641 0.230 -0.073 -0.021
## PEVO_PRIX_cou_13.23 -0.0348 0.159 -0.107 0.117
## EVO_PRIX_cou_13.23 -0.0543 0.220 -0.123 0.083
## EAUprel_21 0.2033 0.192 0.339 0.559
## EAUprelHab_21 0.2479 0.366 0.021 0.157
## Perteseau_Pourc_VIK21 0.1587 0.275 -0.153 -0.036
## Energie.kwh.m3_21ViK -0.4141 -0.014 -0.050 -0.267
## Continuite_21_VIK 0.1007 0.000 -0.159 -0.143
## Avaries100Km_VIK_21 -0.3736 -0.339 0.023 -0.165
## Conso_eau_hab_23 0.2495 0.261 0.135 0.351
## PEVO_POP_01.21 NonBG_21 AgriVA_23 IndusVA_23
## RevHab_24 0.035 -0.220 -0.263 0.07448
## Pop65_24 -0.046 -0.070 -0.029 -0.22527
## TxMig_14.24 0.658 0.345 -0.352 0.19109
## DensPop_24 0.707 0.256 -0.628 0.22161
## PEVO_POP_01.21 1.000 0.410 -0.519 0.10168
## NonBG_21 0.410 1.000 0.391 -0.13858
## AgriVA_23 -0.519 0.391 1.000 -0.43529
## IndusVA_23 0.102 -0.139 -0.435 1.00000
## PIBhab_23 0.234 -0.156 -0.679 0.58608
## PopRural_24 -0.275 0.311 0.546 0.15385
## TxRisq_24 0.165 0.133 -0.034 -0.08669
## REV_VIK_21 0.587 0.142 -0.571 -0.05922
## EVO_REVcon_VIK_pourc_13.21 0.317 -0.168 -0.320 -0.16117
## EVO_REVcon_VIK_31.21 0.529 -0.016 -0.510 -0.15263
## Km_reseau_eau_21 0.384 0.153 -0.438 0.02259
## dens_pop_reseau21 0.719 0.083 -0.548 0.00488
## Rev.VIK_milliers.KM_21 0.678 0.264 -0.471 -0.04640
## POP_VIK_21 0.765 0.165 -0.628 -0.00061
## EVO_POP_VIK_13.21 0.068 -0.060 0.011 -0.02747
## PEVO_POP_VIK_13.21 0.226 -0.034 -0.176 0.03663
## REV_POP_VIK_21 0.018 0.125 -0.158 -0.08791
## EVO_REVcon_POP_13.21 -0.007 -0.108 -0.099 -0.33822
## PRIX_ViK_23 -0.136 0.268 0.082 -0.05372
## PEVO_PRIX_cou_13.23 0.033 0.125 -0.023 -0.06227
## EVO_PRIX_cou_13.23 -0.031 0.183 0.013 -0.10928
## EAUprel_21 0.660 0.132 -0.474 -0.07875
## EAUprelHab_21 0.169 0.056 -0.041 -0.09280
## Perteseau_Pourc_VIK21 -0.036 0.152 0.311 -0.19841
## Energie.kwh.m3_21ViK -0.237 0.474 0.450 -0.09162
## Continuite_21_VIK -0.164 -0.408 -0.011 -0.20696
## Avaries100Km_VIK_21 0.204 -0.234 -0.035 -0.01099
## Conso_eau_hab_23 0.183 -0.470 -0.578 -0.00183
## PIBhab_23 PopRural_24 TxRisq_24 REV_VIK_21
## RevHab_24 -0.0153 -0.150 -0.0067 0.1557
## Pop65_24 -0.0159 -0.479 -0.0702 0.1862
## TxMig_14.24 0.2558 -0.297 -0.1197 0.3901
## DensPop_24 0.3761 -0.474 0.3028 0.6020
## PEVO_POP_01.21 0.2336 -0.275 0.1652 0.5869
## NonBG_21 -0.1557 0.311 0.1331 0.1416
## AgriVA_23 -0.6795 0.546 -0.0336 -0.5708
## IndusVA_23 0.5861 0.154 -0.0867 -0.0592
## PIBhab_23 1.0000 -0.310 -0.0024 0.5153
## PopRural_24 -0.3101 1.000 -0.0537 -0.5525
## TxRisq_24 -0.0024 -0.054 1.0000 0.0122
## REV_VIK_21 0.5153 -0.553 0.0122 1.0000
## EVO_REVcon_VIK_pourc_13.21 -0.1245 -0.219 -0.0446 0.1593
## EVO_REVcon_VIK_31.21 0.2308 -0.494 0.1245 0.6789
## Km_reseau_eau_21 0.6062 -0.425 -0.0549 0.8236
## dens_pop_reseau21 0.0342 -0.348 0.1758 0.4628
## Rev.VIK_milliers.KM_21 0.3272 -0.432 0.2589 0.7711
## POP_VIK_21 0.3700 -0.471 0.1178 0.8578
## EVO_POP_VIK_13.21 -0.3291 0.105 -0.2802 -0.2485
## PEVO_POP_VIK_13.21 -0.1740 -0.098 -0.2002 -0.0519
## REV_POP_VIK_21 0.5165 -0.328 -0.0238 0.5977
## EVO_REVcon_POP_13.21 0.1972 -0.125 0.3138 0.3565
## PRIX_ViK_23 0.2308 0.107 -0.0372 0.0885
## PEVO_PRIX_cou_13.23 0.1746 0.069 0.0672 0.1783
## EVO_PRIX_cou_13.23 0.1923 0.071 0.0458 0.1508
## EAUprel_21 0.3584 -0.547 0.2082 0.8107
## EAUprelHab_21 0.1740 -0.328 0.2027 0.3114
## Perteseau_Pourc_VIK21 -0.2442 -0.020 0.2491 -0.0385
## Energie.kwh.m3_21ViK -0.0220 0.219 -0.0214 0.0156
## Continuite_21_VIK -0.1600 -0.418 0.0745 0.0067
## Avaries100Km_VIK_21 -0.2552 -0.091 0.1160 -0.1337
## Conso_eau_hab_23 0.2822 -0.447 -0.0446 0.3464
## EVO_REVcon_VIK_pourc_13.21 EVO_REVcon_VIK_31.21
## RevHab_24 0.355 0.189
## Pop65_24 0.219 0.220
## TxMig_14.24 0.082 0.237
## DensPop_24 0.335 0.563
## PEVO_POP_01.21 0.317 0.529
## NonBG_21 -0.168 -0.016
## AgriVA_23 -0.320 -0.510
## IndusVA_23 -0.161 -0.153
## PIBhab_23 -0.125 0.231
## PopRural_24 -0.219 -0.494
## TxRisq_24 -0.045 0.125
## REV_VIK_21 0.159 0.679
## EVO_REVcon_VIK_pourc_13.21 1.000 0.747
## EVO_REVcon_VIK_31.21 0.747 1.000
## Km_reseau_eau_21 -0.093 0.488
## dens_pop_reseau21 0.578 0.570
## Rev.VIK_milliers.KM_21 0.330 0.612
## POP_VIK_21 0.326 0.727
## EVO_POP_VIK_13.21 0.520 0.164
## PEVO_POP_VIK_13.21 0.578 0.358
## REV_POP_VIK_21 -0.218 0.173
## EVO_REVcon_POP_13.21 0.236 0.425
## PRIX_ViK_23 -0.261 -0.128
## PEVO_PRIX_cou_13.23 0.130 0.167
## EVO_PRIX_cou_13.23 -0.038 0.031
## EAUprel_21 0.211 0.600
## EAUprelHab_21 0.018 0.177
## Perteseau_Pourc_VIK21 -0.117 -0.123
## Energie.kwh.m3_21ViK -0.717 -0.406
## Continuite_21_VIK 0.056 0.105
## Avaries100Km_VIK_21 -0.078 -0.049
## Conso_eau_hab_23 0.459 0.428
## Km_reseau_eau_21 dens_pop_reseau21
## RevHab_24 -0.056 0.2399
## Pop65_24 0.139 -0.0122
## TxMig_14.24 0.437 0.2747
## DensPop_24 0.347 0.7265
## PEVO_POP_01.21 0.384 0.7188
## NonBG_21 0.153 0.0830
## AgriVA_23 -0.438 -0.5482
## IndusVA_23 0.023 0.0049
## PIBhab_23 0.606 0.0342
## PopRural_24 -0.425 -0.3480
## TxRisq_24 -0.055 0.1758
## REV_VIK_21 0.824 0.4628
## EVO_REVcon_VIK_pourc_13.21 -0.093 0.5775
## EVO_REVcon_VIK_31.21 0.488 0.5702
## Km_reseau_eau_21 1.000 0.0788
## dens_pop_reseau21 0.079 1.0000
## Rev.VIK_milliers.KM_21 0.364 0.7302
## POP_VIK_21 0.678 0.7375
## EVO_POP_VIK_13.21 -0.312 0.2094
## PEVO_POP_VIK_13.21 -0.084 0.3651
## REV_POP_VIK_21 0.513 -0.1569
## EVO_REVcon_POP_13.21 0.253 -0.0372
## PRIX_ViK_23 0.212 -0.2723
## PEVO_PRIX_cou_13.23 0.104 0.0495
## EVO_PRIX_cou_13.23 0.125 -0.0562
## EAUprel_21 0.615 0.5763
## EAUprelHab_21 0.202 0.0702
## Perteseau_Pourc_VIK21 -0.142 0.0055
## Energie.kwh.m3_21ViK 0.263 -0.4605
## Continuite_21_VIK -0.031 -0.0617
## Avaries100Km_VIK_21 -0.181 0.1227
## Conso_eau_hab_23 0.052 0.4493
## Rev.VIK_milliers.KM_21 POP_VIK_21 EVO_POP_VIK_13.21
## RevHab_24 0.2888 0.07509 0.086
## Pop65_24 0.1508 0.03480 0.026
## TxMig_14.24 0.2827 0.44994 0.115
## DensPop_24 0.7662 0.68681 -0.096
## PEVO_POP_01.21 0.6776 0.76458 0.068
## NonBG_21 0.2637 0.16545 -0.060
## AgriVA_23 -0.4713 -0.62759 0.011
## IndusVA_23 -0.0464 -0.00061 -0.027
## PIBhab_23 0.3272 0.36996 -0.329
## PopRural_24 -0.4322 -0.47131 0.105
## TxRisq_24 0.2589 0.11783 -0.280
## REV_VIK_21 0.7711 0.85775 -0.248
## EVO_REVcon_VIK_pourc_13.21 0.3303 0.32601 0.520
## EVO_REVcon_VIK_31.21 0.6117 0.72711 0.164
## Km_reseau_eau_21 0.3639 0.67766 -0.312
## dens_pop_reseau21 0.7302 0.73748 0.209
## Rev.VIK_milliers.KM_21 1.0000 0.72772 -0.222
## POP_VIK_21 0.7277 1.00000 0.024
## EVO_POP_VIK_13.21 -0.2222 0.02381 1.000
## PEVO_POP_VIK_13.21 -0.0702 0.26252 0.915
## REV_POP_VIK_21 0.4939 0.17766 -0.610
## EVO_REVcon_POP_13.21 0.3748 0.10440 -0.511
## PRIX_ViK_23 -0.0073 -0.12759 -0.279
## PEVO_PRIX_cou_13.23 0.2277 0.01099 -0.222
## EVO_PRIX_cou_13.23 0.1807 -0.06166 -0.325
## EAUprel_21 0.7613 0.81502 -0.203
## EAUprelHab_21 0.3840 0.15446 -0.270
## Perteseau_Pourc_VIK21 0.1032 -0.05372 -0.082
## Energie.kwh.m3_21ViK -0.1942 -0.15056 -0.393
## Continuite_21_VIK -0.0830 -0.08974 -0.053
## Avaries100Km_VIK_21 -0.1306 0.03541 -0.046
## Conso_eau_hab_23 0.4362 0.37752 0.168
## PEVO_POP_VIK_13.21 REV_POP_VIK_21
## RevHab_24 0.1319 0.100
## Pop65_24 0.0647 0.338
## TxMig_14.24 0.2802 0.095
## DensPop_24 0.1575 0.218
## PEVO_POP_01.21 0.2256 0.018
## NonBG_21 -0.0336 0.125
## AgriVA_23 -0.1764 -0.158
## IndusVA_23 0.0366 -0.088
## PIBhab_23 -0.1740 0.516
## PopRural_24 -0.0977 -0.328
## TxRisq_24 -0.2002 -0.024
## REV_VIK_21 -0.0519 0.598
## EVO_REVcon_VIK_pourc_13.21 0.5775 -0.218
## EVO_REVcon_VIK_31.21 0.3578 0.173
## Km_reseau_eau_21 -0.0842 0.513
## dens_pop_reseau21 0.3651 -0.157
## Rev.VIK_milliers.KM_21 -0.0702 0.494
## POP_VIK_21 0.2625 0.178
## EVO_POP_VIK_13.21 0.9145 -0.610
## PEVO_POP_VIK_13.21 1.0000 -0.577
## REV_POP_VIK_21 -0.5769 1.000
## EVO_REVcon_POP_13.21 -0.4988 0.574
## PRIX_ViK_23 -0.3425 0.439
## PEVO_PRIX_cou_13.23 -0.2601 0.295
## EVO_PRIX_cou_13.23 -0.3907 0.413
## EAUprel_21 -0.0018 0.370
## EAUprelHab_21 -0.2338 0.419
## Perteseau_Pourc_VIK21 -0.1368 0.076
## Energie.kwh.m3_21ViK -0.4178 0.271
## Continuite_21_VIK -0.0427 -0.088
## Avaries100Km_VIK_21 -0.0934 -0.404
## Conso_eau_hab_23 0.2285 0.186
## EVO_REVcon_POP_13.21 PRIX_ViK_23 PEVO_PRIX_cou_13.23
## RevHab_24 0.226 -0.0641 -0.035
## Pop65_24 0.217 0.2302 0.159
## TxMig_14.24 -0.239 -0.0726 -0.107
## DensPop_24 0.069 -0.0214 0.117
## PEVO_POP_01.21 -0.007 -0.1356 0.033
## NonBG_21 -0.108 0.2680 0.125
## AgriVA_23 -0.099 0.0818 -0.023
## IndusVA_23 -0.338 -0.0537 -0.062
## PIBhab_23 0.197 0.2308 0.175
## PopRural_24 -0.125 0.1068 0.069
## TxRisq_24 0.314 -0.0372 0.067
## REV_VIK_21 0.357 0.0885 0.178
## EVO_REVcon_VIK_pourc_13.21 0.236 -0.2607 0.130
## EVO_REVcon_VIK_31.21 0.425 -0.1282 0.167
## Km_reseau_eau_21 0.253 0.2118 0.104
## dens_pop_reseau21 -0.037 -0.2723 0.049
## Rev.VIK_milliers.KM_21 0.375 -0.0073 0.228
## POP_VIK_21 0.104 -0.1276 0.011
## EVO_POP_VIK_13.21 -0.511 -0.2790 -0.222
## PEVO_POP_VIK_13.21 -0.499 -0.3425 -0.260
## REV_POP_VIK_21 0.574 0.4389 0.295
## EVO_REVcon_POP_13.21 1.000 0.2845 0.430
## PRIX_ViK_23 0.284 1.0000 0.707
## PEVO_PRIX_cou_13.23 0.430 0.7070 1.000
## EVO_PRIX_cou_13.23 0.441 0.8712 0.944
## EAUprel_21 0.255 -0.2570 -0.133
## EAUprelHab_21 0.274 -0.2295 -0.206
## Perteseau_Pourc_VIK21 -0.018 -0.1996 -0.286
## Energie.kwh.m3_21ViK -0.083 0.5518 0.198
## Continuite_21_VIK 0.040 -0.2350 -0.120
## Avaries100Km_VIK_21 -0.175 -0.3242 -0.106
## Conso_eau_hab_23 0.209 -0.3360 -0.124
## EVO_PRIX_cou_13.23 EAUprel_21 EAUprelHab_21
## RevHab_24 -0.054 0.2033 0.248
## Pop65_24 0.220 0.1917 0.366
## TxMig_14.24 -0.123 0.3388 0.021
## DensPop_24 0.083 0.5586 0.157
## PEVO_POP_01.21 -0.031 0.6602 0.169
## NonBG_21 0.183 0.1319 0.056
## AgriVA_23 0.013 -0.4744 -0.041
## IndusVA_23 -0.109 -0.0788 -0.093
## PIBhab_23 0.192 0.3584 0.174
## PopRural_24 0.071 -0.5470 -0.328
## TxRisq_24 0.046 0.2082 0.203
## REV_VIK_21 0.151 0.8107 0.311
## EVO_REVcon_VIK_pourc_13.21 -0.038 0.2112 0.018
## EVO_REVcon_VIK_31.21 0.031 0.6001 0.177
## Km_reseau_eau_21 0.125 0.6154 0.202
## dens_pop_reseau21 -0.056 0.5763 0.070
## Rev.VIK_milliers.KM_21 0.181 0.7613 0.384
## POP_VIK_21 -0.062 0.8150 0.154
## EVO_POP_VIK_13.21 -0.325 -0.2027 -0.270
## PEVO_POP_VIK_13.21 -0.391 -0.0018 -0.234
## REV_POP_VIK_21 0.413 0.3700 0.419
## EVO_REVcon_POP_13.21 0.441 0.2552 0.274
## PRIX_ViK_23 0.871 -0.2570 -0.230
## PEVO_PRIX_cou_13.23 0.944 -0.1331 -0.206
## EVO_PRIX_cou_13.23 1.000 -0.1819 -0.194
## EAUprel_21 -0.182 1.0000 0.665
## EAUprelHab_21 -0.194 0.6654 1.000
## Perteseau_Pourc_VIK21 -0.252 0.3718 0.772
## Energie.kwh.m3_21ViK 0.353 -0.2067 -0.202
## Continuite_21_VIK -0.192 0.1117 0.170
## Avaries100Km_VIK_21 -0.201 0.0116 -0.146
## Conso_eau_hab_23 -0.173 0.3189 0.081
## Perteseau_Pourc_VIK21 Energie.kwh.m3_21ViK
## RevHab_24 0.1587 -0.414
## Pop65_24 0.2753 -0.014
## TxMig_14.24 -0.1526 -0.050
## DensPop_24 -0.0360 -0.267
## PEVO_POP_01.21 -0.0363 -0.237
## NonBG_21 0.1520 0.474
## AgriVA_23 0.3114 0.450
## IndusVA_23 -0.1984 -0.092
## PIBhab_23 -0.2442 -0.022
## PopRural_24 -0.0201 0.219
## TxRisq_24 0.2491 -0.021
## REV_VIK_21 -0.0385 0.016
## EVO_REVcon_VIK_pourc_13.21 -0.1172 -0.717
## EVO_REVcon_VIK_31.21 -0.1233 -0.406
## Km_reseau_eau_21 -0.1422 0.263
## dens_pop_reseau21 0.0055 -0.461
## Rev.VIK_milliers.KM_21 0.1032 -0.194
## POP_VIK_21 -0.0537 -0.151
## EVO_POP_VIK_13.21 -0.0824 -0.393
## PEVO_POP_VIK_13.21 -0.1368 -0.418
## REV_POP_VIK_21 0.0763 0.271
## EVO_REVcon_POP_13.21 -0.0183 -0.083
## PRIX_ViK_23 -0.1996 0.552
## PEVO_PRIX_cou_13.23 -0.2863 0.198
## EVO_PRIX_cou_13.23 -0.2521 0.353
## EAUprel_21 0.3718 -0.207
## EAUprelHab_21 0.7723 -0.202
## Perteseau_Pourc_VIK21 1.0000 -0.077
## Energie.kwh.m3_21ViK -0.0767 1.000
## Continuite_21_VIK 0.1740 -0.201
## Avaries100Km_VIK_21 -0.0177 -0.109
## Conso_eau_hab_23 -0.1255 -0.564
## Continuite_21_VIK Avaries100Km_VIK_21
## RevHab_24 0.1007 -0.374
## Pop65_24 0.0000 -0.339
## TxMig_14.24 -0.1593 0.023
## DensPop_24 -0.1429 -0.165
## PEVO_POP_01.21 -0.1640 0.204
## NonBG_21 -0.4078 -0.234
## AgriVA_23 -0.0110 -0.035
## IndusVA_23 -0.2070 -0.011
## PIBhab_23 -0.1600 -0.255
## PopRural_24 -0.4182 -0.091
## TxRisq_24 0.0745 0.116
## REV_VIK_21 0.0067 -0.134
## EVO_REVcon_VIK_pourc_13.21 0.0562 -0.078
## EVO_REVcon_VIK_31.21 0.1050 -0.049
## Km_reseau_eau_21 -0.0311 -0.181
## dens_pop_reseau21 -0.0617 0.123
## Rev.VIK_milliers.KM_21 -0.0830 -0.131
## POP_VIK_21 -0.0897 0.035
## EVO_POP_VIK_13.21 -0.0531 -0.046
## PEVO_POP_VIK_13.21 -0.0427 -0.093
## REV_POP_VIK_21 -0.0879 -0.404
## EVO_REVcon_POP_13.21 0.0403 -0.175
## PRIX_ViK_23 -0.2350 -0.324
## PEVO_PRIX_cou_13.23 -0.1197 -0.106
## EVO_PRIX_cou_13.23 -0.1917 -0.201
## EAUprel_21 0.1117 0.012
## EAUprelHab_21 0.1697 -0.146
## Perteseau_Pourc_VIK21 0.1740 -0.018
## Energie.kwh.m3_21ViK -0.2013 -0.109
## Continuite_21_VIK 1.0000 0.515
## Avaries100Km_VIK_21 0.5147 1.000
## Conso_eau_hab_23 -0.1093 -0.128
## Conso_eau_hab_23
## RevHab_24 0.2495
## Pop65_24 0.2611
## TxMig_14.24 0.1353
## DensPop_24 0.3513
## PEVO_POP_01.21 0.1835
## NonBG_21 -0.4701
## AgriVA_23 -0.5782
## IndusVA_23 -0.0018
## PIBhab_23 0.2822
## PopRural_24 -0.4475
## TxRisq_24 -0.0446
## REV_VIK_21 0.3464
## EVO_REVcon_VIK_pourc_13.21 0.4588
## EVO_REVcon_VIK_31.21 0.4276
## Km_reseau_eau_21 0.0516
## dens_pop_reseau21 0.4493
## Rev.VIK_milliers.KM_21 0.4362
## POP_VIK_21 0.3775
## EVO_POP_VIK_13.21 0.1677
## PEVO_POP_VIK_13.21 0.2285
## REV_POP_VIK_21 0.1863
## EVO_REVcon_POP_13.21 0.2092
## PRIX_ViK_23 -0.3360
## PEVO_PRIX_cou_13.23 -0.1237
## EVO_PRIX_cou_13.23 -0.1732
## EAUprel_21 0.3189
## EAUprelHab_21 0.0809
## Perteseau_Pourc_VIK21 -0.1255
## Energie.kwh.m3_21ViK -0.5641
## Continuite_21_VIK -0.1093
## Avaries100Km_VIK_21 -0.1283
## Conso_eau_hab_23 1.00001. Analyses uni- et bi-variées
1.A Caractéristiques des ViKs
1.A.a. Démographie des opérateurs de réseau
Population totale des ViKs
Il y a une forte variation entre les ViKs en terme de nombre
d’habitant·es desservie.
Ce chiffre ne recoupe pas totalement la population régionale totale.
D’une part, parce que certains ViKs ne sont pas entièrement régionaux
(il y a des ViKs municipaux au sein de la région). D’autre part, parce
qu’il peut y avoir un léger écart entre les sources.
En 2021, la population minimum que desservait un ViK était de 85 000
habitant·es à Vidin et la population maximum de 1,3 millions
d’habitant·es à Sofia, avec une médiane de 173 000 hab., une moyenne de
243 000 hab. et un coefficient de variation de 102,7%. Le diagramme en
barres montre que Sofia est de très loin le plus gros opérateur en terme
d’abonné·es, suivi des régions où se situent les principales villes du
pays (Plovdiv, Varna, Burgas, Stara Zagora)
summary(ViK$POP_VIK_21)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 84865 111486 173012 243148 238316 1328790sd(ViK$POP_VIK_21)#écart-type## [1] 249801sd(ViK$POP_VIK_21)/mean(ViK$POP_VIK_21)*100#coefficient de variation## [1] 103ggplot(data=ViK, aes(x = reorder(REG_FR, POP_VIK_21), y= POP_VIK_21)) +
geom_bar(stat = "identity", width=0.5, color="black", fill="white")+
geom_text(aes(label=POP_VIK_21), vjust=-0.3, size=2)+
theme(axis.text.x = element_text(face="bold", color="black",size=7, angle=45),
axis.text.y = element_text(face="bold", color="black",size=7, angle=90))+
labs(x="ViK", y="Population des ViKs", title="Population des ViKs")+theme(plot.title=element_text(face="bold", color="black", size=10))+theme(axis.title.x = element_text(face="plain", color="black", size=9))+theme(axis.title.y = element_text(face="plain", color="black", size=9)) #diagramme en barres
Évolution de la population des ViKs (2013-21)
Entre 2013 et 2021, l’évolution globale de la population desservie
par les ViKs est globalement négative, reflétant la tendance à la baisse
démographique très importante dans tout le pays.
Les données sont cependant très dispersées, avec une amplitude
importante, de 136 000 et un maximum et un minimum à presque la même
distance du zéro : -68 201 habitant·es pour le ViK de Vidin et + 67 725
habitant·es pour le ViK de Haskovo.
La moyenne est de -962 hab. entre 2013 et 2021, la médiane de -7238 et
l’écart type de 30 576. Le coefficient de variation est de -3177%.
Parmi les ViKs dont la population augmente la situation est toutefois très différente.
À Sofia, l’évolution de la population positive du ViK correspond à une augmentation réelle de la population dans la ville et région capitale (voir ci-dessous, évolution de la population des oblasts entre 2001 et 2021). Cette augmentation, de +6% s’ajoute à une population importante au départ (>1M abonné·es)
Dans les autres ViKs où l’évolution de la population est positive, la population de départ dans les ViKs n’est pas très importante (< 250 000 hab). L’augmentation de la population des ViKs n’est pas liée à une augmentation régionale réelle de la population (voir graphique). Elle est artificielle, et est liée à des fusions de ViKs municipaux dans le ViK régional. Pour ceux que j’ai pu lister :
Les ViKs de Stambolovo et de Dimitrovgrad ont fusionné avec celui de Haskovo ;
Le ViK de Petrich avec Blagoevgrad ;
Le ViK de Kjustendil avec Dupnica
Les ViKs de Panagjurishte et Strelcha avec Pazardzhik, Svishtov avec Veliko Tarnovo
Le ViK de Sevlievo avec Gabrovo
Le ViK de Botevgrad avec la région de Sofia.
À Pleven et à Razgrad il y a aussi eu des fusions mais la trop forte baisse de la population générale n’a pas permis de compenser la hausse de population dans ces ViKs
summary(ViK$EVO_POP_VIK_13.21)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## -68201 -13065 -7238 -962 19906 67725sd(ViK$EVO_POP_VIK_13.21)## [1] 30576sd(ViK$EVO_POP_VIK_13.21)/mean(ViK$EVO_POP_VIK_13.21)*100## [1] -3177ggplot(data=ViK, aes(x = reorder(REG_FR, EVO_POP_VIK_13.21), y= EVO_POP_VIK_13.21)) +
geom_bar(stat = "identity", width=0.5, color="black", fill="white")+
geom_text(aes(label=EVO_POP_VIK_13.21), vjust=-0.3, size=2)+
theme(axis.text.x = element_text(face="bold", color="black",size=7, angle=45),
axis.text.y = element_text(face="bold", color="black",size=7, angle=90))+
labs(x="ViK", y="Évolution population des ViKs", title="Évolution de la population des ViKs")+theme(plot.title=element_text(face="bold", color="black", size=10))+theme(axis.title.x = element_text(face="plain", color="black", size=9))+theme(axis.title.y = element_text(face="plain", color="black", size=9))
Relation entre l’évolution de la population des ViKs et l’évolution générale de la population des régions
ggplot(data=ViK, aes(y = EVO_POP_VIK_13.21, x= PEVO_POP_01.21))+
geom_point()+
geom_text(aes(label=REG_FR), vjust=-0.3, size=2)+
labs(x="Population au sein des régions : évolution 2001-2021", y="Évolution de la population du ViK entre 2001 et 2021", title="Relation entre l'évolution de la population au sein des ViKs entre 2013 et 2021 et l'évolution démographique régionale entre 2001 et 2021")+theme(plot.title=element_text(face="bold", color="black", size=10))+geom_vline(aes(xintercept=0),color="red", linetype="dashed", size=1)+geom_hline(aes(yintercept=0), color="red", linetype="dashed", size=1)## Warning: Using `size` aesthetic for lines was deprecated in ggplot2 3.4.0.
## ℹ Please use `linewidth` instead.
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.
Évolution de la population et fusion des ViKs
EvoVIK <- data.frame(ViK$REG_FR, ViK$EVO_POP_VIK_13.21, ViK$FUSION_ViK)
EvoVIK <- EvoVIK[c(26,14,24,23,18,16,3,27,1,25), ]
row.names(EvoVIK)<-EvoVIK$ViK.REG_FR
EvoVIK <- EvoVIK[, c(2,3)]
EvoVIK## ViK.EVO_POP_VIK_13.21 ViK.FUSION_ViK
## Kardzhali 9640
## Sofia region 13767 Botevgrad
## Sofia capitale 19156
## Gabrovo 22157 Sevlievo
## Veliko Tarnovo 26029 Svishtov
## Targovishte 26083
## Pazardzhik 32474 Panagjurishte, Strelcha
## Kjustendil 45072 Dupnica
## Blagoevgrad 55768 Petric
## Haskovo 67725 Dimitrovgrad + StambolovoÉvolution de la population dans les oblasts bulgares entre 2001 et 2021
ggplot(data=ViK, aes(x = reorder(REG_FR, PEVO_POP_01.21), y= PEVO_POP_01.21)) +
geom_bar(stat = "identity", width=0.5, color="black", fill="white")+
geom_text(aes(label=PEVO_POP_01.21), vjust=-0.3, size=2)+
theme(axis.text.x = element_text(face="bold", color="black",size=7, angle=45),
axis.text.y = element_text(face="bold", color="black",size=7, angle=90))+
labs(x="ViK", y="Évolution population", title="Évolution de la population des oblasts entre 2001 et 2021")+theme(plot.title=element_text(face="bold", color="black", size=10))+theme(axis.title.x = element_text(face="plain", color="black", size=9))+theme(axis.title.y = element_text(face="plain", color="black", size=9))
1.A.b. Kilomètres de réseau
Kilomètres de réseau
Le nombre de kilomètres de réseau par ViK varie entre 1 150 km (à Pazardzhik) et 4 344 km à Plovdiv, suivi de près par Varna, Burgas. Le réseau de Sofia capitale est moins étendu, donc bien plus dense, que celui des autres régions du fait de la superficie réduite de la région-capitale et de son caractère essentiellement urbain.La médiane est de 2 085 km et le coefficient de variation de 40% avec un écart type de 991 par rapport à une moyenne de 2 470 km
summary(ViK$Km_reseau_eau_21)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 1150 1748 2085 2470 3292 4344sd(ViK$Km_reseau_eau_21)## [1] 991sd(ViK$Km_reseau_eau_21)/mean(ViK$Km_reseau_eau_21)*100## [1] 40ggplot(data=ViK, aes(x = reorder(REG_FR,Km_reseau_eau_21), y= Km_reseau_eau_21)) +
geom_bar(stat = "identity", width=0.5, color="black", fill="white")+
geom_text(aes(label=Km_reseau_eau_21), vjust=-0.3, size=2)+
theme(axis.text.x = element_text(face="bold", color="black",size=7, angle=45),
axis.text.y = element_text(face="bold", color="black",size=7, angle=90))+
labs(x="ViK", y="km réseau", title="Kilomètres de réseau d'eau potable des ViKs en 2021")+theme(plot.title=element_text(face="bold", color="black", size=10))+theme(axis.title.x = element_text(face="plain", color="black", size=9))+theme(axis.title.y = element_text(face="plain", color="black", size=9))
Densité de population sur le réseau
La densité de population sur le réseau mesure le nombre d’abonné·es par km de réseau sur un territoire donné. Une faible densité de population sur le réseau est identifié comme un potentiel indicateur de surdimentionnement des réseaux, en particulier si la population diminue. Plus la population sur un réseau est dense, moins le km de réseau est coûteux à entretenir. Des problèmes techniques peuvent être présents sur un réseau dense et sur un réseau peu dense. Sur un réseau dense, le nombre de cassures pour 100 km de réseau peut augmenter, car beaucoup de cassures ont lieu au niveau des embranchements individuels (joints qui sautent) Sur un réseau peu dense, les fuites sont moins accessibles, plus difficiles à identifier et à réparer. Dans les deux cas, les problèmes techniques ont un impact économique.
La densité d’abonné·es sur le réseau varie entre 50,4 hab/km à Vidin et 359,2 hab/km à Sofia. Sofia représente une valeur extrême. Les autres densités élevées se situent dans les régions de Pazardzhik, de Plovdiv et de Blagoevgrad qui ont environ 150 hab./km de réseau. La moyenne est de 91,9 hab.km et la médiane de 73,6.hab.km
Pazardzhik et Blagoevgrad ont un nombre de km de réseau très faible et un nombre d’abonné·es peu élevé, mais plus important que la médiane des autres ViKs. Plovdiv, lui, a un nombre de km de réseau élevé mais, comme Sofia, un grand nombre d’abonné·es. Globalement, on a une corrélation entre le nombre de km de réseau et la population, même en enlevant la valeur extrême de Sofia. Elle est de 0,68 avec une p.value de moins de 0,05 (Spearman). Un grand nombre de ViK ont une population peu importante (< 200 000 hab) et un réseau de petite taille à entretenir (< 2000 km de canalisations d’eau potable + eaux usées). Les ViKs potentiellement les plus en difficulté pour entretenir leur réseau sont ceux qui ont une faible population et un réseau étendu (notamment Shumen, Vraca et Dobrich).
summary(ViK$dens_pop_reseau21)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 50.4 59.8 73.6 91.9 96.4 359.2sd(ViK$dens_pop_reseau21)## [1] 60sd(ViK$dens_pop_reseau21)/mean(ViK$dens_pop_reseau21)*100## [1] 66ggplot(data=ViK, aes(x = reorder(REG_FR,dens_pop_reseau21), y= dens_pop_reseau21)) + geom_bar(stat = "identity", width=0.5, color="black", fill="white")+ geom_text(aes(label=dens_pop_reseau21), vjust=-0.3, size=2)+ theme(axis.text.x = element_text(face="bold", color="black",size=7, angle=45), axis.text.y = element_text(face="bold", color="black",size=7, angle=90))+ labs(x="ViK", y="hab./km de réseau", title="Densité d'abonné·es sur le réseau d'eau potable en 2021")+theme(plot.title=element_text(face="bold", color="black", size=10))+theme(axis.title.x = element_text(face="plain", color="black", size=9))+theme(axis.title.y = element_text(face="plain", color="black", size=9))
ggplot(data=ViK, aes(y = POP_VIK_21, x= Km_reseau_eau_21))+
geom_point()+
geom_text(aes(label=REG_FR), vjust=-0.3, size=2)+
labs(x="kilomètres de réseau d'eau 2021", y="Population du ViK en 2021", title="Relation entre la population au sein des ViKs et les km de réseau")+theme(plot.title=element_text(face="bold", color="black", size=10))+geom_vline(aes(xintercept=median(Km_reseau_eau_21)),color="red", linetype="dashed", size=1)+geom_hline(aes(yintercept=median(POP_VIK_21), color="red"))
cor.test(ViKnoSOF$Km_reseau_eau_21, ViKnoSOF$POP_VIK_21, method=c("spearman"))##
## Spearman's rank correlation rho
##
## data: ViKnoSOF$Km_reseau_eau_21 and ViKnoSOF$POP_VIK_21
## S = 1056, p-value = 0.0002
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
## rho
## 0.68Densité de population sur le réseau et évolution démographique régionale
On utilise le test de Pearson car on a des exaequos , malgré le fait que la distribution des données ne soit pas tout à fait linéaire (mais elle reste monotone). Selon le test de Pearson, , la corrélation entre les deux variables est de 0.79 avec Sofia et de 0.69 sans Sofia. La p.value est chaque fois < à 0.05 ce qui signifie que la corrélation est statistiquement significative. Cela signifie que les ViK qui ont les densités de population sur les réseaux les plus faibles sont aussi ceux dont la population diminue, en pourcentage annuelle, le plus depuis 2001. La chute démographique a sans surprise un effet direct sur les réseau d’eau et d’assainissement, et l’on peut penser avec peu de doute que les régions qui ont à la fois la plus faible densité et la plus forte baisse de population sont les plus concernés par le surdimentionnement des réseaux. Ici, il s’agit très nettement du ViK de Vidin, la région du pays se dépeuplant le plus et qui a aussi la plus faible densité par habitant.
cor.test(ViK$dens_pop_reseau21, ViK$PEVO_POP_01.21, method=c("pearson"))##
## Pearson's product-moment correlation
##
## data: ViK$dens_pop_reseau21 and ViK$PEVO_POP_01.21
## t = 6, df = 26, p-value = 0.0000007
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.59 0.90
## sample estimates:
## cor
## 0.79cor.test(ViKnoSOF$dens_pop_reseau21, ViKnoSOF$PEVO_POP_01.21, method=c("pearson"))##
## Pearson's product-moment correlation
##
## data: ViKnoSOF$dens_pop_reseau21 and ViKnoSOF$PEVO_POP_01.21
## t = 5, df = 25, p-value = 0.00006
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.42 0.85
## sample estimates:
## cor
## 0.69ggplot(data=ViK, aes(y = dens_pop_reseau21, x= PEVO_POP_01.21))+
geom_point()+
geom_text(aes(label=REG_FR), vjust=-0.3, size=2)+
labs(x="Pourcentage annuel moyen d'évolution de la population entre 2001-2021", y="Densité d'abonné·es sur le réseau, 2021", title="Relation entre l'évolution de la population dans les ViK et la densité d'abonné·es sur le réseau")+theme(plot.title=element_text(face="bold", color="black", size=10))+geom_vline(aes(xintercept=median(PEVO_POP_01.21)),color="red", linetype="dashed", size=1)+ geom_hline(aes(yintercept=median(dens_pop_reseau21)), color="red", linetype="dashed", size=1)
–> Réalisation d’une carte densité d’abonnés sur le réseau et évolution annuelle moyenne de la population
1.A.c. Prélèvements en eau des ViKs
Prélèvements en eau bruts
Les prélèvements en eau sont très inégaux selon les ViKs, de 6,1 millions de m3 par an (Vidin) à 142 millions de m3 par an (Sofia Capitale), avec une médiane à 20 millions, une moyenne à 29,7 millions et un écart type de 28. Le coefficient de variation est de 96%.
summary(ViK$EAUprel_21/1000000)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 6.1 12.4 20.2 29.8 32.7 142.4sd(ViK$EAUprel_21/1000000)## [1] 28sd(ViK$EAUprel_21)/mean(ViK$EAUprel_21)*100## [1] 96ggplot(data=ViK, aes(x = reorder(REG_FR,EAUprel_21/1000000), y= EAUprel_21/1000000)) +
geom_bar(stat = "identity", width=0.5, color="black", fill="white")+
geom_text(aes(label=EAUprel_21/1000000), vjust=-0.3, size=2)+
theme(axis.text.x = element_text(face="bold", color="black",size=7, angle=45),
axis.text.y = element_text(face="bold", color="black",size=7, angle=90))+
labs(x="ViK", y="millions de m3/an", title="Total des prélèvements en eau pour le réseau en 2021")+theme(plot.title=element_text(face="bold", color="black", size=10))+theme(axis.title.x = element_text(face="plain", color="black", size=9))+theme(axis.title.y = element_text(face="plain", color="black", size=9))
Comme le montre les matrices de corrélation avec ou sans Sofia (spearman), les prélèvements en eau potable sont, en valeur brute, fortement corrélés à la population du ViK en nombre d’habitant, au revenu du ViK, à l’évolution de la population sur le territoire, au nombre de km de réseau, à la densité de population desservie ; ils sont négativement corrélés à la population rurale, et à la valeur ajouté agricole régionale. Il y a un lien logique entre le fait que plus un réseau dessert d’habitant·es, plus les prélèvements en eau soient importants. Ceux ci sont aussi corrélés (mais de façon moins importante) au PIB / habitant des régions et probablement de façon très importante au PIB globale, qui n’apparaît pas dans les indicateurs.
Dans la matrice de corrélation sans Sofia, les prélèvements en eau sont aussi corrélés (faiblement : 0.37) aux pertes en eau. Cette corrélation apparaît bien plus forte si l’on compare les prélèvements en eau des ViKs par habitant aux pertes en eau (et non plus les prélèvements en eau bruts).
Prélèvements en eau par habitant desservi
Les prélèvements en eau par habitant vont de 70,5 m3/hab (à Kardzhali) à 235 m3/hab (à Dobrich), avec une médiane de 112 m3/hab et une moyenne de 126 m3/h. La dispersion des données autour de la moyenne est plus faible que pour les prélèvements bruts, avec un écart type de 46 et un coefficient de variation de 37%. Les ViKs qui prélèvent le plus d’eau par habitant ne sont pas les mêmes que ceux qui prélèvent le plus d’eau. Dobrich, Shumen, Sliven, Pernik et Montana sont les plus consommatrices en eau par habitant. Kardzhali, Vidin, Haskovo et Smoljan ont en revanche de plus faibles prélèvements .
summary(ViK$EAUprelHab_21)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 70.5 98.9 112.0 125.5 140.6 235.2sd(ViK$EAUprelHab_21)## [1] 46sd(ViK$EAUprelHab_21)/mean(ViK$EAUprelHab_21)*100## [1] 37ggplot(data=ViK, aes(x = reorder(REG_FR,EAUprelHab_21), y= EAUprelHab_21)) +
geom_bar(stat = "identity", width=0.5, color="black", fill="white")+
geom_text(aes(label=EAUprelHab_21), vjust=-0.3, size=2)+
theme(axis.text.x = element_text(face="bold", color="black",size=7, angle=45),
axis.text.y = element_text(face="bold", color="black",size=7, angle=90))+
labs(x="ViK", y="m3/hab an", title="Total des prélèvements en eau pour le résau par habitant en 2021")+theme(plot.title=element_text(face="bold", color="black", size=10))+theme(axis.title.x = element_text(face="plain", color="black", size=9))+theme(axis.title.y = element_text(face="plain", color="black", size=9))
Eau prélevée par habitant et pertes en eau
Il existe une corrélation entre l’eau prélevée par habitant et les pertes en eau, de 0,75 selon le coefficient de Spearman. Le test de Spearman pour ces deux variables détaillé ci-dessous montre une relation significative puisque la p.value est \< 0,05 (Pearson : coefficient de corrélation = 0.84).
cor.test(ViK$Perteseau_Pourc_VIK21, ViK$EAUprelHab_21, method=c("spearman"))##
## Spearman's rank correlation rho
##
## data: ViK$Perteseau_Pourc_VIK21 and ViK$EAUprelHab_21
## S = 920, p-value = 0.000009
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
## rho
## 0.75cor.test(ViK$Perteseau_Pourc_VIK21, ViK$EAUprelHab_21, method=c("pearson"))##
## Pearson's product-moment correlation
##
## data: ViK$Perteseau_Pourc_VIK21 and ViK$EAUprelHab_21
## t = 8, df = 26, p-value = 0.00000003
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.67 0.92
## sample estimates:
## cor
## 0.84Le nuage de point ci-dessous représente cette relation statistique qui est monotone, plutôt linéaire et positive entre les variables. Ainsi, globalement, plus le pourcentage de pertes en eau est élevé, plus le prélèvement en eau par habitant est lui aussi élevé. Ceci est assez logique, dans la mesure où les pertes en eau représentent des fuites dont l’augmentation des prélèvements vise à limiter les effets.
ggplot(data=ViK, aes(y = EAUprelHab_21, x= Perteseau_Pourc_VIK21))+
geom_point()+
geom_text(aes(label=REG_FR), vjust=-0.3, size=2)+
labs(x="pourcentage de pertes en eau", y="eau prélevée par an (m3/hab)", title="Relation entre les pertes en eau (%) et le volume de prélèvement en eau par habitant en 2021")+theme(plot.title=element_text(face="bold", color="black", size=10))+geom_hline(aes(yintercept=median(EAUprelHab_21), color="red"))+geom_vline(aes(xintercept=median(Perteseau_Pourc_VIK21), color="red"))
Les régions ayant à la fois le plus haut niveau de pertes en eau et le plus haut niveau de prélèvement par habitant sont Sliven, Sumen, Dobrich, Pernik, Montana et Jambol. Les conséquences sont écologiques, mais aussi économiques, les prélèvements en eau coûtant cher, surtout si la région est approvisionnée par pompage et a donc des consommations énergétiques élevées. Ici, le nuage de point suivant permet de voir que Dobrich et Shumen sont particulièrement concernés et, dans une moindre mesure, Jambol et Razgrad
ggplot(data=ViK, aes(y = EAUprelHab_21, x= Energie.kwh.m3_21ViK))+
geom_point()+
geom_text(aes(label=REG_FR), vjust=-0.3, size=2)+
labs(x="Facture énergétique", y="eau prélevée par an (m3/hab)", title="Relation entre les factures énergétiques et le volume de prélèvement en eau par habitant en 2021")+theme(plot.title=element_text(face="bold", color="black", size=10))+geom_hline(aes(yintercept=median(EAUprelHab_21), color="red"))+geom_vline(aes(xintercept=median(Energie.kwh.m3_21ViK), color="red"))
1.A.d. Consommations en eau par habitant (L/hab/J)
Consommations par habitant
Selon les régions, les consommations en eau par habitant varient entre 74 litres d’eau par jour (à Targovishte) et 126 litres par jour (à Sofia), ce qui est important pour ce type d’indicateur. La moyenne de consommation par habitant est de 99 et la médiane est de 98. On observe un lien fort entre tous les indicateurs d’urbanité,de revenu et la consommation en eau par habitant par jour, avec ou sans Sofia, qui représente une valeur extrême.
Globalement donc, les habitants des régions rurales et agricoles tendent à consommer moins d’eau individuellement que les habitants des régions urbaines. Cela peut être lié :
- au mode de vie et à l’équipement en appareil électroménagers
- au revenu globalement plus élevé de la population urbaine qui consomme plus
- à l’existence de davantage sources alternatives à l’eau du robinet en zone rurale (puits, sources naturelles) [outre l’eau en bouteille, qui peut être achetée partout]
- à une moindre confiance en zone rurale pour l’eau du réseau, réputée de moins bonne qualité
summary(ViK$Conso_eau_hab_23)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 74.0 89.5 98.5 99.1 109.2 126.0sd(ViK$Conso_eau_hab_23)## [1] 14sd(ViK$Conso_eau_hab_23)/mean(ViK$Conso_eau_hab_23)*100## [1] 14ggplot(data=ViK, aes(x = reorder(REG_FR,Conso_eau_hab_23), y= Conso_eau_hab_23)) +
geom_bar(stat = "identity", width=0.5, color="black", fill="white")+
geom_text(aes(label=Conso_eau_hab_23), vjust=-0.3, size=2)+
theme(axis.text.x = element_text(face="bold", color="black",size=7, angle=45),
axis.text.y = element_text(face="bold", color="black",size=7, angle=90))+
labs(x="ViK", y="m3/hab an", title="Consommations d'eau annuelles par habitant en 2023, en m3")+theme(plot.title=element_text(face="bold", color="black", size=10))+theme(axis.title.x = element_text(face="plain", color="black", size=9))+theme(axis.title.y = element_text(face="plain", color="black", size=9))
Consommation en eau par habitant et urbanisation
Ce sont en particulier les régions rurales et agricoles qui ont une faible consommation.
Blagoevgrad et Pazardzhik ont une faible VA agricole et ont une consommation en eau et par habitant plus élevée. Cela peut être lié au mode de vie de populations rurales, mais urbanisées par l’industrialisation.
ggplot(data=ViK, aes(y = Conso_eau_hab_23, x= PopRural_24))+
geom_point()+
geom_text(aes(label=REG_FR), vjust=-0.3, size=2)+
labs(x="% pop rurale", y="consommations en eau par habitant L/eau/j", title="Relation entre la ruralité d'une région et les consommations par hab par jour")+theme(plot.title=element_text(face="bold", color="black", size=10))+geom_hline(aes(yintercept=median(Conso_eau_hab_23), color="red"))+geom_vline(aes(xintercept=median(PopRural_24), color="red"))
ggplot(data=ViK, aes(y = Conso_eau_hab_23, x= AgriVA_23))+
geom_point()+
geom_text(aes(label=REG_FR), vjust=-0.3, size=2)+
labs(x="% agriculture dans la VA", y="consommations en eau par habitant par j", title="Relation entre l'agriculture et les consommations en eau par hab")+theme(plot.title=element_text(face="bold", color="black", size=10))+geom_hline(aes(yintercept=median(Conso_eau_hab_23), color="red"))+geom_vline(aes(xintercept=median(AgriVA_23), color="red"))
Consommations par habitant et minorités
Une autre corrélation négative est entre la consommation en eau par jour et par habitant et le pourcentage de population non bulgare dans les régions. Cela est lié à la fois au fait que ces populations vivent dans les régions les plus agricoles, mais aussi peut-être, comme à Sliven, au niveau de pauvreté des populations minoritaires et au fait qu’elles sont mal raccordés.
ggplot(data=ViK, aes(y = Conso_eau_hab_23, x= NonBG_21))+
geom_point()+
geom_text(aes(label=REG_FR), vjust=-0.3, size=2)+
labs(x="% de population se déclarant non bulgare", y="consommations en eau par habitant par j", title="Relation entre le pourcentage de population se déclarant non bulgare et les consommations en eau par hab")+theme(plot.title=element_text(face="bold", color="black", size=10))+geom_hline(aes(yintercept=median(Conso_eau_hab_23), color="red"))+geom_vline(aes(xintercept=median(NonBG_21), color="red"))
1.A.e Les revenus des ViKs
Revenus bruts
Les revenus des ViKs ne sont pas les bénéfices. Ce
sont uniquement les revenus générés par le ViK, dont on a vu dans le
chapitre 6 qu’ils dépendaient essentiellement des factures d’eau, donc
de la consommation et de la capacité de paiement des consommateurs. Un
ViK qui a à la fois des dépenses importantes (en particulier
énergétiques, mais aussi pour les prélèvements) et des revenus faibles
sera donc en difficulté.
Le revenu le plus bas est de 7 millions de leva à Vidin, suivi de près
par Kjustendil et par Smoljan.Le revenu le plus élevé est de 196
millions de leva à Sofia, suivi par Burgas à 93 millions de leva
(amplitude de 189 millions de leva).L’écart type de 39 est plus élevé
que la moyenne à 29 millions de leva par ViK. On a donc une très forte
dispersion des données (coefficient de variation = 132%) et des valeurs
extrêmes puisque 75% des ViKs ont entre 7 et 26 millions de leva par
an.
summary(ViK$REV_VIK_21/1000000)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 7.19 11.21 15.81 29.20 25.97 196.24sd(ViK$REV_VIK_21/1000000)## [1] 39sd(ViK$REV_VIK_21)/mean(ViK$REV_VIK_21)*100## [1] 132ggplot(data=ViK, aes(x = reorder(REG_FR,REV_VIK_21/1000000), y= REV_VIK_21/1000000)) +
geom_bar(stat = "identity", width=0.5, color="black", fill="white")+
geom_text(aes(label=REV_VIK_21/1000000), vjust=-0.3, size=2)+
theme(axis.text.x = element_text(face="bold", color="black",size=7, angle=45),
axis.text.y = element_text(face="bold", color="black",size=7, angle=90))+
labs(x="ViK", y="millions de leva", title="Revenus des ViKs en 2021")+theme(plot.title=element_text(face="bold", color="black", size=10))+theme(axis.title.x = element_text(face="plain", color="black", size=9))+theme(axis.title.y = element_text(face="plain", color="black", size=9))
Revenu et PIB régional : une corrélation notable
cor.test(ViK$REV_VIK_21, ViK$PIBhab_23, method=c("pearson"))##
## Pearson's product-moment correlation
##
## data: ViK$REV_VIK_21 and ViK$PIBhab_23
## t = 8, df = 26, p-value = 0.00000001
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.70 0.93
## sample estimates:
## cor
## 0.85cor.test(ViK$REV_VIK_21, ViK$PIBhab_23, method=c("spearman"))##
## Spearman's rank correlation rho
##
## data: ViK$REV_VIK_21 and ViK$PIBhab_23
## S = 1588, p-value = 0.002
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
## rho
## 0.57Toutefois, les corrélations ne sont plus significatives entre le PIB/hab et le revenu qu’il s’agisse du revenu par habitant des ViK ou du revenu brut des Viks si l’on retire Sofia.
corr.test(ViKnoSOF$PIBhab_23, ViKnoSOF$REV_POP_VIK_21, method=c("spearman"))## Call:corr.test(x = ViKnoSOF$PIBhab_23, y = ViKnoSOF$REV_POP_VIK_21,
## method = c("spearman"))
## Correlation matrix
## [1] 0.52
## Sample Size
## [1] 27
## These are the unadjusted probability values.
## The probability values adjusted for multiple tests are in the p.adj object.
## [1] 0.01
##
## To see confidence intervals of the correlations, print with the short=FALSE optioncorr.test(ViKnoSOF$PIBhab_23, ViKnoSOF$REV_VIK_21, method=c("spearman"))## Call:corr.test(x = ViKnoSOF$PIBhab_23, y = ViKnoSOF$REV_VIK_21, method = c("spearman"))
## Correlation matrix
## [1] 0.52
## Sample Size
## [1] 27
## These are the unadjusted probability values.
## The probability values adjusted for multiple tests are in the p.adj object.
## [1] 0.01
##
## To see confidence intervals of the correlations, print with the short=FALSE optionggplot(data=ViKnoSOF, aes(x = REV_VIK_21 , y= PIBhab_23))+
geom_point()+
geom_text(aes(label=REG_FR), vjust=-0.3, size=2)+
labs(x="Revenu des ViKs (excepté Sofia) en 2021", y="PIB par habitant en 2023 (excepté sofia)", title="Relation entre les revenus des ViKs et le PIB par habitant régional, sans Sofia")+theme(plot.title=element_text(face="bold", color="black", size=10))+geom_vline(aes(xintercept=median(REV_VIK_21)),color="red", linetype="dashed", size=1)+geom_hline(aes(yintercept=median(PIBhab_23)), color="red", linetype="dashed")
Revenu par habitant des ViKs
Les revenus par habitants s’étendent entre 59 leva/hab à Kjustendil et 227 leva/hab à Burgas, suivi de Varna et de la capitale, Sofia. L’écart type est de 33 pour une moyenne de 108 leva/hab. La mediane est elle aussi de 108 leva/hab. On a donc un coefficient de variation bien moins important que pour le revenu brut, de 31%. La matrice des corrélation montre une corrélation négative entre le revenu par habitant des ViKs et :
- l’augmentation de la population d’un ViK : en effet, la plupart des ViKs ayant vu leur population augmenter fortement n’ont pas vu leurs revenus augmenter en parallèle et ont des revenus par habitant plus bas après qu’avant la fusion (-0,5)
- le % de population rurale (-0,4) ce qui signifie que ce sont surtout les régions urbaines qui ont un revenu par habitant plus élevé
…une corrélation positive entre le rev/hab.
des ViKs et :
- le revenu global, le revenu pour 1000 km de réseau, et l’augmentation de ce revenu depuis 2013
- le PIB/hab montrant que les régions avec une plus forte activité économique dépensent aussi plus d’eau ; c’est d’autant plus vrai que les tarifs pour les entreprises sont plus chers que pour les ménages et qu’elles consomment aussi plus (surtout l’industrie)
- la taille du réseau d’eau potable en valeur brute et les prélèvements en eau (ce qui est lié à l’urbanisation surtout ici, comme le % de population ayant fait des études supérieures)
- le prix de l’eau et l’évolution du prix de l’eau : plus il est élevé, plus le revenu par tête d’un ViK est élevé
On a donc un lien entre les ViKs de régions urbaines et peuplées et le revenu par habitant qu’ils génèrent et, de l’autre côté, les ViKs de régions rurales et leur revenu par habitant.
summary(ViK$REV_POP_VIK_21)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 59.1 84.6 108.4 107.5 119.2 226.6sd(ViK$REV_POP_VIK_21)## [1] 33sd(ViK$REV_POP_VIK_21)/mean(ViK$REV_POP_VIK_21)*100## [1] 31ggplot(data=ViK, aes(x = reorder(REG_FR,REV_POP_VIK_21), y= REV_POP_VIK_21)) +
geom_bar(stat = "identity", width=0.5, color="black", fill="white")+
geom_text(aes(label=REV_POP_VIK_21), vjust=-0.3, size=2)+
theme(axis.text.x = element_text(face="bold", color="black",size=7, angle=45),
axis.text.y = element_text(face="bold", color="black",size=7, angle=90))+
labs(x="ViK", y="leva/hab", title="Revenus par habitant desservi des ViKs en 2021")+theme(plot.title=element_text(face="bold", color="black", size=10))+theme(axis.title.x = element_text(face="plain", color="black", size=9))+theme(axis.title.y = element_text(face="plain", color="black", size=9)) 
Revenus par habitant et dépenses énergétiques
Les dépenses énergétiques sont l’une des principales dépenses des ViKs. Il n’y a pas de corrélation entre le revenu par habitant des ViKs et leurs dépenses énergétiques. Toutefois : certains ViK à faible revenu doivent-ils faire face à davantage de dépenses que les autres ?
Au dessus de la médiane, se distingue nettement le cas de Haskovo, avec un très faible revenu par habitant des ViKs et de très fortes dépenses énergétiques ; puis de Razgrad, Silistra, Targovishte et Vidin. D’autres ViKs comme Burgas, Varna et Sofia-capitale, ont inversement de faibles consommations énergétiques et un plus haut revenu des ViKs par habitant
ggplot(data=ViK, aes(y = Energie.kwh.m3_21ViK, x=REV_POP_VIK_21))+
geom_point()+
geom_text(aes(label=REG_FR), vjust=-0.3, size=2)+
labs(y="consommation énergétique en KwH m3", x="revenu par hab. du ViK", title ="Relation entre les consommations énergétiques et la densité de population sur le réseau")+theme(plot.title=element_text(face="bold", color="black", size=10))+geom_vline(aes(xintercept=median(REV_POP_VIK_21)),color="red", linetype="dashed", size=1)+geom_hline(aes(yintercept=median(Energie.kwh.m3_21ViK)), color="red", linetype="dashed", size=1)
Revenu par habitant et prix de l’eau
On voit ici se dessiner une légère relation positive et monotone entre le prix de l’eau et le revenu par habitant des ViKs. Le coefficient de corrélation est assez faible, 0,38 et la p.value est de 0,05 donc l’hypothèse de relation est validée. Les ViKs où le revenu par habitant est le plus élevé sortent de la tendance générale :
- Sofia capitale et Dobrich ont des prix de l’eau faibles par rapport à leur revenu par habitant élevé ; dans une moindre mesure c’est aussi le cas pour Burgas et Varna
- Ici, on voit sans doute le rôle de l’activité économique (industrie, et surtout tourisme dans ces régions situées sur la cote de la Mer Noire ou dans la région capitale).
Cela réhausse les consommations donc le revenu par habitant, sans nécessairement le faire dépendre du prix de l’eau.
cor.test(ViK$PRIX_ViK_23, ViK$REV_POP_VIK_21, method = c("spearman")) ##
## Spearman's rank correlation rho
##
## data: ViK$PRIX_ViK_23 and ViK$REV_POP_VIK_21
## S = 2270, p-value = 0.05
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
## rho
## 0.38ggplot(data=ViK, aes(x = PRIX_ViK_23 , y= REV_POP_VIK_21))+
geom_point()+
geom_text(aes(label=REG_FR), vjust=-0.3, size=2)+
labs(x="prix de l'eau en 2023", y="revenu du ViK par habitant (leva) en 2021", title="Relation entre le prix de l'eau au m3 et le revenu des ViK par habitant")+theme(plot.title=element_text(face="bold", color="black", size=10))+geom_vline(aes(xintercept=median(PRIX_ViK_23)),color="red", linetype="dashed", size=1)+geom_hline(aes(yintercept=median(REV_POP_VIK_21)), color="red", linetype="dashed", size=1)
Évolution des revenus des ViKs
L’évolution des revenus des ViK varie entre -1 millions de leva constants et +45 millions de leva avec un écart type de 10 millions et un coefficient de variation de +169%.
- L’évolution des revenus en valeur brute est négativement corrélée à l’agriculture en valeur ajouté, au % de population rurale et à l’énergie en KWH dépensée pour le réseau.
- Corrélation positive avec la population des ViKs (les ViKs les plus peuplés ont tendance à s’enrichir), avec leur revenu total en 2023 et avec tous les indicateurs désignant des régions urbaines
summary(ViK$EVO_REVcon_VIK_31.21/1000000)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## -1.02 1.95 2.58 6.09 4.07 45.18sd(ViK$EVO_REVcon_VIK_31.21/1000000)## [1] 10sd(ViK$EVO_REVcon_VIK_31.21)/mean(ViK$EVO_REVcon_VIK_31.21)*100## [1] 169ggplot(data=ViK, aes(x = reorder(REG_FR,EVO_REVcon_VIK_31.21/1000000), y= EVO_REVcon_VIK_31.21/1000000)) +
geom_bar(stat = "identity", width=0.5, color="black", fill="white")+
geom_text(aes(label=EVO_REVcon_VIK_31.21/1000000), vjust=-0.3, size=2)+
theme(axis.text.x = element_text(face="bold", color="black",size=7, angle=45),
axis.text.y = element_text(face="bold", color="black",size=7, angle=90))+
labs(x="ViK", y="millions de leva constant", title=" Évolution des revenus des ViKs entre 2013 et 2021")+theme(plot.title=element_text(face="bold", color="black", size=10))+theme(axis.title.x = element_text(face="plain", color="black", size=9))+theme(axis.title.y = element_text(face="plain", color="black", size=9))+geom_hline(yintercept=0, color="red")
Évolution des revenus (leva constants) et consommations énergétiques
Sur la matrice de corrélations, on observe une corrélation négative,
de -0,7, entre les consommations énergétiques d’un ViK et ses
revenus.
Plus un ViK consomme d’énergie, plus l’évolution des revenus est faible
(voire négative). Les ViKs ayant une forte consommation énergétique ont
aussi une faible évolution de leurs revenus en % du revenu de 2013. Il
s’agit surtout de : Haskovo, Razgrad, Silistra, Dobrich, Ruse, Stara
Zagora, Shumen, Targovishte, Lovech, Vidin et Yambol.
Pourtant, les consommations énergétiques n’entrent pas dans le calcul
des revenus. Il s’agit sans doute là d’une corrélation liée au fait que
ce sont les régions les plus rurales et agricoles qui, en Bulgarie,
utilisent le plus d’énergie pour les services d’eau potable, car elles
ont peu de ressource en eau et/ou dépendent du pompage (développement
historique de la Dobrudzha notamment + région de Haskovo).
Les régions plus urbaines et industrielles, qui se sont aussi
historiquement développées par la construction de barrages et autour de
ressources en eau importantes (zones montagneuses, souvent), subissent
moins ce facteur supplémentaire du coût de l’énergie .
p.value \< 0,05
cor.test(ViK$EVO_REVcon_VIK_pourc_13.21, ViK$Energie.kwh.m3_21ViK, method=c("spearman"))## Warning in cor.test.default(ViK$EVO_REVcon_VIK_pourc_13.21,
## ViK$Energie.kwh.m3_21ViK, : Impossible de calculer la p-value exacte avec des
## ex-aequos##
## Spearman's rank correlation rho
##
## data: ViK$EVO_REVcon_VIK_pourc_13.21 and ViK$Energie.kwh.m3_21ViK
## S = 6286, p-value = 0.00002
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
## rho
## -0.72cor.test(ViK$EVO_REVcon_VIK_pourc_13.21, ViK$Energie.kwh.m3_21ViK, method=c("pearson"))##
## Pearson's product-moment correlation
##
## data: ViK$EVO_REVcon_VIK_pourc_13.21 and ViK$Energie.kwh.m3_21ViK
## t = -3, df = 26, p-value = 0.003
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.76 -0.20
## sample estimates:
## cor
## -0.53ggplot(data=ViK, aes(y =EVO_REVcon_VIK_pourc_13.21 , x= Energie.kwh.m3_21ViK))+
geom_point()+
geom_text(aes(label=REG_FR), vjust=-0.3, size=2)+
labs(x="Consommation énergétique kwh/m3", y="% évolution des revenus entre 2013 et 2021 (constants)", title="Relation entre la consommation énergétique des ViKs en 2021 et l'importance de l'évolution de leurs revenus entre 2013 et 2021")+theme(plot.title=element_text(face="bold", color="black", size=10))+
geom_vline(aes(xintercept=median(Energie.kwh.m3_21ViK)),color="red", linetype="dashed", size=1)+
geom_hline(aes(yintercept=median(EVO_REVcon_VIK_pourc_13.21)), color="red", linetype="dashed", size=1)
1.A.f. Le prix de l’eau
Prix de l’eau
-Quels ViK ont les prix de l’eau les plus élevés ?
Et dans les régions où les prix de l’eau sont élevés, le niveau de vie général est-il lui aussi élevé ?
En 2023, les prix de l’eau s’étalent entre 1,4 leva le m3 (Pazardzhik) et 4,5 leva le m3 (Razgrad), soit trois fois plus. La médiane et la moyenne sont proches, entre 3,1 et 3,7 et la variation des prix entre les ViKs est moyen, avec un coefficient de variation de 27% et un écart type de 0,83.
Cela est sans doute lié à la régulation des prix de l’eau, qui ne permet pas aux ViK de les augmenter au delà d’un seuil de tolérance sociale, fixé par rapport au revenu moyen des régions.
summary(ViK$PRIX_ViK_23)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 1.37 2.49 3.24 3.10 3.71 4.49sd(ViK$PRIX_ViK_23)## [1] 0.83sd(ViK$PRIX_ViK_23)/mean(ViK$PRIX_ViK_23)*100## [1] 27ggplot(data=ViK, aes(x = reorder(REG_FR,PRIX_ViK_23), y= PRIX_ViK_23)) +
geom_bar(stat = "identity", width=0.5, color="black", fill="white")+
geom_text(aes(label=PRIX_ViK_23), vjust=-0.3, size=2)+
theme(axis.text.x = element_text(face="bold", color="black",size=7, angle=45),
axis.text.y = element_text(face="bold", color="black",size=7, angle=90))+
labs(x="ViK", y="leva", title="Prix de l'eau au m3 en leva en 2023")+theme(plot.title=element_text(face="bold", color="black", size=10))+theme(axis.title.x = element_text(face="plain", color="black", size=9))+theme(axis.title.y = element_text(face="plain", color="black", size=9))
Évolution des prix
Ces 10 dernières années, les prix de l’eau ont globalement augmenté en Bulgarie, en valeur courante comme constante.
En valeur courante, les prix ont, selon les ViK, varié de 0 leva à +2,42 leva entre 2013 et 2023 dans les ViKs bulgares.
Certains ViK comme Kjustendil et Pazardzhik n’ont eu aucune, ou presque augune, augmentation du prix de l’eau,.
Silistra est la région où les prix ont le plus augmenté.
L’écart type est de 0,71 et le coefficient de variation de 56%, avec une augmentation moyenne des prix de 1,25 (et médiane de +1,27).
En valeur constante (indice des prix décembre 2013-décembre 2023 = 143), les prix ont évolué de 0 leva à +1,70 leva, avec un écart type de 0,49 et un coefficient de variation de 56%.
Moyenne : +0.88 leva par ViK et médiane : +0.89 leva.
summary(ViK$EVO_PRIX_cou_13.23/143*100)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.000 0.558 0.888 0.875 1.280 1.695sd(ViK$EVO_PRIX_cou_13.23/143*100)## [1] 0.49sd(ViK$EVO_PRIX_cou_13.23/143*100)/mean(ViK$EVO_PRIX_cou_13.23/143*100)*100## [1] 56ggplot(data=ViK, aes(x =reorder(REG_FR,EVO_PRIX_cou_13.23), y= EVO_PRIX_cou_13.23)) +
geom_bar(stat = "identity", width=0.5, color="black", fill="white")+
geom_text(aes(label=EVO_PRIX_cou_13.23), vjust=-0.3, size=2)+
theme(axis.text.x = element_text(face="bold", color="black",size=7, angle=45),
axis.text.y = element_text(face="bold", color="black",size=7, angle=90))+
labs(x="ViK", y="leva", title="évolution du prix de l'eau au m3 en leva entre 2013 et 2023")+theme(plot.title=element_text(face="bold", color="black", size=10))+theme(axis.title.x = element_text(face="plain", color="black", size=9))+theme(axis.title.y = element_text(face="plain", color="black", size=9))
Corrélations avec le prix de l’eau
La matrice de corrélation (Spearman) montre **des corrélations positives entre le prix de l’eau et :**
- L’évolution du prix de l’eau (+0,9) = les ViKs où le prix de l’eau a le plus augmenté ces 10 dernières années en valeur sont ceux où il est aujourd’hui le plus élevé
- Les consommations énergétiques
- Le revenu par habitant des ViKs (mais pas leurs revenus globaux!).
Sans doute, la hausse du prix de l’eau joue-t-elle sur les revenus par habitant, mais elle ne compense pas, au niveau des revenus globaux, les diminutions de population et de consommation.
**…et des corrélations négatives entre le prix de l’eau et :**
- L’augmentation de la population sur le territoire d’un ViK
- pour Sofia, seul VIK dont la population augmente réellement et non par l’effet d’une fusion, le prix de l’eau est relativement bas par rapport aux autres ViK
- Dans les autres régions, où il y a eu une fusion, les prix de l’eau sont assez faibles et ont peu augmenté.
En effet : il est politiquement peu populaire d’augmenter le prix de l’eau l’année de la fusion, c’est l’un des points de contestation majeur de cette politique.
- La consommation en eau des habitants.
- Il pourrait y avoir un lien direct : plus le prix de l’eau est fort, plus la consommation des habitants diminue ; plus la consommation en eau des habitants est faible, plus les ViKs compensent en augmentant le prix de l’eau ; ou inversement plus les prix de l’eau sont faibles, plus la population consomme. Mais on sait que les consommations sont peu élastiques par rapport aux prix concernant l’eau potable. Il n’est donc pas certain que ce lien soit si direct.
- Déjà parce que seuls les prix de l’eau, et non l’*évolution* des prix de l’eau, sont corrélés à la consommation en eau par habitant. Cela suggère donc qu’une forte augmentation des prix n’a pas eu d’incidence généralisée sur la consommation en eau potable
- Les ViKs qui ont des prix de l’eau élevés sont pour beaucoup, on l’a vu, des ViK ruraux et agricoles.
On a vu que dans ce type de ViK les consommations en eau par jour et par habitant était aussi plus faibles qu’ailleurs, probablement pour des raisons culturelles et du fait de la disponibilité d’autres sources en eau (voir, ci-dessus). C’est sans doute davantage le facteur de régions rurales et agricoles qui fait le lien entre les variables, bien que des prix de l’eau élevés n’encouragent probablement pas l’usage ou la confiance dans le réseau conventionnel.
ggplot(data=ViK, aes(y =Conso_eau_hab_23 , x= PRIX_ViK_23))+
geom_point()+
geom_text(aes(label=REG_FR), vjust=-0.3, size=2)+
labs(x="Prix de l'eau en 2023", y="Consommation en eau par habitant", title="Relation entre la consommation en eau par hab et le prix de l'eau")+theme(plot.title=element_text(face="bold", color="black", size=10))+
geom_vline(aes(xintercept=median(PRIX_ViK_23)),color="red", linetype="dashed", size=1)+
geom_hline(aes(yintercept=median(Conso_eau_hab_23)), color="red", linetype="dashed", size=1)
1.A.g.Prix de l’eau et revenu moyen des populations
Il convient d’abord de souligner que le prix de l’eau moyen au m3 par mois est faible en Bulgarie au regard du revenu moyen de la population. Pour une consommation moyenne de 36m3 d’eau par an par habitant, en 2023, la facture d’eau maximale moyenne par habitant (prix : 3,1 leva/m3) représente 0,47% du revenu moyen par hab en 2024. Malgré l’inflation entre les deux années, et malgré le fait que la facture d’eau d’un ménage qui compte plusieurs enfants représente nécessairement une part plus élevé du revenu du ménage, on peut estimer que cette part reste faible.
mean(ViK$PRIX_ViK_23) #prix maximum moyen du m3 d'eau en BG = 3,1 leva## [1] 3.1mean(ViK$Conso_eau_hab_23)*365/1000 #consommation moyenne de l'eau = 36m3/an## [1] 36mean(ViK$PRIX_ViK_23)*(mean(ViK$Conso_eau_hab_23)*365/1000) #facture annuelle moyenne des bulgares : 112 leva## [1] 112mean(ViK$RevHab_24) #revenu annuel moyen des bulgares = 24 106 leva## [1] 24106(mean(ViK$PRIX_ViK_23)*(mean(ViK$Conso_eau_hab_23)*365/1000))/mean(ViK$RevHab_24)*100 #la facture d'eau maximale moyenne par habitant pour un niveau de consommation moyen de 99l/j/hab représente 0,47% du revenu moyen par habitant, soit peu## [1] 0.47Eau courante dans le budget des ménages
Création d’une variable : part de l’eau courante dans le budget d’un ménage par région
PourcEauBudget <- (ViK$PRIX_ViK_23)*(ViK$Conso_eau_hab_23*365/1000)/(ViK$RevHab_24)*100
ViK$PourcEauBudget<-PourcEauBudgetLa part du prix de l’eau dans le budget d’un individu varie entre 0,23% et 0,73% au maximum, avec une moyenne de 0,47 et une médiane à 0,44 et un coefficient de variation à 30%
summary(ViK$PourcEauBudget)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.228 0.385 0.443 0.469 0.545 0.726sd(ViK$PourcEauBudget)/mean(ViK$PourcEauBudget)*100## [1] 30ggplot(data=ViK, aes(x =reorder(REG_FR,PourcEauBudget), y= PourcEauBudget)) +
geom_bar(stat = "identity", width=0.5, color="black", fill="white")+
geom_text(aes(label=PourcEauBudget), vjust=-0.3, size=2)+
theme(axis.text.x = element_text(face="bold", color="black",size=7, angle=45),
axis.text.y = element_text(face="bold", color="black",size=7, angle=90))+
labs(x="ViK", y="%", title="Pourcentage du coût de l'eau dans le revenu moyen d'un habitant")+theme(plot.title=element_text(face="bold", color="black", size=10))+theme(axis.title.x = element_text(face="plain", color="black", size=9))+theme(axis.title.y = element_text(face="plain", color="black", size=9))+geom_hline(aes(yintercept=0.54))
À Varna, Lovech, Yambol, Razgrad, Ruse, Silistra, Burgas, Haskovo, Stara Zagora et Veliko Tarnovo, le prix de l’eau dans le budget des ménages se situe dans le dernier quartile (\> 0,54%).Cela reste limité, mais ce sont les régions où ce dernier “pèse” le plus.
Prix de l’eau et revenu moyen de la population sur un territoire
Il n’y a pas de corrélation statistique entre le prix de l’eau (ni son augmentation) et le niveau de vie global dans une région mesuré en revenu par habitant.
Le nuage de point ci-dessous montre une diversité de situations.
Parmi les régions où la population est l
a plus pauvre (\< médiane), on voit que certaines ont des prix de
l’eau plutôt faibles (Région de Sofia, Dobrich, Vidin, Veliko Tarnovo et
Kardzhali) et d’autres des prix de l’eau élevés, voire très élevés :
(Razgrad, Varna, Ruse, Haskovo, Yambol, Lovech).
ggplot(data=ViK, aes(y =RevHab_24 , x= PRIX_ViK_23))+
geom_point()+
geom_text(aes(label=REG_FR), vjust=-0.3, size=2)+
labs(x="Prix de l'eau en 2023", y="Revenu par habitant en 2024", title="Relation entre le revenu par habitant et le prix de l'eau")+theme(plot.title=element_text(face="bold", color="black", size=10))+
geom_vline(aes(xintercept=median(PRIX_ViK_23)),color="red", linetype="dashed", size=1)+
geom_hline(aes(yintercept=median(RevHab_24)), color="red", linetype="dashed", size=1)
De plus, Certaines régions au revenu moyen élevé ont aussi des taux de pauvreté et d’exclusion sociale importants. C’est particulièrement le cas de Pernik, dont une partie de la population travaille à Sofia et une autre - restée sur place - a des revenus bien moins importants ; Certaines régions assez pauvres peuvent aussi avoir des taux élevés, notamment Razgrad, Sliven, Vidin…Pour ces foyers très pauvres le prix de l’eau - même faible - peut représenter un poids très important
ggplot(data=ViK, aes(y =TxRisq_24 , x= RevHab_24))+
geom_point()+
geom_text(aes(label=REG_FR), vjust=-0.3, size=2)+
labs(x="Revenu par habitant", y="Risque de pauvreté et d'exclusion sociale en 2024", title="Relation entre le revenu par habitant et le risque de pauvreté et d'exclusion sociale")+theme(plot.title=element_text(face="bold", color="black", size=10))+
geom_vline(aes(xintercept=median(RevHab_24)),color="red", linetype="dashed", size=1)+
geom_hline(aes(yintercept=median(TxRisq_24)), color="red", linetype="dashed", size=1)
1.A.h Consommations énergétiques des ViKs
consommations énergétiques
Il y a des écarts importants entre les ViKs du point de vue de leur consommation énergétique.
Leurs dépenses vont de 0,04 kwh /m3 d’eau prélevée à 2,34 Kwh avec une moyenne à 0,74 et une médiane à 0,58 .
L’écart type est de 0,49 et le coefficient de variation de 67%.
Le maximum, de 2,34 est à Haskovo et constitue une valeur extrême suivi de Razgrad, de Silistra et de Dobrich.
Tous ces ViKs ont une eau qui dépend majoritairement du pompage et de grands systèmes techniques avec des transferts d’eau importants.
summary(ViK$Energie.kwh.m3_21ViK)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.040 0.445 0.585 0.738 0.970 2.340sd(ViK$Energie.kwh.m3_21ViK)## [1] 0.49sd(ViK$Energie.kwh.m3_21ViK)/mean(ViK$Energie.kwh.m3_21ViK)*100## [1] 67ggplot(data=ViK, aes(x = reorder(REG_FR,Energie.kwh.m3_21ViK), y= Energie.kwh.m3_21ViK)) +
geom_bar(stat = "identity", width=0.5, color="black", fill="white")+
geom_text(aes(label=Energie.kwh.m3_21ViK), vjust=-0.3, size=2)+
theme(axis.text.x = element_text(face="bold", color="black",size=7, angle=45),
axis.text.y = element_text(face="bold", color="black",size=7, angle=90))+
labs(x="ViK", y="KwH/m3", title="Consommation énergétique des VIKs en KwH/m3 en 2021")+theme(plot.title=element_text(face="bold", color="black", size=10))+theme(axis.title.x = element_text(face="plain", color="black", size=9))+theme(axis.title.y = element_text(face="plain", color="black", size=9))
ggplot(data=ViK, aes(x=Energie.kwh.m3_21ViK))+
geom_histogram(color="black", fill="white", binwidth = 0.5)+
geom_vline(aes(xintercept=mean(Energie.kwh.m3_21ViK)),color="blue", linetype="dashed", size=1)+
geom_vline(aes(xintercept=median(Energie.kwh.m3_21ViK)),color="red", linetype="dashed", size=1)+
theme(plot.title=element_text(face="bold", color="black", size=10))+
theme(axis.title.x = element_text(face="plain", color="black", size=9))+
theme(axis.title.y = element_text(face="plain", color="black", size=9))+
labs(x="KwH/m3", y="effectif", title="Dépenses énergétiques des ViKs en 2021")+
theme(plot.title=element_text(face="bold", color="black", size=10))
Prix et consommations énergétiques
cor.test(ViK$Energie.kwh.m3_21ViK, ViK$PRIX_ViK_23,method= c("spearman"))## Warning in cor.test.default(ViK$Energie.kwh.m3_21ViK, ViK$PRIX_ViK_23, method =
## c("spearman")): Impossible de calculer la p-value exacte avec des ex-aequos##
## Spearman's rank correlation rho
##
## data: ViK$Energie.kwh.m3_21ViK and ViK$PRIX_ViK_23
## S = 1533, p-value = 0.001
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
## rho
## 0.58ggplot(data=ViK, aes(x = Energie.kwh.m3_21ViK, y= PRIX_ViK_23))+
geom_point()+
geom_text(aes(label=REG_FR), vjust=-0.3, size=2)+
labs(x="consommation énergétique en KwH m3", y="Prix de l'eau en 2023", title="Relation entre les consommations énergétiques et le prix de l'eau")+theme(plot.title=element_text(face="bold", color="black", size=10))+geom_vline(aes(xintercept=median(Energie.kwh.m3_21ViK)),color="red", linetype="dashed", size=1)+geom_hline(aes(yintercept=median(PRIX_ViK_23)))
Plusieurs profils se détachent avant même la réalisation d’une analyse à composante principale, des caractéristiques générales des ViKs.
- Le groupe des ViKs dans des territoires ruraux et agricoles dont, l’eau est issue majoritairement du pompage et a un prix plus élevé qu’ailleurs
- Groupe de régions plus riches, urbaines, où les ViKs génèrent des revenus (et des revenus par hab.) plus élevés - en premier lieu Sofia.
Les prix de l’eau y sont variables : le revenu est davantage dépendant du nombre d’habitant et de leur consommation
- Groupe de régions plutôt industrialisées (ou industrialisées par le passé), urbaines ou rurales, qui ont un accès relativement facile à l’eau, où les ViKs ne sont pas riches, mais peuvent pratiquer des prix plutôt bas (Sliven, Gabrovo)
Il convient maintenant de regarder si ces caractéristiques générales influencent, ou non, la qualité des services d’eau et d’approvisionnement ?
1.B.Performances des ViKs
1.B.a. Pertes en eau
Les pertes en eau vont de 40% à 85% de l’eau prélevé. C’est, dans tout les cas, très supérieur à la moyenne européenne. L’écart type est de 12 et le coefficient de variation de 19%. La moyenne est de 61% et la médiane de 58%.
La matrice de corrélation montre que les pertes en eau tendent à être plus faibles dans les régions urbaines.
Le pourcentage de pertes en eau sur le réseau est corrélé au nombre de lits d’hôpitaux par habitant (0,41) ; au PIB par habitant (0,3) et au pourcentage de population ayant fait des études supérieures (0,3).
Il s’agit essentiellement de l’influence de Sofia, qui a le plus faible niveau de pertes (40%), suivie par Ruse, Smoljan (45%) et Vidin (47%).
En effet, sans Sofia, l’indicateur du PIB et celui des études supérieurs ne sont plus corrélés au pourcentage de pertes en eau.
Il n’y a donc pas de réel lien entre le fait d’être un territoire riche et urbain et le pourcentage de pertes en eau.
Inversement, les pertes en eau sont corrélées au pourcentage d’agriculture dans la VA régionale.
La corrélation est de 0,38, avec une p.value de 0,05 avec Sofia).
En revanche, si on enlève Sofia, la corrélation n’est plus que de 0,31 et la p.value de 0,1 montre que cette corrélation est plutôt liée au hasard.
Cependant, on peut expliquer cette très légère corrélation par le fait que plusieurs régions qui ont un haut pourcentage de pertes en eau sont les régions qui sont situées dans le Nord Est soit à la fois des régions agricoles et des régions où les systèmes techniques d’eau et d’assainissement sont de très grande taille et où l’eau vient de loin.
Enfin, les pertes en eau sont, on l’a déjà vu, positivement corrélées au niveau de prélèvements en eau par habitant (0,75 avec Sofia / 0,77 sans Sofia, coefficients de corrélation \< 0,05).
En effet, ceci semble logique: plus on perd d’eau, plus il faut prélever pour compenser les pertes en comparaison du nombre d’habitant·es qui consomment.
La corrélation positive avec la population âgée (0,35 avec Sofia) est faible et n’est pas significative (p.value \> 0,05) ; de plus elle n’existe plus si l’on retire la capitale.
summary(ViK$Perteseau_Pourc_VIK21)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 40.0 55.9 58.3 61.4 68.8 84.6sd(ViK$Perteseau_Pourc_VIK21)## [1] 12sd(ViK$Perteseau_Pourc_VIK21)/mean(ViK$Perteseau_Pourc_VIK21)*100## [1] 19ggplot(data=ViK, aes(x = reorder(REG_FR,Perteseau_Pourc_VIK21), y= Perteseau_Pourc_VIK21)) +
geom_bar(stat = "identity", width=0.5, color="black", fill="white")+
geom_text(aes(label=Perteseau_Pourc_VIK21), vjust=-0.3, size=2)+
theme(axis.text.x = element_text(face="bold", color="black",size=7, angle=45),
axis.text.y = element_text(face="bold", color="black",size=7, angle=90))+
labs(x="ViK", y="%", title="Pertes en eau en 2021")+theme(plot.title=element_text(face="bold", color="black", size=10))+theme(axis.title.x = element_text(face="plain", color="black", size=9))+theme(axis.title.y = element_text(face="plain", color="black", size=9))
1.B.b. Discontinuité de service
L’indicateur de dis/continuité de service prend en compte le nombre et la durée des ruptures de service d’eau dans les régions pour l’année 2021.
Il est peu fiable pour un certain nombre de régions, selon le KEVR : Vraca, Plovdiv, Montana, Gabrovo, Sliven, Burgas (6/28).
Son minimum est de 0,06 à Pazardzhik et son maximum de 15,3 à Montana (suivi de Gabrovo à 9), avec une médiane à 1 et une moyenne à 2,5.
L’écart type est de 3,5, soit légèrement supérieur à la moyenne, avec un coefficient de variation de 139%.
Les écarts entre les oblast apparaissent donc importants.
Dans les matrices de corrélation, ce ratio est négativement corrélé à la part de population ne se déclarant pas bulgare (-0,40 avec Sofia ; -0,41 sans Sofia avec des p.value \<0,05) et à la population rurale (-0,39 avec Sofia ; -0,42 sans Sofia avec des p.value \<0,05).
Ces faibles corrélations sont difficiles à interpréter en tant que telles.
Plus la population se déclarerait non bulgare et plus elle vivrait en zone rurale, moins le niveau de discontinuité de service serait important.
Le ratio de continuité est positivement corrélé au nombre de pannes pour 100km de réseau par an ; ceci est compréhensible, puisque les interruptions de services sont généralement liées à des pannes et leur réparation ( coefficient de corrélation de spearman = 0.5 et p.value \<0,05 avec ou sans Sofia).
summary(ViK$Continuite_21_VIK)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.060 0.278 0.960 2.524 2.945 15.282sd(ViK$Continuite_21_VIK)## [1] 3.5sd(ViK$Continuite_21_VIK)/mean(ViK$Continuite_21_VIK)*100## [1] 1391.B.c. Pannes (avaries) pour 100km de canalisations
Le minimum de pannes pour 100km de réseau est de 55, à Burgas puis à Targovishte.
Le maximum à Pernik et à Jambol avec plus de 140 avaries pour 100km de réseau en 2021.
La moyenne est de 84 pannes par an pour 100km de réseau avec des données peu dispersées autour (écart type de 24, coefficient de variation : 29%) .
La médiane est de 78 avaries pour 100km.
Les avaries sont, on l’a vu, positivement et significativement corrélées à la discontinuité de service.
Elles sont négativement et faiblement corrélées à : la part de population âgée sur un territoire (-0,37, p.value \<0,05). Cette corrélation n’est cependant plus significative lorsque l’on retire Sofia. Les avaries sont aussi corrélées (avec Sofia), mais sans significativité, au prix de l’eau (-0,34, mais cela n’est pas significatif : p.value \> 0.05) , au revenu par habitant des VIKs (-0.33, mais cela n’est pas significatif : p.value \> 0.05).
summary(ViK$Avaries100Km_VIK_21)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 55.3 67.7 78.1 84.0 93.7 147.2sd(ViK$Avaries100Km_VIK_21)## [1] 24sd(ViK$Avaries100Km_VIK_21)/mean(ViK$Avaries100Km_VIK_21)*100## [1] 29Au final, si l’on prend les trois indicateurs de performance retenus, il apparaît très difficile d’interpréter les résultats.
On ne peut ni les lier à des caractéristiques des ViKs (population desservie, revenu, km de réseau) ni à des caractéristiques territoriales très claires.
Les ruptures de continuité de service tendent à être plus présentes en milieu urbain.
Cela pourrait être lié au fait que “l*es pertes linéaires (volume perdu divisé par la longueur du réseau) dépendent du nombre d’abonnés par kilomètre de réseau (densité linéaire) car des fuites existent préférentiellement au niveau des branchements. Ainsi, en moyenne, un service urbain dense aura plus de pertes linéaires qu’un service rural (IFEN, 2001). Mais sur le terrain, pour une même densité linéaire, la variabilité autour de la valeur moyenne reste très élevée*.” (Guérin Schneider et Nakhla, 2003).
Cependant, le nombre de pannes pour 100 km de réseau n’est pas corrélé avec le caractère urbain ou rural des régions.
De plus, les pertes en eau (%) sont plutôt lié au caractère agricole des régions sans doute du fait des très longs réseaux qui desservent les zones du nord est.\
On ne peut donc pas dire que la performance des ViKs est liée à des caractéristiques territoriales ou à des caractéristiques des ViKs.
**Il n’y a donc pas de conséquences évidentes en terme de qualité du service des inégalités grandissantes entre les compagnies des eaux et les territoires aujourd’hui.**
Si rien n’est automatique et évident, il est cependant possible d’établir des profils des compagnies des eaux sur leurs territoires.
Certains territoires cumulent-t-ils difficultés régionales , difficultés économiques des ViKs et problème de qualité du service.
C’est l’objectif de l’analyse multivariée ci-dessous.
2. Analyse multivariée : l’ACP
Objectif : dresser des profils territorialisés des ViKs et de la qualité des services d’eau
2.A. Préparation de la base de données
Il faut retirer de la BD ViK les variables qualitatives (colonnes 1,2,3,17) ainsi qu’un certain nombre de variables que j’ai choisi de ne pas utiliser pour l’ACP.
Variables retenues : pourcentage de l’évolution de la population entre 2001 et 2021 (8), part de l’agriculture dans la Valeur ajouté régionale en 2023 (10), PIB/hab (12), pourcentage de population vivant en zone rurale (13), revenu des ViKs par habitant en 2023 (25), prix de l’eau en 2023 (27), consommations énergétiques (33), consommation d’eau par habitant en L/J/Hab (36). Dans les mesures où les indicateurs de performances ne font pas l’objet de corrélations claires je les représenterai après sur une carte
ACP<-ViK[, c(8, 10,12, 14,18, 25, 27, 33, 36)]
ACP## PEVO_POP_01.21 AgriVA_23 PIBhab_23 TxRisq_24
## Blagoevgrad -11.4 5.534 16449 26
## Montana -31.6 11.619 15127 24
## Pazardzhik -17.8 5.953 16057 32
## Pernik -23.0 3.486 15608 37
## Pleven -30.5 5.763 17703 32
## Plovdiv -11.2 2.427 22209 35
## Razgrad -29.7 11.602 17081 48
## Ruse -27.5 5.187 21678 23
## Shumen -25.5 10.169 17129 19
## Silistra -30.6 16.487 13114 27
## Sliven -20.5 6.727 13493 40
## Burgas -10.0 2.801 22169 30
## Smoljan -30.6 5.659 17016 16
## Sofia region -14.2 3.293 24251 23
## Stara Zagora -19.6 2.083 36470 37
## Targovishte -30.3 10.452 18802 33
## Varna -6.3 1.904 26946 29
## Veliko Tarnovo -28.8 5.545 20187 34
## Vidin -41.1 11.928 15071 37
## Vraca -31.6 4.561 27113 30
## Jambol -28.9 8.961 17865 27
## Dobrich -29.6 12.360 17075 27
## Gabrovo -31.1 3.843 24297 25
## Sofia capitale 8.1 0.092 61833 25
## Haskovo -23.1 7.257 13012 25
## Kardzhali -13.6 9.823 15593 30
## Kjustendil -33.2 6.365 15030 25
## Lovech -31.7 6.698 18024 24
## EVO_REVcon_VIK_31.21 REV_POP_VIK_21 PRIX_ViK_23
## Blagoevgrad 9118111 93 2.2
## Montana 2664439 111 3.4
## Pazardzhik 2501905 60 1.4
## Pernik 2036535 100 2.0
## Pleven 4835506 114 3.7
## Plovdiv 13395872 97 2.5
## Razgrad 820830 115 4.5
## Ruse 2275191 129 3.7
## Shumen -1016848 119 2.9
## Silistra 1855493 110 4.4
## Sliven 2794877 84 2.0
## Burgas 33216685 227 3.5
## Smoljan 1987513 83 3.9
## Sofia region 2117623 71 2.1
## Stara Zagora 2362934 119 3.8
## Targovishte -589925 90 3.8
## Varna 17642787 148 4.2
## Veliko Tarnovo 6871973 107 3.1
## Vidin 521055 85 2.9
## Vraca 2654986 120 3.3
## Jambol 1511370 99 3.7
## Dobrich 3400699 136 2.5
## Gabrovo 3814481 110 2.9
## Sofia capitale 45179858 148 2.4
## Haskovo 3278131 82 3.6
## Kardzhali 3016240 77 3.2
## Kjustendil 2015503 59 1.7
## Lovech 269743 119 3.4
## Energie.kwh.m3_21ViK Conso_eau_hab_23
## Blagoevgrad 0.06 125
## Montana 0.50 90
## Pazardzhik 0.38 115
## Pernik 0.29 111
## Pleven 0.81 100
## Plovdiv 0.60 102
## Razgrad 1.53 88
## Ruse 1.08 109
## Shumen 1.03 88
## Silistra 1.48 96
## Sliven 0.52 75
## Burgas 0.40 115
## Smoljan 0.59 84
## Sofia region 0.56 98
## Stara Zagora 1.08 95
## Targovishte 0.95 74
## Varna 0.46 110
## Veliko Tarnovo 0.56 105
## Vidin 0.78 92
## Vraca 0.58 103
## Jambol 0.67 91
## Dobrich 1.33 96
## Gabrovo 0.31 105
## Sofia capitale 0.04 126
## Haskovo 2.34 83
## Kardzhali 0.53 81
## Kjustendil 0.38 120
## Lovech 0.82 99ACPnoSOF <- ViKnoSOF [, c(8, 10,12, 14,18, 25, 27, 33, 36)]
ACPnoSOF## PEVO_POP_01.21 AgriVA_23 PIBhab_23 TxRisq_24
## Blagoevgrad -11.4 5.5 16449 26
## Montana -31.6 11.6 15127 24
## Pazardzhik -17.8 6.0 16057 32
## Pernik -23.0 3.5 15608 37
## Pleven -30.5 5.8 17703 32
## Plovdiv -11.2 2.4 22209 35
## Razgrad -29.7 11.6 17081 48
## Ruse -27.5 5.2 21678 23
## Shumen -25.5 10.2 17129 19
## Silistra -30.6 16.5 13114 27
## Sliven -20.5 6.7 13493 40
## Burgas -10.0 2.8 22169 30
## Smoljan -30.6 5.7 17016 16
## Sofia region -14.2 3.3 24251 23
## Stara Zagora -19.6 2.1 36470 37
## Targovishte -30.3 10.5 18802 33
## Varna -6.3 1.9 26946 29
## Veliko Tarnovo -28.8 5.5 20187 34
## Vidin -41.1 11.9 15071 37
## Vraca -31.6 4.6 27113 30
## Jambol -28.9 9.0 17865 27
## Dobrich -29.6 12.4 17075 27
## Gabrovo -31.1 3.8 24297 25
## Haskovo -23.1 7.3 13012 25
## Kardzhali -13.6 9.8 15593 30
## Kjustendil -33.2 6.4 15030 25
## Lovech -31.7 6.7 18024 24
## EVO_REVcon_VIK_31.21 REV_POP_VIK_21 PRIX_ViK_23
## Blagoevgrad 9118111 93 2.2
## Montana 2664439 111 3.4
## Pazardzhik 2501905 60 1.4
## Pernik 2036535 100 2.0
## Pleven 4835506 114 3.7
## Plovdiv 13395872 97 2.5
## Razgrad 820830 115 4.5
## Ruse 2275191 129 3.7
## Shumen -1016848 119 2.9
## Silistra 1855493 110 4.4
## Sliven 2794877 84 2.0
## Burgas 33216685 227 3.5
## Smoljan 1987513 83 3.9
## Sofia region 2117623 71 2.1
## Stara Zagora 2362934 119 3.8
## Targovishte -589925 90 3.8
## Varna 17642787 148 4.2
## Veliko Tarnovo 6871973 107 3.1
## Vidin 521055 85 2.9
## Vraca 2654986 120 3.3
## Jambol 1511370 99 3.7
## Dobrich 3400699 136 2.5
## Gabrovo 3814481 110 2.9
## Haskovo 3278131 82 3.6
## Kardzhali 3016240 77 3.2
## Kjustendil 2015503 59 1.7
## Lovech 269743 119 3.4
## Energie.kwh.m3_21ViK Conso_eau_hab_23
## Blagoevgrad 0.06 125
## Montana 0.50 90
## Pazardzhik 0.38 115
## Pernik 0.29 111
## Pleven 0.81 100
## Plovdiv 0.60 102
## Razgrad 1.53 88
## Ruse 1.08 109
## Shumen 1.03 88
## Silistra 1.48 96
## Sliven 0.52 75
## Burgas 0.40 115
## Smoljan 0.59 84
## Sofia region 0.56 98
## Stara Zagora 1.08 95
## Targovishte 0.95 74
## Varna 0.46 110
## Veliko Tarnovo 0.56 105
## Vidin 0.78 92
## Vraca 0.58 103
## Jambol 0.67 91
## Dobrich 1.33 96
## Gabrovo 0.31 105
## Haskovo 2.34 83
## Kardzhali 0.53 81
## Kjustendil 0.38 120
## Lovech 0.82 992.B. Pertinence des variables retenues pour l’ACP
2.B.a. matrice de corrélations
Création d’une matrice de corrélation pour étudier la pertinence de réaliser une ACP à partir de ces variables.
Ici, on calcule les coefficients de corrélation avec la fonction cor.
C’est par défaut le test de pearson qui est utilisé.
On ne peut en effet effectuer le test de Bartlett sur une matrice de corrélation calculée à partir de spearman sous R
MatrixACP<-cor(ACP)
MatrixACP<-as.data.frame(MatrixACP)
MatrixACP## PEVO_POP_01.21 AgriVA_23 PIBhab_23 TxRisq_24
## PEVO_POP_01.21 1.000 -0.61 0.622 -0.0367
## AgriVA_23 -0.610 1.00 -0.601 0.0502
## PIBhab_23 0.622 -0.60 1.000 -0.0661
## TxRisq_24 -0.037 0.05 -0.066 1.0000
## EVO_REVcon_VIK_31.21 0.776 -0.54 0.723 -0.0487
## REV_POP_VIK_21 0.308 -0.23 0.392 -0.0077
## PRIX_ViK_23 -0.232 0.29 -0.017 0.0376
## Energie.kwh.m3_21ViK -0.375 0.51 -0.314 0.0753
## Conso_eau_hab_23 0.427 -0.55 0.425 -0.1503
## EVO_REVcon_VIK_31.21 REV_POP_VIK_21 PRIX_ViK_23
## PEVO_POP_01.21 0.776 0.3077 -0.232
## AgriVA_23 -0.539 -0.2285 0.285
## PIBhab_23 0.723 0.3920 -0.017
## TxRisq_24 -0.049 -0.0077 0.038
## EVO_REVcon_VIK_31.21 1.000 0.6386 -0.064
## REV_POP_VIK_21 0.639 1.0000 0.381
## PRIX_ViK_23 -0.064 0.3811 1.000
## Energie.kwh.m3_21ViK -0.390 -0.0385 0.520
## Conso_eau_hab_23 0.553 0.2930 -0.395
## Energie.kwh.m3_21ViK Conso_eau_hab_23
## PEVO_POP_01.21 -0.375 0.43
## AgriVA_23 0.509 -0.55
## PIBhab_23 -0.314 0.43
## TxRisq_24 0.075 -0.15
## EVO_REVcon_VIK_31.21 -0.390 0.55
## REV_POP_VIK_21 -0.038 0.29
## PRIX_ViK_23 0.520 -0.39
## Energie.kwh.m3_21ViK 1.000 -0.53
## Conso_eau_hab_23 -0.533 1.00write.csv2(MatrixACP, "C:/Users/Raph/Documents/thesebulgarie/STATISTIQUES/analysesR/VIK_INEG/MatriceCorr.csv")2.B.b. Test de Bartlett
L’ACP est une compression de l’information. Elle n’est possible que si les données présentent une certaines redondances. Si les variables sont toutes parfaitement corrélées alors un axe factoriel suffit et restiuera 100% de l’information. Si les variables sont toutes indépendantes, alors il y a autant d’axes que de variables.
Le test de sphéricité de Bartlett permet de voir dans quelle mesure la matrice de corrélation calculée avec nos données diverge de la matrice unité, càd de la matrice théorique sous l’hypothèse nulle (H=0). II faut que l’hypothèse nulle soit démentie par les données pour pouvoir faire un résumé.
Pour mesurer le lien entre les variables, il faut calculer le déterminant [R] de la matrice de corrélation
H0 : R=1
colinéarité parfaite : R=0
Lorsque R \< 0,00001 on considère qu’il y a une forte redondance des données (elles ne recèlent qu’un seul type d’information).
Et si R se rapproche de 1 alors on est proche de H0.
**Ici p.value (R) =0.00000000000000044**\
L’hypothèse nulle est invalidée et il y a une redondance importante entre les données.
cortest.bartlett(MatrixACP, n=28)## $chisq
## [1] 113
##
## $p.value
## [1] 0.00000000083
##
## $df
## [1] 362.B.c. Test KMO solution factorielle
Ce test permet de voir la corrélation nette entre deux variables, c’est à dire en retranchant l’influence des autres variables dans les corrélations obtenues sur la matrice. Si la corrélation nette est inférieure à la corrélation brute alors les liaisons sont déterminées par d’autres variables. Cela veut dire qu’il y a de la redondance et donc la possibilité de faire une réduction efficace de l’information
Si la corrélation brute > corrélation nette ou équivalente alors la relation directe entre les variables est réelle et sera prise en compte par l’ACP. Ces deux variables détermineront donc souvent un axe factoriel à elles seules
L’indice KMO global donne une idée de la redondance générale des données.
Il varie entre 0 et 1.
Si il est proche de 0 : corrélations partielles identiques aux corrélations brutes : variables indépendantes.
Si proche de 1 : excellent résumé de l’information sur les premiers axes factoriels.
**Ici il est de 0,65.**
KMO(MatrixACP)## Kaiser-Meyer-Olkin factor adequacy
## Call: KMO(r = MatrixACP)
## Overall MSA = 0.73
## MSA for each item =
## PEVO_POP_01.21 AgriVA_23 PIBhab_23
## 0.74 0.80 0.81
## TxRisq_24 EVO_REVcon_VIK_31.21 REV_POP_VIK_21
## 0.58 0.68 0.57
## PRIX_ViK_23 Energie.kwh.m3_21ViK Conso_eau_hab_23
## 0.58 0.84 0.822.B.d. Tests sans Sofia
MatrixACPNoSof<-cor(ACPnoSOF)
MatrixACPNoSof<-as.data.frame(MatrixACPNoSof)
MatrixACPNoSof## PEVO_POP_01.21 AgriVA_23 PIBhab_23 TxRisq_24
## PEVO_POP_01.21 1.000 -0.54 0.299 0.042
## AgriVA_23 -0.541 1.00 -0.629 0.010
## PIBhab_23 0.299 -0.63 1.000 0.068
## TxRisq_24 0.042 0.01 0.068 1.000
## EVO_REVcon_VIK_31.21 0.630 -0.46 0.268 0.063
## REV_POP_VIK_21 0.214 -0.16 0.365 0.022
## PRIX_ViK_23 -0.177 0.25 0.212 0.019
## Energie.kwh.m3_21ViK -0.273 0.46 -0.154 0.044
## Conso_eau_hab_23 0.276 -0.48 0.217 -0.114
## EVO_REVcon_VIK_31.21 REV_POP_VIK_21 PRIX_ViK_23
## PEVO_POP_01.21 0.630 0.214 -0.177
## AgriVA_23 -0.460 -0.162 0.251
## PIBhab_23 0.268 0.365 0.212
## TxRisq_24 0.063 0.022 0.019
## EVO_REVcon_VIK_31.21 1.000 0.711 0.077
## REV_POP_VIK_21 0.711 1.000 0.436
## PRIX_ViK_23 0.077 0.436 1.000
## Energie.kwh.m3_21ViK -0.285 0.030 0.503
## Conso_eau_hab_23 0.442 0.226 -0.368
## Energie.kwh.m3_21ViK Conso_eau_hab_23
## PEVO_POP_01.21 -0.273 0.28
## AgriVA_23 0.459 -0.48
## PIBhab_23 -0.154 0.22
## TxRisq_24 0.044 -0.11
## EVO_REVcon_VIK_31.21 -0.285 0.44
## REV_POP_VIK_21 0.030 0.23
## PRIX_ViK_23 0.503 -0.37
## Energie.kwh.m3_21ViK 1.000 -0.48
## Conso_eau_hab_23 -0.482 1.00write.csv2(MatrixACPNoSof, "C:/Users/Raph/Documents/thesebulgarie/STATISTIQUES/analysesR/VIK_INEG/MatriceCorrNoSof.csv")cortest.bartlett(MatrixACPNoSof, n=27)## $chisq
## [1] 98
##
## $p.value
## [1] 0.00000011
##
## $df
## [1] 36KMO(MatrixACPNoSof)## Kaiser-Meyer-Olkin factor adequacy
## Call: KMO(r = MatrixACPNoSof)
## Overall MSA = 0.57
## MSA for each item =
## PEVO_POP_01.21 AgriVA_23 PIBhab_23
## 0.59 0.66 0.43
## TxRisq_24 EVO_REVcon_VIK_31.21 REV_POP_VIK_21
## 0.17 0.50 0.48
## PRIX_ViK_23 Energie.kwh.m3_21ViK Conso_eau_hab_23
## 0.56 0.76 0.78BARTLETT : p.value =0.00000000003\
KMO : 0,64
2.C. Réalisation ACP
2.C.a. Préparation des données
Centrer réduire
ACP_scale<-scale(ACP)Réalisation ACP avec 6 composantes
ACPViK <- PCA(ACP_scale, ncp = 6, graph = TRUE)## Warning: ggrepel: 4 unlabeled data points (too many overlaps). Consider
## increasing max.overlaps
Choix du nombre de dimensions pertinentes (3 : 78% données)
get_eigenvalue(ACPViK)## eigenvalue variance.percent cumulative.variance.percent
## Dim.1 3.982 44.2 44
## Dim.2 1.744 19.4 64
## Dim.3 1.004 11.2 75
## Dim.4 0.686 7.6 82
## Dim.5 0.481 5.3 88
## Dim.6 0.423 4.7 92
## Dim.7 0.344 3.8 96
## Dim.8 0.240 2.7 99
## Dim.9 0.097 1.1 1002.C.b. Explor : analyse de l’ACP
explor(ACPViK)Sans Sofia :
scaleACPnoSOF<-scale(ACPnoSOF)
ACPViKnoSOF <- PCA(scaleACPnoSOF, ncp = 6, graph = FALSE)
explor(ACPViKnoSOF)Réalisation d’un tableau excel en annexe résumant l’analyse des principaux axes. Voir annexes
3. La CAH, classification ascendante hiérarchique
3.A. préparation
Pour cela, on utilise la fonction AGNES du package “Cluster”.
On utilise une métrique euclidienne pour la distance entre les individus et la méthode de ward qui donne une perte d’inertie minimale
CAH <- agnes(ACP_scale, metric="euclidean", method="ward")3.B. Dendogrammes et analyse des groupes de région
On travaille surtout à partir des 2 premiers axes de l’ACP.
À partir du dendogramme ci-dessous, j’identifie 6 classes
Régions urbaines, où le PIB/hab sont élevés, où le revenu des ViK et le revenu/hab sont importants et où les consommations en eau sont fortes
1) Sofia : prix de l’eau plutôt bas et conso énergétique très basses
2) Burgas et Varna : consommations énergétiques plutôt faibles mais prix de l’eau élevé. Tourisme
Régions où les ViKs sont davantage en difficulté économique
a.Au prix de l’eau bas
3) Blagoevgrad, Pazardzhik, Kjustendil, Pernik, Veliko Tarnovo, Gabrovo, Plovdiv, Région de Sofia= ViK s’étant consolidées, prix de l’eau assez faible, plutôt en zone montagneuse. Le revenu par habitant du ViK est assez faible mais les dépenses énergétiques aussi. Plovdiv est un cas limite car au niveau du PIB/hab il se situe plutôt dans une moyenne haute. Les consommations énergétiques sont faibles
b. Au prix de l’eau élevé
4) Agricoles, où la population diminue et a de faibles revenus et où les prix de l’énergie sont très élevés : Haskovo, Razgrad, Silistra
5) Plus urbaines et industrielles, : Stara Zagora, Ruse, Vraca, Pleven, Lovech. Les consommations énergétiques et les prix de l’eau sont importants notamment par rapport au groupe 3 avec qui équivalent au niveau de l’urbanisation.
6) “moyennes” Moins agricoles, où l**a diminution de la population est forte** et où le prix de l’eau m’est pas bas est variable : Sliven, Kardzhali, Shumen, Dobrich, Smoljan, Targovishte, Vidin, Montana, Jambol. Les consommations énergétiques sont plutot importantes ou moyennes. Beaucoup des régions rurales ex. industrielles en crise
CAH_Dist <- dist(ACP_scale, method="euclidian")
CAH2<-hclust(CAH_Dist, method="ward.D2")
plot(CAH2)
Je choisis donc 6 classes.
Selon le graphique suivant, ces 6 classes résument 48% de l’information de la CAH
HeightCAH <- sort(CAH$height, decreasing=TRUE)
relHeightCAH <- HeightCAH/sum(HeightCAH)*100
cumHeightCAH <- cumsum(relHeightCAH)
barplot(relHeightCAH [1:28], names.arg =seq (1,28,1), col = "lightblue4", border= "lightblue4", xlab="Nombre de groupes", ylab="Part de l'inertie totale (%)")
cumHeightCAH## [1] 13 22 29 35 41 46 50 54 58 61 65 68 71 74 77 79 82 84 87
## [20] 89 91 93 95 96 98 99 100Dendogramme 3 D avec 6 classes :
CAH2_dendo <- as.data.frame(ACP_scale)
CAH2_den <- HCPC(CAH2_dendo, nb.clust = 6, graph=FALSE)
plot.HCPC(CAH2_den, choice = "3D.map")
Dendogramme en arbre avec 6 classes
plot.HCPC(CAH2_den, choice = "tree")
3.C. Finalisation de la CAH : les groupes de clusters
#### Création d’une nouvelle variable indiquant les groupes de clusters
clusCAH_ViK <- cutree(CAH, k=6) #on le fait pour 6 classes
ViK_Clusters <- as.data.frame(ViK)
ViK_Clusters$clusCAH_ViK <- factor(clusCAH_ViK, levels = 1:6, labels = paste("CLUS", 1:6))
ViK_Clusters## REG_FR CODE clusCAH RevHab_24 Pop65_24 TxMig_14.24
## Blagoevgrad Blagoevgrad BGR.1_1 CLUS 4 25997 23 -20.0
## Montana Montana BGR.10_1 CLUS 1 26217 28 -28.2
## Pazardzhik Pazardzhik BGR.11_1 CLUS 1 24736 26 -17.9
## Pernik Pernik BGR.12_1 CLUS 1 30570 30 37.5
## Pleven Pleven BGR.13_1 CLUS 1 26815 31 -36.1
## Plovdiv Plovdiv BGR.14_1 CLUS 2 24742 26 49.8
## Razgrad Razgrad BGR.15_1 CLUS 4 21316 27 -23.3
## Ruse Ruse BGR.16_1 CLUS 1 22413 27 -5.2
## Shumen Shumen BGR.17_1 CLUS 4 27406 26 38.3
## Silistra Silistra BGR.18_1 CLUS 4 25220 29 -14.6
## Sliven Sliven BGR.19_1 CLUS 4 23765 23 -44.4
## Burgas Burgas BGR.2_1 CLUS 2 26769 29 69.0
## Smoljan Smoljan BGR.20_1 CLUS 1 29927 26 -87.8
## Sofia region Sofia region BGR.21_1 CLUS 3 18820 24 59.1
## Stara Zagora Stara Zagora BGR.22_1 CLUS 3 24704 26 6.7
## Targovishte Targovishte BGR.23_1 CLUS 4 23476 25 -3.5
## Varna Varna BGR.24_1 CLUS 2 23076 28 54.1
## Veliko Tarnovo Veliko Tarnovo BGR.25_1 CLUS 1 22549 25 -20.4
## Vidin Vidin BGR.26_1 CLUS 4 21539 23 -46.7
## Vraca Vraca BGR.27_1 CLUS 3 23540 23 -48.2
## Jambol Jambol BGR.28_1 CLUS 1 17821 27 -32.2
## Dobrich Dobrich BGR.3_1 CLUS 1 20742 28 -11.0
## Gabrovo Gabrovo BGR.4_1 CLUS 1 28608 31 -10.0
## Sofia capitale Sofia capitale BGR.5_1 CLUS 5 31829 19 41.2
## Haskovo Haskovo BGR.6_1 CLUS 1 20292 27 8.5
## Kardzhali Kardzhali BGR.7_1 CLUS 6 21024 25 188.4
## Kjustendil Kjustendil BGR.8_1 CLUS 1 23491 28 -18.2
## Lovech Lovech BGR.9_1 CLUS 4 17555 22 -19.6
## DensPop_24 PEVO_POP_01.21 NonBG_21 AgriVA_23 IndusVA_23
## Blagoevgrad 44 -11.4 18.9 5.534 25
## Montana 31 -31.6 17.0 11.619 24
## Pazardzhik 50 -17.8 19.6 5.953 39
## Pernik 46 -23.0 5.9 3.486 31
## Pleven 46 -30.5 11.5 5.763 28
## Plovdiv 107 -11.2 19.1 2.427 35
## Razgrad 41 -29.7 62.3 11.602 29
## Ruse 65 -27.5 23.1 5.187 31
## Shumen 44 -25.5 45.9 10.169 32
## Silistra 33 -30.6 47.2 16.487 16
## Sliven 48 -20.5 33.1 6.727 24
## Burgas 50 -10.0 26.3 2.801 23
## Smoljan 28 -30.6 9.8 5.659 36
## Sofia region 32 -14.2 11.8 3.293 54
## Stara Zagora 56 -19.6 19.1 2.083 66
## Targovishte 35 -30.3 52.7 10.452 33
## Varna 115 -6.3 18.4 1.904 26
## Veliko Tarnovo 43 -28.8 13.9 5.545 28
## Vidin 23 -41.1 9.6 11.928 13
## Vraca 40 -31.6 10.0 4.561 60
## Jambol 31 -28.9 17.3 8.961 29
## Dobrich 31 -29.6 27.4 12.360 25
## Gabrovo 46 -31.1 9.1 3.843 44
## Sofia capitale 965 8.1 16.9 0.092 13
## Haskovo 37 -23.1 27.2 7.257 22
## Kardzhali 47 -13.6 73.5 9.823 42
## Kjustendil 35 -33.2 9.0 6.365 27
## Lovech 27 -31.7 11.1 6.698 32
## PIBhab_23 PopRural_24 TxRisq_24 FUSION_ViK
## Blagoevgrad 16449 39.8 26 Petric
## Montana 15127 36.3 24
## Pazardzhik 16057 38.3 32 Panagjurishte, Strelcha
## Pernik 15608 21.7 37
## Pleven 17703 32.5 32 Knezha
## Plovdiv 22209 25.2 35
## Razgrad 17081 52.6 48 Kurbat
## Ruse 21678 22.0 23
## Shumen 17129 38.1 19
## Silistra 13114 55.5 27
## Sliven 13493 35.0 40
## Burgas 22169 23.4 30
## Smoljan 17016 42.5 16
## Sofia region 24251 39.8 23 Botevgrad
## Stara Zagora 36470 28.6 37
## Targovishte 18802 43.9 33
## Varna 26946 15.5 29
## Veliko Tarnovo 20187 29.9 34 Svishtov
## Vidin 15071 36.8 37
## Vraca 27113 40.7 30
## Jambol 17865 30.1 27
## Dobrich 17075 30.2 27
## Gabrovo 24297 20.6 25 Sevlievo
## Sofia capitale 61833 4.9 25
## Haskovo 13012 28.5 25 Dimitrovgrad + Stambolovo
## Kardzhali 15593 59.3 30
## Kjustendil 15030 31.3 25 Dupnica
## Lovech 18024 38.3 24
## REV_VIK_21 EVO_REVcon_VIK_pourc_13.21 EVO_REVcon_VIK_31.21
## Blagoevgrad 24782754 75.4 9118111
## Montana 12392993 33.5 2664439
## Pazardzhik 10761841 37.3 2501905
## Pernik 12066065 24.6 2036535
## Pleven 27415887 25.9 4835506
## Plovdiv 64806330 31.8 13395872
## Razgrad 9907945 10.7 820830
## Ruse 27851183 10.5 2275191
## Shumen 20580338 -5.5 -1016848
## Silistra 12022800 22.0 1855493
## Sliven 15631200 26.4 2794877
## Burgas 92973600 71.6 33216685
## Smoljan 8781623 35.9 1987513
## Sofia region 15990648 18.3 2117623
## Stara Zagora 37340883 8.0 2362934
## Targovishte 9997000 -6.4 -589925
## Varna 69827463 41.9 17642787
## Veliko Tarnovo 25490400 46.0 6871973
## Vidin 7185600 9.3 521055
## Vraca 19227600 19.2 2654986
## Jambol 11711968 17.7 1511370
## Dobrich 23638800 20.2 3400699
## Gabrovo 11879122 60.0 3814481
## Sofia capitale 196243690 36.8 45179858
## Haskovo 19172873 25.0 3278131
## Kardzhali 11143000 46.2 3016240
## Kjustendil 7603100 44.8 2015503
## Lovech 11235860 2.9 269743
## Km_reseau_eau_21 dens_pop_reseau21 Rev.VIK_milliers.KM_21
## Blagoevgrad 1770 150 14.0
## Montana 1824 61 6.8
## Pazardzhik 1150 155 9.4
## Pernik 1481 82 8.1
## Pleven 3255 74 8.4
## Plovdiv 4344 154 14.9
## Razgrad 1428 60 6.9
## Ruse 2629 82 10.6
## Shumen 2367 73 8.7
## Silistra 1861 59 6.5
## Sliven 1857 100 8.4
## Burgas 4304 95 21.6
## Smoljan 1512 70 5.8
## Sofia region 3632 62 4.4
## Stara Zagora 3475 90 10.7
## Targovishte 2212 50 4.5
## Varna 4334 109 16.1
## Veliko Tarnovo 2299 103 11.1
## Vidin 1684 50 4.3
## Vraca 2865 56 6.7
## Jambol 1785 67 6.6
## Dobrich 3404 51 6.9
## Gabrovo 1958 55 6.1
## Sofia capitale 3699 359 53.1
## Haskovo 3193 73 6.0
## Kardzhali 1628 89 6.8
## Kjustendil 1413 91 5.4
## Lovech 1803 52 6.2
## POP_VIK_21 EVO_POP_VIK_13.21 PEVO_POP_VIK_13.21 REV_POP_VIK_21
## Blagoevgrad 265892 55768 26.5 93
## Montana 111151 -6552 -5.6 111
## Pazardzhik 178564 32474 22.2 60
## Pernik 120880 -5771 -4.6 100
## Pleven 240380 -9165 -3.7 114
## Plovdiv 668334 -14693 -2.2 97
## Razgrad 85825 -21508 -20.0 115
## Ruse 215865 -11650 -5.1 129
## Shumen 172367 -27459 -13.7 119
## Silistra 109271 -7925 -6.8 110
## Sliven 185191 -10503 -5.4 84
## Burgas 410331 -4732 -1.1 227
## Smoljan 105421 -3384 -3.1 83
## Sofia region 226385 13767 6.5 71
## Stara Zagora 313396 -50287 -13.8 119
## Targovishte 111597 26083 30.5 90
## Varna 471252 -38558 -7.6 148
## Veliko Tarnovo 237628 26029 12.3 107
## Vidin 84865 -68201 -44.6 85
## Vraca 160829 -32143 -16.7 120
## Jambol 118897 -9535 -7.4 99
## Dobrich 173657 -12522 -6.7 136
## Gabrovo 108481 22157 25.7 110
## Sofia capitale 1328790 19156 1.5 148
## Haskovo 234354 67725 40.6 82
## Kardzhali 145218 9640 7.1 77
## Kjustendil 128673 45072 53.9 59
## Lovech 94650 -10230 -9.8 119
## EVO_REVcon_POP_13.21 PRIX_ViK_23 PEVO_PRIX_cou_13.23
## Blagoevgrad 22.2 2.2 38.5
## Montana 28.0 3.4 42.8
## Pazardzhik 5.7 1.4 5.1
## Pernik 20.0 2.0 25.9
## Pleven 23.0 3.7 53.5
## Plovdiv 21.4 2.5 46.4
## Razgrad 27.5 4.5 50.1
## Ruse 15.7 3.7 49.7
## Shumen 9.0 2.9 22.9
## Silistra 22.2 4.4 54.7
## Sliven 18.2 2.0 14.2
## Burgas 82.2 3.5 48.3
## Smoljan 20.5 3.9 51.0
## Sofia region 6.0 2.1 12.1
## Stara Zagora 20.6 3.8 40.9
## Targovishte -30.3 3.8 29.5
## Varna 44.2 4.2 52.3
## Veliko Tarnovo 21.2 3.1 52.7
## Vidin 35.7 2.9 39.8
## Vraca 30.8 3.3 41.5
## Jambol 18.0 3.7 55.7
## Dobrich 26.1 2.5 5.5
## Gabrovo 20.1 2.9 37.2
## Sofia capitale 32.6 2.4 39.9
## Haskovo -8.8 3.6 31.4
## Kardzhali 17.6 3.2 57.0
## Kjustendil -3.2 1.7 0.0
## Lovech 12.5 3.4 40.8
## EVO_PRIX_cou_13.23 EAUprel_21 EAUprelHab_21
## Blagoevgrad 0.85 27368811 103
## Montana 1.48 19510132 176
## Pazardzhik 0.07 19844284 111
## Pernik 0.52 24727882 205
## Pleven 1.99 27109307 113
## Plovdiv 1.18 80112632 120
## Razgrad 2.25 11059338 129
## Ruse 1.84 18441714 85
## Shumen 0.66 36739201 213
## Silistra 2.42 10280450 94
## Sliven 0.28 39440815 213
## Burgas 1.69 60083418 146
## Smoljan 1.96 7924761 75
## Sofia region 0.26 23889253 106
## Stara Zagora 1.57 31325423 100
## Targovishte 1.11 12009848 108
## Varna 2.18 65314429 139
## Veliko Tarnovo 1.64 27424738 115
## Vidin 1.16 6102371 72
## Vraca 1.36 18876730 117
## Jambol 2.05 20653219 174
## Dobrich 0.14 40835687 235
## Gabrovo 1.08 12564286 116
## Sofia capitale 0.98 142388515 107
## Haskovo 1.14 17299143 74
## Kardzhali 1.83 10237822 70
## Kjustendil 0.00 13299419 103
## Lovech 1.37 9049870 96
## Perteseau_Pourc_VIK21 Energie.kwh.m3_21ViK Continuite_21_VIK
## Blagoevgrad 58 0.06 1.85
## Montana 76 0.50 15.28
## Pazardzhik 61 0.38 0.06
## Pernik 80 0.29 5.43
## Pleven 62 0.81 1.27
## Plovdiv 58 0.60 2.41
## Razgrad 69 1.53 0.16
## Ruse 45 1.08 0.18
## Shumen 81 1.03 0.73
## Silistra 58 1.48 0.20
## Sliven 85 0.52 8.28
## Burgas 54 0.40 0.18
## Smoljan 45 0.59 0.59
## Sofia region 56 0.56 0.29
## Stara Zagora 57 1.08 3.53
## Targovishte 69 0.95 0.24
## Varna 64 0.46 0.78
## Veliko Tarnovo 58 0.56 1.38
## Vidin 47 0.78 4.79
## Vraca 60 0.58 0.12
## Jambol 76 0.67 2.75
## Dobrich 74 1.33 1.01
## Gabrovo 57 0.31 9.21
## Sofia capitale 40 0.04 0.83
## Haskovo 57 2.34 2.33
## Kardzhali 49 0.53 0.46
## Kjustendil 66 0.38 0.91
## Lovech 55 0.82 5.40
## Avaries100Km_VIK_21 Conso_eau_hab_23 PourcEauBudget clusCAH_ViK
## Blagoevgrad 86 125 0.39 CLUS 1
## Montana 76 90 0.43 CLUS 2
## Pazardzhik 79 115 0.23 CLUS 1
## Pernik 143 111 0.27 CLUS 1
## Pleven 71 100 0.51 CLUS 3
## Plovdiv 81 102 0.38 CLUS 1
## Razgrad 65 88 0.68 CLUS 4
## Ruse 63 109 0.66 CLUS 3
## Shumen 68 88 0.33 CLUS 4
## Silistra 60 96 0.62 CLUS 4
## Sliven 84 75 0.23 CLUS 2
## Burgas 55 115 0.55 CLUS 5
## Smoljan 72 84 0.39 CLUS 2
## Sofia region 93 98 0.41 CLUS 1
## Stara Zagora 87 95 0.54 CLUS 3
## Targovishte 55 74 0.43 CLUS 2
## Varna 79 110 0.73 CLUS 5
## Veliko Tarnovo 69 105 0.53 CLUS 3
## Vidin 127 92 0.45 CLUS 2
## Vraca 67 103 0.52 CLUS 3
## Jambol 147 91 0.69 CLUS 2
## Dobrich 78 96 0.42 CLUS 4
## Gabrovo 63 105 0.39 CLUS 3
## Sofia capitale 95 126 0.35 CLUS 6
## Haskovo 114 83 0.54 CLUS 4
## Kardzhali 109 81 0.45 CLUS 2
## Kjustendil 72 120 0.31 CLUS 1
## Lovech 96 99 0.69 CLUS 3Export du nouveau tableau
write.csv2(ViK_Clusters, "C:/Users/Raph/Documents/thesebulgarie/STATISTIQUES/analysesR/VIK_INEG/ViK_Clusters_DEF.csv")Visualisation du poids des variables de l’ACP dans chacun des clusters
Scale_clusCAH_ViK <- as.data.frame(ACP_scale)
Scale_clusCAH_ViK$clusCAH_ViK <- factor(clusCAH_ViK, levels = 1:6, labels = paste("CLUS", 1:6))
clusProfil6cl_ViK <- aggregate(Scale_clusCAH_ViK[,1:9], by=list(Scale_clusCAH_ViK$clusCAH_ViK), mean)
colnames(clusProfil6cl_ViK) [1] <- "GROUPES"
GROUPES=c("Groupe 3", "Groupe 2", "Groupe 4","Groupe 1","Groupe 5","Groupe 6")
clusProfil6cl_ViK$GROUPES<-GROUPES
clusLongViK<-melt(clusProfil6cl_ViK, id.vars="GROUPES")
ggplot(clusLongViK)+geom_bar(aes(x=variable, y=value, fill=GROUPES), stat="identity")+facet_wrap(~GROUPES)+coord_flip()+scale_fill_manual(values=c("#75bc86", "#bcd4cf", "#b9dffd", "#72ccf9", "#0199ff", "#0020ae"))
Cartographie associée : sur QGIS