Regresión lineal

y no lineal

Sesión 9
Logro:

El alumno aprende las técnicas de regresión a través de métodos

numéricos, particularmente se desarrolla problemas orientados a las
diferentes ramas de la ingeniería.

Datos/Observaciones
 REGRESIÓN LINEAL
FORMULACIÓN DEL PROBLEMA DE LA REGRESIÓN
Dados n puntos (x1, y1), (x2, y2), … , (xn, yn) hallar la función
y= f(x), que se ajuste mejor a ellos.
El criterio que normalmente se usa para establecer el mejor
ajuste es el de mínimos cuadrados que consiste en minimizar
la suma de los cuadrados de los errores residuales.
El error residual para el punto i, se define como:

i = yi − f ( x i )

La suma de los cuadrados de los errores residuales se

expresa como:
n n 2

Sr =  =  (yi − a0 − a1xi )
2
i
i =1 i =1
2
 REGRESIÓN LINEAL

Dados n puntos (x1, y1), (x2, y2), … , (xn, yn) hallar la

función y = a 0 + a1 x que se ajuste mejor a ellos.
Para este caso el error residual para el dato i se puede
expresar como:
εi = yi − (a0 + a 1 xi )

La suma de los cuadrados de los errores residuales, en

este caso, sería:
n n 2

Sr =  i
2
=  (yi − a0 − a1xi )
i =1 i =1

3
 REGRESIÓN LINEAL
Las constantes a0 y a1, del modelo lineal serán aquellas que
minimicen a Sr, lo cual requieres derivar Sr con respecto a cada uno de
los coeficientes:
S r n
= −2 (yi − a0 − a1 xi )(− 1) = 0
a0 i =1

S r
= −2 (yi − a0 − a1 xi )(− xi ) = 0
n

a1 i=1

Lo cual puede escribirse como:

n n n

a +a x =y
0 1 i i
i=1 i=1 i=1
n n n

a x + a x =  yi xi
2
0 i 1 i
i=1 i=1 i=1
 REGRESIÓN LINEAL
n

Teniendo en cuenta  a 0 = na 0 se obtiene:

i=1
n n
na0 +  a1 xi =  yi
i=1 i=1

n n n

a x + a x =  yi xi
2
0 i 1 i
i=1 i=1 i=1

Que se puede expresar matricialmente como;

 n
  n 
 n  i  a0    i 
x y
 n i=1
n
   =  ni=1 
2   a1 
 xi
i=1 
i=1
x i

 yi xi 
 i =1 

5
 REGRESIÓN LINEAL
Aplicando la regla de Crámer, para resolver este
sistema, se ontiene
n n

 yi  xi
i=1 i=1
n n n n n n

i=1
y ix i i=1
x i
2
 y i x i
2
−  yi xi  xi
a0 = n
= i =1 i =1
n
i =1
n n
i =1

n  xi n  xi −
i=1
2
x i xi
i=1 i=1 i=1
n n


i=1
xi 
i=1
x i2

n n n n

x y −x x y
i
2
i i i i
a0 = i=1 i=1 i=1 i=1
2
n
 n 
n x −  xi  i
2

i=1  i=1 

6
 REGRESIÓN LINEAL

n n

 yi i =1
xi
i=1
n n

n n n n

i=1
yixi 
i=1
x i
2
 y i x i
2
−  yi xi  xi
a0 = n
= i=1 i=1
n
i=1
n n
i=1

n  xi n x i
2
−  xi  xi
i=1 i=1 i=1 i=1
n n


i=1
xi 
i=1
x i2

n n n n

x y −x x yi
2
i i i i
a0 = i=1 i=1 i=1 i=1
2
n
 n 
n x −  xi 
i
2

i=1  i=1 

7
 REGRESIÓN LINEAL

n
n  i=1
yi
n n n n n
 xi  yi xi n  yi xi −  xi  yi
i=1 i=1
a1 = n
= i=1
n
i=1
n
i=1
n
n  xi n x i
2
−  xi  xi
i=1 i=1 i=1 i=1
n n

i=1
xi i=1
x i2

n n n
n xi yi − xi  yi
a1 = i =1 i =1 i =1
2
2 
n n

n xi −  xi 
i=1  i=1 
 EJEMPLO
Tomado del texto: Chapra, Steven C y Canale, Raymond. Métodos Numéricos para
ingeniero. Quinta edición
Se sabe que el esfuerzo a la
tensión de un plástico se Contador Tiempo Esfuerzo a la
incrementa como función del (min) tensión(N/cm2)
tiempo que recibe tratamiento i xi yi
a base de calor. Se obtuvieron
1 10 5
los datos que se muestran en
2 15 20
la tabla.
3 20 18
Ajuste una línea recta a estos 4 25 40
datos y utilice una ecuación 5 40 33
para determinar el esfuerzo a
6 50 54
la tensión en un tiempo de 32
7 55 70
minutos
8 60 60
9 75 78
 EJEMPLO
