source("http://192.38.117.59/~linearpredictors/datafiles/readVitaminD.R")

main.model <- lm(vitd ~ bmi, data = irlwomen)

outliers <- data.frame(bmi = c(25, 50, 50, 15.9375), vitd = c(130, 80, 10, 39.2))

## fitting models 1 to 4, adding different outliers
for (i in 1:4){
  ## generating a model name
  m <- paste("model", i, sep="")

  ## fitting the model
  assign(m,
         lm(vitd ~ bmi,
            data = data.frame(vitd = c(irlwomen$vitd, outliers$vitd[i]),
              bmi = c(irlwomen$bmi, outliers$bmi[i]))
            )
         )
  print(summary(get(m)))
  
  ## prediction of vitaminD at bmi=25
  print(predict(get(m), data.frame(bmi=25), se.fit=T))

}