sink('NO2_EX.sink') a=read.table('NO2_EX.txt', col.names=c("logNO2","lognbcar","temp","vitvent","difftemp","dirvent","heure","jour")) attach(a) m1.lm=lm(logNO2~lognbcar+temp+vitvent+difftemp) print(summary(m1.lm)) ####### print(anova(m1.lm)) sink()