Comprender los parámetros de la función de base gaussiana que se utilizarán en la regresión lineal


12

Me gustaría aplicar la función de base gaussiana en una implementación de regresión lineal. Lamentablemente, me está costando entender un par de parámetros en la función base. Específicamente y .μσ

Mi conjunto de datos es una matriz de 10,000 x 31. 10.000 muestras y 31 características. He leído que "Cada función base convierte el vector de entrada x en un valor escalar". Así que supongo que x es 1 muestra, por lo que es un vector 1 x 31. Desde aquí estoy confundido. ¿Qué es exactamente el parámetro ? He leído que esto gobierna las ubicaciones de las funciones básicas. Entonces, ¿no es esto el significado de algo? También estoy desilusionado por el subíndice j ( y ), esto me hace pensar en la fila j. Pero eso no parece tener sentido. ¿Es un vector? Ahora para laμjμϕμjσque "gobierna la escala espacial". ¿Qué es eso exactamente? He visto algunas implementaciones que prueban valores como .1, .5, 2.5 para este parámetro. ¿Cómo se calculan estos valores? He estado investigando y buscando ejemplos para aprender, pero hasta ahora no he podido encontrar ninguno. Cualquier ayuda o dirección es muy apreciada! Gracias.

Respuestas:


11

Como está confundido, permítame comenzar señalando el problema y respondiendo sus preguntas una por una. Tiene un tamaño de muestra de 10,000 y cada muestra se describe mediante un vector de características . Si desea realizar una regresión utilizando funciones de base radial gaussianas, entonces está buscando una función de la forma donde son sus funciones . Específicamente, necesita encontrar los pesos para que para los parámetros dados y minimice el error entre y la predicción correspondiente = f ( x ) = Σ j w j * g j ( x ; μ j , σ j ) , j = 1 .. m g i m w j μ j σ j y y f ( x )xR31

f(x)=jwjgj(x;μj,σj),j=1..m
gimwjμjσjyy^f(x^) - normalmente minimizará el error de mínimos cuadrados.

¿Qué es exactamente el parámetro Mu subíndice j?

Necesita encontrar funciones . (Aún debe determinar el número ) Cada función básica tendrá un y un (también desconocido). El subíndice varía de a .g j m μ j σ j j 1 mmgjmμjσjj1m

¿Es un vector?μj

Sí, es un punto en . En otras palabras, es un punto en algún lugar de su espacio de características y se debe determinar a para cada una de las funciones básicas de . μmR31μm

He leído que esto gobierna las ubicaciones de las funciones básicas. Entonces, ¿no es esto el significado de algo?

La función de base se centra en . Tendrá que decidir dónde están estas ubicaciones. Entonces, no, no es necesariamente el significado de nada (pero vea más abajo para conocer las formas de determinarlo) μ jjthμj

Ahora para la sigma que "gobierna la escala espacial". ¿Qué es eso exactamente?

σ es más fácil de entender si nos centramos en las funciones básicas.

Es útil pensar en las funciones de la base radial gaussiana en las dimensiones inferiores, digamos o . En la función de base radial gaussiana es solo la curva de campana conocida. La campana puede, por supuesto, ser estrecha o ancha. El ancho está determinado por : cuanto mayor sea más estrecha será la forma de la campana. En otras palabras, escala el ancho de la forma de la campana. Entonces, para = 1 no tenemos escala. Para grandes tenemos una escala sustancial.R 2 R 1 σσσσσR1R2R1σσσσσ

Puede preguntar cuál es el propósito de esto. Si piensa en la campana que cubre una parte del espacio (una línea en ), una campana estrecha solo cubrirá una pequeña parte de la línea *. Los puntos cerca del centro de la campana tendrán un valor mayor de . Los puntos alejados del centro tendrán un valor menor de . El escalado tiene el efecto de empujar los puntos más lejos del centro, ya que la campana estrecha los puntos se ubicarán más lejos del centro, reduciendo el valor de x g j (x) g j (x) g j (x)R1xgj(x)gj(x)gj(x)

Cada función base convierte el vector de entrada x en un valor escalar

Sí, está evaluando las funciones básicas en algún momento .xR31

exp(xμj222σj2)

Obtienes un escalar como resultado. El resultado escalar depende de la distancia del punto desde el centro dado pory el escalar .xx - μ jσ jμjxμjσj

