# traitement des données issues du fichier culture_bio1.txt
# Orsay, 28/11/2008
# Objet du script : réaliser un histogramme des données réparties en classes
# Recupérer les données
cult<-read.table("./culture_bio1.txt",h=TRUE,sep=";")
attach(cult)
names(cult)
cult
summary(cult)
# définir l'intervalle de classe
# (remarque : l'unité considérée est : qt/ha)
interclasse=5
# limites en X du graphique (la variable s'appelle ici Rendement)
min=80
max=150
limclasses=seq(min,max,interclasse)
limclasses
# dessin de l'histogramme
hist(Rendement,breaks=limclasses,main="Nouvelle variété céréalière",ylim=c(0,15))
# dessin de la loi normale pouvant approcher la distribution 
# utilisant la fonction dnorm(etendue, moy, sigma) de R
# A/ choix de moy et sigma : ceux de la ditribution dans un premier temps
moy=mean(Rendement)
moy
sigma=sd(Rendement)
sigma
etendue=min:max
dist_theo=dnorm(etendue,moy,sigma)
# B/ dessin de la normale
par(new=TRUE)
plot(dist_theo~etendue,col="red",type="l")