Factor models, VARMA processes and parameter instability with applications in macroeconomics

Avec les avancements de la technologie de l'information, les données temporelles économiques et financières sont de plus en plus disponibles. Par contre, si les techniques standard de l'analyse des séries temporelles sont utilisées, une grande quantité d'information est accompagnée du problème de dimensionnalité. Puisque la majorité des séries d'intérêt sont hautement corrélées, leur dimension peut être réduite en utilisant l'analyse factorielle. Cette technique est de plus en plus populaire en sciences économiques depuis les années 90. Étant donnée la disponibilité des données et des avancements computationnels, plusieurs nouvelles questions se posent. Quels sont les effets et la transmission des chocs structurels dans un environnement riche en données? Est-ce que l'information contenue dans un grand ensemble d'indicateurs économiques peut aider à mieux identifier les chocs de politique monétaire, à l'égard des problèmes rencontrés dans les applications utilisant des modèles standards? Peut-on identifier les chocs financiers et mesurer leurs effets sur l'économie réelle? Peut-on améliorer la méthode factorielle existante et y incorporer une autre technique de réduction de dimension comme l'analyse VARMA? Est-ce que cela produit de meilleures prévisions des grands agrégats macroéconomiques et aide au niveau de l'analyse par fonctions de réponse impulsionnelles? Finalement, est-ce qu'on peut appliquer l'analyse factorielle au niveau des paramètres aléatoires? Par exemple, est-ce qu'il existe seulement un petit nombre de sources de l'instabilité temporelle des coefficients dans les modèles macroéconomiques empiriques? Ma thèse, en utilisant l'analyse factorielle structurelle et la modélisation VARMA, répond à ces questions à travers cinq articles. Les deux premiers chapitres étudient les effets des chocs monétaire et financier dans un environnement riche en données. Le troisième article propose une nouvelle méthode en combinant les modèles à facteurs et VARMA. Cette approche est appliquée dans le quatrième article pour mesurer les effets des chocs de crédit au Canada. La contribution du dernier chapitre est d'imposer la structure à facteurs sur les paramètres variant dans le temps et de montrer qu'il existe un petit nombre de sources de cette instabilité. Le premier article analyse la transmission de la politique monétaire au Canada en utilisant le modèle vectoriel autorégressif augmenté par facteurs (FAVAR). Les études antérieures basées sur les modèles VAR ont trouvé plusieurs anomalies empiriques suite à un choc de la politique monétaire. Nous estimons le modèle FAVAR en utilisant un grand nombre de séries macroéconomiques mensuelles et trimestrielles. Nous trouvons que l'information contenue dans les facteurs est importante pour bien identifier la transmission de la politique monétaire et elle aide à corriger les anomalies empiriques standards. Finalement, le cadre d'analyse FAVAR permet d'obtenir les fonctions de réponse impulsionnelles pour tous les indicateurs dans l'ensemble de données, produisant ainsi l'analyse la plus complète à ce jour des effets de la politique monétaire au Canada. Motivée par la dernière crise économique, la recherche sur le rôle du secteur financier a repris de l'importance. Dans le deuxième article nous examinons les effets et la propagation des chocs de crédit sur l'économie réelle en utilisant un grand ensemble d'indicateurs économiques et financiers dans le cadre d'un modèle à facteurs structurel. Nous trouvons qu'un choc de crédit augmente immédiatement les diffusions de crédit (credit spreads), diminue la valeur des bons de Trésor et cause une récession. Ces chocs ont un effet important sur des mesures d'activité réelle, indices de prix, indicateurs avancés et financiers. Contrairement aux autres études, notre procédure d'identification du choc structurel ne requiert pas de restrictions temporelles entre facteurs financiers et macroéconomiques. De plus, elle donne une interprétation des facteurs sans restreindre l'estimation de ceux-ci. Dans le troisième article nous étudions la relation entre les représentations VARMA et factorielle des processus vectoriels stochastiques, et proposons une nouvelle classe de modèles VARMA augmentés par facteurs (FAVARMA). Notre point de départ est de constater qu'en général les séries multivariées et facteurs associés ne peuvent simultanément suivre un processus VAR d'ordre fini. Nous montrons que le processus dynamique des facteurs, extraits comme combinaison linéaire des variables observées, est en général un VARMA et non pas un VAR comme c'est supposé ailleurs dans la littérature. Deuxièmement, nous montrons que même si les facteurs suivent un VAR d'ordre fini, cela implique une représentation VARMA pour les séries observées. Alors, nous proposons le cadre d'analyse FAVARMA combinant ces deux méthodes de réduction du nombre de paramètres. Le modèle est appliqué dans deux exercices de prévision en utilisant des données américaines et canadiennes de Boivin, Giannoni et Stevanovic (2010, 2009) respectivement. Les résultats montrent que la partie VARMA aide à mieux prévoir les importants agrégats macroéconomiques relativement aux modèles standards. Finalement, nous estimons les effets de choc monétaire en utilisant les données et le schéma d'identification de Bernanke, Boivin et Eliasz (2005). Notre modèle FAVARMA(2,1) avec six facteurs donne les résultats cohérents et précis des effets et de la transmission monétaire aux États-Unis. Contrairement au modèle FAVAR employé dans l'étude ultérieure où 510 coefficients VAR devaient être estimés, nous produisons les résultats semblables avec seulement 84 paramètres du processus dynamique des facteurs. L'objectif du quatrième article est d'identifier et mesurer les effets des chocs de crédit au Canada dans un environnement riche en données et en utilisant le modèle FAVARMA structurel. Dans le cadre théorique de l'accélérateur financier développé par Bernanke, Gertler et Gilchrist (1999), nous approximons la prime de financement extérieur par les credit spreads. D'un côté, nous trouvons qu'une augmentation non-anticipée de la prime de financement extérieur aux États-Unis génère une récession significative et persistante au Canada, accompagnée d'une hausse immédiate des credit spreads et taux d'intérêt canadiens. La composante commune semble capturer les dimensions importantes des fluctuations cycliques de l'économie canadienne. L'analyse par décomposition de la variance révèle que ce choc de crédit a un effet important sur différents secteurs d'activité réelle, indices de prix, indicateurs avancés et credit spreads. De l'autre côté, une hausse inattendue de la prime canadienne de financement extérieur ne cause pas d'effet significatif au Canada. Nous montrons que les effets des chocs de crédit au Canada sont essentiellement causés par les conditions globales, approximées ici par le marché américain. Finalement, étant donnée la procédure d'identification des chocs structurels, nous trouvons des facteurs interprétables économiquement. Le comportement des agents et de l'environnement économiques peut varier à travers le temps (ex. changements de stratégies de la politique monétaire, volatilité de chocs) induisant de l'instabilité des paramètres dans les modèles en forme réduite. Les modèles à paramètres variant dans le temps (TVP) standards supposent traditionnellement les processus stochastiques indépendants pour tous les TVPs. Dans cet article nous montrons que le nombre de sources de variabilité temporelle des coefficients est probablement très petit, et nous produisons la première évidence empirique connue dans les modèles macroéconomiques empiriques. L'approche Factor-TVP, proposée dans Stevanovic (2010), est appliquée dans le cadre d'un modèle VAR standard avec coefficients aléatoires (TVP-VAR). Nous trouvons qu'un seul facteur explique la majorité de la variabilité des coefficients VAR, tandis que les paramètres de la volatilité des chocs varient d'une façon indépendante. Le facteur commun est positivement corrélé avec le taux de chômage. La même analyse est faite avec les données incluant la récente crise financière. La procédure suggère maintenant deux facteurs et le comportement des coefficients présente un changement important depuis 2007. Finalement, la méthode est appliquée à un modèle TVP-FAVAR. Nous trouvons que seulement 5 facteurs dynamiques gouvernent l'instabilité temporelle dans presque 700 coefficients.

