Resumen1. INTRODUCCIÓNEsta nota técnica propone una forma alternativa de tratar el componente estacional presente en algunas series de tiempo que forman parte de los modelos econométricos. Sugiere incorporar variables instrumentales (dummy) estacionales para capturar este fenómeno, como una opción más apropiada en la estimación de modelos que utilizan la usual transformación de las variables para eliminar ese fenómeno, tal como series desestacionalizadas, en tendencia-ciclo, etc.
Además, se describen los principales comandos en los paquetes econométricos SHAZAM y EVIEWS necesarios para incluir esas variables ficticias y ejemplifica su uso con la estimación de un modelo de demanda monetaria para Costa Rica.
También, compara los resultados del enfoque tradicional (modelos con variables desestacionalizadas) con el propuesto (modelos con variables estacionales).
Concluye, en línea con las nuevas tendencias en el área de la econometría, la conveniencia técnica de emplear variables dummy estacionales por sobre el uso de variables transformadas en los modelos de regresión, debido a que brindan resultados más atractivos, abundantes y valiosos para el análisis.
El estudio de ciertos fenómenos económicos en periodos subanuales suele dificultarse cuando las observaciones cuantitativas sobre éstos están dominadas por una marcada estacionalidad. Esta se refiere a un comportamiento o patrón en las cifras que se repite cada año, con cierta regularidad, pero que también puede evolucionar, y que es causado generalmente por fuerzas no económicas, de tipo exógeno, que los tomadores de decisión no pueden controlar.
Si se deja este componente estacional en los modelos de regresión, se distorsionarían las relaciones de largo plazo y con ello la estimación de los coeficientes de los parámetros. Para reducir este problema el investigador tiende a obviar esta situación, que le impide observar en forma clara y definida el comportamiento de largo plazo del fenómeno, utilizando transformaciones estadísticas de los datos. Esto posiblemente influido por la idea de que la variación estacional es simplemente un tipo de "ruido" que contamina los datos económicos y no una parte fundamental de la serie que debería tratar de explicarse1/. Al parecer, este enfoque ha propiciado la utilización de relaciones funcionales con variables no directamente observables, tales como aquellas desestacionalizadas e incluso en tendencia-ciclo.
No obstante lo anterior, las nuevas tendencias de la teoría econométrica han abandonado esa práctica para concentrarse en relaciones funcionales que utilizan variables observables, grupo al que pertenecen las estacionales. Algunas razones por las cuales se recomienda descartar el uso de ese tipo de transformaciones, son las siguientes:
Un aspecto importante es que las variables estacionales se han utilizado en diversas investigaciones del FMI, por ejemplo en los modelos VAR y en el Banco de España en los modelos de ecuaciones de demanda de dinero.8/
En este trabajo se presenta este nuevo enfoque econométrico, mediante el cual se incorporan variables instrumentales estacionales en los modelos de regresión. Su contenido es el siguiente: en la segunda sección se explica el marco conceptual de las variables estacionales, en la tercera se mencionan los principales comandos de SHAZAM y EVIEWS necesarios para aplicarlas, en la cuarta se presenta una ilustración de su aplicación a una ecuación de demanda de dinero seleccionada, en la quinta se comparan los resultados de un modelo que utiliza variables desestacionalizadas con uno con variables estacionales y, por último, se incluyen algunas observaciones finales.
2. EL CONCEPTO DE LAS VARIABLES ESTACIONALES
Las variables dummy para ajuste estacional son variables artificiales que asumen valores discretos, generalmente de 0 y 1. Estas fueron originalmente aplicadas por Lovell a inicios de los años 60 y sirven para "explicar" la estacionalidad en las series de tiempo, la cual, como se adelantó, es un patrón de comportamiento regular de una serie a lo largo de cada año, que obedecería a factores tales como costumbres, días festivos decretados, vacaciones de verano, época de navidad y otros factores similares que ocasionan incrementos o disminuciones en las magnitudes de ciertas variables, como por ejemplo la producción, las ventas, etc.
Si se trabaja con datos trimestrales, la primera idea que surge es la de utilizar variables artificiales; una para cada trimestre, que pueden definirse como: q1, q2, q3 y q4; cuya representación para dos años cualesquiera sería:

