Analyzing multiple vector autoregressions through matrix-variate normal distribution with two covariance matrices

Abstract

This article proposes a new approach to analyze multiple vector autoregressive (VAR) models that render us a newly constructed matrix autoregressive (MtAR) model based on a matrix-variate normal distribution with two covariance matrices. The MtAR is a generalization of VAR models where the two covariance matrices allow the extension of MtAR to a structural MtAR analysis. The proposed MtAR can also incorporate different lag orders across VAR systems that provide more flexibility to the model. The estimation results from a simulation study and an empirical study on macroeconomic application show favorable performance of our proposed models and method.     

Bayesian analysis of threshold autoregressions

A nonasymptotic Bayesian approach is developed for analysis of data from threshold autoregressive processes with two regimes. Using the conditional likelihood function, the marginal posterior distribution for each of the parameters is derived along with posterior means and variances. A test for linear functions of the autoregressive coefficients is presented. The approach presented uses a posterior p-value averaged over the values of the threshold. The one-step ahead predictive distribution is derived along with the predictive mean and variance. In addition, equivalent results are derived conditional upon a value of the threshold. A numerical example is presented to illustrate the approach.     

Recursive online EM estimation of mixture autoregressions

This paper develops a recursive expectation–maximization (REM) algorithm for estimating a mixture autoregression (MAR) with an independent and identically distributed regime transition process. The proposed method, which is useful for long time series as well as for data available in real time, follows a recursive predictor error-type scheme. Based on a slightly modified system to the expectation–maximization (EM) equations for an MAR model, the REM algorithm consists of two steps at each iteration: the expectation step, in which the current unobserved regime transition is estimated from new data using previous recursive estimates, and the minimization step, in which the MAR parameter estimates are recursively updated following a minimization direction. Details of implementation of the REM algorithm are given and its finite-sample performance is shown via simulation experiments. In particular, the EM and REM provide roughly similar estimates, especially for moderate and long time series.     

On trends and constants in periodic autoregressions

Periodic autoregressions are characterised by autoregressive structures that vary with the season. If a time series is periodically integrated, one needs a seasonally varying differencing filter to remove the stochastic trend. When the periodic regression model contains constants and trends with unrestricted parameters, the data can show diverging seasonal deterministic trends. In this paper we derive explicit expressions for parameter restrictions that result in common deterministic trends under periodic trend stationarity and periodic integration.     

On trends and constants in periodic autoregressions

Periodic autoregressions are characterised by autoregressive structures that vary with the season. If a time series is periodically integrated, one needs a seasonally varying differencing filter to remove the stochastic trend. When the periodic regression model contains constants and trends with unrestricted parameters, the data can show diverging seasonal deterministic trends. In this paper we derive explicit expressions for parameter restrictions that result in common deterministic trends under periodic trend stationarity and periodic integration.     

Estimation bias and bias correction in reduced rank autoregressions

This paper characterizes the finite-sample bias of the maximum likelihood estimator (MLE) in a reduced rank vector autoregression and suggests two simulation-based bias corrections. One is a simple bootstrap implementation that approximates the bias at the MLE. The other is an iterative root-finding algorithm implemented using stochastic approximation methods. Both algorithms are shown to be improvements over the MLE, measured in terms of mean square error and mean absolute deviation. An illustration to US macroeconomic time series is given.     

Bootstrap confidence regions computed from autoregressions of arbitrary order

Given a linear time series, e.g. an autoregression of infinite order, we may construct a finite order approximation and use that as the basis for confidence regions. The sieve or autoregressive bootstrap, as this method is often called, is generally seen as a competitor with the better-understood block bootstrap approach. However, in the present paper we argue that, for linear time series, the sieve bootstrap has significantly better performance than blocking methods and offers a wider range of opportunities. In particular, since it does not corrupt second-order properties then it may be used in a double-bootstrap form, with the second bootstrap application being employed to calibrate a basic percentile method confidence interval. This approach confers second-order accuracy without the need to estimate variance. That offers substantial benefits, since variances of statistics based on time series can be difficult to estimate reliably, and—partly because of the relatively small amount of information contained in a dependent process—are notorious for causing problems when used to Studentize. Other advantages of the sieve bootstrap include considerably greater robustness against variations in the choice of the tuning parameter, here equal to the autoregressive order, and the fact that, in contradistinction to the case of the block bootstrap, the percentile t version of the sieve bootstrap may be based on the 'raw' estimator of standard error. In the process of establishing these properties we show that the sieve bootstrap is second order correct.     

Third and fourth moments of vector autoregressions with regime switching

