SciELO - Scientific Electronic Library Online

 
vol.79 número314Tasa impositiva única y cumplimiento fiscalLos efectos del incremento del salario mínimo interprofesional en el empleo de los trabajadores inmigrantes en España índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay artículos similaresSimilares en SciELO

Compartir


El trimestre económico

versión On-line ISSN 2448-718Xversión impresa ISSN 0041-3011

El trimestre econ vol.79 no.314 Ciudad de México abr./jun. 2012

 

Artículos

Estacionariedad, cambios estructurales y crecimiento económico en México (1895-2008) *

Antonio E. Noriega** 

Cid Alonso Rodríguez-Pérez*** 

** Banco de México y Universidad de Guanajuato, México.

*** Banco de México, México.


Resumen:

Este artículo estima que la evolución del producto real y per capita de México entre 1895 y 2008 se puede describir adecuadamente por medio de un modelo estacionario en tendencia, afectado por cuatro cambios estructurales, cuya ocurrencia coincide con cambios institucionales internos, guerras y crisis económicas y financieras. Estos cambios son modelados por funciones logísticas de transición suave, en las que los periodos de transición son estimados endógenamente. En términos de tasas de crecimiento, nuestros resultados indican que para el PIB real y PIB real per capita de México existen cuatro etapas de crecimiento estacionario separadas por tres periodos de transición. Por ejemplo, para el PIB real se identifican las siguientes etapas de crecimiento estacionario: 1895-1924, 1935-1952, 1956-1978 y 1989-2008, separadas por tres épocas de transición (1925-1934, 1953-1955 y 1979-1988).

Palabras clave: producto interno bruto; crecimiento económico; estacionariedad; raíz unitaria; cambio estructural; función logística; transición suave

Abstract:

This paper shows that the evolution of the level of Mexico real and real per capita output between 1895 and 2008 can be adequately described through a trend stationary model, affected by 4 structural breaks, which ocurred at dates that seem to coincide with domestic institutional arrangements, wars, and financial and economic crises. These changes are modeled through logistic smooth transition functions, in which transition periods are endogenously estimated. In terms of growth rates, our results indicate that for Mexico real and real per capita GDP, there are four stationary growth paths, separated by three transition periods. For instance, for real GDP we identify the following stationary growth paths: 1895-1924, 1935-1952, 1956-1978, and 1989-2008, separated by three transition periods: 1925-1934, 1953-1955, and 1979-1988.

Clasificación JEL: C12, C13, C15, C22, O47, O54.

Introducción

Este trabajo tiene por objetivo analizar el comportamiento de baja frecuencia del producto real y per capita de México de 1895 a 2008. En el contexto del análisis univariado de series de tiempo se busca determinar si la evolución de la producción en México se puede caracterizar por una tendencia estocástica -esto es, por una raíz unitaria- o por una tendencia lineal determinista, sujeta a cambios estructurales modelados por medio de funciones de transición suave. En otras palabras, este estudio retoma la vieja pregunta de si la tendencia que domina la evolución de largo plazo del producto es estocástica o determinista.

Considérese el siguiente modelo de tendencia lineal sin cambios estructurales en los parámetros,

en los que yt es, digamos, el logaritmo del PIB real en niveles, vt son las desviaciones respecto a la tendencia lineal, α, β y ρ son parámetros, y εt denota un proceso de ruido blanco. Si ρ=1 entonces yt tiene una raíz unitaria o en otras palabras es estacionario en diferencias, por lo que para obtener la tasa de crecimiento del producto sería necesario aplicar la transformación de primeras diferencias y obtener una estimación del parámetro β. Por lo contrario, en el caso de ρ ˂ 1, yt es estacionario en tendencia y, por tanto, la pendiente de la función de tendencia lineal de la serie (β) representa una estimación de la tasa de crecimiento de yt .

Considérese la evolución del producto real de los Estados Unidos, que se presenta en la parte A de la Gráfica 1. Resultados en Noriega (2004), así como cálculos propios utilizando las pruebas de Kejriwal y Perron (2010), indican que el producto real estadunidense ha fluctuado de manera estacionaria alrededor de una tendencia lineal sin cambios estructurales. En otras palabras, el modelo de arriba con parámetros α y β constantes y ρ ˂ 1 parecen describir adecuadamente la evolución de la serie.1

GRÁFICA 1 PIB real (logs) y tendencias quebradas 

Sin embargo, esta representación lineal puede no ser adecuada para todos los casos, debido a que la función de tendencia puede estar sujeta a cambios estructurales, que usualmente están relacionados con guerras, arreglos institucionales, crisis económicas, etc. Como ejemplo de esto la Gráfica 1 presenta la evolución del producto real del Reino Unido (parte B) y de Australia (parte C), junto con una tendencia lineal quebrada, basada en los resultados de Noriega, Soria y Velázquez (2008). En el primer caso se identifica un cambio estructural tanto en α, como en la pendiente β de la función de tendencia en 1918, que alteran permanentemente el nivel y la tasa de crecimiento del producto. Este hallazgo ha sido relacionado por la bibliografía con la primera Guerra Mundial (véase por ejemplo Mills, 1991, y Duck, 1992). En la parte C, el PIB real de Australia presenta tres cambios estructurales en su nivel (tres cambios en el parámetro sin cambios en el parámetro de la pendiente β). En ambos casos, pruebas de raíz unitaria aplicadas a las series indican que el producto evoluciona de una manera estacionaria, es decir que ρ ˂1 una vez que se consideran estos cambios estructurales en la función de tendencia.2

En este artículo, ponemos particular atención a la cuestión de si los parámetros de la función de tendencia (α y β) se han mantenido constantes a lo largo del periodo en estudio, tanto para el producto como para el producto per capita de México.

Existen diversos estudios internacionales que han explorado este cuestionamiento. Por ejemplo, Ben-David y Papell (1995) plantean cómo las pruebas de cambio estructural en series de tiempo univariadas proporcionan un marco natural para investigar el hecho (estilizado) de una tasa de crecimiento del producto constante. Así, Ben-David y Papell (1995) buscan determinar la existencia de un cambio estructural en la serie del producto y producto per capita de un conjunto de países seleccionados, así como sus consecuencias en el crecimiento económico. Siguiendo con esta línea de análisis, y dentro del contexto de modelos que permiten cambios estructurales múltiples, Ben-David, Lumsdaine y Papell (2003) encuentran pruebas de tasas de crecimiento que aumentan en el tiempo.

Es importante señalar que las evidencias encontradas en esta bibliografía tienen como característica común suponer que los cambios estructurales ocurren instantáneamente, al ser modelados por medio de variables dicotómicas. Sin embargo, tal especificación no forzosamente puede ser la más adecuada.

Este artículo propone modelar los potenciales cambios estructurales, con otra hipótesis, por medio de funciones de transición suave entre regímenes (o tendencias lineales) más que como cambios estructurales instantáneos. La idea básica de las transiciones suaves es simple: el cambio estructural no se asume como un cambio instantáneo, sino como un periodo de transición suave entre diferentes regímenes a lo largo del tiempo. Además, como se analizará más adelante, un cambio estructural instantáneo representa un caso particular de la función de transición suave. Así, nuestro enfoque consiste en probar la hipótesis nula de una raíz unitaria contra la opción de estacionariedad alrededor de una tendencia lineal determinista afectada por transiciones suaves múltiples.3

Basados en el trabajo de Bacon y Watts (1971), Maddala (1977), Granger y Teräsvirta (1993), y de Lin y Teräsvirta (1994), Greenaway, Leybourne y Sapsford (1997) modelan un cambio estructural para series del producto de un conjunto de países seleccionados, como una transición suave entre regímenes. En particular, suponen que cualquier cambio en el desempeño económico después de la liberación comercial puede ser modelado más apropiadamente como una transición suave más que como un cambio discreto en la trayectoria de la serie de tiempo. Según esta línea para modelar los cambios estructurales en la tendencia determinista de una serie de tiempo, Leybourne, Newbold y Vougas (1998) (LNV, en adelante) proponen una prueba tipo Dickey-Fuller, en la cual plantean como otra hipótesis diferente de la de raíz unitaria, la de un cambio estructural en la tendencia lineal de una serie por medio de una transición suave entre dos diferentes tendencias.4 Esta especificación, propuesta por LNV en el contexto de pruebas de raíz unitaria, es bastante flexible en el sentido de que implica como casos límite: uno en el que no existe cambio estructural, y otro en el que hay un cambio y éste es instantáneo. Para esto utilizan la función logística de transición suave que estima endógenamente, de una manera no lineal, la ubicación y duración de la transición (cambio estructural).5

A pesar de sus virtudes, la función empleada por Greenaway, Leybourne y Sapsford (1997) y LNV posee la limitante de que la senda de transición es monótona y simétrica alrededor del punto medio de la transición. Sollis, Leybourne y Newbold (1999) (SLN, en adelante), también dentro del contexto de pruebas de raíz unitaria, extienden tal especificación para modelar un cambio estructural en la tendencia de la serie permitiendo asimetría en la función logística de transición suave.6

Por su parte, Harvey y Mills (2002) desarrollan una extensión a la prueba de LNV al permitir como otra hipótesis diferente de la de raíz unitaria, la de estacionariedad alrededor de una tendencia lineal con transiciones suaves dobles, especificadas por la función logística simétrica. Harvey y Mills (2002) también analizan la posibilidad de que las transiciones sean potencialmente asimétricas. Sin embargo, destacan que aunque tal acercamiento resulta teóricamente deseable, el número de optimaciones numéricas necesarias para obtener simulaciones de los valores críticos para la prueba, lo torna un problema poco tratable desde un punto de vista computacional, por lo que sólo presentan valores críticos para el caso de transiciones suaves dobles en el caso simétrico.

En este artículo se extiende los planteamientos previamente analizados para permitir, como opción a la hipótesis de raíz unitaria, cambios estructurales múltiples, modelados por medio de funciones de transición suave (tanto con la función logística simétrica como con la asimétrica) en la tendencia determinista de las series. En particular, se consideran especificaciones que permiten desde una hasta cuatro transiciones suaves. Hay que recalcar que nuestro planteamiento se diferencia del de Harvey y Mills (2002), quienes permiten transiciones suaves dobles, en el hecho de que aquí se impone la restricción de que las transiciones estimadas no se mezclen. Es decir, no se permite que existan transiciones suaves (cambios estructurales) en maneras consecutiva o sobrepuesta. Esto con la finalidad de que sea posible separar debidamente entre regímenes o tendencias lineales y periodos de transición hacia una nueva tendencia.7

