```
set.seed(1) # so can regenerate results
x1 <- rnorm(100)
x2 <- rnorm(100)
y <- exp(x1+x2+rnorm(100))
f <- lm(log(y) ~ x1 + poly(x2,3), y=TRUE)
abs.error.pred(lp=exp(fitted(f)), y=y)
rm(x1,x2,y,f)
```

Run the code above in your browser using DataLab