oldfaith <-read.table(file = data1,header = T,sep = ',')
plot(oldfaith$DURATION,oldfaith$INTERVAL,ylab="Interval (mins)",xlab="Duration (mins)", main="Old Faithful Data")
fit<-lm(oldfaith$INTERVAL~oldfaith$DURATION)
abline(fit)