Para México puede encontrarse evidencia de cambios estructurales en el PIB real y real per capita en Noriega y Ramírez-Zamora (1999), Noriega y de Alba (2001), y Astorga, Bergés y Fitzgerald (2005), mientras que Ben-David y Papell (1998) sólo presentan evidencia para el producto per capita, y Castillo y Díaz-Bautista (2002) y Noriega, Soria y Velázquez (2008) sólo para el producto real. La característica común de estos estudios es que presuponen que el cambio estructural en las variables es instantáneo, de tal manera que no existen antecedentes en la bibliografía de estimar cambios estructurales modelados como transiciones suaves para México. Por otra parte, estos trabajos consideran diferentes periodos muestrales a lo largo del siglo XX, en el que 1921 es el periodo más lejano, y ninguno incluye datos del presente siglo. Aquí se utiliza una base de datos que no ha sido utilizada antes para estudiar la trayectoria de largo plazo del producto y producto per capita de México. A saber, se actualiza hasta 2008 la base de datos recién desarrollada por Barro y Ursúa (2008, 2009), con lo cual se obtiene información para el PIB real y real per capita de México que abarca el periodo 1895-2008.

En síntesis, el uso de funciones de transición suave para estimar cambios estructurales múltiples nos ayudará a distinguir entre etapas de crecimiento sostenido y fases de transición en la trayectoria del producto y producto per capita de México. Nuestros resultados muestran que para el PIB real (PIB real per capita) es posible rechazar la hipótesis nula de raíz unitaria, en favor de la otra hipótesis que modela cuatro etapas de crecimiento estacionario alrededor de una tendencia determinista: 1895-1924, 1935-1952, 1956-1978 y 1989-2008 (1895-1908, 1916-1951, 1954-1978 y 1990-2008), separadas por tres épocas de transición, 1925-1934, 1953-1955 y 1979-1988 (1909-1915, 1952-1953 y 1979-1989). Esto es, parecen existir tres cambios estadísticamente significativos en la tasa de crecimiento de largo plazo para cada una de las series.

Este artículo está organizado de la siguiente manera. La sección I explica la metodología propuesta para estimar transiciones suaves múltiples, y probar la hipótesis de raíz unitaria en las series en estudio. En la sección II se presenta la base de datos utilizada y se realiza un análisis exploratorio preliminar de la evolución de las variables. La presentación y análisis de las estimaciones se realiza en la sección III, mientras que el análisis de las consecuencias de éstas para el crecimiento económico se aborda en la sección IV. Al final se presenta las conclusiones.

I. Metodología

El procedimiento que se utiliza para estimar la especificación que permite cambios estructurales múltiples, modelados por funciones de transición suave, y realizar la prueba de raíz unitaria, se basa en el espíritu de la metodología propuesta por LNV y SLN, quienes permiten la posibilidad de una transición suave entre tendencias lineales. Sin embargo, a diferencia de LNV y SLN, en este estudio se propone una extensión de estos métodos, al permitir la presencia de transiciones suaves múltiples.

Específicamente, se obtendrá resultados permitiendo desde una hasta cuatro transiciones suaves, tanto para el caso estrictamente simétrico (basados en LNV) como permitiendo la posibilidad de asimetría en las transiciones (basados en SLN). Así, se estiman ocho distintas especificaciones tanto para el PIB real como para el PIB real per capita. La decisión de permitir un máximo de cuatro cambios estructurales (transiciones suaves) se justifica por motivos de carácter técnico, y es argumentada líneas abajo, en el que se explica con detalle el procedimiento de estimación empleado en este estudio.

El método propuesto se efectúa en dos etapas: i) se estima el modelo que permite m = 1,...,4 cambios estructurales (transiciones suaves), del cual se extraen los residuales, y ii) se aplica una prueba tipo Dickey-Fuller (1979) aumentada (DFA) en estos residuales. Este procedimiento es aplicado al PIB real y PIB real per capita, para m = 1,..., 4. A continuación se describe cada etapa.

Etapa 1

El modelo para estimar m = 1,...,4 transiciones suaves entre regímenes utiliza el método de mínimos cuadrados no lineales (MCNL), y se especifica como sigue:

(1)

En esta ecuación, yt representa el logaritmo natural del PIB real o el logaritmo natural del PIB real per capita, en la que t = 1, 2,..., T , T es el tamaño de la muestra, y m denota el número de cambios estructurales (transiciones suaves) por estimar. Nótese que la especificación en (1) permite m cambios simultáneos de nivel y pendiente. Con m = 1 esta ecuación se reduce al denominado modelo C planteado por SLN que permite estimar un cambio tanto de nivel como de pendiente en una tendencia lineal. Gt (τ, γ, θ) es la función logística generalizada, basada en una muestra de tamaño T, la cual controla la transición entre dos regímenes (tendencias lineales), y se define como (véase SLN):

(2)

La interpretación de los parámetros de la función es como sigue. El parámetro γ gobierna la velocidad de transición. Con un γ pequeño, la función Gt (τ, γ, θ) toma un largo periodo en atravesar el intervalo (0, 1), y en el caso límite de γ = 0, Gt (τ, γ, θ) = 2-θ para todo t. Por lo contrario, para valores grandes de γ, Gt(τ, γ, θ) atraviesa el intervalo (0, 1) de manera rápida. Específicamente, conforme γ se acerca a +∞ la función cambia de 0 a 1 de modo instantáneo, mientras que al acercarse a a - ∞ cambia de 1 a 0 de manera instantánea. El parámetro θ introduce la posibilidad de asimetría en la función Gt , o sea, permite la posibilidad de que la velocidad de la transición sea distinta antes y después del punto de inflexión, para 0 ˂ θ ˂1. El caso de θ=1, corresponde a la función logística simétrica, St , introducida por LNV [St (τ, γ) = Gt (τ, γ 1)], en el contexto de pruebas de raíz unitaria, en el que la trayectoria de transición es monótona y simétrica alrededor del punto de inflexión. En el otro extremo, conforme θ se acerca a 0 una asimetría extrema es generada. Finalmente, el parámetro τ se refiere a la ubicación del cambio estructural. En el caso de una transición simétrica, θ=1, la interpretación de τ es relativamente sencilla, puesto que τT es el punto de inflexión de Gt (τ, γ 1) y también puede interpretarse como su punto medio, ya que GtT (τ, γ 1) = 0.5. Sin embargo, con algún grado de asimetría (0 < θ < 1) el punto de inflexión de Gt (τ, γ θ) está dado por t = τT + γ -1 θInθ, lo cual hace que la interpretación del parámetro τ no sea tan directa.

En la Gráfica 2 se presenta a función logística generalizada, con T = 100, dejando fijo al parámetro de ubicación τ = 0.5, tal que el cambio estructural se encuentre alrededor de la mitad de la muestra. Con la finalidad de tener una mejor perspectiva visual de las distintas formas que puede adquirir la función Gt (τ, γ θ), en la Gráfica 2 se realiza un acercamiento al eliminar los primeros y últimos 30 periodos. Primeramente obsérvese el caso de St1τ=0.5,γ=9=Gt1τ=0.5,γ9,θ=1, que con una γ suficientemente grande (γ = 9) la función logística cambia de 0 a 1 en aproximadamente un periodo, tal que imita a una variable dicotómica tradicional, por lo que puede estimar un cambio estructural instantáneo. En contraparte, con un parámetro relativamente bajo (γ = 0.25), ejemplificado por St2=Gt2 y Gt3, la función implica un largo periodo de transición. En el caso simétrico, θ = 1 véase St2 en la Gráfica 2, la función tarda más de 20 periodos en moverse de 0 a 1, teniendo el punto de inflexión, τT, en el periodo 50. Por otra parte, en el caso en que θ < 1, la naturaleza de la asimetría es determinada por el signo del parámetro γ. Para un γ ˃ 0, la transición es más lenta al inicio que al final, como lo muestra Gt3 en la Gráfica 2. Para estos valores de los parámetros, la transición entre regímenes es más lenta antes del punto de inflexión, [τT + γ -1 θInθ] = 49.1, que después de éste. Por último, Gt4 muestra la situación inversa, con γ ˂ 0, en que la transición es más lenta hacia el final de la transición.

GRÁFICA 2 Función logística simétrica St(τ, γ ) y generalizada Gt(τ, γ , θ) 

Como se observa, la estimación de las funciones de transición suave puede resultar en periodos de transición muy largos con una γ pequeña. Entonces, al momento de estimar simultáneamente m funciones de transición suave, una posibilidad es que las funciones estimadas impliquen periodos de transición muy largos que se sobreponen o bien que sean consecutivos. Por tanto, es necesario introducir una variable, h > 0, que denota el tamaño de segmento mínimo permitido entre la ocurrencia de dos cambios estructurales en la estimación. Esta h desempeña el mismo papel que en el caso estándar en que los cambios estructurales son modelados con variables dicotómicas, lo que evita que existan dos periodos de manera consecutiva con cambio estructural.8 Así, a diferencia de Harvey y Mills (2002), quienes desarrollan una prueba de raíz unitaria extendiendo el procedimiento de LNV y permitiendo transiciones suaves dobles, aquí imponemos la restricción de que las transiciones estén separadas por un segmento mínimo de tamaño h. De esta manera, la estimación de la tendencia determinista de las series puede ser interpretada como una tendencia lineal quebrada o separada por cambios estructurales graduales (transiciones suaves) o instantáneos, dependiendo del valor final que tomen los parámetros de la función logística.

Para instrumentar esta restricción, es necesario identificar adecuadamente el periodo en que inicia y termina cada una de las transiciones. Partimos de que en general es posible encontrar el periodo en que la función logística generalizada toma algún valor predeterminado entre 0 y 1, al dejar fijo el valor de la función y despejar para t.9 Esto es, déjese Gtτ,γ,θ=G¯, con 0<G¯<1, entonces de la ecuación (2) tenemos G-=1/1+exp-γt-τT/θθ. Tras algunas manipulaciones algebraicas es posible obtener una expresión que nos dé el periodo en que la función logística toma el valor G¯. Esto es,

(3)

Sean Pi y Pf los periodos en que inicia y termina, respectivamente, la transición. Dados los valores de los parámetros γ, τ y θ es posible calcular Pi y Pf de acuerdo con la ecuación (3), fijando un valor predeterminado para G¯. Por ejemplo, considérese la función Gt3τ=0.5,γ=0.25, θ=0.1 ilustrada en la Gráfica 2. Entonces, dado que γ = 0 sabemos que la función se mueve de 0 a 1, por lo cual podríamos considerar, por ejemplo, que la transición inicia a partir del periodo en que Gt3=0.01 y finaliza donde Gt3=0.99. Así, se tendría que PiG¯=0.01=31.6 y PfG¯=0.99=50.9, con lo cual tenemos una transición entre regímenes que dura aproximadamente Pf - Pi = 19.3 periodos. Nótese que con γ = 0 la función logística se mueve de 1 a 0, por ejemplo Gt4 en la Gráfica 2, por lo cual el cálculo de los periodos en que inicia y termina la transición debe efectuarse cambiando los valores de G¯ entre Pi y Pf , esto es, PiG¯=0.99=49.1 y PfG¯=0.01=68.4. De esta manera puede calcularse la duración de transición y conocer los periodos en que ésta principia y culmina. En este artículo, al igual que en los ejemplos mencionados, se deja que G¯ tome los valores de 0.01 y 0.99, con lo cual podemos tener el cálculo de aproximadamente 98% de la transición.10

