1
2
3
4
5
6
7
8
9
10
11
12
## Buscamos un estimador de la pendiente de un modelo lineal simple de la
##   forma y=1+bx, probamos 50 posibles valores de b, y calculamos la suma
##   de cuadrados del error (lo que llamamos en clases "RSS") del modelo
##   para cada uno de esos valores de b.
library(animation)
library(MASS)
par(mar = c(4, 4, 0.5, 0.1), mgp = c(2, 0.5, 0), tcl = -0.3)
least.squares(y = y, x = x, ani.type = "slope", a = b0.hat,
    b.range = b.range, las = 1)
## R version 2.15.1 (2012-06-22)
## Platform: x86_64-pc-linux-gnu (64-bit)
## Other packages: animation 2.1, MASS 7.3-18