Aspects of second‐order analysis of structured inhomogeneous spatio‐temporal point processes

Statistical methodology for spatio‐temporal point processes is in its infancy. We consider second‐order analysis based on pair correlation functions and K‐functions for general inhomogeneous spatio‐temporal point processes and for inhomogeneous spatio‐temporal Cox processes. Assuming spatio‐temporal separability of the intensity function, we clarify different meanings of second‐order spatio‐temporal separability. One is second‐order spatio‐temporal independence and relates to log‐Gaussian Cox processes with an additive covariance structure of the underlying spatio‐temporal Gaussian process. Another concerns shot‐noise Cox processes with a separable spatio‐temporal covariance density. We propose diagnostic procedures for checking hypotheses of second‐order spatio‐temporal separability, which we apply on simulated and real data.     

Socio‐economic distance and spatial patterns in unemployment

This paper examines the spatial patterns of unemployment in Chicago between 1980 and 1990. We study unemployment clustering with respect to different social and economic distance metrics that reflect the structure of agents' social networks. Specifically, we use physical distance, travel time, and differences in ethnic and occupational distribution between locations. Our goal is to determine whether our estimates of spatial dependence are consistent with models in which agents' employment status is affected by information exchanged locally within their social networks. We present non‐parametric estimates of correlation across Census tracts as a function of each distance metric as well as pairs of metrics, both for unemployment rate itself and after conditioning on a set of tract characteristics. Our results indicate that there is a strong positive and statistically significant degree of spatial dependence in the distribution of raw unemployment rates, for all our metrics. However, once we condition on a set of covariates, most of the spatial autocorrelation is eliminated, with the exception of physical and occupational distance. Racial and ethnic composition variables are the single most important factor in explaining the observed correlation patterns. Copyright © 2002 John Wiley & Sons, Ltd.     

Portfolio selection based on semivariance and distance correlation under minimum variance framework

In the minimum variance model, the covariance matrix plays an important role because it measures the risk and relationship of asset returns simultaneously under the normality assumption. However, in practice, the distribution of asset returns is nonnormal and has an obvious fat‐tail nature. In addition, the risk is one‐sided. In this paper, the main objective is to propose a better tool to replace the covariance matrix. The covariance matrix can be decomposed into two parts: a diagonal variance matrix and a square matrix with its elements being the Pearson correlation coefficient. A substitution of the covariance matrix is presented by replacing the variance and Pearson correlation coefficient in the decomposition of the covariance matrix with a semivariance and distance correlation coefficient, respectively. The proposed portfolio optimization strategy is applied to empirical data, and the numerical studies show the strategy performs well.     

Vector Stochastic Processes with Pólya‐Type Correlation Structure

This paper introduces a simple method to construct a stationary process on the real line with a Pólya‐type covariance function and with any infinitely divisible marginal distribution, by randomising the timescale of the increment of a second‐order Lévy process with an appropriate positive random variable. With the construction method extended to the multivariate case, we construct vector stochastic processes with Pólya‐type direct covariance functions and with any specified infinitely divisible marginal distributions. This makes available a new class of non‐Gaussian vector stochastic processes with flexible correlation structure for use in modelling and simulation.     

Moment-based estimation of extendible Marshall-Olkin copulas

Associated with any parametric family of Lévy subordinators there is a parametric family of extendible Marshall-Olkin copulas, which shares the dependence structure with the vector of first passage times of the Lévy subordinator across i.i.d. exponential threshold levels. The present article derives a strongly consistent and asymptotically normal estimator for the parameters in such models. The estimation strategy is to minimize the Euclidean distance between certain empirical and theoretical functionals of the distribution. As a byproduct, the covariance structure of the order statistics of a d-dimensional extendible Marshall-Olkin distribution is computed.     

Testing block‐diagonal covariance structure for high‐dimensional data

A test statistic is developed for making inference about a block‐diagonal structure of the covariance matrix when the dimensionality p exceeds n, where n = N ? 1 and N denotes the sample size. The suggested procedure extends the complete independence results. Because the classical hypothesis testing methods based on the likelihood ratio degenerate when p > n, the main idea is to turn instead to a distance function between the null and alternative hypotheses. The test statistic is then constructed using a consistent estimator of this function, where consistency is considered in an asymptotic framework that allows p to grow together with n. The suggested statistic is also shown to have an asymptotic normality under the null hypothesis. Some auxiliary results on the moments of products of multivariate normal random vectors and higher‐order moments of the Wishart matrices, which are important for our evaluation of the test statistic, are derived. We perform empirical power analysis for a number of alternative covariance structures.     

