###########################
# Nos premiers pas avec R #
# Orsay le 12/10/2007     #
###########################
#
####################################
# 1/ importation des données       #
####################################
### Attention l'importation du fichier demande une gestion parfaite des répertoire et du raccourci R #####
#
croissance<-read.table("./datafiles/bacteria.txt",h=TRUE)
attach(croissance)
#
#///////////////////////////////////
#Contrôle des variables :
#summary(croissance)
#temps
#DO
#summary(DO)
#min(DO)
#length(DO)
#sum(DO)/length(DO)
#mean(DO)
#sort(DO)
#croissance[DO=="4.25",]
#croissance[DO==4.25 & temps==210,]
########################################
# 2/ Graphe ds données expérimentales  #
########################################
#
plot(x=temps,y=DO,xlab="temps en min", ylab="Abs-600nm",main="croissance bactérienne-Graph#1-2007/10/12",col="red",las=1,xlim=c(0,250),ylim=c(0,5),pch=20)
################################################
# (zone exponentielle : une première approche) #
################################################
#phase_expo=croissance[temps<=170,]
#phase_expo
#attach(phase_expo)
#time=temps
#abs=DO
#par(new=TRUE)
#plot(x=time,y=abs,col="darkgreen",las=1,xlim=c(0,250),ylim=c(0,5),pch=2)
#
####################################################################
## 3/ Génération d’une courbe de tendance exponentielle théorique  #
####################################################################
# on commence par générer suffisamment d’abscisses
t=seq(0,185,.1)
# on génère ensuite toutes les valeurs théoriques à partir de l’équation récupérée 
# dans EXCEL
DO_lissage_excel=0.0055*exp(0.0354*t)
# Et hop, voici la courbe tracée avec la fonction lines
# rq le paramètre lwd fixe la taille des lignes
lines(t,DO_lissage_excel,col="lightblue",lwd=2)
#
####################
# FIN DU PROGRAMME #
####################