[R-es] RMSE

genaro llusco gellusco en gmail.com
Mar Mayo 25 18:13:30 CEST 2010


Saludos:

bueno estoy trabajando con modelos FAVAR con series mensuales de N=5,
T=172, y bueno
1) determine el orden de rezago en la cual obtuve un rezago de 2 por
los criterios AIC y BIC, pero los errores no cumplen los supuestos, de
No-autocorrelacion u NO-heterocedasticidad,
pero claro es donde el error de pronostico es minimo, pero por otro
lado obtuve un modelo menos eficiente un VAR(12) basado simplemente en
el MSE y tambien mediante el test LR
 H_{0} :  A_{i}=0, pero con el VAR(12)(si cumple con los supuestos de
los residuos) se desminuye el sesgo asi para recompensar el
MSE(obviamente menos parsimonioso). Por otro lado, los PRMSE de ambos
modelos son aproximadamente iguales(te envio los datos calulados tipo
.txt).

CONSULTA:
1) si el error de pronostico se hace minimo en rezago 2, me parece
demasiado extender el rezago a orden 12.  Pero la verdad no se como
justificar el rezago de orden 12?
2) para hallar el RMSE se proyecta para 3 periodos adelante ej:se toma
datos hasta diciembre hacia atras para generar pronosticos de enero,
feb, marzo, luego para genera los siguientes pronosticos (feb,
marzo,abril) se actualizan las variables endogenas con los datos mas
recientes. Por que se pronostica para 3 periodos tambien puede ser 2
periodos o 4 periodos, es decir por que se toma tres periodos?
3) del enciso 2) cuantas de estos agrupaciones o proyeciones se debe hacer?
------------ próxima parte ------------
1)$selection
AIC(n)  HQ(n)  SC(n) FPE(n) 
     2      2      1      2 

$criteria
               1         2         3         4         5
AIC(n)  3.787968  3.527462  3.548320  3.559319  3.568965
HQ(n)   4.015307  3.944251  4.154559  4.355009  4.554104
SC(n)   4.348086  4.554346  5.041970  5.519735  5.996145
FPE(n) 44.173379 34.069986 34.857476 35.375205 35.935649
 

2) test de LR para H_{0}:A_{i}=0

      Lags    Log-Det    Chi^2      p-value
---------------------------------------------
 [1,]   14 -1.3613699 37.92518 4.704504e-02
 [2,]   13 -0.9254484 32.04680 1.566625e-01
 [3,]   12 -0.5771136 46.15732 6.153525e-03
 [4,]   11 -0.1012649 60.08856 1.016567e-04
 [5,]   10  0.4878386 44.15681 1.041534e-02
 [6,]    9  0.9005190 28.67643 2.776677e-01
 [7,]    8  1.1565586 21.54100 6.620897e-01
 [8,]    7  1.3406697 15.73518 9.224872e-01
 [9,]    6  1.4696466 37.73240 4.911767e-02
[10,]    5  1.7667521 39.79437 3.061412e-02
[11,]    4  2.0682246 48.31043 3.422998e-03
[12,]    3  2.4208554 40.15324 2.812448e-02
[13,]    2  2.7036247 91.65500 1.533389e-09
[14,]    1  3.3271281  0.00000 0.000000e+00


3)diagnost
------------------------------------------!-----------------------------
                        VAR(2)            !            VAR(12)
------------------------------------------!-----------------------------
            Statistics  g.l.    p-value   ! Statistics   g.l.   p-value
JB--VAR      946.52952   10   0.000000e+00!  1225.88757   10   0.0000000
ARCH--VAR   1245.00604 1125   6.963536e-03!  1097.64719   1125 0.7147465
Portmanteau 1151.74768  950   6.768695e-06!  746.14900    700  0.1102878
Skewness      24.71555    5   1.581043e-04!  96.09931     5    0.0000000
Kurtosis     921.81397    5   0.000000e+00!  1129.78826   5    0.0000000




4)result
 RMSE-MAPE     VAR(2)     VAR(2)     VAR(12)    VAR(12)
------------------------------------------------------------
Julio-Sept    2.3985697 0.020953702 2.6715388 0.023358440
Agos-Oct      2.3165999 0.020145472 1.6148832 0.013643679
Sept-Nov      1.8195352 0.015167225 1.8401584 0.015789247
Oct-Dic       0.8662240 0.007147114 1.5242311 0.012220599
Nov-Ene       0.8750749 0.007541200 0.9955535 0.006848410
Dic-Feb       0.5855205 0.004192831 0.6467009 0.003635713
Enero-Marzo   0.7298595 0.005223897 0.7037542 0.005024562
Febrero-Abril 1.0830057 0.008033597 1.3054747 0.010073539
Marzo-Mayo    1.4240308 0.010539730 1.2359335 0.008183120
Abril-Junio   0.9721559 0.007794409 0.6947928 0.005305666
------------------------------------------------------------
Promedio      1.3070576 0.010673918 1.3233021 0.010408297
               PRMSE     PMAPE       PRMSE        PMAPE


Más información sobre la lista de distribución R-help-es