Covariance measurement in the presence of non-synchronous trading and market microstructure noise

This paper studies the problem of covariance estimation when prices are observed non-synchronously and contaminated by i.i.d. microstructure noise. We derive closed form expressions for the bias and variance of three popular covariance estimators, namely realised covariance, realised covariance plus lead and lag adjustments, and the Hayashi and Yoshida estimator, and present a comprehensive investigation into their properties and relative efficiency. Our main finding is that the ordering of the covariance estimators in terms of efficiency crucially depends on the level of microstructure noise, as well as the level of correlation. In fact, for sufficiently high levels of noise, the standard realised covariance estimator (without any corrections for non-synchronous trading) can be most efficient. We also propose a sparse sampling implementation of the Hayashi and Yoshida estimator, study the robustness of our findings using simulations with stochastic volatility and correlation, and highlight some important practical considerations.     

Graphical Diagnostics for Modeling Unstructured Covariance Matrices

Prudent statistical analysis of correlated data requires accounting for the correlation among the measurements. Specifying a form for the covariance matrix of the data could reduce the high number of parameters of the covariance and increase efficiency of the inferences about the regression parameters. Motivated by the success of ordinary, partial and inverse correlograms in identifying parsimonious models for stationary time series, we introduce generalizations of these plots for nonstationary data. Their roles in detecting heterogeneity and correlation of the data and identifying parsimonious models for the covariance matrix are illuminated using a longitudinal dataset. Decomposition of a covariance matrix into "variance" and "dependence" components provides the necessary ingredients for the proposed graphs. This amounts to replacing a 3-D correlation plot by a pair of 2-D plots, providing complementary information about dependence and heterogeneity. Models identified and fitted using the variance-correlation decomposition of a covariance matrix are not guaranteed to be positive definite, but those using the modified Cholesky decomposition are. Limitations of our graphical diagnostics for general multivariate data where the measurements are not (time-) ordered are discussed.     

On Some Characteristics of Gaussian Covariance Functions

The concepts of isotropy/anisotropy and separability/non‐separability of a covariance function are strictly related. If a covariance function is separable, it cannot be isotropic or geometrically anisotropic, except for the Gaussian covariance function, which is the only model both separable and isotropic. In this paper, some interesting results concerning the Gaussian covariance model and its properties related to isotropy and separability are given, and moreover, some examples are provided. Finally, a discussion on asymmetric models, with Gaussian marginals, is furnished and the strictly positive definiteness condition is discussed.     

Asymptotic theory for nonparametric regression with spatial data

Nonparametric regression with spatial, or spatio-temporal, data is considered. The conditional mean of a dependent variable, given explanatory ones, is a nonparametric function, while the conditional covariance reflects spatial correlation. Conditional heteroscedasticity is also allowed, as well as non-identically distributed observations. Instead of mixing conditions, a (possibly non-stationary) linear process is assumed for disturbances, allowing for long range, as well as short-range, dependence, while decay in dependence in explanatory variables is described using a measure based on the departure of the joint density from the product of marginal densities. A basic triangular array setting is employed, with the aim of covering various patterns of spatial observation. Sufficient conditions are established for consistency and asymptotic normality of kernel regression estimates. When the cross-sectional dependence is sufficiently mild, the asymptotic variance in the central limit theorem is the same as when observations are independent; otherwise, the rate of convergence is slower. We discuss the application of our conditions to spatial autoregressive models, and models defined on a regular lattice.     

Modelling and forecasting multivariate realized volatility

This paper proposes a methodology for dynamic modelling and forecasting of realized covariance matrices based on fractionally integrated processes. The approach allows for flexible dependence patterns and automatically guarantees positive definiteness of the forecast. We provide an empirical application of the model, which shows that it outperforms other approaches in the extant literature, both in terms of statistical precision as well as in terms of providing a superior mean‐variance trade‐off in a classical investment decision setting. Copyright © 2010 John Wiley & Sons, Ltd.     

Uncertainty Estimation for Pseudo‐Bayesian Inference Under Complex Sampling