Una vez definida la manera en que se efectúa el cálculo de los periodos de inicio y fin de una transición, volvemos a la cuestión de estimar un modelo que permita m transiciones suaves separadas por un segmento mínimo, denotado por h. Definamos como Pij y Pfj a los periodos en que inicia y termina la j-ésima transición, respectivamente. Entonces, se busca minimizar la suma de residuales al cuadrado de la ecuación (1), sujeta a que Pfj - 1 + h < Pij , para todo j = 2,..., m. Esto es, se busca los parámetros que resuelven el siguiente problema11

(4)

Nótese que Pij y Pfj-1 dependen de los parámetros estimados de la correspondiente función logística y están dados por las siguientes expresiones, que a su vez están basadas en (3):

En síntesis, en esta etapa 1, tanto para el PIB real como para el PIB real per capita se estima la tendencia determinista modelada como una tendencia lineal que permite hasta m cambios estructurales (transiciones suaves), para m = 1, 2, 3 y 4.

Como se estableció anteriormente, dado que la especificación de la ecuación (4), para m ≥ 1, es no lineal en los parámetros (τj, γj, θj ), j = 1,..., m, es necesario emplear la técnica de mínimos cuadrados no lineales para obtener una estimación de éstos. En este caso, la función objetivo (suma de residuales al cuadrado) además del mínimo global puede tener mínimos locales, por lo cual es necesario probar un conjunto amplío de valores iniciales con la finalidad de asegurar que se ha alcanzado el mínimo global. Entonces, para la estimación de la ecuación (4) se considera un conjunto amplío de valores iniciales para los parámetros (τj, γj, θj ), j = 1,..., m, denotados como τj0,γj0,θj0, mientras que para el resto de parámetros, αi y βi, i = 0,..., m, dada su cualidad lineal, podemos emplear como valores iniciales sus correspondientes estimadores de mínimos cuadrados ordinarios (MCO) sujetos a Gtj0τj0,γj0,θj0. La rutina para la estimación de la ecuación (4) es instrumentada en EViews 7, que en la estimación por MCNL usa el algoritmo (iterativo) de optimación Gauss-Newton (también conocido como el algoritmo de Berndt, Hall, Hall y Hausman).12

El conjunto de valores iniciales que se prueba para cada m, no forzosamente es el mismo debido a dos razones. La primera se debe a la implicación de los valores iniciales seleccionados. Por ejemplo, con un m = 1 es posible probar con valores iniciales de la velocidad de transición muy lenta, digamos γ 0 = 0.1 que implica una transición de aproximadamente 60 periodos. Sin embargo, para el caso de m = 4 no tiene mucho sentido probar estos valores iniciales ya que implicarían que las cuatro transiciones modeladas se podrían sobreponer dado el tamaño de la muestra. Por tanto, conforme m se incrementa, se reduce el conjunto de valores iniciales considerados para la velocidad de transición, γj . La segunda razón es producto de los límites de la capacidad computacional. Conforme el número de transiciones por estimar crece, el conjunto de valores iniciales también se hace más grande, al punto en que la memoria virtual de la computadora es agotada.13

Considérese, por ejemplo, el siguiente conjunto de valores iniciales utilizado en la sección III para estimar el modelo con m = 4 transiciones simétricas. Para cada función logística se prueba diversos valores iniciales para γj y τ j . Específicamente, se deja que la velocidad de transición tome los siguientes valores γj0=1,2,3,4,5,6, mientras que para el parámetro de ubicación se utilizaron los valores τj0=.08,.16,.24,.32,.39,.47,.55,.63,.71,.79,.87,.95, lo que implica que cada nueve periodos se ubica el punto medio de las transiciones. Así, se tiene un total de 124*64=641520 combinaciones de valores iniciales. Al aumentar el número de transiciones a m = 5, utilizando el mismo conjunto de valores iniciales para γj y τ j que en el caso m = 4, encontramos que la rutina agota la memoria disponible de la computadora antes de obtener los resultados. Esto se debe a que en el caso m = 5, deben de procesarse 125*65=6158592 de combinaciones.14 Aunque reducir el conjunto de valores iniciales podría en principio resolver este problema, en la práctica observamos que, para el caso m = 5, los resultados se tornan extremadamente sensibles a los valores iniciales escogidos, de manera que no es posible asegurar que se ha alcanzado el mínimo global de la ecuación (4). Esto no ocurre para las estimaciones con m = 1, 2, 3 y 4 presentadas en la sección III, en la que los resultados son robustos a varios conjuntos de valores iniciales. Debido a este problema técnico se decide dejar en cuatro el número máximo de cambios estructurales por estimar.

Etapa 2

Se efectúa una regresión tipo DFA y se calcula el estadístico t asociado con ρ, el parámetro autorregresivo en la regresión:

(5)

en el que ϑ^t son los residuales calculados de la estimación de la ecuación (4).

Así, en esta segunda etapa se extrae los residuales de la especificación estimada en la etapa 1, y se utilizan para estimar por MCO la ecuación (5), para obtener el estadístico t del parámetro autorregresivo ρ, denotado como t^ρ para probar la hipótesis nula de raíz unitaria. Para realizar la prueba de raíz unitaria, con base en Rudebush (1992) y Diebold y Senhadji (1996), simulamos la distribución y obtenemos la densidad empírica del estadístico t^ρ tanto con la hipótesis de que el verdadero modelo es estacionario en tendencia con m transiciones [I (0) + m cambios, denotado como ETm ], como con la hipótesis de que el verdadero modelo es estacionario en diferencias [I (1), denotado como ED]. Una vez simuladas las dos distribuciones, se compara la posición en que el estimador muestral del estadístico t para probar la hipótesis de raíz unitaria, cae relativo a las densidades empíricas simuladas (ETm y ED ). Específicamente, las hipótesis que se consideran son:15

(6)

en la que

(6.1)

(7)

Con el modelo ETm, yt fluctúa estacionariamente alrededor de una tendencia lineal determinista perturbada por m cambios estructurales (transiciones suaves), mientras que con el modelo ED , yt sigue un proceso de raíz unitaria con deriva (o constante).16

Para discriminar entre los modelos ETm y ED , para cada valor de m = 1, 2, 3, 4, se simula la distribución de con las hipótesis de que el verdadero modelo es ETm y ED ; en ambos casos las distribuciones son estimadas de los datos, utilizando el método de remuestreo.

Para simular la distribución del estadístico t^ρ con el modelo ETm se siguen los siguientes pasos.

A. Estimación de los parámetros para el modelo ETm : A 1. Se estima la ecuación (6) por mínimos cuadrados no lineales, y se calcula la serie ajustada, y^t y los residuales de la regresión, ϑ^t. A2. Utilizando ϑ^t se estima (5) por MCO, se obtienen los parámetros estimados (ρ^ y δ^t) se calculan los residuales de esta regresión, ε^t.

B. Simulación de la densidad empírica: B 1. Usando las primeras k + 1 observaciones de ϑ^t como condiciones iniciales, los parámetros estimados en el paso (A2.), y con residuales escogidos aleatoriamente con remplazo [ε^t de la ecuación (5)] se genera una muestra de ϑ^t como en la ecuación (6.1), que llamaremos ϑ^t+B 2. Se genera una muestra de yt como:y^t+=y^t+ϑ^t+. B3. Utilizando la serie generada yt+ se estima por MCO la ecuación (6) (no se estima por MCNL, pues ya se conoce la localización y duración de las transiciones), y se calcula los residuales de esta regresión. B4. Sobre estos residuales se realiza la regresión tipo DFA, ecuación (5), y se obtiene el valor del estadístico de prueba t^ρ, para probar la hipótesis nula ρ = 0. B 5. Los pasos B1-B4 se repiten 10 mil veces, y los correspondientes 10 mil valores del estadístico de prueba t^ρ son usados para construir la función de densidad empírica de este estadístico con el modelo ETm , llamada fETmtρ, m = 1, 2, 3, 4.

La simulación de la distribución del estadístico t^ρ con el modelo ED sigue los siguientes pasos.

C. Estimación de los parámetros para el modelo ED: C1. Se estima la ecuación (7) por MCO, se obtiene los parámetros estimados (K^ y φ^i) se calcula los residuales de esta regresión, η^t.

D. Simulación de la densidad empírica: D 1. Usando las primeras k + 1 observaciones de yt como condiciones iniciales, los parámetros estimados en el paso C1, y utilizando residuales escogidos aleatoriamente con remplazo [η^t de la ecuación (7)] se genera una muestra de yt con base en la ecuación (7), que llamaremos yt*. Se estima por MCO la ecuación (6) utilizando la serie yt* (no se estima por mínimos cuadrados no lineales, pues ya se conoce la localización y duración de las transiciones), y se calcula los residuales de esta regresión. D3. Sobre estos residuales se realiza la regresión tipo DFA, ecuación (5), y se obtiene el valor del estadístico de prueba t^ρ, para probar la hipótesis nula ρ = 0. D4. Los pasos de D1-D3 se repiten 10 mil veces, y los correspondientes 10 mil valores del estadístico de prueba son usados para construir la función de densidad empírica de este estadístico con el modelo ED, llamada fEDt^ρ.

Por último, se obtiene la posición en la que el estimador muestral del estadístico tt^ρ muestral para probar la raíz unitaria, cae relativo a las densidades empíricas simuladas. Estas posiciones son calculadas como la masa de probabilidad a la izquierda de tt^ρ muestral, denotadas como PETmPrt^pt^ρ muestralfETmt^ρ y PEDPrt^pt^ρ muestralfETmt^ρ. Se concluye en favor de alguna especificación ETm cuando se tiene que pED0.10 y 0.10<pETm<0.90. 17

II. Datos

Las series de datos anuales del PIB real y real per capita son tomadas de Barro y Ursúa (2008, 2009), quienes con base en las tasas de crecimiento de los cinco mayores sectores de la economía (agricultura, minería, energía, transformación y servicios) construyen la información para el periodo de la Revolución Mexicana entre 1911 y 1920. Para este estudio, se actualiza esta base de datos (originalmente para 1895-2006) hasta el año 2008. En el Apéndice 1 se muestra los pormenores de la actualización.

En la Gráfica 3 se presenta los datos tanto en logaritmos (parte A) como en primeras diferencias (parte B). Como se puede apreciar en esta gráfica (en particular en la parte B), las series muestran al menos tres pautas de crecimiento diferenciadas. En la primera, que comprende desde el comienzo de la muestra hasta principios de 1930, se observa una senda de crecimiento muy errático en ambas series, alternando entre tasas de crecimiento positivas y negativas. El segundo periodo, que dura aproximadamente 50 años, abarca de mediados del decenio de los treinta hasta principios de los ochenta, y es caracterizado por tasas de crecimiento invariantemente positivas, que en algunos años incluso llegaron a alcanzar tasas superiores a 10% en el caso del PIB real, y presenta una tasa promedio de aproximadamente 6 y 3%, para el PIB real y real per capita, respectivamente. Finalmente, de principios de los años ochenta y hasta el final de la muestra se observa una pauta de crecimiento baja, afectado principalmente por el estancamiento de la mayor parte de los años ochenta, la crisis de 1994-1995 y la desaceleración de principios del decenio del 2000. En estos últimos años, en promedio se observan tasas de crecimiento para el PIB real y real per capita de 3 y 1%, respectivamente.

