# Cours02_16  Densité gaussienne
#Lecture des données de l'exemple des prix du caoutchouc dans rubber.csv
rubber<-read.csv2(file.choose())

rubber.cc<-rubber[complete.cases(rubber),]

hist(rubber$Return)
x<-rubber$Prices
length(x)
hx=hist(x)
names(hx)
hx$counts
hx$breaks
y<-rubber.cc$Return
length(y)
hy=hist(y)
names(hy)
hy$counts
min(y)
max(y)
bornes=seq(-0.18,+0.20,0.01)
bornes
hist(y,breaks=bornes,col="red")
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,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=10000
z=rnorm(N)
plot(z,type="h")
hist(z)
hist(z,freq=F)
gaussienne=function(x){return(dnorm(x,0,1))}
plot(gaussienne,min(z),max(z),col="red",add=TRUE) # tracé du graphe de la fonction entre  xmin et xmax