A multiple-building diffraction attenuation function expressed in terms of UTD coefficients for microcellular communications

This work presents an explicit formulation for multiple- edge diffraction for mobile radiowave propagation in terms of uniform theory of diffraction (UTD) coefficients when a spherical incident wave is considered. This solution can be used in an UTD context and sharply reduces the computing time over existing formulation. Results can be applied in the planning of microcellular systems

Evanescent wave sensor to monitor the etching rate of an optical fibre

The design and development of an evanescent wave sensor to determine the etching rate of the core of an optical fibre is discussed in this paper. The working of the device is based on the principle of propagation and loss of the evanescent wave in the cladding region of the fibre. The fraction of light intensity creeping out of the core of an uncladded fibre is a function of the core radius. As this radius decreases, the evanescent wave coupling to the medium surrounding the core enhances. This results in a decrease of the transmitted light intensity through the fibre. This technique is useful to design and fabricate optical fibres with different core geometries.

Spatial and temporal characteristics of rain intensity in the peninsular Malaysia using TRMM rain rate

The present study is focused on the intensity distribution of rainfall in different classes and their contribution to the total seasonal rainfall. In addition, we studied the spatial and diurnal variation of the rainfall in the study areas. For the present study, we retrieved data from TRMM (Tropical Rain Measuring Mission) rain rate available in every 3 h temporal and 25 km spatial resolutions. Moreover, station rainfall data is used to validate the TRMM rain rate and found significant correlation between them (linear correlation coefficients are 0.96, 0.85, 0.75 and 0.63 for the stations Kota Bharu, Senai, Cameron highlands and KLIA, respectively). We selected four areas in the Peninsular Malaysia and they are south coastal, east coastal, west coastal and highland regions. Diurnal variation of frequency of rain occurrence is different for different locations. We noticed bimodal variation in the coastal areas in most of the seasons and unimodal variation in the highland/inland area. During the southwest monsoon period in the west coastal stations, there is no distinct diurnal variation. The distribution of different intensity classes during different seasons are explained in detail in the results