Social and economic studies are often implemented as complex survey designs. For example, multistage, unequal probability sampling designs utilised by federal statistical agencies are typically constructed to maximise the efficiency of the target domain level estimator (e.g. indexed by geographic area) within cost constraints for survey administration. Such designs may induce dependence between the sampled units; for example, with employment of a sampling step that selects geographically indexed clusters of units. A sampling‐weighted pseudo‐posterior distribution may be used to estimate the population model on the observed sample. The dependence induced between coclustered units inflates the scale of the resulting pseudo‐posterior covariance matrix that has been shown to induce under coverage of the credibility sets. By bridging results across Bayesian model misspecification and survey sampling, we demonstrate that the scale and shape of the asymptotic distributions are different between each of the pseudo‐maximum likelihood estimate (MLE), the pseudo‐posterior and the MLE under simple random sampling. Through insights from survey‐sampling variance estimation and recent advances in computational methods, we devise a correction applied as a simple and fast postprocessing step to Markov chain Monte Carlo draws of the pseudo‐posterior distribution. This adjustment projects the pseudo‐posterior covariance matrix such that the nominal coverage is approximately achieved. We make an application to the National Survey on Drug Use and Health as a motivating example and we demonstrate the efficacy of our scale and shape projection procedure on synthetic data on several common archetypes of survey designs.     

Trivariate non‐Gaussian copulas to analyze the spatial behavior of rat sightings

The brown rat lives with man in a wide variety of environmental contexts and adversely affects public health by transmission of diseases, bites, and allergies. Understanding behavioral and spatial correlation aspects of pest species can contribute to their effective management and control. Rat sightings can be described by spatial coordinates in a particular region of interest defining a spatial point pattern. In this paper, we investigate the spatial structure of rat sightings in the Latina district of Madrid (Spain) and its relation to a number of distance‐based covariates that relate to the proliferation of rats. Given a number of locations, biologically considered as attractor points, the spatial dependence is modeled by distance‐based covariates and angular orientations through copula functions. We build a particular spatial trivariate distribution using univariate margins coming from the covariate information and provide predictive distributions for such distances and angular orientations.     

Missing in Asynchronicity: A Kalman‐em Approach for Multivariate Realized Covariance Estimation

Motivated by the need for a positive‐semidefinite estimator of multivariate realized covariance matrices, we model noisy and asynchronous ultra‐high‐frequency asset prices in a state‐space framework with missing data. We then estimate the covariance matrix of the latent states through a Kalman smoother and expectation maximization (KEM) algorithm. Iterating between the two EM steps, we obtain a covariance matrix estimate which is robust to both asynchronicity and microstructure noise, and positive‐semidefinite by construction. We show the performance of the KEM estimator using extensive Monte Carlo simulations that mimic the liquidity and market microstructure characteristics of the S&P 500 universe as well as in a high‐dimensional application on US stocks. KEM provides very accurate covariance matrix estimates and significantly outperforms alternative approaches recently introduced in the literature. Copyright © 2014 John Wiley & Sons, Ltd.     

A measure of mutual complete dependence

Two random variables X and Y on a common probability space are mutually completely dependent (m.c.d.) if each one is a function of the other with probability one. For continuous X and Y, a natural approach to constructing a measure of dependence is via the distance between the copula of X and Y and the independence copula. We show that this approach depends crucially on the choice of the distance function. For example, the L ^p -distances, suggested by Schweizer and Wolff, cannot generate a measure of (mutual complete) dependence, since every copula is the uniform limit of copulas linking m.c.d. variables. Instead, we propose to use a modified Sobolev norm, with respect to which mutual complete dependence cannot approximate any other kind of dependence. This Sobolev norm yields the first nonparametric measure of dependence which, among other things, captures precisely the two extremes of dependence, i.e., it equals 0 if and only if X and Y are independent, and 1 if and only if X and Y are m.c.d. Examples are given to illustrate the difference to the Schweizer–Wolff measure.     

Efficient estimation of factor models with time and cross‐sectional dependence

This paper studies the efficient estimation of large‐dimensional factor models with both time and cross‐sectional dependence assuming (N,T) separability of the covariance matrix. The asymptotic distribution of the estimator of the factor and factor‐loading space under factor stationarity is derived and compared to that of the principal component (PC) estimator. The paper also considers the case when factors exhibit a unit root. We provide feasible estimators and show in a simulation study that they are more efficient than the PC estimator in finite samples. In application, the estimation procedure is employed to estimate the Lee–Carter model and life expectancy is forecast. The Dutch gender gap is explored and the relationship between life expectancy and the level of economic development is examined in a cross‐country comparison.     

