```r
residuals<-rnorm(30,mean=0,sd=5)
results<-fitted+residuals
d<-data.frame(predictor,results)
boxplot(results~predictor,data=d)
```

<img src="Quantitative_and_Spatial_files/figure-html/unnamed-chunk-129-1.png" width="672" />