An Improved Technique for Evolving Wavelet Coefficients for Fingerprint Image Compression

In this paper, an improved technique for evolving wavelet coefficients refined for compression and reconstruction of fingerprint images is presented. The FBI fingerprint compression standard [1, 2] uses the cdf 9/7 wavelet filter coefficients. Lifting scheme is an efficient way to represent classical wavelets with fewer filter coefficients [3, 4]. Here Genetic algorithm (GA) is used to evolve better lifting filter coefficients for cdf 9/7 wavelet to compress and reconstruct fingerprint images with better quality. Since the lifting filter coefficients are few in numbers compared to the corresponding classical wavelet filter coefficients, they are evolved at a faster rate using GA. A better reconstructed image quality in terms of Peak-Signal-to-Noise-Ratio (PSNR) is achieved with the best lifting filter coefficients evolved for a compression ratio 16:1. These evolved coefficients perform well for other compression ratios also.

The counter-propagating Rossby-wave perspective on baroclinic instability. I: Mathematical basis

It is shown that Bretherton's view of baroclinic instability as the interaction of two counter-propagating Rossby waves (CRWs) can be extended to a general zonal flow and to a general dynamical system based on material conservation of potential vorticity (PV). The two CRWs have zero tilt with both altitude and latitude and are constructed from a pair of growing and decaying normal modes. One CRW has generally large amplitude in regions of positive meridional PV gradient and propagates westwards relative to the flow in such regions. Conversely, the other CRW has large amplitude in regions of negative PV gradient and propagates eastward relative to the zonal flow there. Two methods of construction are described. In the first, more heuristic, method a ‘home-base’ is chosen for each CRW and the other CRW is defined to have zero PV there. Consideration of the PV equation at the two home-bases gives ‘CRW equations’ quantifying the evolution of the amplitudes and phases of both CRWs. They involve only three coefficients describing the mutual interaction of the waves and their self-propagation speeds. These coefficients relate to PV anomalies formed by meridional fluid displacements and the wind induced by these anomalies at the home-bases. In the second method, the CRWs are defined by orthogonality constraints with respect to wave activity and energy growth, avoiding the subjective choice of home-bases. Using these constraints, the same form of CRW equations are obtained from global integrals of the PV equation, but the three coefficients are global integrals that are not so readily described by ‘PV-thinking’ arguments. Each CRW could not continue to exist alone, but together they can describe the time development of any flow whose initial conditions can be described by the pair of growing and decaying normal modes, including the possibility of a super-modal growth rate for a short period. A phase-locking configuration (and normal-mode growth) is possible only if the PV gradient takes opposite signs and the mean zonal wind and the PV gradient are positively correlated in the two distinct regions where the wave activity of each CRW is concentrated. These are easily interpreted local versions of the integral conditions for instability given by Charney and Stern and by Fjørtoft.