Multivariate residual‐based finite‐sample tests for serial dependence and ARCH effects with applications to asset pricing models

In this paper, we propose several finite‐sample specification tests for multivariate linear regressions (MLR). We focus on tests for serial dependence and ARCH effects with possibly non‐Gaussian errors. The tests are based on properly standardized multivariate residuals to ensure invariance to error covariances. The procedures proposed provide: (i) exact variants of standard multivariate portmanteau tests for serial correlation as well as ARCH effects, and (ii) exact versions of the diagnostics presented by Shanken ( 1990 ) which are based on combining univariate specification tests. Specifically, we combine tests across equations using a Monte Carlo (MC) test method so that Bonferroni‐type bounds can be avoided. The procedures considered are evaluated in a simulation experiment: the latter shows that standard asymptotic procedures suffer from serious size problems, while the MC tests suggested display excellent size and power properties, even when the sample size is small relative to the number of equations, with normal or Student‐t errors. The tests proposed are applied to the Fama–French three‐factor model. Our findings suggest that the i.i.d. error assumption provides an acceptable working framework once we allow for non‐Gaussian errors within 5‐year sub‐periods, whereas temporal instabilities clearly plague the full‐sample dataset. Copyright © 2009 John Wiley & Sons, Ltd.     

Dynamic stochastic copula models: estimation,inference and applications

We propose a new dynamic copula model in which the parameter characterizing dependence follows an autoregressive process. As this model class includes the Gaussian copula with stochastic correlation process, it can be viewed as a generalization of multivariate stochastic volatility models. Despite the complexity of the model, the decoupling of marginals and dependence parameters facilitates estimation. We propose estimation in two steps, where first the parameters of the marginal distributions are estimated, and then those of the copula. Parameters of the latent processes (volatilities and dependence) are estimated using efficient importance sampling. We discuss goodness‐of‐fit tests and ways to forecast the dependence parameter. For two bivariate stock index series, we show that the proposed model outperforms standard competing models. Copyright © 2010 John Wiley & Sons, Ltd.     

Heteroskedasticity, autocorrelation, and spatial correlation robust inference in linear panel models with fixed-effects

This paper develops an asymptotic theory for test statistics in linear panel models that are robust to heteroskedasticity, autocorrelation and/or spatial correlation. Two classes of standard errors are analyzed. Both are based on nonparametric heteroskedasticity autocorrelation (HAC) covariance matrix estimators. The first class is based on averages of HAC estimators across individuals in the cross-section, i.e. “averages of HACs”. This class includes the well known cluster standard errors analyzed by Arellano (1987) as a special case. The second class is based on the HAC of cross-section averages and was proposed by Driscoll and Kraay (1998). The ”HAC of averages” standard errors are robust to heteroskedasticity, serial correlation and spatial correlation but weak dependence in the time dimension is required. The “averages of HACs” standard errors are robust to heteroskedasticity and serial correlation including the nonstationary case but they are not valid in the presence of spatial correlation. The main contribution of the paper is to develop a fixed-b asymptotic theory for statistics based on both classes of standard errors in models with individual and possibly time fixed-effects dummy variables. The asymptotics is carried out for large time sample sizes for both fixed and large cross-section sample sizes. Extensive simulations show that the fixed-b approximation is usually much better than the traditional normal or chi-square approximation especially for the Driscoll-Kraay standard errors. The use of fixed-b critical values will lead to more reliable inference in practice especially for tests of joint hypotheses.     

Telling tales from the tails: High‐dimensional tail interdependence

We propose a simple and flexible framework that allows for a comprehensive analysis of tail interdependence in high dimensions. We use co‐exceedances to capture the structure of the dependence in the tails and, relying on the concept of multi‐information, define the coefficient of tail interdependence. Within this framework, we develop statistical tests of (i) independence in the tails, (ii) goodness‐of‐fit of the tail interdependence structure of a hypothesized model with the one observed in the data, and (iii) dependence symmetry between any two tails. We present an analysis of tail interdependence among 250 constituents of the S&P 250 index.