esfuerzo a la tensión(N /cm 2)
90

80

70

60

50
tensión(N/cm2

40
Esfuerzoala

30

20

10

0
0 10 20 30 40 50 60 70 80
Tie m p o (m in )

GRÁFICO DE DISPERSIÓN PARA: Esfuerzo a la tensión vs. Tiempo

 EJEMPLO
La solución básicamente consiste en hallar los coeficientes de la
ecuación de la recta de mejor ajuste. Para el problema planteado, la
forma de la ecuación de la recta sería:
y = a0 + a1 x
Donde y= esfuerzo a la tensión(N/cm2)
x= tiempo(min)

Para determinar los coeficientes a0 y a1, del modelo lineal, se utiliza

el procedimiento de regresión lineal, para lo cual se deben calcular las
siguientes cantidades:

n n n n
 xi yi
i =1
 x i2 x i y i
i =1 i =1 i=1
 EJEMPLO
Teniendo en cuenta que para este problema, n=9 datos, entonces:
9

 x i y i = x1 y 1 + x 2 y 2 + x 3 y 3 + x 4 y 4 + x 5 y 5 + x 6 y 6 + x 7 y 7 + x 8 y 8 + x 9 y 9
i =1
9

 x y =105+1520 + 2018+ 2540 + 4033+5054 +5570 + 6060 + 7578

i i
i =1
9

 x y =19030
i i
i =1

x = x +x i 1 2 + x3 + x4 + x5 + x6 + x7 + x8 + x9
i =1

 x =10 +15+ 20 + 25+ 40 + 50 + 55 + 60 + 75

i
i =1

 x = 350 i
i =1
 EJEMPLO
9

y =y +y i 1 2 + y3 + y4 + y5 + y6 + y7 + y8 + y9
i =1

 y = 5+ 20 +18+ 40 + 33+ 54 + 70 + 60 + 78
i
i =1

 x = 378 i
i =1

x
i =1
2
i
= x12 + x22 + x32 + x42 + x52 + x26 + x27 + x28 + x29

x
i =1
2
i
= 102 +152 + 202 + 252 + 402 + 502 + 552 + 602 +
752
9

xi =1
2
i =17700
 EJEMPLO
Con estas cantidades calculadas, se pueden utilizar las
fórmulas para los coeficientes, obtenidas anteriormente:
n n n n

x y −x x y
i
2
i i i i
a0 = i=1 i=1 i=1 i=1
2
n
 n 
n x −  xi 
i
2

i=1  i=1 

n n n
n xi yi − xi  yi
a1 = i =1 i =1 i =1
2

n
 n
n xi2 −  xi 
i=1  i=1 
 EJEMPLO
Reemplazando:
9

 xi = 378
9 9 9

 x y =19030
i i  x = 350 i
i =1
x 2
i =17700
i =1 i =1 i =1

Se obtiene:

17700  378 − 350 19030

a0 = = 0.81793478
9 17700 −(350 )
2

9 19030 −350  378

a1 = = 1.05896739
9 17700 − (350 )
2
 EJEMPLO
Reemplazando estos valores en el modelo lineal, se obtiene
la ecuación de la recta de mejor ajuste:
y =a0 + a1 x

y = 0.81793478+1.05896739x
Remplazando y, x por los nombre reales de las
variables, quedaría:

Esfuerzo a la tensión(N / cm2 ) = 0.81793478+1.05896739  tiempo (min )

En la siguiente diapositiva se presenta una gráfica de los

datos iniciales con la recta de regresión
 EJEMPLO
esfuerzo a la tensión(N /cm 2 )
90

80
y = 1 ,0 5 9 x + 0 ,8 1 7 9
70

60

50
tensión(N/cm2

40
Esfuerzoala

30

20

10

0
0 10 20 30 40 50 60 70 80
T ie m p o ( m in )

Datos originales con la recta de regresión obtenida

 EJEMPLO

Para estimar el valor del esfuerzo a la tensión para un

tiempo igual a 32 minutos, se utiliza la ecuación
obtenida:

Esfuerzo a la tensión(32) = 0.81793478+1.05896739  32

Esfuerzo a la tensión(32) = 34.7048913

En la siguiente figura se ilustra gráficamente la

obtención del esfuerzo para tiempo=32 minutos
 EJEMPLO
esfuerzo a la tensión(N/cm 2)
90

80
y = 1 ,0 5 9 x + 0 ,8 1 7 9
70

60

50
tensión(N/cm2

40 x = 34.7048913
Esfuerzoala

30

20

10
x = 32
0
0 10 20 30 40 50 60 70 80
Tie m p o (m in )