#..
CartoGraphie
# virer tous les objets (mais pas les packahges)
rm(list = ls())
# sans succès: use Ctrl+Shift+F10
install.packages("geojsonio")
install.packages("sp")
library(sp)
library(geojsonio)
# fond france
france <- geojson_read("https://raw.githubusercontent.com/gregoiredavid/france-geojson/master/communes.geojson", what = "sp")
# Since it is a bit to much data, I select only a subset of it:
loiret <- france[ substr(france@data$code,1,2) %in% c("45") , ]
#carte blanche
plot(loiret)
# we recommend installing the entire tidyverse
# modeling set, which includes broom:
install.packages("tidymodels")
library(broom)
loiret_fortified <- tidy(loiret, region = "code")
install.packages("mapproj")
library(mapproj)
library(ggplot2)
ggplot() +
geom_polygon(data = loiret_fortified, aes( x = long, y = lat, group = group), fill="white", color="grey") +
theme_void() +
coord_map()
# we recommend installing the entire tidyverse
# modeling set, which includes broom:
install.packages("tidymodels")
#package ‘tidymodel’ is not available (for R version 4.0.2)
# si broom erreur backport
library(broom)
# GAFFE: i faut refuser la compilation
install.packages("tidyverse")
library(tidyr)
# loiret
loiret <- france[ substr(france@data$code,1,2) %in% c("45") , ]
head(loiret)
##je ne sais la signification du "fortified", conservé de l'exemple
loiret_fortified <- tidy(loiret, region = "code")
head (loiret_fortified)
#lecture des données
# dans la vraie vie, proviendra probablement d'une table excel en .csv
# pour mémoire:
#id,nb_equip
#45270,10
#45234,1000
#lecture depuis table .csv
# d1<-read.csv("loiret.csv",colClasses=c("character","numeric"))
#pour obtenir un exemple "stand-alone" qui ne doive rien à personne, ici les données sont
##forgées dans le programe:
valeurs <- data.frame(id = c("45001","45010","45270", "45234"),
nb_equip = c(6,5,10,1))
head(valeurs)
#creation de la table fusion de u fond et des valeurs:
install.packages("dplyr")
library(dplyr)
fusion<-left_join(loiret_fortified,valeurs)
head(fusion)
#Ok:
ggplot() + geom_polygon(data = fusion, aes(fill = nb_equip, x = long, y = lat, group = group)) + theme_void()+ coord_map()