CAPÍTULO 7. AJUSTE CON MODELOS ESTADÍSTICOS
179
l
l l
l
l
l l
l
l
l
l
l
0 2 4 6 8 10
0
5
10 15
x
y
Figura 7.3: Medidas del lanzamiento de un proyectil
modelo.tiro
<-
lm
(formula.tiro,
data
=datos.tiro)
# Veamos la información básica del modelo:
modelo.tiro
##
## Call:
## lm(formula = formula.tiro, data = datos.tiro)
##
## Coefficients:
## (Intercept)
x
I(x^2)
##
14.902
1.069
-0.224
De esta manera, la ecuación de la trayectoria del proyectil ajustada a los
datos estaría dada por:
y
=
14.902
+
1.069
x
−
0.224
x
2
(7.5)
Tanto para el uso del modelo, como para graficar la curva junto con los
datos que lo originaron, es conveniente hacer una función basada en la función
predict()
del lenguaje, de la manera siguiente:
func.tiro
<-
function
(
x
) {
# El argumento x puede ser un solo número o un vector de
# números
predict
(modelo.tiro,
newdata
=
data.frame
(
x
=x),
type
=
"response"
)