Tissue culture propagation alters plant-microbe interactions in tobacco rhizosphere

We have compared properties of roots from different lines (genotypes) of tobacco raised either in tissue culture or grown from seed. The different lines included unmodified plants and plants modified to express reduced activity of the enzyme cinnamoyl-CoA reductase, which has a pivotal role in lignin biosynthesis. The size and structure of the rhizosphere microbial community, characterized by adenosine triphosphate and phospholipid fatty acid analyses, were related to root chemistry (specifically the soluble carbohydrate concentration) and decomposition rate of the roots. The root material from unmodified plants decomposed faster following tissue culture compared with seed culture, and the faster decomposing material had significantly higher soluble carbohydrate concentrations. These observations are linked to the larger microbial biomass and greater diversity of the rhizosphere communities of tissue culture propagated plants.

Tissue culture propagation alters plant-microbe interactions in tobacco rhizosphere

We have compared properties of roots from different lines (genotypes) of tobacco raised either in tissue culture or grown from seed. The different lines included unmodified plants and plants modified to express reduced activity of the enzyme cinnamoyl-CoA reductase, which has a pivotal role in lignin biosynthesis. The size and structure of the rhizosphere microbial community, characterized by adenosine triphosphate and phospholipid fatty acid analyses, were related to root chemistry (specifically the soluble carbohydrate concentration) and decomposition rate of the roots. The root material from unmodified plants decomposed faster following tissue culture compared with seed culture, and the faster decomposing material had significantly higher soluble carbohydrate concentrations. These observations are linked to the larger microbial biomass and greater diversity of the rhizosphere communities of tissue culture propagated plants.

Reactive uptake coefficients for heterogeneous reaction of N2O5 with submicron aerosols of NaCl and natural sea salt

The kinetics of uptake of gaseous N2O5 on submicron aerosols containing NaCl and natural sea salt have been investigated in a flow reactor as a function of relative humidity (RH) in the range 30-80% at 295±2K and a total pressure of 1bar. The measured uptake coefficients, γ, were larger on the aerosols containing sea salt compared to those of pure NaCl, and in both cases increased with increasing RH. These observations are explained in terms of the variation in the size of the salt droplets, which leads to a limitation in the uptake rate into small particles. After correction for this effect the uptake coefficients are independent of relative humidity, and agree with those measured previously on larger droplets. A value of γ=0.025 is recommended for the reactive uptake coefficient for N2O5 on deliquesced sea salt droplets at 298K and RH>40%.

A block coprocessor for user data rate improvements to GPRS coding scheme 4