GRÁFICA 3 PIB real y PIB real per capita (1895-2008) 

Uno de los objetivos del presente trabajo es precisamente darle mayor precisión a la identificación de las posibles sendas de crecimiento de las series del producto y producto per capita de México a lo largo de la muestra.

III. Resultados empíricos

Como punto de partida para el análisis empírico comenzamos aplicando pruebas de la hipótesis nula de raíz unitaria, que se contrasta contra la otra hipótesis de que las series siguen una tendencia lineal sin cambios estructurales. Estas pruebas servirán como marco de referencia (benchmark), contra los cuales compararemos los resultados (presentados líneas abajo) de la aplicación de la prueba de raíz unitaria permitiendo con la otra hipótesis que la tendencia se modele utilizando cambios estructurales múltiples en la forma de transiciones suaves.

1. Pruebas de raíz unitaria sin cambios estructurales

Se presenta los resultados de las pruebas de Ng y Perron (2001) que tienen la hipótesis nula de una raíz unitaria (los nombres de estas pruebas son MZa , MZt , MSB y MPT ). Como sugieren Perron y Qu (2007), la prueba se realiza eligiendo el número de rezago óptimo por incluir usando el criterio modificado de Akaike (MAIC), ya que al emplear este criterio de información modificado se mejora el poder de la prueba en muestras finitas.

El Cuadro 1 presenta los resultados de la aplicación de las pruebas de raíz unitaria de Ng-Perron (utilizando el MAIC), en el contexto de la metodología propuesta por Pantula (1989), para determinar el número de raíces unitarias presentes en cada serie. Ésta consiste en un procedimiento secuencial asintóticamente congruente para probar la hipótesis nula de Hr : exactamente r raíces unitarias, contra la opción de Hr - 1: exactamente (r - 1) raíces unitarias, con r = n, ..., d + 1, d, en que n(≥ d) es un número máximo de raíces unitarias que se supone que está presente en los datos, y d es el verdadero número de raíces unitarias en los datos. Pantula (1989) sugiere que la hipótesis debe ser probada de manera secuencial en el orden Hn, Hn-1 ,...,Hd .18 Aquí suponemos que el máximo número posible de raíces unitarias presentes en los datos es de dos. Basados en los resultados de Pantula (1989), las hipótesis deben ser probadas secuencialmente en el orden H2 y H1 . El componente determinista se modifica según sea el número de raíces que están siendo probadas. Así, en el caso H2 contra H1 , se incluye sólo intercepto, mientras que en el caso H1 contra H0 se incluye intercepto y pendiente de la tendencia lineal.19 Estos resultados son resumidos en el Cuadro 1, en el que se deja en cinco el rezago máximo empleado, kmax .

a k* es elegido de acuerdo con Perron y Qu (2007).

***, ** y * significan rechazo de la hipótesis nula al nivel de 1, 5 y 10%, respectivamente.

CUADRO 1 Resultados de las pruebas de Ng-Perron a 

Para ambas series la inferencia que se desprende de estas pruebas es la misma: se rechaza la presencia de dos raíces unitarias (al nivel de 1%), y no se rechaza la presencia de una raíz unitaria. De esta manera, podemos concluir que las series del PIB real y del PIB real per capita son no estacionarias y a lo más tienen una raíz unitaria. Esta es la conclusión a la que se llegaría cuando la otra hipótesis es modelada como una simple tendencia lineal. Sin embargo, como mostramos líneas abajo, las conclusiones cambian diametralmente cuando permitimos la presencia de posibles cambios estructurales múltiples en la función de tendencia determinista para ambas series.

2. Raíces unitarias contra transiciones suaves múltiples

En esta subsección se presenta los resultados de las pruebas de raíz unitaria cuando permitimos la estimación de transiciones suaves múltiples para el PIB real y PIB real per capita de México. En el Cuadro 2 se resume los principales resultados, que nos permiten discriminar entre los modelos estacionario en diferencias (ED) y estacionario en tendencia con m = 0, 1, 2, 3 y 4 transiciones suaves (ETm ). La primera columna indica el número de transiciones permitidas en la función de tendencia lineal, m. La columna dos muestra el rezago elegido, k *, en la estimación de la ecuación (5). Según Ng y Perron (1995), el rezago óptimo se elige de acuerdo con su significancia. Específicamente, se empieza por fijar un número máximo arbitrario de rezagos, kmax ; cuando el coeficiente estimado del kmax -ésimo rezago es no significativo, se estima un modelo con kmax - 1 rezago; si el último vuelve a ser no significativo, se estima un nuevo modelo con kmax - 2 rezagos. Así se continúa hasta encontrar el primer rezago significativo, k *. Como en la subsección anterior se deja kmax = 5.20 En la tercera columna se presenta el estadístico t para probar la hipótesis nula de raíz unitaria, t^ρ muestral. La masa de probabilidad a la izquierda de t^ρ muestral, en cada una de las especificaciones simuladas, pETmETm y pEDED, es presentada en las últimas dos columnas del cuadro. Estos resultados son ordenados en seis bloques. Primeramente se muestran los resultados para el PIB real con m = 0 (parte A), después cuando se permite la asimetría en las funciones de transición (parte B), y para finalizar cuando se usan las funciones simétricas para estimar las transiciones (parte C). Similarmente las partes D, E y F resumen los resultados para el PIB real per capita.

CUADRO 2 Modelos ED contra modelos ET que permiten hasta cuatro transiciones suaves, simétricas y asimétricas 

En la parte A del Cuadro 2 se complementa las pruebas presentadas en el Cuadro 1: en el caso en que la otra hipótesis se modela como una tendencia lineal, no se puede discriminar entre las dos hipótesis puesto que no se puede rechazar la nula, pED > 0.10, ni la opción, 0.10 < pET0 < 0.90. Al mismo tiempo, el Cuadro 2 muestra que para el PIB real, el modelo ED (de raíz unitaria) es rechazado (pED < 0.10) para cualquier número positivo de transiciones estimadas, m , y sin importar si las transiciones son modeladas de manera simétrica o permitiendo asimetría. De hecho la especificación que permite cuatro transiciones (cambios estructurales) modeladas con la función simétrica, es la que rechaza el modelo ED con más fuerza (pED = 0.000). Además, para este mismo modelo, pET4 es la más cercana a la mitad de la distribución empírica. En la parte izquierda de la Gráfica 4 se muestra las distribuciones empíricas para esta especificación, junto con el correspondiente valor estimado del estadístico de prueba t^ρ muestral = -8.153. Como se observa, el valor del estadístico, t^ρ muestral se encuentra muy cerca de la mitad de la distribución fET4t^ρ, mientras que éste se encuentra muy alejado en la cola izquierda de la distribución muestral con la hipótesis de raíz unitaria, que corresponde al modelo ED. Por tanto, estos resultados sugieren que se puede rechazar la hipótesis de una raíz unitaria para el PIB real, en favor de un modelo estacionario en tendencias con cuatro cambios estructurales modelados como transiciones suaves simétricas.

GRÁFICA 4 Funciones de distribución empírica modelos ET4 y ED, m = 4 

Los resultados para el PIB real per capita son prácticamente iguales, por lo que omitimos su análisis. La representación gráfica de estos resultados se muestra en la parte derecha de la Gráfica 4. Al igual que en el caso del PIB real, la especificación con cuatro transiciones simétricas es la que rechaza con más fuerza el modelo ED, y en la cual pET4 es la más cercana a la mitad de la distribución empírica.

La Gráfica 4 muestra cómo, para ambas variables, es mucho más probable que el valor muestral del estadístico t^ρ muestral haya sido generado por un proceso estacionario en tendencias con cuatro transiciones suaves, que por un proceso de raíz unitaria.

En el Cuadro 3 se presenta los valores estimados de los parámetros de las funciones de transición para las ocho diferentes especificaciones consideradas para cada serie. La columna uno hace referencia a la variable considerada. La segunda columna indica el número de funciones de transición introducidas en la ecuación (1), mientras que la tercera columna enumera y distingue entre la función logística de transición simétrica Stjτj,γj=Gtjτj,γj,θj=1 y la función logística generalizada que permite asimetría Gtjτj,γj,θj. Los valores estimados de las funciones de transición se presentan de la cuarta a la sexta columna. Las últimas tres columnas muestran las fechas de los periodos (redondeados) de inicio, inflexión y finalización de la transición implicada por las funciones de transición estimadas, calculados como se explica en la sección I.21

***, ** y * significa rechazo de la hipotesis nula al nivel de 1, 5 y 10%, respectivamente.

CUADRO 3 Estimación de las funciones de transición 

Primero, nótese que cuando se usa la función logística generalizada, en ninguno de los casos el parámetro de asimetría, θ^j, es significativamente distinto de 0 al menos a un nivel de 10%. La velocidad de transición γ^j, con excepción de un par de casos tampoco es estadísticamente distinta de 0 en las distintas estimaciones. Nótese, sin embargo, que esta ausencia de significancia de la velocidad de transición se elimina cuando las transiciones son modeladas con la función logística simétrica, Stj (·). De hecho, sólo cuando las estimaciones son casi instantáneas el parámetro γ^j no puede rechazar estadísticamente que es 0. Al respecto, Greenaway, Leybourne y Sapsford (1997), pp. 801, explican que "este fenómeno es porque un gran número de valores estimados de γ llevan a valores muy similares de St (·). Las consecuencias en la práctica de esto es que los errores estándar de la estimación de MCNL de γ pueden aparecer artificialmente grandes y por tanto no deben de ser de manera foroza tomados como indicador de la insignificación de la estimación".22 El único parámetro que con regularidad resulta significativo en las especificaciones de las transiciones con la función logística generalizada es el parámetro de ubicación, τ^j.

Por otra parte, es de destacar que cuando se incluye un número reducido de cambios estructurales modelados por las funciones de transición suave, las transiciones son extraordinariamente largas para ambas variables. Para el PIB real y real per capita, con m = 1, 2 y 3 con Gtj (·) y m = 1 con Stj (·) se tienen valores estimados de la velocidad de transición, γj muy pequeños. Esto provoca larguísimos periodos de transición, que simplemente parecen carecer de sentido económico alguno.23 Por ejemplo, para el PIB real con una transición modelada por Gtj (·), la estimación implica una transición asimétrica θ^1=0.079 que inicia en 1920 y termina en 2462. En contraparte, para el PIB real con m = 4, con cualquiera de las funciones de transición [Gtj (·) y Stj (·)], y para el PIB real per capita con m = 4 con Gtj (·), y m = 2, 3 y 4, utilizando Stj (·) los resultados tienen más sentido con el desempeño de la economía mexicana a lo largo de la muestra considerada.

