Homework5_2

library(faraway)
data("faithful", package = "faraway")

## Warning in data("faithful", package = "faraway"): data set 'faithful' not found

head(faithful)

##   eruptions waiting
## 1     3.600      79
## 2     1.800      54
## 3     3.333      74
## 4     2.283      62
## 5     4.533      85
## 6     2.883      55

lmod<-lm(eruptions~waiting, data=faithful)
plot(eruptions~waiting, faithful, ylab="distance")
abline(lmod)

lmod1<-lm(eruptions~I(waiting+rnorm(272)),faithful)
coef(lmod1)

##             (Intercept) I(waiting + rnorm(272)) 
##             -1.86071287              0.07550788

require(simex)

## Loading required package: simex

set.seed(123)
lmod <- lm(eruptions ~ waiting, faithful, x=TRUE)
simout <- simex(lmod,"waiting",45, B=1000)
simout

## 
## Naive model:
## lm(formula = eruptions ~ waiting, data = faithful, x = TRUE)
## 
## SIMEX-Variables: waiting
## Number of Simulations: 1000
## 
## Coefficients:
## (Intercept)      waiting  
##    -11.4245       0.2104