CAPÍTULO 7. AJUSTE CON MODELOS ESTADÍSTICOS
180
}
Y el problema de saber la altura del proyectil a 5.5 unidades, junto con otros
valores, se resuelve así:
# Altura a 5.5 unidades horizontales:
func.tiro
(
5.5
)
## 1
## 14
# Si se quiere conocer las alturas predichas por el modelo
# para las distancias 3.5, 4.3 y 8.2:
func.tiro
(
c
(
3.5
,
4.3
,
8.2
))
##
1
2
3
## 15.900 15.357 8.603
Finalmente, si se quiere graficar la curva correspondiente al modelo, junto
con los datos, se hace agregando al gráfico de la Fig. 7.3, la salida del código
que se presenta a continuación. En este código se está añadiendo, además, sólo
para comparar, lo que sería el resultado si los datos se hubieran ajustado a una
línea recta.
# La curva ajustada:
curve
(func.tiro,
lwd
=
2
,
col
=
"blue"
,
add
=T)
# Hagamos el modelo de una recta:
mod.recta
<-
lm
( y
~
x,
data
=datos.tiro)
# Dibujémosla en el gráfico también:
abline
(mod.recta)
# va en negro, el default
# ---------------
# Pondremos una leyenda apropiada para identificar
# los objetos en el gráfico:
legend
(
"bottomleft"
,
# Donde va la leyenda
legend
=
c
(
"modelo parabola"
,
"modelo recta"
),
# textos
lwd
=
c
(
2
,
1
),
# anchos de línea
col
=
c
(
"blue"
,
"black"
)
# colores para cada modelo
)
El resultado de este código, añadido a la figura anterior, se puede apreciar
en la Fig. 7.4.
7.2. Modelos lineales generalizados
En los modelos que se han revisado en la sección anterior, se parte de la su-
posición de que la distribución de los residuos tiene una distribución normal,