La inclusión de los coeficientes de estas variables y de la constante en un modelo de regresión produciría una matriz bianual X de la siguiente forma:

No obstante que lo anterior parece natural, las columnas correspondientes a las variables estacionales darían lugar a una combinación lineal exacta con la constante,9/ lo cual produciría que el determinante de la matriz X'X fuera igual a cero y, por tanto, singular (no invertible), lo cual no permitiría estimar los coeficientes del modelo de regresión.
Para evitar este inconveniente deben utilizarse únicamente tres de las cuatro variables dummy y por supuesto la constante. Así, por ejemplo, se puede excluir q4 en la matriz X y en este caso esa variable omitida estaría implícitamente recogida con la columna de la constante. En efecto, al utilizar sólo esas tres variables estacionales (junto con la constante y las restantes variables que sean de interés) en la matriz, se evitaría la colinealidad perfecta con la columna de la constante asociada con el intercepto. Esto por cuanto los cuatro trimestres estarían indicados con sólo las tres dummy: q1, q2, q3 y la constante, con lo cual la matriz sería:

Otra forma muy utilizada recientemente consiste en expresar las variables artificiales estacionales como desviaciones con respecto a la que corresponde al cuarto trimestre. Estas nuevas variables, que podrían denominarse S1, S2 y S3, corresponderían a las siguientes diferencias vectoriales:
S1 = q1 - q4Una vez efectuadas las operaciones anteriores e incorporado el vector de la constante, la nueva matriz X queda definida de la siguiente manera:S2 = q2 - q4
S3 = q3 - q4

