########################################
#TELECOM Saint-Etienne - année 2010/2011
#L. Carraro - O. Roustant
########################################

#fonction de répartition de la loi uniforme sur {1,...,6}
x <- 1:6
Fn <- ecdf(x)
plot(Fn, verticals=FALSE, do.points=TRUE, xlim=c(-2,10), col.hor="blue", col.ver="blue", col.points="blue",lwd=2, main="Fonction de répartition de la loi uniforme sur {1,...,6}", xlab="", ylab="")

#fonction de répartition de la loi uniforme sur [0,1]
m <- -0.5
M <- 1.5
plot(c(0,1), c(0,1), xlim=c(m,M), col="blue", type="l", lwd=2, main="Fonction de répartition de la loi uniforme sur [0,1]", xlab="", ylab="")
lines(c(m,0), c(0,0), col="blue",lwd=2)
lines(c(1,M), c(1,1), col="blue",lwd=2)


#bien se placer dans le répertoire des données avant de les charger
oxyde<-read.table(file="oxyde.txt",header=TRUE,sep=";",dec=".",skip=1)
plot(ecdf(oxyde$thickness),xlab="épaisseur en nm",main="fonction de répartition empirique")
op<-par(mfrow=c(1,3))
plot(ecdf(oxyde$thickness[oxyde$site==1]),xlab="épaisseur en nm",main=c("fonction de répartition empirique","plot 1"))
plot(ecdf(oxyde$thickness[oxyde$site==2]),xlab="épaisseur en nm",main=c("fonction de répartition empirique","plot 2"))
plot(ecdf(oxyde$thickness[oxyde$site==3]),xlab="épaisseur en nm",main=c("fonction de répartition empirique","plot 3"))
par(op)

#loi binomiale
#diagramme en batons
x<-0:20
plot(x,dbinom(x,size=20,prob=0.15),type="h",ylab="probabilité",main="diagramme en bâtons de la loi binomiale B(20,0.15)")
points(x,dbinom(x,size=20,prob=0.15),pch="o")

#fonction de repartition
f<-stepfun(0:20,c(0,pbinom(0:20,size=20,prob=0.15)))
plot(f,verticals=F,lwd=2,pch=19,col.hor="blue",do.points=TRUE, col.points="blue", ylab="F(x)",main="fonction de répartition de la loi binomiale B(20,0.15)")

#loi uniforme sur [0,1]
#fonction de répartition et densité
op<-par(mfrow=c(2,1))
m <- -0.5
M <- 1.5
plot(c(0,1), c(0,1), xlim=c(m,M), col="blue", type="l", lwd=2, main="Fonction de répartition de la loi uniforme sur [0,1]", xlab="", ylab="")
lines(c(m,0), c(0,0), col="blue",lwd=2)
lines(c(1,M), c(1,1), col="blue",lwd=2)
plot(c(0,1), c(1,1), xlim=c(m,M), ylim=c(0,1),col="blue", type="l", lwd=2, main="Densité de la loi uniforme sur [0,1]", xlab="", ylab="")
lines(c(m,0), c(0,0), col="blue",lwd=2)
lines(c(1,M), c(0,0), col="blue",lwd=2)
par(op)