Changes between Version 1 and Version 2 of Piecewise

Timestamp:

Jun 13, 2014, 1:35:30 PM (11 years ago)

Author:

Pedro Gea

Comment:

--

Piecewise

@@ -6 +6 @@
 Un ejemplo muy conocido es la función valor absoluto, que puede describirse como la función identidad para los valores positivos o cero y como la función opuesto para los valores negativos:
 {{{
-         { x ,, x>=0
-abs(x) = {
-         { -x ,, x<0
+          | x ,, x>=0
+abs(x) = <|
+          | -x ,, x<0
 }}}
 
@@ -15 +15 @@
 En una regresión describimos una variable {{{Y}}} (el ''output'') como combinación lineal de otras variables {{{X_i}}} (los ''inputs''):
 {{{
-Y = Sum_i(beta_i * X_i) + N
+Y = Sum_i(beta_i *X_i) + N
 }}}
 más una parte no explicada {{{N}}} que es el ruido o residuo de la regresión.
@@ -25 +25 @@
 En otras ocasiones, esta forma funcional es conocida, pero depende de algún parámetro no lineal que puede ser estimado por los algoritmos correspondientes.
 
-En el caso que nos ocupa, se propone describir dicha forma no-lineal como una función lineal a trozos, de modo que el término {{{E_i = beta_i * X_i}}} se pueda reescribir como:
+En el caso que nos ocupa, se propone describir dicha forma no-lineal como una función lineal a trozos, de modo que un término {{{E(X) = beta * X}}} se pueda reescribir como:
 {{{
-E_i = beta_i,k * X_i ,,  X_i <: Omega_i,k
+E(X) = alpha_k + beta_k * X ,,  X <: [X_k, X_{k+1}]
 }}}
+cumpliendo dos condiciones:
+ * continuidad: {{{alpha_k + beta_k * X_{k+1} ==  alpha_{k+1} + beta_{k+1} * X_{k+1} }}}
+ * efecto nulo para input nulo: {{{E(0) == 0}}}
+
+== Particionamiento del input ==
+
+Para incorporar dicha función a trozos a la regresión, utilizaremos un particionamiento del input {{{P_k(X)}}} de modo que la función a trozos se pueda escribir como combinación lineal de dicho particionamiento:
+{{{
+E(X) = Sum_k(beta_k * P_k(X))
+}}}
+
+=== Particionamiento ===
+
+Se admite que las marcas utilizadas para la partición del input recorren todo el dominio del input. Ya que, de lo contrario, el input tendría que truncarse al intervalo definido por las marcas usadas.
+
+Sean {{{n+1}}} marcas: {{{x_k}}} con {{{k<:[1,n+1]}}} y un input {{{x}}} que será modelado con diferentes parámetros según sus valores estén en un intervalo {{{[x_k, x_{k+1}]}}} u otro:
+{{{
+        |  a_1 + b_1 * x  ,,  x<:[x_1, x_2]
+        |  ...
+E(x) = <|  a_m + b_m * x  ,,  x<:[x_m, x_{m+1}]
+        |  ...
+        |  a_n + b_n * x  ,,  x<:[x_n, x_{n+1}]
+}}}
+ 
+Se busca que el efecto ({{{E}}}) sea continuo de modo que el valor para una marca: {{{x_k}}} sea el mismo por la izquierda o por la derecha: 
+{{{
+a_{k-1} + b_{k-1} * x_k = a_k + b_k * x_k
+}}}
+y se quiere que la el efecto sea nulo para el input nulo (esté libre de constante u ordenada en el origen):
+{{{
+E(x=0) = 0
+}}}
+
+Si reescribimos la expresión anterior teniendo en cuenta cual es el intervalo que contiene al {{{0}}}: {{{0<:[x_q, x_{q+1}]}}}
+{{{
+        |  ... (+)
+        |  ...
+E(x) = <|  b_q * x                        ,,  x<:[x_q, x_(q+1)]
+        |  b_{q+1}*(x-x_{q+1})
+        |    + b_q*x_{q+1}                ,,  x<:[x_{q+1}, x_{q+2}]
+        |  b_{q+2}*(x-x_{q+2}) 
+        |    + b_{q+1}*(x_{q+2}-x_{q+1}) 
+        |    + b_q*x_{q+1}                ,,  x<:[x_{q+2}, x_{q+3}]
+        |  ...
+        |  b_n*(x-x_n) 
+        |    + b_{n-1}*(x_n-x_{n-1}) 
+        |    + ... + b_q*x_{q+1}          ,,  x<:[x_n, x_{n+1}]
+}}}
+Y para los valores negativos de x hacemos lo mismo pero usando el extremo superior del intervalo: {{{(+)}}}
+{{{
+        |  b_{q-1}*(x-x_q)
+        |    + b_q*x_q                    ,,  x<:[x_{q-1}, x_q]
+        |  b_{q-2}*(x-x_{q-1}) 
+        |    + b_{q-1}*(x_{q-1}-x_q) 
+        |    + b_q*x_q                    ,,  x<:[x_{q-2}, x_{q-1}]
+        |  ...
+        |  b_1*(x-x_2) 
+        |    + b_2*(x_2-x_3) 
+        |    + ... + b_q*x_q              ,,  x<:[x_1, x_2]
+}}}
+Podemos reagrupar los términos en función de los parámetros {{{b_k}}} del siguiente modo:
+{{{
+        |         / (x_k-x_{k+1})  ,,  x<x_k
+        |  b_k * |  (x-x_{k+1})    ,,  x<:[x_k, x_{k+1}]  ,,  k<q
+        |         \ 0              ,,  x>x_{k+1}
+        |         / x_q      ,,  x<x_q
+E(x) = <|  b_q * |  x        ,,  x<:[x_q, x_{q+1}]
+        |         \ x_{q+1}  ,,  x>x_{q+1}
+        |         / 0              ,,  x<x_m
+        |  b_m * |  (x-x_m)        ,,  x<:[x_m, x_{m+1}]  ,,  m>q
+        |         \ (x_{m+1}-x_m)  ,,  x>x_{m+1}
+}}} 
+A estos nuevos inputs los denominamos ''inputs-piecewise''.