We derive matrix formulae in closed form for the unconditional third and fourth moments of a broad class of vector autoregressive time series with regime switching. First and second moments are well known. New measures of multivariate skewness and kurtosis are introduced and basic properties are investigated. The knowledge of series level, variation, co-movements, skewness, and kurtosis is useful to support model interpretation in real data application. Numerical examples complete the paper.     

Estimating multivariate autoregressive moving average models by fitting long autoregressions

Some simple methods for the estimation of mixed multivariate autoregressive moving average time series models are introduced. The methods require the fitting of a long autoregression to the data and the computation of consistent initial estimates for the parameters of the model. After these preliminaries the estimators of the paper are obtained by applying weighted least squares to a multivariate auxiliary regression model. Two types of weight matrices are considered. Both of them yield estimators which are strongly consistent and asymptotically normally distributed. The first estimators are also asymptotically efficient while the second ones are not fully efficient but computationally simple. A simulation study is performed to illustrate the behaviour of the estimators in finite samples.     

Asymptotic robustness of least median of squares for autoregressions with additive outliers

This paper considers the robustness properties in the time series context of the least median of squares (LMS) estimator. The influence function of the LMS estimator is derived under additive outlier contamination. This influence function is redescending and bounded for fixed values of the AR parameters. The gross-error sensitivity, however, is an unbounded function of the AR parameters. In order to asses the global robustness behavior of the LMS estimator, we consider several notions of breakdown. The breakdown points of the LMS estimator depend on the value of the underlying AR parameter. Generally, the breakdown point is below one half for high values of the AR parameter. The bias curves of the LMS estimator reveal, however, that the magnitude of outliers has to be considerable in order to cause breakdown.     

Intrinsic Priors for Testing Two Normal Means with Intrinsic Bayes Factors

ABSTRACT

In this article we consider the problem of comparing two normal means with unknown common variance using a Bayesian approach. Conventional Bayes factors with improper non informative priors are not well defined. The intrinsic Bayes factors are used to overcome such a difficulty. We derive intrinsic priors whose Bayes factors are asymptotically equivalent to the corresponding intrinsic Bayes factors. We illustrate our results with numerical examples.     

Continuous-time threshold autoregressions with jumps: Properties,estimation, and application to electricity markets

Continuous-time autoregressive processes have been applied successfully in many fields and are particularly advantageous in the modeling of irregularly spaced or high-frequency time series data. A convenient nonlinear extension of this model are continuous-time threshold autoregressions (CTAR). CTAR allow for greater flexibility in model parameters and can represent a regime switching behavior. However, so far only Gaussian CTAR processes have been defined, so that this model class could not be used for data with jumps, as frequently observed in financial applications. Hence, as a novelty, we construct CTAR processes with jumps in this paper. Existence of a unique weak solution and weak consistency of an Euler approximation scheme is proven. As a closed form expression of the likelihood is not available, we use kernel-based particle filtering for estimation. We fit our model to the Physical Electricity Index and show that it describes the data better than other comparable approaches.     

Modelling daily multivariate pollutant data at multiple sites

Summary. This paper considers the spatiotemporal modelling of four pollutants measured daily at eight monitoring sites in London over a 4-year period. Such multiple-pollutant data sets measured over time at multiple sites within a region of interest are typical. Here, the modelling was carried out to provide the exposure for a study investigating the health effects of air pollution. Alternative objectives include the design problem of the positioning of a new monitoring site, or for regulatory purposes to determine whether environmental standards are being met. In general, analyses are hampered by missing data due, for example, to a particular pollutant not being measured at a site, a monitor being inactive by design (e.g. a 6-day monitoring schedule) or because of an unreliable or faulty monitor. Data of this type are modelled here within a dynamic linear modelling framework, in which the dependences across time, space and pollutants are exploited. Throughout the approach is Bayesian, with implementation via Markov chain Monte Carlo sampling.     

如何开展对非正规部门的统计——介绍国际劳工组织的有关决议

在很多国家 ,特别是发展中国家 ,非正规部门在经济和社会发展的各个领域起着非常重要的作用 ,特别是在创造就业岗位、形成国民收入方面的作用不容忽视。但由于非正规部门的特殊性 ,如生产规模小 ,没有在工商等部门进行注册或登记等 ,而难以界定和统计其活动 ,不少国家的国民经济核算及就业统计对此都存在比较严重的遗漏 ,影响了国民经济核算及就业统计的数据质量。为了真实地反映各国的就业状况 ,确保国民经济核算的全面性 ,国际劳工组织 (ＩＬＯ)在 1 993年第 1 5届国际劳工统计大会做出了《关于非正规部门就业统计的决议》 ,这个决议涉及…     

Joint modeling of multivariate hearing thresholds measured longitudinally at multiple frequencies