The general packet radio service (GPRS) has been developed to allow packet data to be transported efficiently over an existing circuit-switched radio network, such as GSM. The main application of GPRS are in transporting Internet protocol (IP) datagrams from web servers (for telemetry or for mobile Internet browsers). Four GPRS baseband coding schemes are defined to offer a trade-off in requested data rates versus propagation channel conditions. However, data rates in the order of > 100 kbits/s are only achievable if the simplest coding scheme is used (CS-4) which offers little error detection and correction (EDC) (requiring excellent SNR) and the receiver hardware is capable of full duplex which is not currently available in the consumer market. A simple EDC scheme to improve the GPRS block error rate (BLER) performance is presented, particularly for CS-4, however gains in other coding schemes are seen. For every GPRS radio block that is corrected by the EDC scheme, the block does not need to be retransmitted releasing bandwidth in the channel and improving the user's application data rate. As GPRS requires intensive processing in the baseband, a viable field programmable gate array (FPGA) solution is presented in this paper.

Multispacecraft observation of magnetic cloud erosion by magnetic reconnection during propagation

During propagation, Magnetic Clouds (MC) interact with their environment and, in particular, may reconnect with the solar wind around it, eroding away part of its initial magnetic flux. Here we quantitatively analyze such an interaction using combined, multipoint observations of the same MC flux rope by STEREO A, B, ACE, WIND and THEMIS on November 19–20, 2007. Observation of azimuthal magnetic flux imbalance inside a MC flux rope has been argued to stem from erosion due to magnetic reconnection at its front boundary. The present study adds to such analysis a large set of signatures expected from this erosion process. (1) Comparison of azimuthal flux imbalance for the same MC at widely separated points precludes the crossing of the MC leg as a source of bias in flux imbalance estimates. (2) The use of different methods, associated errors and parametric analyses show that only an unexpectedly large error in MC axis orientation could explain the azimuthal flux imbalance. (3) Reconnection signatures are observed at the MC front at all spacecraft, consistent with an ongoing erosion process. (4) Signatures in suprathermal electrons suggest that the trailing part of the MC has a different large-scale magnetic topology, as expected. The azimuthal magnetic flux erosion estimated at ACE and STEREO A corresponds respectively to 44% and 49% of the inferred initial azimuthal magnetic flux before MC erosion upon propagation. The corresponding average reconnection rate during transit is estimated to be in the range 0.12–0.22 mV/m, suggesting most of the erosion occurs in the inner parts of the heliosphere. Future studies ought to quantify the influence of such an erosion process on geo-effectiveness.

Computational aspects of harmonic wavelet Galerkin methods and an application to a precipitation front propagation model

This article is dedicated to harmonic wavelet Galerkin methods for the solution of partial differential equations. Several variants of the method are proposed and analyzed, using the Burgers equation as a test model. The computational complexity can be reduced when the localization properties of the wavelets and restricted interactions between different scales are exploited. The resulting variants of the method have computational complexities ranging from O(N(3)) to O(N) (N being the space dimension) per time step. A pseudo-spectral wavelet scheme is also described and compared to the methods based on connection coefficients. The harmonic wavelet Galerkin scheme is applied to a nonlinear model for the propagation of precipitation fronts, with the front locations being exposed in the sizes of the localized wavelet coefficients. (C) 2011 Elsevier Ltd. All rights reserved.

Inflammation, oxidative stress, glomerular filtration rate, and albuminuria in elderly men : a cross-sectional study

