# Cours01_16  Histogrammes
#Lecture des données de l'exemple des prix du caoutchouc dans rubber.csv
rubber<-read.csv2(file.choose())
head(rubber)
summary(rubber)
rubber.cc<-rubber[complete.cases(rubber),]
head(rubber.cc)
#les prix sont en THB/kg
plot(rubber$Prices,type="l")
plot(rubber$Return,type="l")
hist(rubber$Prices)
hist(rubber$Return)
x<-rubber$Prices
length(x)
y<-rubber.cc$Return
length(y)
min(y)
max(y)
bornes=seq(-0.18,+0.20,0.01)
bornes
hist(y,breaks=bornes,col="blue",freq=F)

mu=mean(y);mu #notez qu'on peut remplacer  <-  par  = 
sigma=sd(y);sigma  # sd vient de "standard deviation"
magaussienne=function(x){return(dnorm(x,mu,sigma))} 
# par(new=TRUE) #passe magique pour ajouter un dessin à la fenêtre courante
# mais sans ajustage des coordonnées.
plot(magaussienne,min(y),max(y),col="red",add=TRUE) # tracé du graphe de la fonction entre  xmin et xmax
magaussienne(mu) # maximum de la fonction "magaussienne"



hist(y,breaks=bornes,add=TRUE,col="blue",freq=T)

length(which(y==0))
length(which(y==0))/length(y)
hist(y,breaks=bornes,col="green",freq=F,right=F)

# N=1000
# z=rnorm(N)
# plot(z,type="h")