En suma, dada la evidencia mostrada en el Cuadro 3 acerca de las transiciones estimadas para las distintas especificaciones, y la conclusión en favor de los modelos estacionarios en tendencia extraída del Cuadro 2, se decide en favor de la especificación con cuatro transiciones simétricas tanto para el PIB real como para el PIB real per capita. Por tanto, de aquí en adelante nos enfocaremos en los resultados de estas especificaciones.24 Sin embargo, notemos que los resultados en el Apéndice 2 muestran que para el PIB real el primer cambio de pendiente en la [β 1 ecuación (1)], no resulta significativo al nivel de significancia de 10%. Lo mismo ocurre para el segundo cambio de pendiente [β 2 en la ecuación (1)] en el caso del PIB real per capita. Esto significa que tanto el primer cambio para el producto real, como el segundo para el per capita, son solamente de nivel, y no de nivel y pendiente. Los otros tres cambios en cada variable parecen ser cambios tanto de nivel como de pendiente.

Dado este resultado, decidimos estimar una versión restringida de estas especificaciones. Así, se estima la ecuación (4) con m = 4, empleando la función logística simétrica Stjτj,γj para modelar cada transición, restringiendo β 1 = 0 para el PIB real, y β 2 = 0 para el PIB real per capita.25

Los resultados de la estimación de estas especificaciones restringidas son resumidos en los Cuadros 4 y 5. El Cuadro 4 presenta la estimación por MCNL de la ecuación (1) con las restricciones mencionadas para cada variable. En la parte superior se presenta las principales medidas de evaluación de la ecuación estimada, mientras que en la inferior se presenta los valores estimados de los parámetros y su significancia.

a EE, error estandar de la regresion; SRC, suma de residuales al cuadrado; AIC, BIC y HQC, se refieren a los criterios de informacion de Akaike, Schwarz y Hannan-Quinn, respectivamente.

***, ** y * rechazo de la hipotesis nula al nivel de 1, 5 y 10% respectivamente.

CUADRO 4 Resultados de la estimación restringida con m = 4a 

***, ** y * Rechazo de la hipotesis nula al nivel de 1, 5 y 10%, respectivamente.

CUADRO 5 Estimación de las funciones de transición 

Nótese que al comparar los resultados de la parte superior del Cuadro 4 con los presentados en el Cuadro A2 del Apéndice, la especificación restringida para el PIB real y real per capita, resulta mejor que el resto de las especificaciones estimadas en términos de los criterios de información AIC, BIC y HQC. Además todos los parámetros resultan estadísticamente distintos de 0 al nivel del 1 por ciento.

En el Cuadro 5 se presenta los valores estimados de los parámetros de las funciones de transición y los años correspondientes de las transiciones para cada una de las series. Comparando estos resultados con los obtenidos de las estimaciones del modelo no restringido (véase Cuadro 3), podemos apreciar que no existen modificaciones importantes en la ubicación y duración de las transiciones. Los valores estimados del parámetro de ubicación, χ j sólo difieren en milésimas, mientras que en el caso de la duración de la transición las diferencias son mínimas entre los modelos restringidos y no restringidos. Las diferencias observadas de a lo más un año en las fechas de inicio y final de la transición, se deben a una diferencia de decimales en el cálculo de los periodos. Para observar esto, considérese la tercera transición del modelo no restringido del PIB real en el que el periodo (y no año) final de la transición es 60.451, mientras que el periodo final del restringido ocurre en el periodo 60.504, por lo que las diferencias observadas en las fechas en el Cuadro 5 se deben básicamente al redondeo de las cifras. Realizamos también la estimación de la matriz de covarianza utilizando errores estándar congruentes con autocorrelación y heteroscedasticidad (llamados errores estándar HAC por sus siglas en inglés), sin encontrar ningún cambio cualitativo. Esto es, la inferencia que se desprende de utilizar estos errores estándar es idéntica a la presentada en los Cuadros 3, 4, 5 y A2.

Las Gráficas 5 y 6 presentan los resultados de la estimación del PIB real y PIB real per capita de México, respectivamente. En la parte superior de las gráficas se presenta la serie original (en logaritmos) junto a la serie ajustada con cuatro transiciones simétricas, en las que las áreas sombreadas corresponden a las etapas de transición presentadas en el Cuadro 5.26 En estas gráficas es posible apreciar el buen ajuste que realizan las especificaciones de la tendencia lineal afectada por cuatro transiciones simétricas. Los residuales extraídos de la estimación, y sobre los cuales es aplicada la prueba tipo DFA en la etapa 2 de la metodología, se muestran en la parte inferior de las gráficas. En ambos casos, visualmente los residuales lucen como procesos estacionarios. Estos residuales pueden interpretarse como el ciclo económico o ciclo de negocios en México.

GRÁFICA 5 Estimación (restringida) de cuatro transiciones simétricas: PIB real 

GRÁFICA 6 Estimación (restringida) de cuatro transiciones simétricas: PIB real per capita 

En torno de este asunto de la estacionariedad, en el Cuadro 6 y Gráfica 7 se presenta los resultados de aplicar la segunda etapa de la metodología a las especificaciones restringidas. En el Cuadro 6 observamos cómo al imponer la restricción de β 1 = 0 para la especificación del PIB real, y de β 2 = 0 para la del real per capita, se mantiene el rechazo del modelo estacionario en diferencias (pED < 0.10), mientras que la probabilidad de que la especificación sea generada por un modelo estacionario en tendencia con cuatro transiciones suaves, se mueve aún más hacia el centro de la distribución empírica, como se observa en la Gráfica 7.

CUADRO 6 Modulos ED contra modelos ET restringidos que permiten cuatro transicciones suaves simétricas 

GRÁFICA 7 Funciones de distribución empírica, modelos ET4 y ED, m = 4 

Ahora, adviértase que aunque hemos mostrado que estos residuales fluctúan de manera estacionaria alrededor de una tendencia afectada por transiciones suaves, ciertamente no esperamos que estas fluctuaciones sean por completo aleatorias, o en otras palabras, que sigan un proceso de ruido blanco. Por lo contrario, dada la naturaleza del ciclo de producción agregado que representan, es de esperar que estos residuales observen cierto grado de correlación temporal o autocorrelación. De hecho, el correlograma correspondiente a cada serie de residuales muestra correlaciones significativas durante los primeros cuatro periodos, que se tornan no significativas para rezagos mayores. Este hecho se refleja en los resultados de las pruebas de raíz unitaria presentados en los Cuadros 2 y 6, en en los que el modelo autorregresivo utilizado para llevar a cabo la prueba (ecuación 5) requiere la inclusión de tres rezagos de la variable dependiente en primeras diferencias (lo que equivale a cuatro rezagos de la variable en niveles), precisamente para tomar en cuenta la autocorrelación contenida en el ciclo económico. Estos resultados sugieren, pues, que las fluctuaciones de corto plazo del producto están autocorrelacionadas por periodos de cuatro años.27

Así, los resultados indican la presencia de cuatro cambios estructurales durante el periodo 1895-2008, que afectan tanto al producto real como al per capita. Nótese que la forma y tiempo de ocurrencia y duración de estos cambios son muy similares para ambas series. Las únicas diferencias son las siguientes: i) la Revolución Mexicana redujo el nivel de ambas series (en un proceso que duró aproximadamente siete años), pero afectó solamente la tasa de crecimiento del producto per capita, y ii) el segundo cambio estructural, normalmente relacionado con la Gran Depresión y varios arreglos financieros e institucionales internos (como la fundación del banco central, de la Comisión Nacional Bancaria, del Partido Nacional Revolucionario, la entrada en vigor del impuesto sobre la renta, del código civil, entre otros), da paso a una tasa de crecimiento mayor sólo para el producto real (una respuesta de dimensión equivalente de la población provocó que el producto per capita se mantuviera constante antes y después del cambio).

Nótese además que tanto el cambio estructural de principios de los años cincuenta como el cambio asociado a la crisis de la deuda de los ochenta tuvieron formas y tiempos de ocurrencia/duración prácticamente idénticas en ambas series.

Por último y con el fin de presentar un marco comparativo a nuestros resultados, en los Cuadros 7 y 8 se resume la evidencia existente en la bibliografía de cambios estructurales instantáneos para las series del PIB real y real per capita de México, respectivamente; mientras que en el Apéndice 3 se muestran los pormenores relativos a estos estudios. Estos cuadros representan esta evidencia como un segmento de línea fragmentado por las fechas en que son estimados los cambios estructurales, y en la parte superior de cada fragmento de línea se presenta la tasa de crecimiento de la variable en la correspondiente etapa.

a En estos estudios se utiliza dos tecnicas en la deteccion de los cambios estructurales como se detalla en el Apéndice 3.

CUADRO 7 Evidencia de cambios estructurales para el (ln) PIB real 

a En estos estudios se utiliza dos tecnicas en la deteccion de los cambios estructurales como se detalla en el Apendice 3.

CUADRO 8 Evidencia de cambios estructurales para el (ln) PIB real per capita 

Como se puede observar, todos los estudios en los que la muestra de estimación inicia en 1921, encuentran un cambio estructural a principios de los años treinta (entre 1930 y 1932). Esta fecha es cercana al periodo en que termina la segunda transición estimada en este artículo para ambas variables (1934, véase Cuadro 5), por lo que, por un lado, se podría decir que la conclusión de ambos acercamientos es casi la misma: existe una nueva tendencia en la trayectoria de las variables que inicia entre 1933 y 1935. Pero, por otra parte, nuestras estimaciones indican que no existe una tendencia lineal definida sino una etapa de transición desde 1925 en el caso del PIB real, y 1923 en el caso del per capita, lo cual marca una considerable diferencia en cuanto a las consecuencias de los resultados del crecimiento de las variables.

Una conclusión análoga puede extraerse de los Cuadros 7 y 8 al comparar la ubicación de la tercera transición estimada con el cambio estructural encontrado a principios de los años ochenta por estos trabajos. Sin embargo, las consecuencias son diferentes, ya que nuestros resultados indican que tras esta perturbación en la economía existe un periodo de transición a lo largo de casi todo el decenio de los ochenta, mientras que la evidencia de cambios estructurales instantáneos implica que la economía pasó a una nueva tendencia inmediatamente después del cambio, ocurrido alrededor de 1982.

Finalmente, el cambio estructural o transición que estimamos en ambas variables en la primera mitad de los años cincuenta es el único que resulta instantáneo, por lo que nuestros resultados confirman el cambio estructural de mediados de siglo que encuentran Noriega, Soria y Velázquez (2008) y Noriega y Ramírez-Zamora (1999), representados en los Cuadros 7 y 8.

IV. Implicaciones para el crecimiento económico

Como se analiza en la Introducción, el uso de las funciones de transición para estimar cambios estructurales múltiples nos ayuda a distinguir entre etapas de crecimiento estacionario y fases de transición en la trayectoria del producto y producto per capita.

De acuerdo con el planteamiento de este estudio, los cambios en la pendiente de la tendencia lineal de la serie, corresponden a cambios en la tasa de crecimiento de largo plazo de la economía. En otras palabras, la tasa de crecimiento de yt está medida por los coeficientes relacionados con la variable de tendencia t, por lo que, según la especificación en (1), la tasa de crecimiento cambia de un valor inicial β0 a β0+j=1mβj. Así, cualquier evidencia estadística de algún β^j0 j = 1, ..., m, es interpretada como un cambio en la tasa de crecimiento de largo plazo de la economía. Como se observa en el Cuadro 4, existen tres cambios de pendiente estadísticamente significativos en cada una de las series analizadas.

