Regresión no lineal y causalidad
En Santa Cruz de la Sierra, en los meses de calor se venden más helados y también hay más accidentes en piscinas y balnearios. Los datos "van juntos". ¿Quiere decir que comer helado provoca accidentes? Y otra pregunta: cuando los datos no siguen una recta, ¿cómo los modelamos? En este apartado vas a responder las dos cosas. Al final volvemos al helado.
Cuándo una recta no es el modelo adecuado
No toda nube de puntos sigue una recta. A veces sube y después baja, a veces se acelera cada vez más. Mirar la forma de la nube es el primer paso para elegir el modelo.
Las tres nubes de puntos de la figura corresponden a tres estudios distintos. Para cada una de ellas, indica qué tipo de modelo parece más adecuado —lineal (…
Las tres nubes de puntos de la figura corresponden a tres estudios distintos. Para cada una de ellas, indica qué tipo de modelo parece más adecuado —lineal ($y = a + bx$), cuadrático ($y = a + bx + cx^2$) o exponencial ($y = a \cdot b^x$)— y justifica brevemente tu elección.
Un observatorio de la AEMET en una capital del interior peninsular registra las temperaturas en un caluroso día de julio desde las 08:00 ($x = 0$) hasta las 20:00 (…
Un observatorio de la AEMET en una capital del interior peninsular registra las temperaturas en un caluroso día de julio desde las 08:00 ($x = 0$) hasta las 20:00 ($x = 12$). Los termómetros marcan $21\,^\circ\text{C}$ a primera hora, alcanzan un máximo de $41\,^\circ\text{C}$ sobre las 16:00 ($x = 8$) y descienden a $33\,^\circ\text{C}$ a las 20:00. Un equipo de estudiantes duda entre tres modelos para ajustar los datos: el lineal $y = a + bx$, el exponencial $y = a \cdot b^x$ ($a > 0$, $b > 0$) y el cuadrático $y = ax^2 + bx + c$.
- a) Compara las propiedades gráficas y de crecimiento de los tres modelos y descarta razonadamente dos de ellos explicando por qué su comportamiento es físicamente incompatible con la evolución de la temperatura en ese intervalo.
- b) Para el modelo que has considerado idóneo, indica justificadamente qué signo debe tener el parámetro $a$ para reflejar la existencia de una temperatura máxima diaria.
- c) Si al calcular el coeficiente de correlación lineal se obtiene un valor bajo de $r = 0,38$, un estudiante afirma: «Como $r$ es tan bajo, queda demostrado que la hora del día apenas influye en la temperatura». Explica qué error conceptual comete al interpretar este valor de $r$.
Un grupo de biólogos monitoriza en un humedal de Doñana la propagación de una especie vegetal invasora. Se mide la superficie afectada $y$ (en $\text{m}^2$) en función…
Un grupo de biólogos monitoriza en un humedal de Doñana la propagación de una especie vegetal invasora. Se mide la superficie afectada $y$ (en $\text{m}^2$) en función del tiempo transcurrido $x$ (en semanas) desde que comenzó el seguimiento, obteniendo las siguientes mediciones:
-
A las $x = 0$ semanas: $y = 10\text{ m}^2$.
-
A las $x = 2$ semanas: $y = 40\text{ m}^2$.
-
A las $x = 4$ semanas: $y = 160\text{ m}^2$.
- a) Observa la pauta de variación de los datos y realiza una estimación intuitiva de qué superficie ocupará la planta invasora a las $x = 5$ semanas y a las $x = 6$ semanas, explicando el razonamiento que has utilizado.
- b) Para comprobar matemáticamente tu estimación, aplica la transformación logarítmica $z = \ln y$ a los tres pares de datos. Demuestra que los puntos transformados $(x, z)$ están perfectamente alineados en una recta $z = \ln a + x \cdot \ln b$ y calcula los valores exactos de los parámetros $a$ y $b$ del modelo exponencial $y = a \cdot b^x$.
- c) Sustituye $x = 5$ en el modelo analítico obtenido, comprueba si coincide con tu estimación previa y valora si este modelo exponencial sigue siendo verosímil para predecir la superficie ocupada tras un año completo ($x = 52$ semanas), teniendo en cuenta los límites físicos del ecosistema.
Qué tan bueno es un ajuste: el coeficiente de determinación
Ya sabes que el coeficiente de correlación lineal $r$ está entre $-1$ y $1$. Para saber cuánto "explica" el modelo, se usa su cuadrado.
Para comparar un modelo lineal con uno cuadrático se usa el coeficiente de determinación de cada uno (el $r^2$ de la recta y el $R^2$ de la parábola): mayor valor, más variación explicada de esos datos. Pero agregar parámetros puede mejorar el ajuste sin mejorar las predicciones, así que conviene mirar también el sentido del fenómeno.
En un ajuste lineal por mínimos cuadrados se ha obtenido un coeficiente de correlación $r = -0,9$. Calcula el coeficiente de determinación $r^2$ e interpreta su…
En un ajuste lineal por mínimos cuadrados se ha obtenido un coeficiente de correlación $r = -0,9$. Calcula el coeficiente de determinación $r^2$ e interpreta su significado: ¿qué porcentaje de la variación de la variable $y$ queda explicado por la variación de $x$?
Para los mismos datos se han probado dos modelos: un ajuste lineal con $r = 0,84$ y un ajuste cuadrático con coeficiente de determinación $r^2 = 0,93$. a) Calcula el…
Para los mismos datos se han probado dos modelos: un ajuste lineal con $r = 0,84$ y un ajuste cuadrático con coeficiente de determinación $r^2 = 0,93$.
- a) Calcula el coeficiente de determinación del modelo lineal.
- b) ¿Qué porcentaje de la variación de $y$ explica cada modelo?
- c) ¿Cuál elegirías para hacer predicciones? Justifica tu respuesta.
En una distribución bidimensional, la recta de regresión de $y$ sobre $x$ tiene pendiente $m_{yx} = 1,2$ y la recta de regresión de $x$ sobre $y$ tiene pendiente…
En una distribución bidimensional, la recta de regresión de $y$ sobre $x$ tiene pendiente $m_{yx} = 1,2$ y la recta de regresión de $x$ sobre $y$ tiene pendiente $m_{xy} = 0,75$.
- a) Utilizando la relación $m_{yx} \cdot m_{xy} = r^2$, calcula el coeficiente de determinación.
- b) Deduce el valor del coeficiente de correlación $r$, razonando su signo.
- c) Interpreta el valor de $r^2$ como porcentaje de variación explicada y valora la fiabilidad de las predicciones hechas con esas rectas.
Linealizar un modelo exponencial con logaritmos
Si la nube se acelera, puede ser exponencial, como un cultivo de bacterias en el laboratorio. En el laboratorio se ajusta con Excel o GeoGebra (línea de tendencia exponencial); aquí, como ampliación, ves cómo lo hace el programa por dentro: un cambio de variable que convierte la curva en una recta.
Lucía analiza la depreciación del valor de mercado de un vehículo eléctrico en el mercado de segunda mano español en función de su antigüedad $x$ (en años). Para ajustar…
Lucía analiza la depreciación del valor de mercado de un vehículo eléctrico en el mercado de segunda mano español en función de su antigüedad $x$ (en años). Para ajustar un modelo exponencial decreciente $y = a \cdot b^x$ ($y > 0$), aplica la transformación logarítmica neperiana $z = \ln y$ y obtiene mediante mínimos cuadrados la recta:
$$z = 10,2 - 0,15x$$En las conclusiones de su trabajo de fin de curso, Lucía escribe:
«El precio de salida del vehículo nuevo es de $a = 10,2$ bolivianos y cada año que pasa su valor se multiplica por $b = -0,15$. Como el coeficiente de correlación de la recta transformada es $r = -0,98$, el paso del tiempo es la causa directa y única de toda la pérdida de valor del automóvil».
- a) Localiza los dos errores algebraicos que comete Lucía al deducir los parámetros $a$ y $b$ a partir de la recta linealizada y calcula sus valores correctos (puedes dejarlos expresados como potencias del número $e$ y dar una aproximación con dos decimales).
- b) Detecta el error conceptual en su afirmación sobre la causalidad y propón al menos una variable latente o ajena a los años del vehículo que también influya decisivamente en su valor de reventa.
El número de usuarios de una aplicación crece de forma aproximadamente exponencial: $y = a \cdot b^x$, donde $x$ son los meses transcurridos. Tomando logaritmos…
El número de usuarios de una aplicación crece de forma aproximadamente exponencial: $y = a \cdot b^x$, donde $x$ son los meses transcurridos. Tomando logaritmos neperianos, $\ln y = \ln a + x \cdot \ln b$, de modo que los pares $(x, \ln y)$ se ajustan a una recta. Los datos linealizados son:
| $x$ (meses) | $0$ | $1$ | $2$ | $3$ |
|---|---|---|---|---|
| $\ln y$ | $1,0$ | $1,7$ | $2,4$ | $3,1$ |
- a) Comprueba que los cuatro puntos linealizados están exactamente sobre una recta y halla su ecuación.
- b) Deduce los valores de $a$ y $b$ (expresa $a$ y $b$ como potencias de $e$ y aproxímalas con dos decimales).
- c) Estima el número de usuarios al cabo de $6$ meses.
La superficie de un lago cubierta por una planta invasora se ha medido durante $5$ meses: | $x$: mes | $1$ | $2$ | $3$ | $4$ | $5$ | |---|---|---|---|---|---| | $y$:…
La superficie de un lago cubierta por una planta invasora se ha medido durante $5$ meses:
| $x$: mes | $1$ | $2$ | $3$ | $4$ | $5$ |
|---|---|---|---|---|---|
| $y$: superficie ($\text{m}^2$) | $3,5$ | $7,1$ | $14,0$ | $28,2$ | $56,5$ |
- a) Representa los pares $(x, \ln y)$, comprueba que se alinean razonablemente y ajusta una recta a los datos linealizados (tómate $\ln y$ con dos decimales).
- b) Escribe el modelo exponencial $y = a \cdot b^x$ correspondiente y calcula el ritmo de crecimiento mensual en porcentaje.
- c) Predice la superficie cubierta en el mes $8$ y, sabiendo que el lago tiene $1.500\ \text{m}^2$, comenta qué problema presenta extrapolar este modelo a largo plazo.
Linealizar un modelo potencial
En un laboratorio de Física se mide el periodo $T$ de un péndulo en función de su longitud $L$, y se sospecha una relación potencial $T = a \cdot L^b$. Tomando…
En un laboratorio de Física se mide el periodo $T$ de un péndulo en función de su longitud $L$, y se sospecha una relación potencial $T = a \cdot L^b$. Tomando logaritmos, $\ln T = \ln a + b \cdot \ln L$. Los datos linealizados son:
| $\ln L$ | $-1,386$ | $-0,693$ | $0$ | $0,693$ |
|---|---|---|---|---|
| $\ln T$ | $0$ | $0,347$ | $0,693$ | $1,040$ |
- a) Ajusta una recta a los datos linealizados y escribe su ecuación.
- b) Deduce los valores de $a$ y $b$.
- c) Compara el valor obtenido de $b$ con el que predice la Física, $b = \frac{1}{2}$, y comenta el resultado.
La Dirección General de Tráfico (DGT) realiza ensayos de seguridad vial para modelizar la distancia total de frenado de emergencia $y$ (en metros) de un automóvil sobre…
La Dirección General de Tráfico (DGT) realiza ensayos de seguridad vial para modelizar la distancia total de frenado de emergencia $y$ (en metros) de un automóvil sobre pavimento mojado según la velocidad inicial $x$ (medida en decenas de km/h, es decir, $x = 2$ representa $20\text{ km/h}$). Los resultados experimentales son los siguientes:
-
Para $x = 2$: $y = 4\text{ m}$.
-
Para $x = 4$: $y = 16\text{ m}$.
-
Para $x = 8$: $y = 64\text{ m}$.
- a) Se sospecha que la relación responde a un modelo potencial $y = a \cdot x^b$. Aplica la doble transformación logarítmica neperiana ($w = \ln x$ y $z = \ln y$) y expresa las coordenadas transformadas $(w_i, z_i)$ en función de $\ln 2$.
- b) Comprueba que los tres puntos transformados pertenecen a una misma recta $z = \ln a + b \cdot w$, calcula algebraicamente su pendiente $b$ y su ordenada en el origen $\ln a$, y escribe la ecuación definitiva del modelo potencial $y = f(x)$.
- c) Interpreta el valor numérico del exponente $b$ apoyándote en la física (fórmula de la energía cinética, $E_c = \frac{1}{2}m v^2$). ¿Por qué este modelo evidencia que duplicar la velocidad de conducción implica un aumento mucho más drástico del riesgo de impacto?
Un estudio del Instituto Nacional de Estadística (INE) sobre el desarrollo socioeconómico en municipios del área rural de Bolivia analiza la relación entre la cobertura…
Un estudio del Instituto Nacional de Estadística (INE) sobre el desarrollo socioeconómico en municipios del área rural de Bolivia analiza la relación entre la cobertura de fibra óptica ($x$, en $\%$ de hogares conectados) y la productividad agraria media ($y$, en miles de bolivianos por unidad productiva al año). Para ajustar los datos de 50 municipios, un equipo econométrico ensaya dos linealizaciones logarítmicas:
-
Opción 1 (regresión de $\ln y$ sobre $x$): obtiene la recta $z = 1,38 + 0,02 x$ con coeficiente de determinación $R_1^2 = 0,82$.
-
Opción 2 (regresión de $\ln y$ sobre $\ln x$): obtiene la recta $z = 0,69 + 0,55 w$ (donde $w = \ln x$) con coeficiente de determinación $R_2^2 = 0,95$.
- a) Deshaz algebraicamente las transformaciones y escribe la expresión explícita $y = f(x)$ que corresponde a cada uno de los dos modelos (puedes aproximar los coeficientes que involucren potencias de $e$ a dos cifras decimales).
- b) Analiza el comportamiento de la tasa de crecimiento de la productividad en ambos modelos cuando la cobertura de fibra $x$ se hace muy grande. Justifica cuál de los dos refleja mejor el principio económico de rendimientos decrecientes (donde cada punto porcentual extra de cobertura aporta un incremento de productividad cada vez menor) y compáralo con el criterio estadístico que proporciona el coeficiente $R^2$.
- c) El alcalde de una localidad concluye en una sesión del concejo municipal: «Dado que el modelo potencial tiene un ajuste excelente de $R^2 = 0,95$, invertir todo el presupuesto disponible en alcanzar el $100\,\%$ de fibra óptica garantizará directamente que la productividad de nuestros agricultores alcance el máximo previsto». Realiza una crítica metodológica de esta conclusión analizando por qué una correlación no implica causalidad y señalando al menos dos variables ocultas (o factores de confusión) que probablemente expliquen simultáneamente una alta conectividad y una alta productividad.
Ajuste cuadrático por mínimos cuadrados
Cuando la nube tiene un máximo o un mínimo, buscamos la parábola $y = ax^2 + bx + c$ que minimiza la suma de los cuadrados de los errores, igual que en la recta. En el laboratorio la obtienes con Excel o GeoGebra (línea de tendencia polinómica de grado 2).
Se sospecha que la variable $y$ depende cuadráticamente de $x$. Con los cinco datos de la tabla | $x$ | $0$ | $1$ | $2$ | $3$ | $4$ | |---|---|---|---|---|---| | $y$ |…
Se sospecha que la variable $y$ depende cuadráticamente de $x$. Con los cinco datos de la tabla
| $x$ | $0$ | $1$ | $2$ | $3$ | $4$ |
|---|---|---|---|---|---|
| $y$ | $5$ | $3$ | $3$ | $5$ | $9$ |
se han calculado las sumas $\sum x_i = 10$, $\sum x_i^2 = 30$, $\sum x_i^3 = 100$, $\sum x_i^4 = 354$, $\sum y_i = 25$, $\sum x_i y_i = 60$ y $\sum x_i^2 y_i = 204$.
- a) Plantea el sistema de tres ecuaciones que deben cumplir los coeficientes de la parábola de regresión $y = a x^2 + bx + c$.
- b) Resuélvelo y escribe la ecuación del ajuste cuadrático.
- c) Estima el valor de $y$ para $x = 5$ y comenta la fiabilidad de esa predicción.
El consumo de un auto eléctrico, en $\text{kWh}$ cada $100\ \text{km}$, se ha medido a distintas velocidades y los datos aparecen en la figura: | $x$: velocidad…
El consumo de un auto eléctrico, en $\text{kWh}$ cada $100\ \text{km}$, se ha medido a distintas velocidades y los datos aparecen en la figura:
| $x$: velocidad ($\text{km/h}$) | $40$ | $60$ | $80$ | $100$ | $120$ |
|---|---|---|---|---|---|
| $y$: consumo ($\text{kWh}/100\ \text{km}$) | $13,1$ | $12,4$ | $12,5$ | $13,4$ | $15,1$ |
- a) Representa mentalmente (o bosqueja) la nube de puntos y razona por qué un modelo cuadrático es más adecuado que uno lineal.
- b) Calcula las sumas necesarias $\sum x_i$, $\sum x_i^2$, $\sum x_i^3$, $\sum x_i^4$, $\sum y_i$, $\sum x_i y_i$, $\sum x_i^2 y_i$ y plantea el sistema del ajuste cuadrático.
- c) Resuelve el sistema y escribe la ecuación de la parábola de regresión (redondea los coeficientes a cuatro decimales).
- d) Halla la velocidad de consumo mínimo según el modelo y predice el consumo a $140\ \text{km/h}$, valorando la fiabilidad de esa predicción.
Para modelizar la trayectoria de un proyectil en una prueba balística se toma una muestra de $N = 5$ puntos con abscisas simétricas respecto al origen: $x_1 = -2$,…
Para modelizar la trayectoria de un proyectil en una prueba balística se toma una muestra de $N = 5$ puntos con abscisas simétricas respecto al origen: $x_1 = -2$, $x_2 = -1$, $x_3 = 0$, $x_4 = 1$ y $x_5 = 2$. Se desea ajustar una parábola por mínimos cuadrados $y = ax^2 + bx + c$, cuyo sistema general de ecuaciones normales es:
$$\begin{cases} a \sum x_i^4 + b \sum x_i^3 + c \sum x_i^2 = \sum x_i^2 y_i \\ a \sum x_i^3 + b \sum x_i^2 + c \sum x_i = \sum x_i y_i \\ a \sum x_i^2 + b \sum x_i + c \cdot N = \sum y_i \end{cases}$$- a) Calcula las sumas $\sum x_i$, $\sum x_i^2$, $\sum x_i^3$ y $\sum x_i^4$. Explica qué propiedad de la simetría de las $x_i$ hace que el parámetro $b$ quede completamente desacoplado de $a$ y $c$, simplificando el sistema de 3 ecuaciones en dos bloques independientes.
- b) Plantea el problema inverso: se pretende que la curva ajustada tenga vértice en $(0, 4)$ con apertura hacia abajo, resultando exactamente el modelo $y = -0,5x^2 + 4$. Si la suma de las ordenadas observadas es $\sum y_i = 15$, deduce qué valor exacto debe tomar la suma ponderada $\sum x_i^2 y_i$ para que el modelo satisfaga el sistema de mínimos cuadrados.
- c) Construye una tabla con cinco valores concretos $(x_i, y_i)$ que cumplan simultáneamente la simetría $y_1 = y_5$, $y_2 = y_4$, $\sum y_i = 15$ y que encajen con dicha parábola. Si calculases el coeficiente de correlación lineal tradicional $r$ para esos 5 puntos, ¿qué valor exacto obtendrías? Razona la aparente paradoja entre el valor de $r$ y la fuerza de la dependencia funcional existente.
Predecir con cuidado: interpolar y extrapolar
Un modelo se ajusta con datos de cierto rango. Usarlo dentro de ese rango es interpolar; usarlo fuera es extrapolar.
Hay un caso especial: variables con una cota natural. Un porcentaje no puede superar el $100\,\%$. Si una adopción de tecnología va creciendo, una recta la hace pasar de $100\,\%$ en algún momento, lo que no tiene sentido.
Para una variable acotada conviene un modelo que se "aplane" cerca de la cota, y no una recta.
Con los datos de una distribución bidimensional en la que $x$ varía entre $10$ y $30$, con $\bar{x} = 20$ y $r = 0,97$, se ha obtenido una recta de regresión. Dos…
Con los datos de una distribución bidimensional en la que $x$ varía entre $10$ y $30$, con $\bar{x} = 20$ y $r = 0,97$, se ha obtenido una recta de regresión. Dos alumnos la usan para estimar $y$: uno para $x = 22$ y otro para $x = 55$.
- a) ¿Qué estimación es más fiable? Razona la respuesta usando la posición de cada valor respecto del rango de los datos y de $\bar{x}$.
- b) ¿Cómo se llama el tipo de predicción que se hace fuera del rango de los datos y por qué es arriesgada?
La figura muestra el porcentaje de hogares bolivianos con conexión de fibra óptica en función de los años transcurridos desde 2015, según datos del INE. a) Ajusta una…
La figura muestra el porcentaje de hogares bolivianos con conexión de fibra óptica en función de los años transcurridos desde 2015, según datos del INE.
- a) Ajusta una recta de regresión a los cinco primeros datos usando las medias $\bar{x} = 4$, $\bar{y} = 64,6$ y los valores $s_{xy} = 55,2$, $s_x^2 = 8$.
- b) Predice el porcentaje para $x = 12$ (año 2027) con esa recta.
- c) Observando la forma de la nube de puntos completa, critica la validez de esa predicción a largo plazo: ¿qué cota natural tiene la variable y qué tipo de modelo sería más sensato?
Correlación no es causalidad
Volvamos al helado. Una correlación alta indica que dos variables se mueven juntas, no que una provoque a la otra.
Para desenmascarar una relación espuria se puede controlar la variable oculta: comparar solo datos del mismo nivel de esa variable (por ejemplo, solo meses de verano) y ver si la correlación desaparece. Lo ideal es un experimento donde se asigne al azar el factor que se quiere estudiar.
Vuelve al principio: el helado no provoca accidentes; los dos suben con el calor y la cantidad de gente en las piscinas y balnearios.
Razona si la siguiente afirmación es verdadera o falsa: «El coeficiente de correlación entre las ventas de helados y el número de ahogamientos en piscinas es alto y…
Razona si la siguiente afirmación es verdadera o falsa: «El coeficiente de correlación entre las ventas de helados y el número de ahogamientos en piscinas es alto y positivo; por tanto, vender más helados provoca más ahogamientos». Si es falsa, explica el error de interpretación y nombra una variable que podría explicar ambas a la vez.
En una ciudad se observa una correlación $r = 0,92$ entre el número de camiones de bomberos que acuden a un incendio y los daños económicos que este provoca. Un político…
En una ciudad se observa una correlación $r = 0,92$ entre el número de camiones de bomberos que acuden a un incendio y los daños económicos que este provoca. Un político afirma: «Para reducir los daños de los incendios, deberíamos enviar menos camiones de bomberos». Analiza la falacia de su razonamiento, identifica la variable que explica la correlación y formula una conclusión correcta a partir del mismo dato.
En varios distritos de una región europea se ha observado una correlación $r = 0,85$ entre el número de cigüeñas censadas y el número de nacimientos registrados. a)…
En varios distritos de una región europea se ha observado una correlación $r = 0,85$ entre el número de cigüeñas censadas y el número de nacimientos registrados.
- a) Explica por qué sería absurdo concluir que las cigüeñas «traen» los bebés, identificando el tipo de correlación que presenta este ejemplo.
- b) Propón una variable oculta plausible que explique la correlación y razona cómo actúa sobre ambas variables.
- c) Diseña brevemente qué datos recogerías y qué cálculo harías para desenmascarar la relación espuria.