# ouvrir la table #
V<-read.delim("displacement.txt",header=TRUE)
# définir les nom de colonne (ici P1, P2, etc.) et mettre les données en numérique (et NA)#
fix(V)
# voir la répartition des données et le range #
# ainsi que la moyenne, et la déviation standard pour les différentes variable (ici période) #
# et faire le test de shapiro-wilk pour tester la normalité de la distribution (accepté si p-value > 0.1)#
summary(V)
mean(V, na.rm=TRUE)
sd (V, na.rm=TRUE)
shapiro.test(V$P1)

# faire le boxplot pour les 4 périodes et rajouter une ligne de référence horizontale #
boxplot(V$P1,V$P2,V$P3,V$P4,V$P5,ylim=c(0,0.2),main="Vitesse 2006-... (points in flat area)",lab=c(1,6,2),names=c("2006-2007","2007-2008","2008-2009","2009-2010","2010-2011"), notch=TRUE, ylab="déplacements [m/mois]", border="red")
#  fonction de distribution cumulée pour chaque année sur un même graphe#
#  mais avant il faut loader les package "graphics" et "gplot"  #
#  les 2 premier chiffre dans legend sont le point d'ancrage dans le graphe # 
library(gdata,gtools,graphics,bitops,KernSmooth,caTools,gplots)
plot.ecdf (V$P1,xlim=c(0,0.2),main="Distribution de déplacements annuels 2006-2011",col.01line="grey", col="lightblue")
plot.ecdf (V$P2, xlim=c(0,0.2), col.01line="grey", add=TRUE, col="blue", abline(h=0.5))  
legend (0.16,1,legend=c("2006-2007","2007-2008"), col=c("lightblue","blue"), lty=1,lwd=2, bty="n",bg="white",cex=0.8)

#  fonction de distribution pour chaque année sur 5 graphes#
par(mfrow=c(2,3))
plot.ecdf (V$P1,xlim=c(0,0.16), main="2006-2007",col.01line="grey",col="lightblue")
plot.ecdf (V$P2, xlim=c(0,0.16), main="2007-2008",col.01line="grey", col="blue",abline(h=0.5))
legend (0.16,1,legend=c("2006-2007","2007-2008"), col=c("lightblue","blue"), lty=1,lwd=2, bty="n",bg="white",cex=0.8)

# La même chose, mais en non cumulé#
plot.density(V$P1,xlim=c(0,2),main="2006-2007",col="lightblue")
plot(density(V$P2, xlim=c(0,2), col.01line="grey", add=TRUE, col="blue"))  
legend ("topright",legend=c("2006-2007","2007-2008"), col=c("lightblue","blue"), lty=1,lwd=2, bty="n",bg="white",cex=0.8)