123456890123456890123456890123456890123456890123456890123456890123456890
#Initialisation
#--------------
rm(list=ls(all=TRUE))
par(mfrow=c(1,2))

#
# Recupération données
#---------------------
dataset=read.table(file.choose(), header = TRUE, sep = ";", dec=".")
# la commande suivante permet d'utiliser les entêtes de colonnes  
# comme noms de variables utilisables séparément
nrow(dataset)
ncol(dataset)
attach(dataset)
names(dataset)

# Plots
#------
plot(metal~dist,col="blue",pch=16,las=1,frame.plot=FALSE,xaxs="i",yaxs="i",
xlim=c(-10,150),ylim=c(-10,300))
# Comparaison 2 modèles
#----------------------
mod2=lm(metal~dist)
summary(mod2)
abline(mod2,col="red",lwd=2)
par(new=TRUE)
curve(exp(5.4789-0.01965705*x),from=0,to=150,las=1,xlim=c(-10,150),
ylim=c(-10,300),xaxs="i",yaxs="i",xlab="",ylab="")
mod=lm(log(metal)~dist)
summary(mod)
plot(log(metal)~dist,col="blue",pch=16,las=1,frame.plot=FALSE,xaxs="i",yaxs="i",
xlim=c(-10,150),ylim=c(0,6))
abline(mod,col="red",lwd=2)
mod$coef