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 (13 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 |