{"word_count": 282, "character_count": 2809}
Sur http://opendata.stif.info/, on peut avoir accès à des trafics journaliers sur des lignes de bus en région parisienne.&#62; base1=read.csv("/home/charpentier/Téléchargements/validations-sur-le-reseau-de-surface-nombre-de-validations-par-jour-1er-sem (1).csv", sep=";")&#62; base1=base1[-which(base1$LIBELLE_LIGNE%in%c("Inconnu","NON DEFINI")),]&#62; dim(base1)[1] 991842      8Il faut faire un peu attention avec le trafic, qui est ici un facteur (à cause de réponses du type 'moins de 5')&#62; base1$NB_VALD=as.numeric(+ as.character(base1$NB_VALD))On peut considerer une ligne de bus particulière, par  exemple entre la porte de St-Cloud et la Bibliothèque Nationale,&#62; id="PARIS-16 (Porte de St-Cloud)  - PARIS-13 (Bibliothèque F. Mitterrand)"ce qui pourrait être la ligne 62,On peut spécifier ensuite la catégorie du titre de transport, comme une carte navigo, ou une carte imagine'r (étudiants et scolaires)&#62; B_navigo=base1[(base1$LIBELLE_LIGNE==id)&#38;(base1$CATEGORIE_TITRE=="NAVIGO"),]&#62; B_ir=base1[(base1$LIBELLE_LIGNE==id)&#38;(base1$CATEGORIE_TITRE=="IMAGINE R"),]Comme on va visualiser une série temporelle, on doit trier en fonction de la date&#62; B_navigo$date=as.Date(B_navigo$JOUR,"%d/%m/%Y")&#62; B_navigo=B_navigo[order(B_navigo$date),]Sur les 6 premiers mois de 2015, on observe la série temporelle suivante (avec une carte navigo)&#62; plot(B_navigo$date,B_navigo$NB_VALD,type="l")Si on regarde en fonction du jour de la semaine,  du dimanche au samedi,&#62; jour=as.POSIXlt(B_navigo$date)$wday&#62; x=B_navigo$NB_VALD&#62; plot(jour,x,axes=FALSE,ylim=c(-1000,max(x)))&#62; Z_m=tapply(as.numeric(x),as.factor(jour),mean)&#62; lines(as.numeric(names(Z_m)),Z_m,col="red")&#62; axis(2)&#62; text(0:6,rep(-1000,7),+ c("D","L","M","M","J","V","S"))on obtientavec une tendance qui ne présente pas trop de surprise : beaucoup de monde les jours de semaine, moins le samedi, et encore moins le dimanche. Pour la carte étudiante, on a une tendance un peu différente&#62; B_ir$date=as.Date(B_ir$JOUR,"%d/%m/%Y")&#62; B_ir=B_ir[order(B_ir$date),]&#62; plot(B_ir$date,B_ir$NB_VALD,type="l")avec sur les premiers mois de 2016,On retrouve les deux périodes de vacances, et une étrange tendance: une hausse jusqu'en avril, puis une décroissance (sur les jours de semaine). Si on regarde en fonction du jour de la semaine&#62; jour=as.POSIXlt(B_ir$date)$wday&#62; x=B_ir$NB_VALD&#62; plot(jour,x,axes=FALSE,ylim=c(-1000,max(x)))&#62; Z_m=tapply(as.numeric(x),as.factor(jour),mean)&#62; lines(as.numeric(names(Z_m)),Z_m,col="red")&#62; axis(2)&#62; text(0:6,rep(-1000,7),+ c("D","L","M","M","J","V","S"))On a ici, là encore, beaucoup moins de trafic le dimanche que les jours de la semaine,On verra - si on a le temps - pour utiliser ces données dans le cadre du cours de séries temporelles.