He visto algunas implementaciones que prueban valores como .1, .5, 2.5 para este parámetro. ¿Cómo se calculan estos valores?

Por supuesto, este es uno de los aspectos interesantes y difíciles del uso de funciones de base radial gaussianas. si busca en la web, encontrará muchas sugerencias sobre cómo se determinan estos parámetros. Esbozaré en términos muy simples una posibilidad basada en la agrupación. Puede encontrar esta y otras sugerencias en línea.

Comience agrupando sus 10000 muestras (primero puede usar PCA para reducir las dimensiones seguidas de la agrupación de k-medias). Puede dejar que sea ​​el número de clústeres que encuentre (generalmente, empleando validación cruzada para determinar el mejor ). Ahora, cree una función de base radial para cada grupo. Para cada función de base radial, sea el centro (p. Ej., La media, el centroide, etc.) del grupo. Deje que refleje el ancho del clúster (p. Ej., Radio ...) Ahora continúe y realice su regresión (esta descripción simple es solo una descripción general: ¡necesita mucho trabajo en cada paso!)m g j μ j σ jmmgjμjσj

* Por supuesto, la curva de campana se define de - a por lo que tendrá un valor en todas partes en la línea. Sin embargo, los valores lejos del centro son insignificantes.


¡Buena respuesta! Sin embargo, al buscar , ¿no terminamos con la regresión de máquina de vectores de soporte (con núcleo gaussiano)? μ
O_Devinyak

@ O_Devinyak: muchos métodos de expansión base requerirán algún tipo de estimación de parámetros. Hay muchas formas de encontrar así que no creo que esto necesariamente signifique que estamos reduciendo el problema a SVR. Para ser honesto, no soy un experto en SVR, pero la función de pérdida que se minimiza es ciertamente diferente y estoy seguro de que muchas de las características se ignoran, esa es la forma de Vector de soporte. Con las funciones básicas usamos todas las funciones para la evaluación, pero afortunadamente el soporte compacto significa que muchas de las funciones básicas devuelven valores insignificantes o cero. De todos modos, sería una buena pregunta en este foroμ
martino

¿Por qué necesitamos una escala lugar de una matriz de covarianza que haga que la función base se vea como la parte exponencial de un gaussiano multivariado? σj
stackunderflow

1

Déjame intentar dar una explicación simple. En dicha notación, puede ser un número de fila, pero también puede ser un número de característica. Si escribimos entonces denota el número de característica, es columna-vector, es escalar y es una columna -vector. Si escribimos entonces denota el número de fila, es escalar, es vector de columna y es un vector de fila. La notación donde denota fila y denota columna es más común, así que usemos la primera variante.jy=β0+j=1:31βjϕj(x)jyβjϕj(x)yj=βϕj(x)jyjβϕj(x)ij

Al introducir la función de base gaussiana en la regresión lineal, (escalar) ahora no depende de los valores numéricos de las (vector), sino de las distancias entre y el centro de todos los demás puntos . En tal modo no depende de si valor de la característica -ésimo de -ésima observación es alta o pequeña, pero depende de si valor de la característica-ésimo está cerca o lejos de la media para que -feature . Entonces no es un parámetro, ya que no se puede ajustar. Es solo una propiedad de un conjunto de datos. El parámetrox i x i μ i y i j i j j μ i j μ j σ 2 y y σ 2yixixiμiyijijjμijμjσ2es un valor escalar, controla la suavidad y se puede ajustar. Si es pequeño, los pequeños cambios en la distancia tendrán un gran efecto (recuerde el fuerte gaussiano: todos los puntos ubicados ya a poca distancia del centro tienen pequeños valores ). Si es grande, los pequeños cambios en la distancia tendrán un efecto bajo (recuerde gaussiano plano: la disminución de con el aumento de la distancia desde el centro es lenta). Se debe buscar el valor óptimo de (generalmente se encuentra con validación cruzada).yyσ2


0

Las funciones de base gaussianas en las configuraciones multivariadas tienen centros multivariados. Suponiendo que su , entonces también. El gaussiano tiene que ser multivariado, es decir, donde es Una matriz de covarianza. El índice no es un componente de un vector, es solo el vector . Del mismo modo, es la ésima matriz. μ jR 31 e ( x - μ j ) Σ - 1 j ( x - μ j ) Σ jR 31 × 31 j j Σ j jxR31μjR31e(xμj)Σj1(xμj)ΣjR31×31jjΣjj

Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.