close
Warning:
Can't synchronize with repository "(default)" (/var/svn/tolp does not appear to be a Subversion repository.). Look in the Trac log for more information.
- Timestamp:
-
Jun 9, 2011, 6:15:40 PM (14 years ago)
- Author:
-
Víctor de Buen Remiro
- Comment:
-
--
Legend:
- Unmodified
- Added
- Removed
- Modified
-
v8
|
v9
|
|
26 | 26 | y las transformaciones de Box-Cox no son una aproximación razonable. |
27 | 27 | |
28 | | En este sentido, la forma ás natural de extensión de la clase de modelos es partir de la |
| 28 | En este sentido, la forma más natural de extensión de la clase de modelos es partir de la |
29 | 29 | [http://en.wikipedia.org/wiki/Generalized_linear_model regresión lineal generalizada] como base sobre |
30 | 30 | la que añadir el resto de características que se precisan para la modelación de problemas reales. |
… |
… |
|
33 | 33 | atisbar la clase general de modelos que queremos ser capaces de estimar: |
34 | 34 | |
35 | | * Distribución de los residuos |
| 35 | * '''Diseño por bloques de variables''': El mayor salto cualitativo de BSR sobre otros intentos anteriores |
| 36 | de simulación bayesiana es el diseño por bloques de Gibbs. Para poder acometer muchas de las mejoras |
| 37 | que se pretenden hay que ir a un enfoque de tipo ''Metropolis within Gibbs'' |
| 38 | * '''Diseño por segmentos de regresión''': Está claro que uno de los grandes logros de BSR es la capacidad para |
| 39 | definir diferentes segmentos de regresión, cada uno de los cuales tiene definida una distribución de |
| 40 | probabilidad propia independiente de la del resto de segmentos. Un segmento de regresión puede estar |
| 41 | ligado a un conjunto de datos de contraste, a una relación entre variables latentes y primarias o a |
| 42 | información a priori sobre cualquiera de ellas. |
| 43 | * '''Distribución de un segmento''': En BSR hasta ahora sólo es posible definir distribuciones normales |
| 44 | para cada uno de los segmentos pero eso podría cambiar en un futuro. |
36 | 45 | * caso normal |
37 | | * caso generalizado ([http://en.wikipedia.org/wiki/Generalized_linear_model#Probability_distribution overdispersed exponential family]). Enumero los casos que han sido solicitados alguna vez |
| 46 | * caso generalizado ([http://en.wikipedia.org/wiki/Generalized_linear_model#Probability_distribution overdispersed exponential family]). |
| 47 | Enumero los casos que han sido solicitados alguna vez |
38 | 48 | * regresión cualitativa |
39 | 49 | * regresión booleana |
… |
… |
|
49 | 59 | * distribuciones arbitrarias: definidas por el usuario, bien entre una gama parametrizable bien implementando la log-densidad. |
50 | 60 | * distribuciones distintas en cada segmento mediante bloques de Gibbs |
51 | | * Matriz de Varianzas-Covarianzas del ruido |
| 61 | * '''Matriz de Varianzas-Covarianzas''': aunque esto forma parte de la distribución es conveniente tratarla aparte pues |
| 62 | la teoría de regresión suele funcionar así, primero se parte del modelo incorrelado y homocedástico y luego se aplica la covarianza. |
52 | 63 | * Diagonal: Regresión ponderada o heterocedástica |
53 | 64 | * ARIMA: |
… |
… |
|
55 | 66 | * Caso generalizado |
56 | 67 | * Modelos multivariantes |
57 | | * Tratamiento de omitidos en los datos (input y output) |
| 68 | * '''Tratamiento de omitidos''' en los datos (input y output) |
58 | 69 | * Caso normal |
59 | 70 | * Caso generalizado |
… |
… |
|
61 | 72 | * paramétricos |
62 | 73 | * no paramétricos |
63 | | * Inclusión de restricciones para todas las variables del problema: |
| 74 | * Inclusión de '''restricciones''' para todas las variables del problema: |
64 | 75 | * de igualdad |
65 | | * de punto fijo |
66 | | * lineal |
67 | | * no lineal |
| 76 | * de punto fijo: (éstas ya están en BSR) [[BR]][[LatexEquation( \alpha_i = h_i )]] |
| 77 | * lineal: [[BR]][[LatexEquation( H \alpha = h )]] |
| 78 | * no lineal: [[BR]] [[LatexEquation( H\left( \alpha \right) = 0 )]] |
68 | 79 | * de desigualdad |
69 | | * de dominio |
70 | | * lineal |
71 | | * no lineal |
72 | | * Información a priori para todas las variables del problema: |
| 80 | * de dominio: [[BR]][[LatexEquation( l_i \le \alpha_i \le u_i )]] |
| 81 | * lineal: [[BR]][[LatexEquation( G \alpha = g )]] |
| 82 | * no lineal: [[BR]][[LatexEquation( G \left( \alpha \right) = 0 )]] |
| 83 | * '''Información a priori''' para todas las variables del problema: |
73 | 84 | * escalar |
74 | 85 | * normal |
… |
… |
|
76 | 87 | * otros |
77 | 88 | * vectorial |
| 89 | * variables no independientes |
78 | 90 | * modelos jerárquicos (paralelización) |
79 | 91 | * redes bayesianas |
… |
… |
|
84 | 96 | estimación, diagnosis, mantenimiento y explotación, y hacerlo de la forma más general posible: |
85 | 97 | |
86 | | * Chequeo previo: Comprobación de que el diseño del modelo es coherente con la clase establecida |
| 98 | * '''Chequeo previo''': Comprobación de que el diseño del modelo es coherente con la clase establecida |
87 | 99 | * Ausencia de colinealidades |
88 | 100 | * Existencia de soluciones factibles |
89 | 101 | * Compatibilidad de la región factible con la verosimilitud del modelo |
90 | | * Estimación: |
| 102 | * '''Estimación''': |
91 | 103 | * Simulación MCMC: Generación de cadenas de Markov de Montecarlo. |
92 | 104 | * División en bloques de Gibbs: Debe haber un proceso maestro (master) que se encargue de disparar |
… |
… |
|
118 | 130 | asintótica de forma que la salida sea análoga a la de la simulación MCMC para poder ser utilizada |
119 | 131 | por los métodos de diagnosis e inferencia baysianas. |
120 | | * Diagnosis: |
| 132 | * '''Diagnosis''': |
121 | 133 | * Chequeo de convergencia |
122 | 134 | * Raftery: Autocorrelación de la muestra |