En la Gráfica 8 se presenta la tasa de crecimiento observada y estimada para nuestra especificación preferida para el PIB real y real per capita de México. Específicamente, graficamos la tasa de crecimiento observada como la primera diferencia logarítmica de las series (línea discontinua en la Gráfica 8), mientras que la tasa de crecimiento estimada (línea continua en la Gráfica 8) es resultado de la suma de coeficientes relacionados con la pendiente de la tendencia lineal de la ecuación (1), β^0+j=1mβ^jStjτ^j,γ^j que son estadísticamente distintos de 0. Las áreas sombreadas muestran las etapas de transición estimadas que corresponden sólo a los cambios significativos en la pendiente. En la Gráfica 8 no se registran las transiciones de 1908-1916 en el caso del PIB real y 1923-1934 en el del PIB real per capita, puesto que estas transiciones implican sólo un cambio de nivel en la serie correspondiente, y no afectan la tasa de crecimiento de largo plazo. Por último, en la parte inferior de las gráficas se presenta en términos porcentuales la tasa de crecimiento de la etapa correspondiente junto al intervalo de confianza de dos desviaciones estándar.28

GRÁFICA 8 Transiciones y crecimiento económico 

A continuación procedemos a análizar con más detalle los periodos de transición y los cambios en la tasa de crecimiento económico encontrados para México. Primero, la Revolución no parece haber alterado la tasa de crecimiento del producto real, pero sí parece haber aumentado la del per capita. Este resultado ha sido analizado en varios estudios. Haber (1978), p. 416, señala: "La mayor parte de la planta manufacturera de México emergió de la lucha intacta; no fue destruida." Por su parte, en lo referente al trabajo en el periodo 1910-1920, Womack (1978) afirma que "Por muerte, emigración, conscripción o repliegue en zonas de refugio que garantizaban la subsistencia, la mano de obra disponible se redujo en general, pero no en los grandes centros industriales donde aumentó la oferta y demanda"; mientras que del sector industrial, dice, "Las plantas industriales sufrieron pocos daños físicos: la lucha le cerraba el camino a mercados lejanos, pero ampliaba los mercados inmediatos. La dificultad del transporte provocó una grave recesión en casi todas las zonas industriales del centro y del norte en 1913, y en la ciudad de México durante los dos años siguientes. Pero desde 1916 los principales centros fabriles empezaron a llevar constantemente su producción y para fines del decenio la mayoría había alcanzado nuevamente el nivel de 1910". Por su parte, Solís (1985), p. 72, señala que "En 1907 se inició una nueva crisis internacional, que se prolongó hasta 1910. Bajó el precio del henequén, se restringió el crédito, disminuyeron los préstamos y se aceleraron los cobros de vencimiento". Solís continúa: "El año de 1908 señaló una crisis de fuerte magnitud y el principio de una rápida decadencia del sistema bancario porfirista. A partir de entonces, sólo pudieron sostenerse a base de emisiones de billetes, fomentando una inflación peligrosa y agudizando la crisis económica por la cual atravesaba el país... [Además], Las escasas y poco confiables estadísticas [...] indican una sustancial caída de la vida económica de 1910 a 1915, en todos sus aspectos... [y sigue] Sin embargo, una vez derrotado el ejército federal, y con el auge de la demanda externa ocasionado por la primera Guerra Mundial, el valor en dólares de las exportaciones totales se recuperó rápidamente en 1916 y 1917." Estos hechos pueden ser visualizados en la Gráfica 8.

Segundo, a diferencia de la primera, la segunda transición estimada (1925-1934) afecta de manera positiva sólo la tasa de crecimiento del PIB real, incrementándose de 3.05% en el periodo 1895-1924 a 5.28% en 1935-1952, después de la transición. Ese periodo de transición hacia una tasa de crecimiento del producto superior coincide con la etapa que, por ejemplo, Villareal (1997) y Solís (1985) llaman "periodo de reconstrucción" o "formación institucional", ya que en ese periodo inmediato posterior a la Revolución Mexicana ocurrieron una serie de cambios institucionales que resultaron ser fundamentales para el proceso de crecimiento económico posterior. Por ejemplo, Solís (1985), p. 76, escribe: "El año de 1910 marca el inicio de una fuerte caída de la actividad económica a la que le sigue una etapa que podríamos llamar de formación institucional que se extiende de 1921 a 1935."

A partir de 1925 se emprende un ambicioso programa económico en el que se establecen las bases institucionales y materiales para el futuro crecimiento. Krauze (1977) analiza cómo en el periodo 1924-1928, en la presidencia de Plutarco Elías Calles, se conforma la reconstrucción económica de México. Las primeras acciones fueron la creación de la Comisión Nacional Bancaria (el 2 de enero de 1925), la Comisión Nacional de Caminos (el 30 de marzo de 1925), el Banco de México, S.A. (1 de septiembre de 1925) con el privilegio de ser el único banco emisor, entre otras atribuciones de banca central, y el 3 de diciembre la Comisión Nacional de Irrigación. Además, entra en vigor la ley del ISR en 1925, y en 1926 se expiden la Ley General de Instituciones de Crédito y Establecimientos Bancarios, la ley de Impuestos a la Minería y restituciones de tierra y aguas (1927), y el código civil (1928); además, se abre el Banco de Crédito Agrícola. En 1929 se crea el PNR (Partido Nacional Revolucionario) de la unión de la mayoría de los partidos existentes, anunciando el fin del caudillismo y el surgimiento de estabilidad política en el país; véase más detalles en Solís (1985) y Villareal (1997). Pero, por otra parte, cuando estaba en marcha este proceso de reconstrucción o formación sobrevino la gran depresión de 1929, la cual al parecer se manifiesta en las estimaciones como un cambio de nivel negativo en ambas series.29 De acuerdo a Solís (1985), la gran depresión afectó de una manera considerable la economía mexicana: provocó que el valor de las exportaciones mexicanas en 1932 fueran sólo la tercera parte de las de 1929, el ingreso del gobierno disminuyó 25% de 1930 a 1933, el gasto del gobierno decreció en 25% de 1930 a 1932; además, en 1931 México abandona el patrón oro y el tipo de cambio fluctúa libremente y en 1933 el peso sufre una devaluación de 80 por ciento.

Finalmente, los dos últimos cambios estructurales se presentan tanto en el PIB real como en el per capita. El segundo cambio de tasa de crecimiento en las series es el único instantáneo, que ocurre a principios de los años noventa. La tasa de crecimiento de largo plazo cambia positivamente alcanzando una tasa de 6.58% en el caso del PIB real y de 3.54% en el del per capita. En ambos casos esta etapa de crecimiento económico se extiende por un poco más de 20 años para finalizar en 1978. Al respecto, Solís (1985) documenta que este cambio marca el paso de un periodo de crecimiento con inflación (1935-1956) a un periodo de crecimiento con estabilidad de precios y de tipo de cambio (1956-1972) y, finalmente, un lapso con inflación y menor crecimiento. La tercera transición, 1979-1988 para el PIB real, y 1979-1989 para el PIB real per capita, conduce a una tasa de crecimiento menor que la obtenida en la etapa anterior: la tasa de crecimiento del PIB real se redujo de 6.58 a 3.01% y la del real per capita de 3.54 a 1.63%. La transición que precede a esta etapa de crecimiento claramente refleja el periodo conocido como "la década perdida" (véase por ejemplo Cárdenas, 2010) caracterizado principalmente por una alta inflación y estancamiento económico. La fecha de inicio de esta transición bien puede estar asociada con el auge petrolero y la subsecuente crisis de deuda de 1982.30

En síntesis, las pruebas presentadas sugieren que para el PIB real (PIB real per capita) parecen existir al menos cuatro etapas de crecimiento estacionario, 1895-1924, 1935-1952, 1956-1978 y 1989-2008 (1895-1908, 1916-1951, 1954-1978 y 1990-2008), separadas por tres épocas de transición, 1925-1934, 1953-1955 y 1979-1988 (1909-1915, 1952-1953 y 1979-1989).

Conclusiones

Este artículo examina la senda de crecimiento económico en México a lo largo del periodo 1895-2008. Los resultados parecen favorecer la hipótesis de que la tendencia de la producción en México se caracteriza por medio de una función lineal determinista, perturbada por cambios estructurales que coinciden con arreglos institucionales, guerras y crisis económicas y financieras. Los datos parecen no favorecer, por otra parte, la hipótesis de una tendencia estocástica.

Al modelar los cambios estructurales por medio de funciones de transición suave entre regímenes, somos capaces de distinguir empíricamente entre una fase de transición y una etapa de crecimiento sostenido. Nuestros resultados implican que para el PIB real y PIB real per capita de México existen cuatro etapas de crecimiento estacionario separadas por tres épocas de transición.

Uno de los hallazgos de este artículo es que la Revolución no afectó la tasa de crecimiento de largo plazo del producto. El descenso en la población, sin embargo, provocó el incremento de la producción per capita a partir de 1916. El cambio encontrado en la etapa de reconstrucción económica posterior a la Revolución Mexicana lleva a una tasa de crecimiento del producto que es compensada por el crecimiento de la población, de tal manera que no se manifiesta ningún efecto en el crecimiento del PIB per capita.

Consideramos que sería interesante estudiar con mayor detalle las fuentes del incremento en la tasa de crecimiento del producto per capita obtenido a principios de los años cincuenta, ya que, a diferencia del incremento en la tasa de crecimiento del producto per capita durante la Revolución, que está asociado a un decenso en la población, el aumento del producto per capita a mediados de siglo podría estar vinculado a un incremento en la productividad. Esta hipótesis amerita por supuesto un análisis minucioso, que rebasa las intenciones de este estudio, ya que dicho aumento llevó a la economía por cerca de 25 años a tener las mayores tasas de crecimiento del producto y producto per capita dentro de los 114 años considerados.

Por otra parte, el último cambio ocurrido en las tasas de crecimiento tanto del producto como del producto per capita es negativo y transcurre en un periodo de transición que refleja la llamada "década perdida". El periodo posterior a este cambio retoma una senda de crecimiento estacionario alrededor de una tendencia determinista que va de 1990 hasta el final de la muestra.

Referencias bibliográficas

Aguiar, M., y G. Gopinath (2007), "Emerging Market Business Cycles: The Cycle is the Trend", Journal of Political Economy, 115, 1, pp. 69-102. [ Links ]

Astorga, P., A. R. Bergés y V. Fitzgerald (2005), "Endogenous Growth and Exogenous Shocks in Latin America during the Twentieth Century", Discussion Papers in Economic and Social History, University of Oxford, núm. 57. [ Links ]

Bacon, D. W., y D. G. Watts (1971), "Estimating the Transition between Two Intersecting Straight Lines", Biometrika, 58, 3, pp. 25-34. [ Links ]

Bai, J., y P. Perron (1998), "Estimating and Testing Linear Models with Multiple Structural Changes", Econometrica, 66, 1, pp. 47-78. [ Links ]

______, y ______ (2003), "Computation and Analysis of Multiple Structural Change Models", Journal of Applied Econometrics, 18, pp. 1-22. [ Links ]