BACKGROUND: The role of inflammation and oxidative stress in mild renal impairment in the elderly is not well studied. Accordingly, we aimed at investigating the associations between estimated glomerular filtration rate (eGFR), albumin/creatinine ratio (ACR), and markers of different inflammatory pathways and oxidative stress in a community based cohort of elderly men. FINDINGS: Cystatin C-based GFR, ACR, and biomarkers of cytokine-mediated inflammation (interleukin-6, high-sensitivity C-reactive protein[CRP], serum amyloid A[SAA]), cyclooxygenase-mediated inflammation (urinary prostaglandin F2alpha [PGF2alpha]), and oxidative stress (urinary F2 isoprostanes) were assessed in the Uppsala Longitudinal Study of Adult Men(n = 647, mean age 77 years). RESULTS: In linear regression models adjusting for age, BMI, smoking, blood pressure, LDL-cholesterol, HDL-cholesterol, triglycerides, and treatment with statins, ACE-inhibitors, ASA, and anti-inflammatory agents, eGFR was inversely associated with CRP, interleukin-6, and SAA (beta-coefficient -0.13 to -0.19, p < 0.001 for all), and positively associated with urinary F2-isoprostanes (beta-coefficient 0.09, p = 0.02). In line with this, ACR was positively associated with CRP, interleukin-6, and SAA (beta- coefficient 0.09-0.12, p < 0.02 for all), and negatively associated with urinary F2-isoprostanes (beta-coefficient -0.12, p = 0.002). The associations were similar but with lower regression coefficients in a sub-sample with normal eGFR (>60 ml/min/1.73 m2, n = 514), with the exception that F2-isoprostane and SAA were no longer associated with eGFR. CONCLUSION: Our data indicate that cytokine-mediated inflammation is involved in the early stages of impaired kidney function in the elderly, but that cyclooxygenase-mediated inflammation does not play a role at this stage. The unexpected association between higher eGFR/lower albuminuria and increased F2-isoprostanes in urine merits further studies.

Job destruction and propagation of shocks

We develop and quantitatively implement a dynamic general equilibrium model with labor market matching and endogenous deterllÚnation of the job destruction rate. The mo deI produces a elose match with data on job creation and destruction. Cyelical fluctuations in the job destruction rate serve to magnify the effects of productivity shocks on output; as well as making the effects much more persistent. Interactions between the labor and capital markets, mediated by the rental rate of capital, play the central role in propagating shocks.

In search of exchange rate predictability : a study about accuracy, consistency, and granger causality of forecasts generated by a Taylor Rule Model

Este estudo investiga o poder preditivo fora da amostra, um mês à frente, de um modelo baseado na regra de Taylor para previsão de taxas de câmbio. Revisamos trabalhos relevantes que concluem que modelos macroeconômicos podem explicar a taxa de câmbio de curto prazo. Também apresentamos estudos que são céticos em relação à capacidade de variáveis macroeconômicas preverem as variações cambiais. Para contribuir com o tema, este trabalho apresenta sua própria evidência através da implementação do modelo que demonstrou o melhor resultado preditivo descrito por Molodtsova e Papell (2009), o “symmetric Taylor rule model with heterogeneous coefficients, smoothing, and a constant”. Para isso, utilizamos uma amostra de 14 moedas em relação ao dólar norte-americano que permitiu a geração de previsões mensais fora da amostra de janeiro de 2000 até março de 2014. Assim como o critério adotado por Galimberti e Moura (2012), focamos em países que adotaram o regime de câmbio flutuante e metas de inflação, porém escolhemos moedas de países desenvolvidos e em desenvolvimento. Os resultados da nossa pesquisa corroboram o estudo de Rogoff e Stavrakeva (2008), ao constatar que a conclusão da previsibilidade da taxa de câmbio depende do teste estatístico adotado, sendo necessária a adoção de testes robustos e rigorosos para adequada avaliação do modelo. Após constatar não ser possível afirmar que o modelo implementado provém previsões mais precisas do que as de um passeio aleatório, avaliamos se, pelo menos, o modelo é capaz de gerar previsões “racionais”, ou “consistentes”. Para isso, usamos o arcabouço teórico e instrumental definido e implementado por Cheung e Chinn (1998) e concluímos que as previsões oriundas do modelo de regra de Taylor são “inconsistentes”. Finalmente, realizamos testes de causalidade de Granger com o intuito de verificar se os valores defasados dos retornos previstos pelo modelo estrutural explicam os valores contemporâneos observados. Apuramos que o modelo fundamental é incapaz de antecipar os retornos realizados.