Como se observa en la matriz anterior, los vectores de las variables dummy estacionales han sido definidos de forma tal que su suma sea cero en cada año, por lo que este sistema permite que el efecto estacional se anule en el año y que se obvie el problema de singularidad de la matriz.
A manera de ejemplo, considérese un modelo de regresión con cifras trimestrales, en donde la variable Y depende de la variable Z y en el que se incorporan tres variables dummy trimestrales (Si, para todo i = 1, 2, 3) y un término de error (e). Este modelo estaría representado de la siguiente manera:
Es importante mencionar que en el caso de variables con periodicidad mensual, se crearían únicamente once variables estacionales, en forma equivalente a lo explicado en esta sección. Sin embargo, en este caso se presenta el inconveniente de que se requiere gran cantidad de observaciones.
Por otra parte, no se debe dejar de lado que el uso de las variables estacionales presenta problemas cuando la estacionalidad es móvil, es decir, cuando varía año con año. En este caso, es difícil que capturen de una forma adecuada la estacionalidad.10/
3. PROCEDIMIENTO EN SHAZAM Y EVIEWS
Si bien las variables estacionales se pueden generar en una hoja electrónica, o pueden ser creadas manualmente en el archivo de trabajo, algunos paquetes econométricos incorporan comandos (automáticos) especiales para estos efectos.
El paquete econométrico SHAZAM versión 7.0 tiene el comando SEAS(X), el cual permite crear una variable dummy estacional con periodicidad X.
Por ejemplo, para generar variables dummy estacionales trimestrales, se puede utilizar el siguiente conjunto de comandos:
MATRIX QD=SEAS(#,4)Donde el signo # corresponde al número de observaciones de la muestra con que se está trabajando. Como se observa, y de acuerdo con lo mencionado en la sección anterior, solamente se crean tres variables estacionales.
GENR QD1=QD:1
GENR QD2=QD:2
GENR QD3=QD:3
En EVIEWS versión 2.0 existe el siguiente comando, el cual forma parte de los argumentos de la ecuación de regresión: @SEAS(d). El comando anterior permite generar una dummy estacional igual a uno cuando el trimestre o mes es igual a d y cero en otro caso.
Por otra parte, una ventaja que presenta SHAZAM con respecto a EVIEWS es que permite incorporar las variables estacionales con el procedimiento stepwise, de una forma muy fácil, así:
4. APLICACIÓN
En esta sección se presenta un ejemplo utilizando variables dummy estacionales q1, q2 y q3. Luego se repite el ejercicio empleando S1, S2 y S3 y finalmente se comparan ambos resultados.
Como aplicación del uso de las variables estacionales en el análisis de regresión, se optó por utilizar una de las ecuaciones de demanda de dinero trimestral, estimada recientemente y que produjo resultados satisfactorios. En esa ecuación se trató de una forma un tanto diferente la presencia de estacionalidad en el modelo de regresión.11/
Para llevar a cabo la aplicación, se utilizó el paquete EVIEWS. La nueva especificación se ejecutó de la siguiente manera:

donde:
LNEMIR: Logaritmo natural de los promedios geométricos trimestrales de la emisión monetaria nominal mensual, en saldos a fin de mes, deflactada por el IPC base enero 1995=100. C: ConstanteEl periodo analizado comprende desde el segundo trimestre de 1987 al cuarto de 1997 (43 observaciones).12/ En el siguiente gráfico se observa la marcada estacionalidad de la emisión monetaria y la necesidad de incorporar variables estacionales en el análisis de regresión.LNPIBR: Logaritmo natural del PIB real trimestral original, a precios de productor, año base 1966=100.
IPA3: Tasa de interés pasiva bruta anual en colones, pagada por los bancos comerciales, a 3 meses plazo.
LNEMIR(-1): Variable dependiente (emisión) rezagada un periodo.
q1, q2, q3 y q4: Variables estacionales.

Cuando se ajustó esa regresión se produjo un mensaje de error, tal como se indica en la siguiente pantalla, el cual se originó, como se explicó anteriormente, en la multicolinealidad perfecta de las cuatro variables estacionales con la constante de regresión, lo que impediría estimar los coeficientes del modelo:

Posteriormente, se incorporaron a la ecuación de regresión sólo las tres primeras variables estacionales (q1, q2 y q3), de la siguiente forma:

Una vez eliminada la singularidad de la matriz de coeficientes, se obtuvieron los siguientes resultados principales:
| Cuadro 1 AJUSTE DE REGRESIÓN CON VARIABLES ESTACIONALES q1, q2, q3 VARIABLE DEPENDIENTE: LNEMIR | ||
| Nombre de la variable | Coeficiente | Estadístico t |
| C | 0,8963 | 2,72 |
| LNPIBR | 0,2436 | 5,40 |
| IPA3 | -0,0039 | -5,60 |
| LNEMIR(-1) | 0,7596 | 14,22 |
| Q1 | -0,1909 | -18,28 |
| Q2 | -0,1704 | -17,41 |
| Q3 | -0,1591 | -17,87 |
| R2 ajustado = 0,9787 Durbin-h = 1,350313/ | ||
Los coeficientes de las variables estacionales presentan signos negativos, que acumulan entre los tres -0,52, según se puede apreciar en el cuadro anterior. Esto indica que la demanda de emisión es en realidad menor de lo que sería sin estacionalidad (vía fundamentos o tendencia).14/
En este caso, q1, q2 y q3 tienden a compensar a q4; en donde la estacionalidad es más alta. Para estimar el coeficiente de q4 se debe restar de la unidad el acumulado de -0,52, lo que daría un coeficiente de 1,52.
Por otra parte, una forma adicional de "detectar" la necesidad de incorporar variables estacionales consiste en realizar el ajuste sin incorporarlas, observar el gráfico de los residuos y, si en estos aparece un "marcado" comportamiento estacional, entonces se decide agregarlas.
A continuación se incluyen dos gráficos con los residuos "antes" de incorporar las variables estacionales y otro con los residuos "después" de agregarlas.
Gráfico 3
RESIDUOS DEL MODELO CON VARIABLES DUMMY ESTACIONALES
También se hizo la aplicación utilizando las variables estacionales en diferencias (S1, S2 y S3) como se muestra a continuación:

Los resultados de esta nueva corrida se muestran en el siguiente arreglo estadístico:
| Cuadro 2 AJUSTE DE REGRESIÓN CON LAS VARIABLES ESTACIONALES S1, S2, S3 VARIABLE DEPENDIENTE: LNEMIR | ||
| Nombre de la variable | Coeficiente | Estadístico t |
| C | 0,8963 | 2,72 |
| LNPIBR | 0,2436 | 5,40 |
| IPA3 | -0,0039 | -5,60 |
| LNEMIR(-1) | 0,7596 | 14,22 |
| Q1 | -0,1909 | -18,28 |
| Q2 | -0,1704 | -17,41 |
| Q3 | -0,1591 | -17,87 |
| R2 ajustado = 0,9787 Durbin-h = 1,350315/ | ||
En resumen, la aplicación de las variables estacionales es muy sencilla. Además, las variables resultaron estadísticamente significativas, con los signos teóricos esperados y ajustes satisfactorios.
5. COMPARACIÓN DE RESULTADOS
Se consideró importante comparar los modelos con dummy estacionales y los que utilizan variables desestacionalizadas. Para ello, se siguió trabajando con la misma aplicación de la sección anterior. En el cuadro 3 se incluyen nuevamente los resultados del cuadro 1, correspondientes al caso en que se usan variables estacionales, así como la evidencia para el caso en que se emplea la variable dependiente desestacionalizada (LNEMIRD), que es la forma que tradicionalmente se utiliza en algunos de los estudios económicos.
| Cuadro 3 DEMANDA DE EMISIÓN MONETARIA PARA COSTA RICA: RESULTADOS COMPARATIVOS DE LOS MODELOS CON DUMMY ESTACIONALES Y CON VARIABLES DESESTACIONALIZADAS 1987:2 - 1997:4 | ||||
| Variables | Variable dependiente LNEMIR | Variable dependiente LNEMIRD | ||
| Coeficiente | Estadístico t | Coeficiente | Estadístico t | |
| C | 0,8963 | 2,72 | 0,6495 | 2,00 |
| LNPIBR | 0,2436 | 5,40 | 0,1998 | 4,96 |
| IPA3 | -0,0039 | -5,60 | -0,0037 | -5,39 |
| LNEMIR(-1) | 0,7596 | 14,22 | --- | --- |
| LNEMIRD(-1) | --- | --- | 0,8020 | 16,31 |
| q1 | -0,1909 | -18,28 | --- | --- |
| q2 | -0,1704 | -17,41 | --- | --- |
| q3 | -0,1591 | -17,87 | --- | --- |
| R2 aj =0,9787 | R2 aj = 0,9732 | |||
| Durbin-h = 1,3503 | Durbin-h =0,8981 | |||
Los resultados estadísticos no son esencialmente diferentes, como se observa en el cuadro anterior. En particular hay igualdad en los signos y analogía en magnitud y significancia de los coeficientes de regresión. Por otra parte, se decidió analizar las series ajustadas por regresión de ambas funciones. Con respecto a la evidencia del cuadro 3, la Y ajustada (emisión real), se obtiene en términos desestacionalizados por lo que fue necesario aplicarle el factor de estacionalidad para regresar a los valores originales,16/ mientras que los resultados del cuadro 1 se obtienen directamente en forma original. Posteriormente, se analizaron ambas series, llegándose a la conclusión de que casi no difieren. A continuación se presenta un gráfico con ambas series:

Resalta el alejamiento que muestran ambas líneas en el año 92, el cual probablemente se deba a que la apertura de la cuenta de capitales que se dio a partir de ese periodo, pudo tener algún impacto sobre el fenómeno estacional, produciendo un cambio en el mismo. Sin embargo, la serie ajustada por ambos métodos presenta consistentemente un comportamiento similar.
Para confirmar la afirmación anterior, se realizó una prueba estadística denominada la prueba de las observaciones aparejadas. Consiste en detectar si existe diferencia significativa entre las variaciones porcentuales de la emisión monetaria ajustada según el método de las variables desestacionalizadas y las dummy estacionales.
Se compara el t calculado (tc) con el t tabular (tt)17/. En este caso, el t calculado es igual a -0,62 y el t tabulado con un nivel de significancia del 5% (a = 5%) y 42 grados de libertad es de 2,02 y como | tc| < | tt | se rechaza la hipótesis de que existe diferencia significativa entre ambos métodos con un nivel de significancia del 5%.18/
En síntesis, se nota el proceso más largo que se requiere cuando se trabaja con series desestacionalizadas en una regresión, en tanto que si se usan las dummy estacionales no es necesario realizar este procedimiento.
6. CONSIDERACIONES FINALES
No obstante, las nuevas orientaciones de la teoría econométrica indican que se comete un error al utilizar variables desestacionalizadas o en tendencia-ciclo en las ecuaciones de regresión. Debido a que los coeficientes de los parámetros determinarían una recta diferente a la que se obtendría con las variables estacionales, que es la que se considera como la mejor estimación.
Agénor, Pierre-Richard; Hoffmaister, Alexander y Medeiros, Carlos (1997) "Cyclical Fluctuations in Brazil's Real Exchange Rate: The Role of Domestic and External Factors". Research Department. Fondo Monetario Internacional. Cabrero, A. y Delrieu, J.C. (1996) "Elaboración de un índice sintético para predecir la inflación en España". Documento de trabajo No. 9619. Banco de España. Davidson, R. y MacKinnon, J. (1993) "Estimation and Inference in Econometrics". Oxford University Press.
Granger, C.W.J. (1990) "Modelling Economic Series". Advanced Texts in Econometrics. Oxford University Press.
Hylleberg, Svend (1992) "Modelling Seasonality". Advanced Texts in Econometrics. Oxford University Press.
Maravall, Agustín (1994). "Use and Misuse of Unobseved Components in Economic Forecasting", Journal of Forecasting, Vol. 13, págs. 157-178.
Sánchez C., Carmen (1995) "Ecuación de demanda de ALP2: Propiedades de corto plazo y utilización en el análisis monetario". EC/1995/33. Servicio de Estudios. Banco de España.
Solano, Ivannia y Torres, Carlos (1998) "Estimación de funciones de demanda de dinero trimestrales y anuales". Mimeo. Banco Central de Costa Rica.
Direcciones consultadas en INTERNET:
http://www.ech.bris.ac.uk/compute/shazam/intro/dummy.htmhttp://shazam.econ.ubc.ca/intro/dumseas.htm


Anexo 3
PRUEBA DE LAS OBSERVACIONES APAREJADAS

donde:
tc = estadístico t calculado.Se compara con t tabular:
d = promedio de las diferencias de las variaciones porcentuales de ambos métodos.Sd = desviación estándar de las diferencias porcentuales de ambos métodos.
n = número de diferencias de las variaciones porcentuales de ambos métodos.
donde:
a = nivel de significancia.

1/ Por medio de los coeficientes de variables dummy y no con el uso de variables explicativas económicas.
2/ Véase, Davidson y Mackinnon (1993), pág. 688.
3/ Véase, Cabrero y Delrieu (1996), pág. 9.
4/ Véase, Davidson y Mackinnon (1993), pág. 688.
5/ Véase Davidson y Mackinnon (1993), pág. 696.
6/ Véase, Cabrero y Delrieu (1996), pág. 8 y Maravall (1994), citado por estos autores.
7/ Véase, Davidson y MacKinnon (1993), pág. 687.
8/ Véase, Agénor, Hoffmaister y Medeiros (1997) y Sánchez C., Carmen (1995). 9/ Se trata de un caso de multicolinealidad perfecta.
10/ No obstante, aún en este caso se sugiere utilizarlas adicionando variables dummy estacionales que han sido interactuadas con potencias de una tendencia de tiempo lineal que se incrementa anualmente o usar métodos de dominio de frecuencias. Profundizar en este tema se sale de los objetivos de esta nota técnica; una referencia detallada se podría consultar en Davidson y MacKinnon (1994), pág. 189-190.
11/ Para mayor detalle, véase Solano y Torres (1998).
12/ En el Anexo 1 se incluye la base de datos utilizada en este trabajo.
13/ El Durbin-h estaría indicando que no se presentan problemas de autocorrelación. Asimismo, se detecta multicolinealidad entre la variable rezagada y las estacionales.
14/ Si no se rechaza la hipótesis de que, estadísticamente, los coeficientes son iguales a cero, no habría estacionalidad y la variable dependiente estaría de acuerdo con los fundamentos o la tendencia.
15/ El Durbin-h estaría indicando que no se presentan problemas de autocorrelación. Asimismo, se detecta multicolinealidad entre la variable rezagada y las estacionales.
16/ En el Anexo 2 se incluyen estos cálculos.
17/ Véase Anexo 3.
18/ En el Anexo 4 se encuentran los resultados de esta prueba. Kikut y Torres: "Variables Estacionales en los Modelos de Regresión: Una Aplicación a la Demanda de Dinero en Costa Rica"