Barro, R. J., y J. F. Ursúa (2008), "Macroeconomic Crises since 1870", Brookings Papers on Economic Activity, 39, 1, pp. 255-350. [ Links ]

______, y ______ (2009), "Stock-Market Crashes and Depressions", Working Paper 14760, NBER Working paper series. [ Links ]

Ben-David, D., y D. H. Papell (1995), "The Great Wars, the Great Crash, and Steady State Growth: Some New Evidence about an Old Stylized Fact", Journal of Monetary Economics, 36, 3, pp. 453-475. [ Links ]

______, y ______ (1998), "Slowdowns and Meltdowns: Postwar Growth Evidence from 74 Countries", The Review of Economics and Statistics, 80, 4, pp. 561-571. [ Links ]

______, y ______ (2000), "Some Evidence on the Continuity of the Growth Process among the G-7 Countries", Economic Inquiry, 38, 2, pp. 320-330. [ Links ]

______, R. L. Lumsdaine, y D. H. Papell (2003), "Unit Roots, Postwar Slowdowns and Long-Run Growth: Evidence from Two Structural Breaks", Empirical Economics, 28, 2, pp. 303-319. [ Links ]

Caner, M., y B. Hansen (2001), "Threshold Autoregression with a Unit Root", Econometrica, 69, 6, pp. 1555-1596. [ Links ]

Cárdenas, E. (2010), "La reestructuración económica de 1982 a 1994", Clara García Ayluardo, Ignacio Marván Laborde y Erika Pani (coords.), Historia crítica de las modernizaciones de México, México, Fondo de Cultura Económica. [ Links ]

Castillo, R., y A. Díaz-Bautista (2002), "Testing for Unit Roots in Mexico's GDP", Momento Económico, 124, pp. 2-10. [ Links ]

Chang, T., C. C. Nieh y C. C. Wei (2005), "Is Per Capita Real GDP Stationary? Evidence from Selected African Countries Based On more Powerful Nonlinear (logistic) Unit Root Tests", Economics Bulletin, 3, 24, pp. 1-9. [ Links ]

Darne, O. (2009), "The Uncertain Unit Root in Real GNP: A Re-Examination", Journal of Macroeconomics, 31, 1, pp. 153-166. [ Links ]

Dickey, D. A., y S. G. Pantula (1987), "Determining the Order of Differencing in Autoregressive Processes", Journal of Business and Economic Statistics, 5, 4, pp. 455-461. [ Links ]

______, y W. A. Fuller (1979), "Distribution of the Estimators for Autoregressive Time Series with a Unit Root", Journal of the American Statistical Association, 74, pági- nas 427-431. [ Links ]

Diebold, F. X., y A. S. Senhadji (1996), "The Uncertain Unit Root in Real GNP: Comment", The American Economic Review, 86, 5, pp. 1291-1298. [ Links ]

Duck, N. W. (1992), "UK Evidence on Breaking Trend Functions", Oxford Economic Papers, 44, 3, pp. 426-439. [ Links ]

Foster, N., y R. Stehrer (2007), "Modelling Transformation in CEECs using Smooth Transitions", Journal of Comparative Economics, 35, 1, pp. 57-86. [ Links ]

Granger, C. W. J., y T. Teräsvirta (1993), Modelling Nonlinear Economic Relationships, Oxford University Press. [ Links ]

Greenaway, D., S. Leybourne y D. Sapsford (1997), "Modeling Growth (and Liberalization) using Smooth Transitions Analysis", Economic Inquiry, 35, 4, pp. 798-814. [ Links ]

______, ______, y ______ (2000), "Smooth Transitions and GDP Growth in the European Union", Manchester School, 68, 2, pp. 145-65. [ Links ]

Haber, S. (1978), "La revolución y la industria manufacturera mexicana, 1910-1925", E. Cárdenas, Historia Económica de México, vol. 3, Serie de Lecturas 64 de El Trimestre Económico, México, Fondo de Cultura Económica, pp. 415-446. [ Links ]

Harvey, D. I., y T. C. Mills (2000), "Detecting Climate Change using Double Smooth Transitions", Loughborough University Economic Research Paper, núm. 00/21. [ Links ]

______, y ______ (2002), "Unit Roots and Double Smooth Transitions", Journal of Applied Statistics, 29, 5, pp. 675-683. [ Links ]

______, y ______ (2004), "Tests for Stationarity in Series with Endogenously Determined Structural Change", Oxford Bulletin of Economics & Statistics, 66, 5, pp. 863-894. [ Links ]

Kapetanios, G. (2005), "Unit-Root Testing Against the Alternative Hypothesis of up to m Structural Breaks", Journal of Time Series Analysis, 26, 1, pp. 123-133. [ Links ]

______, y Y. Shin (2006), "Unit Root Tests in Three-Regime SETAR Models", The Econometrics Journal, 9, pp. 252-278. [ Links ]

______, y Snell, A. y Y. Shin (2003), "Testing for a Unit Root in the Nonlinear STAR Framework", Journal of Econometrics, 112, pp. 359-379. [ Links ]

Kejriwal, M., y P. Perron (2010), "A Sequential Procedure to Determine the Number of Breaks in Trend with an Integrated or Stationary Noise Component", Journal of Time Series Analysis, 31, pp. 305-328. [ Links ]

Krauze, E. (1977), Historia de la Revolución Mexicana. La reconstrucción económica, México, El Colegio de México. [ Links ]

Leybourne, S., y P. Mizen (1999), "Understanding the Disinflations in Australia, Canada and New Zealand using Evidence from Smooth Transition Analysis", Journal of International Money and Finance, 18, 5, pp. 799-816. [ Links ]

______, P. Newbold y D. Vougas (1998), "Unit Roots and Smooth Transitions", Journal of Time Series Analysis, 19, 1, pp. 83-97. [ Links ]

Lin, C. J., y T. Teräsvirta (1994), "Testing the Constancy of Regression Parameters Against Continuous Structural Change", Journal of Econometrics, 62, páginas 211-228. [ Links ]

Maddala, G. S. (1977), Econometrics Nueva York, McGraw-Hill. [ Links ]

Mills, T. C (1991), "Are Fluctuations in U.K. Ouput Transitory or Permanent?", The Manchester School, 59, 1, pp. 1-11. [ Links ]

Murray, C., y C. Nelson (2000), "The Uncertain Trend in U.S. GDP", Journal of Monetary Economics, 46, pp. 79-95. [ Links ]

Newbold, P., S. Leybourne, R. Sollis y M. E. Wohar (2001), "U.S. and U.K. Interest Rates 1890 - 1934: New Evidence on Structural Breaks", Journal of Money, Credit and Banking, 33, 2, pp. 235-250. [ Links ]

Ng, S., y P. Perron (1995), "Unit Root Tests in ARMA Models with Data-Dependent Methods for the Selection of the Truncation Lag", Journal of the American Statistical Association, 90, pp. 268-281. [ Links ]

______, y ______ (2001), "Lag Length Selection and the Construction of Unit Root Tests with Good Size and Power", Econometrica, 69, 6, pp. 1519-1554. [ Links ]

Noriega, A. (2004), "Long-run Monetary Neutrality and the Unit Root Hypothesis: Further International Evidence", North American Journal of Economics and Finance, 15, 2, pp. 179-197. [ Links ]

______, y A. Ramírez-Zamora (1999), "Unit Roots and Multiple Structural Breaks in Real Output: How Long Does an Economy Remain Stationary?", Estudios Económicos, 14, 2, pp. 163-188. [ Links ]

______, y E. de Alba (2001), "Stationarity and Structural Breaks: Evidence from Classical and Bayesian Approaches", Economic Modelling, 18, 4, pp. 503-524. [ Links ]

______, L. M. Soria y R. Velázquez (2008), "International Evidence on Stochastic and Deterministic Monetary Neutrality", Economic Modelling, 25, 6, pp. 1261-1275. [ Links ]

Pantula, S. G. (1989), "Testing for Unit Roots in Time Series", Econometric Theory, 5, pp. 56-271. [ Links ]

Papell, D., y R. Prodan (2004), "The Uncertain Unit Root in US real GDP: Evidence with Restricted and Unrestricted Structural Change", Journal of Money, Credit and Banking, 36, pp. 423-427. [ Links ]

Perron, P. (1989), "The Great Crash, the Oil Price Shock, and the Unit Root Hypothesis", Econometrica, 57, 6, pp. 1361-1401. [ Links ]

_______, y Z. Qu (2007), "A Simple Modification to Improve the Finite Sample Properties of Ng and Perron's Unit Root Tests", Economics Letters, 94, 1, pp. 12-19. [ Links ]

Rudebush, G. D. (1992), "Trends and Random Walks in Macroeconomic Time Series: a Re-Examination", International Economic Review, 33, 3, pp. 661-680. [ Links ]

______ (1993), "The Uncertain Unit Root in Real GNP", American Economic Review, 83, 1, pp. 264-272. [ Links ]

Solís, L. (1985), La realidad económica mexicana: Retrovisión y perspectivas, Méxco, Siglo XXI Editores, vigésima edición. [ Links ]

Sollis, R. (2004), "Asymmetric Adjustment and Smooth Transitions: A Combination of some Unit Root Tests", Journal of Time Series Analysis, 25, 3, pp. 409-417. [ Links ]

Sollis, R., S. Leybourne y P. Newbold (1999), "Unit Roots and Asymmetric Smooth Transitions", Journal of Time Series Analysis, 20, 6, pp. 671-677. [ Links ]

Villareal, R. (1997), Industrialización, Deuda y desequilibrio externo en México. Un enfoque neoestructuralista (1929-1997), México, Fondo de Cultura Económica, 3a edición. [ Links ]

Vougas, D. (2006), "On Unit Root Testing with Smooth Transitions", Computational Statistics and Data Analysis, 51, pp. 797-800. [ Links ]

______ (2007), "Is the Trend in Post-WW II U.S. Real GDP Uncertain or Nonlinear?", Economics Letters, 94, pp. 348-355. [ Links ]

Womack, J. (1978), "La economía en la Revolución (1910-1920). Historiografía y análisis", E. Cárdenas, Historia Económica de México, vol. 3, Serie de Lecturas 64 de El Trimestre Económico, México, Fondo de Cultura Económica, pp. 391-414. [ Links ]

*Agradecemos los comentarios y sugerencias de Carlos Capistrán, Gonzalo Rangel y un dictaminador anónimo de El Trimestre Económico, así como de los participantes en el seminario del Banco de México. Todos los errores son responsabilidad nuestra. Las opiniones en este artículo son de los autores y no forzosamente reflejan el punto de vista del Banco de México.

1En relación con la cuestión de si la función de tendencia de largo plazo de la producción real de los Estados Unidos debería ser modelada como estacionaria en tendencia o estacionaria en diferencias véase Rudebusch (1993), Diebold y Senhadji (1996), Murray y Nelson (2000), Papell y Prodan (2004), Vougas (2007) y Darne (2009).

2Véase más pormenores de estos resultados en Noriega, Soria y Velázquez (2008).

