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:
-
Nov 8, 2012, 10:50:47 AM (12 years ago)
- Author:
-
Víctor de Buen Remiro
- Comment:
-
--
Legend:
- Unmodified
- Added
- Removed
- Modified
-
v4
|
v5
|
|
94 | 94 | |
95 | 95 | == Identificación de los grados ARMA == |
| 96 | |
| 97 | Cada par de polinomios regulares AR y MA lleva asociado una ACF teórica única que es posible comparar |
| 98 | con las ACF muestrales para calcular la verosimilitud de un modelo ARMA regular con respecto a esas |
| 99 | ACF muestrales. Integrando dicha verosimilitud por un prior no informativo o el que se considere |
| 100 | oportuno, en la región de estacionariedad asociada al par de grados {{{(p,q)}}} que definen una estrutura |
| 101 | ARMA regular, se puede calcular la probabilidad a posteriori de dicha estructura. |
| 102 | |
| 103 | Dentro del paquete ArimaTools se ofrecen distintos métodos para el cálculo de esta integral, todos |
| 104 | ellso basados en simulación de modelos ARMA estacionarios. |
| 105 | |
| 106 | Para cada polinomio de diferencias con probabilidad significativa se tomará la ACF muestral del ruido |
| 107 | diferenciado, y para cada factor regular o estacional, se extraerán las autocorrelaciones de órdenes |
| 108 | múltiplos de la periodicidad {{{s}}} correspondiente que no sean satélites de las periodicidades superiores |
| 109 | para minimizar la contaminación entre los distintos factores. Para cada uno de estos factores se |
| 110 | calculará la probabilidad de cada par de grados {{{(p*s,q*s)}}} como si se tratara de un modelo regular |
| 111 | {{{(p,q)}}}. |
| 112 | |
| 113 | Finalmente se obtendrán las probabilidades del producto cartesiano de los pares {{{(p*s,q*s)}}} |
| 114 | correspondientes a cada polinomio de diferencias multiplicando las probabilidades de dichos pares |
| 115 | por la probabilidad de las raíces unitarias calculada en el apartado anterior. De esta forma |
| 116 | se obtiene un ''ranking'' de modelos ARIMA ordenados por probabilidad de mayor a menor. |
| 117 | |
96 | 118 | {{{ |
97 | 119 | #!java |
… |
… |
|
117 | 139 | |
118 | 140 | == Estimación y diagnosis de la estructura ARIMA completa == |
| 141 | El último paso es estimar los modelos ordenados por probabilidad sobre la serie |
| 142 | homocedástica original incorporando los inputs de usuario y el AIA encontrados en |
| 143 | el filtro inicial. Una vez estimados hay que efectuar la diagnosis de calidad del |
| 144 | mismo que consiste en comprobar qeu efectivamente cumple los requisitos de un modelo |
| 145 | ARIMA |
| 146 | * hipótesis sobre los residuos |
| 147 | * siguen una distribución normal de media nula |
| 148 | * son independendientes entre sí |
| 149 | * hipótesis de los parámetros |
| 150 | * deben ser significativos |
| 151 | * no pueden ser colineales dos a dos ni estar demasiado cerca (correlación alta) |
| 152 | * no pueden ser multicolineales |
| 153 | * los polinomios AR y MA resultantes deben tener todas las raíces fuera del círculo unidad |
| 154 | |
| 155 | Esto se traduce en una batería de tests para los que se obtiene la potencia del mismo |
| 156 | que es siempre un número entre 0 y 1. Para cada uno se establecen dos límites, aceptación |
| 157 | y rechazo, mediante los cuales se califica el test con el siguiente criterio |
| 158 | * Si la potencia es menor que el límite de aceptación se clasifica el test como bueno. |
| 159 | * Si se sitúa entre el límite de aceptación y el de rechazo el tets es aceptable o regular. |
| 160 | * Si sobrepasa el límite de rechazo se dice que el resultado del test es malo. |
| 161 | El modelo se calificará como el peor de los tests de la batería de diagnosis. |
| 162 | Entre los modelos de igual calificación se utilizará el criterio de Swartz para establecer |
| 163 | las prioridades de ordenación. |
| 164 | |
119 | 165 | {{{ |
120 | 166 | #!java |