Pure-tone thresholds are used to estimate hearing acuity and, when measured longitudinally, can characterize age-related changes in hearing. Measured at multiple-frequencies, multiple-irregular time points, for right and left ears, these longitudinal studies of age-related hearing loss produce data of inherent complexity due to: 1) multivariate outcomes at different frequencies; 2) longitudinal measurements taken at subject-specific time intervals; and 3) inter-ear correlations due to clustering and nesting. To address limitations in existing methods, we propose a multivariate generalized linear mixed model (mGLMM) and assess its performance. We demonstrate its application using a unique dataset from a cohort study of age-related hearing loss.     

Demographic policy in the light of the resolutions of the Twenty-sixth Congress of the Communist Party

The USSR is moving from extensive to intensive type of population reproduction, not because of any governmental measures but as a reaction to objective circumstances. At present there is a reduction in the population cohort of working age. This has an adverse effect on production. And the number of the elderly is increasing; we must involve them more in the production process, in monitoring the service sphere, for instance. The infant mortality rate has declined sharply in the USSR, especially since World War II. The birthrate has also dropped, but it is still higher than in the US, France, or West Germany. Some authorities think that this (the ZPG) is a good idea, but this author thinks it is a bad idea, meaning a drop in society's productive forces. We need an active demographic policy, meaning one which ensures optimal reproduction of the population with an improvement in its health, culture, and living standard. The policy must be differentiated, as conditions vary from one part of the country to another. An example is the 1981 decree on assistance to families with children which provides for payments of 35-50 rubles/month to mothers who must stay home to look after children. The number of children per family to ensure reproduction must be 2.6 on the average. In the Central Asian republics family sizes of 4-5 children are anticipated, but in the RSFSR, Ukraine, Belorussia, and the Baltic republics the anticipated size of families is rather low. Measures must be taken to alter this by providing assistance to families with children.     

流动人口社会融合指标体系内在关系研究

利用2014年流动人口社会融合专题调查数据,分别以经济融合、社区融合、文化接纳、自我认同中任何一个变量为因变量,剩余三个变量为自变量,构建由经济模型、社会模型、文化模型和自我认同模型组成的联立方程组,并采用多元线性回归和logistic回归为主要方法,引入流动性别、年龄、受教育程度、户口性质、流动区域、民族等六个变量为控制变量,对流动人口社会融合的四个维度的内在关系进行了分析。研究发现,经济收入与社区融合统计关系不显著,与文化接纳和身份认同统计关系显著,表明收入越高文化接纳越好,身份认同却越差;经济收入对心理文化及社区参与的影响要略微强于心理文化及社区参与对经济收入的影响,经济融合有一定的独立性特征,而社区参与、文化接纳与身份认同存在更高的一致性关系。     

Using multiple imputation to estimate cumulative distribution functions in longitudinal data analysis with data missing at random

In longitudinal clinical studies, after randomization at baseline, subjects are followed for a period of time for development of symptoms. The interested inference could be the mean change from baseline to a particular visit in some lab values, the proportion of responders to some threshold category at a particular visit post baseline, or the time to some important event. However, in some applications, the interest may be in estimating the cumulative distribution function (CDF) at a fixed time point post baseline. When the data are fully observed, the CDF can be estimated by the empirical CDF. When patients discontinue prematurely during the course of the study, the empirical CDF cannot be directly used. In this paper, we use multiple imputation as a way to estimate the CDF in longitudinal studies when data are missing at random. The validity of the method is assessed on the basis of the bias and the Kolmogorov–Smirnov distance. The results suggest that multiple imputation yields less bias and less variability than the often used last observation carried forward method. Copyright © 2013 John Wiley & Sons, Ltd.     

Intrinsic Bayes factor approach to a test for the power law process

The versatile new criterion called the intrinsic Bayes factor (IBF), introduced by Berger and Pericchi [J. Amer. Statist. Assoc. 91 (1996) 109–122], has made it possible to perform model selection and hypotheses testing using standard (improper) noninformative priors in a variety of situations. In this paper, we use their methodology to test several hypotheses regarding the shape parameter of the power law process, which has been widely used to model failure times of repairable systems. Assuming that we have data from the process according to the time-truncation sampling scheme, we derive the arithmetic IBFs using four default priors, including the reference and Jeffreys priors. We establish the frequentist probability matching properties of these priors. We also identify two priors that are justifiable under both time-truncation and failure-truncation schemes, so that the IBFs for both schemes can be unified. Deducing the intrinsic priors of a certain canonical form, as the time of truncation tends to infinity, we show that the arithmetic IBFs correspond asymptotically to actual Bayes factors. We also discuss the expected IBFs, which are useful with small samples. We then use these results to analyze an actual data set on the interruption times of a transmission line, summarizing our results under the default priors.