3Existen en la bibliografía otros enfoques que prueban la hipótesis nula de raíz unitaria contra opciones no lineales, en las cuales la serie de tiempo sigue un proceso globalmente estacionario, tales como las de los modelos threshold autoregressive (TAR) en Caner y Hansen (2001), smooth transition autoregressive (STAR) en Kapetanios, Snell y Shin (2003) y self-exciting threshold autoregressive (SETAR) en Kapetanios y Shin (2006). Incluso existen pruebas que combinan estos dos tipos de hipótesis como en Sollis (2004), en el que se prueba la nula raíz unitaria contra una opción que combina al modelo TAR con las transiciones suaves entre tendencias lineales deterministas. Finalmente, nótese que aquí usamos un enfoque univariado para determinar el orden de integración del producto. Aguiar y Gopinath (2007), por ejemplo, desde una perspectiva multivariada, analizan cómo la magnitud de la reacción del consumo a choques en el ingreso pueden utilizarse para identificar la persistencia de dichos choques.

4LNV explican la idea de una transición suave en una serie de tiempo económica con el argumento de que en muchas situaciones los cambios en variables económicas agregadas, como el PIB, son influidos por los cambios en el comportamiento de un número grande de agentes económicos, y es poco probable que todos los agentes individuales reaccionen simultáneamente a cierto estímulo económico: mientras que unos están listos para reaccionar instantáneamente otros son propensos a diferentes grados de inercia institucional (dependiendo, por ejemplo, de la eficiencia de los mercados en que éstos actúan), y por consecuencia reaccionarán en diferentes momentos.

5Otros trabajos en los cuales se emplea la función logística de transición suave, en que la variable de transición es el tiempo, son: Leybourne y Mizen (1999), Harvey y Mills (2000), Greenaway, Leybourne y Sapsford (2000), Newbold et al (2001), Sollis (2004), Chang, Nieh y Wei (2005), Vougas (2006), Vougas (2007), y Foster y Stehrer (2007).

6 Harvey y Mills (2004) plantean y desarrollan la prueba de estacionariedad contra la opción de raíz unitaria. Con la hipótesis nula se supone estacionariedad respecto a una tendencia lineal con una transición suave simétrica.

7Los resultados empíricos que presentan Harvey y Mills (2002), resultan precisamente en este tipo de estimaciones, en las cuales los periodos de transición se sobreponen, lo cual tiene la consecuencia de que una etapa de transición empiece antes de que haya terminado la anterior.

8Para los resultados empíricos de este artículo utilizamos h = 5. Otros trabajos que utilizan este valor (o valores similares) son Bai y Perron (1998), Noriega y Ramírez-Zamora (1999), Ben-David y Papell (2000), Bai y Perron (2003), Kapetanios (2005), y Noriega, Soria y Velázquez (2008).

9Por sencillez en la exposición, este análisis se presenta sólo para el caso de la función logística generalizada, Gtj , pero también es válido para el caso de la función simétrica, Stj .

10Por ejemplo, Newbold et al (2001) en su análisis consideran una duración aproximadamente de 90% de la transición.

11Por sencillez en la exposición, el problema de optimación de la ecuación (4) se presenta sólo para el caso de la función logística generalizada, Gtj . El mismo problema se debe resolver cuando se utiliza la función simétrica, Stj .

12Se deja el máximo número de iteraciones posibles en 5 mil para que el resultado converja; el criterio de convergencia es que la norma de la variación de los parámetros relativa a la norma de los valores de los parámetros de la última iteración sea inferior a 0.001, que son los valores predeterminados en Eviews 7.

13La rutina es ejecutada en una computadora con un procesador Intel(R) Core(TM)2 CPU 6400 de 2.13 GHz y 1.99 GB RAM.

14Aunado a esto, además del número de funciones de transición por estimar, el costo computacional también depende de la rapidez de convergencia de la estimación no lineal de los parámetros.

15El modelo ETm se presenta para la función logística generalizada Gtj . Para el caso de la función simétrica, simplemente se sustituye Gtj por Stj .

16Nótese que el estadístico estimado a partir de (5) es invariante respecto a componentes deterministas, por lo que la inclusión o no de una constante en (7) no afecta el resultado de la inferencia. De hecho, los resultados que se presentan más adelante en el Cuadro 2, respecto a la inferencia en la presencia de una raíz unitaria, no cambian si eliminamos de (7) el parámetro. Sobre este tema véase Perron (1989), p. 1393.

17Este procedimiento se utiliza también en Noriega, Soria y Velázquez (2008).

18Este procedimiento está basado en los resultados de Dickey y Pantula (1987), en el que se observa empíricamente que la probabilidad de rechazar la hipótesis nula de una raíz unitaria (H1 ) contra la opción de estacionariedad (H0 ) se incrementa con el número de raíces unitarias presentes.

19Mediante una inspección gráfica y estimaciones previas se verificó si las pruebas deberían incluir como componente determinista constante o constante y tendencia.

20Ninguno de los resultados en el Cuadro 2 se modifica, si en vez de elegir el rezago óptimo, k*, de acuerdo con el nivel de significancia de 10%, se elige utilizando algunos de los criterios de información de Akaike, Schwarz o Hannan-Quinn.

21Hay que advertir que la medición de los valores en que inicia, finaliza y cambia de concavidad la función de transición estimada, no forzosamente arroja valores enteros. Con fines de una mejor interpretación de las transiciones en el Cuadro 3, se presenta los años correspondientes a los periodos después de ser redondeados.

22Traducción nuestra.

23Esto puede ser resultado de un sesgo debido a una mala especificación del modelo, es decir, de permitir menos cambios estructurales de los que realmente se encuentran presentes en las series, o sea, sesgo por omisión de variable.

24En el Apéndice 2 se presenta todos los resultados de la estimación de la ecuación (4), para las ocho distintas especificaciones de cada variable. Adviértase en el Cuadro A2 que los tres criterios de información presentados (AIC, BIC y HQC) coinciden en "escoger" la especificación con cuatro transiciones suaves simétricas, tanto para el PIB real como para el per capita. Tanto estos resultados como los analizados en el párrafo anterior nos permiten concluir que no debería de existir el problema de sobreajuste (overfitting).

25Otra posible opción sería plantear una especificación que combine a las dos funciones de transición Gt y St tal que algunas transiciones sean simétricas mientras otras permitan la posibilidad de asimetría. Sin embargo, esta posibilidad se relega para una investigación futura debido a los costos computacionales que son analizados en la sección de metodología, puesto que en el caso de considerar cuatro transiciones se tendrían que estimar 24 especificaciones para poder elegir con algún criterio la mejor especificación.

26Note que, por ejemplo, a pesar de que la transición presentada para ambas series entre los decenios de los veinte y treinta, pudiera parecer excesivamente larga, y de hecho ser interpretada como un régimen de crecimiento negativo en sí mismo, como se argumenta líneas abajo, el periodo 1921-1935 está marcado por un conjunto de reformas institucionales orientadas a la reconstrucción económica de México. Más bien, fue precisamente este periodo de transición (reconstrucción) el que permitió alcanzar una tasa de crecimiento significativamente mayor en el producto real a partir de mediados de los años treinta.

27Por otra parte, las series de los residuales, ilustrados en las Gráficas 5 y 6, no presentan pruebas de problemas de heteroscedasticidad. Los resultados de las pruebas de diagnóstico a las cuales fueron sometidos los residuales (autocorrelación, normalidad, heteroscedasticidad condicional e incondicional) pueden ser solicitados a los autores.

28Estas tasas podrían interpretarse como las tasas de crecimiento del estado estacionario: véase más pormenores en Ben-David y Papell (1995).

29Nótese que en la Gráfica 8 la transición estimada del PIB real para el periodo 1925-1934 no parece corresponder a lo que se observa en los datos. Esto es debido precisamente a la caída de esta serie durante ese periodo, como se observa claramente en la Gráfica 5. Nótese que el mismo argumento se aplica a la transición entre 1909 y 1915 para el PIB per capita.

30Con base en los descubrimientos de yacimientos de 1978 y el aumento de los precios del petróleo, entre 1978 y 1981, el PIB real tuvo un crecimiento promedio superior a 9%. Sin embargo, esta estrategia basada en el petróleo finalizó con una caída en el precio del petróleo en 1981, dejando al país con una enorme deuda externa, la cual había sido contratada para desarrollar la industria petrolera. En 1982 el aumento de las tasas de interés internacionales dejó a México en una grave crisis de deuda. El PIB real sufrió un decremento de 0.6% y el PIB real per capita de 2.5%, mientras que la inflación anual llegó a 98.8%. Véase más detalles en, por ejemplo, Solís (1985) y Villareal (1997).

31La base de datos construida por Barro y Ursúa (2008) hasta 2006 se encuentra disponible en línea http://www.economics.harvard.edu/facultary/barro/data_sets_barro. En Barro y Úrzua (2008) pueden encontrarse todos los detalles de cómo se calculan los datos del PIB para el periodo 1911-1920.

Apéndice

1. Datos

En Barro y Ursúa (2008) se explica de manera pormenorizada la forma en que se construye la serie del PIB real para el periodo 1895-2006.31 Para actualizar la información hasta 2008, y siendo congruentes con las fuentes presentadas por Barro y Ursúa, se toman los datos del PIB trimestral registrados por el Instituto Nacional de Estadística Geografía e Informática (a precios de 2003) como el promedio del periodo. Para la serie de población los datos de 2007 y 2008 son tomados de la base de datos del World Development Indicators del Banco Mundial. En el siguiente cuadro se presentan los datos en logaritmos.

CUADRO A1 PIB real y real per capita (en logaritmos) de México (1895-2008) 

2. Resultados de la estimación de transiciones múltiples

En el siguiente cuadro se presenta los resultados de las ocho estimaciones realizadas tanto para el PIB real como para el PIB real per capita.

a EE, error estándar de la regresión, SRC, suma de residuales al cuadrado; AIC, BIC, HQC, se refieren a los criterios de información de Akaike, Schwarz y Hannan-Quinn, respectivamente.

***, ** y * rechazo de la hipótesis nula al nivel de 1, 5 y 10%, respectivamente.

CUADRO A2 Resultados de la estimación para el PIB real y real per capita a 

3. Evidencia de cambios estructurales en México

En el Cuadro A3 se presenta la evidencia existente en la bibliografía de cambios estructurales instantáneos para las series del PIB real y real per capita de México, respectivamente.

a Noriega y Ramirez-Zamora emplean dos metodologías para determinar el número de cambios estructurales en las series: (PC) parámetros constantes y (RRU) rechazo de raíz unitaria.

b Noriega y de Alba comparan los resultados de estimar un único cambio estructural a través de métodos clásico (C) y Bayesiano (B).

c Para probar la estabilidad de los parámetros en un modelo autorregresivo de las series, Astorga, Bergés y Fitzgerald (2005) aplican la prueba de Chow y la prueba de Chow de pronóstico (esta última se presenta entre paréntesis).

CUADRO A3 Evidencia de cambios estructurales (endógenos) para el PIB real y PIB real per capita 

Recibido: 02 de Septiembre de 2010; Aprobado: 14 de Marzo de 2011

Creative Commons License Este es un artículo publicado en acceso abierto bajo una licencia Creative Commons