A general central limit theorem and a subsampling variance estimator for α‐mixing point processes

We establish a central limit theorem for multivariate summary statistics of nonstationary α‐mixing spatial point processes and a subsampling estimator of the covariance matrix of such statistics. The central limit theorem is crucial for establishing asymptotic properties of estimators in statistics for spatial point processes. The covariance matrix subsampling estimator is flexible and model free. It is needed, for example, to construct confidence intervals and ellipsoids based on asymptotic normality of estimators. We also provide a simulation study investigating an application of our results to estimating functions.     

Cocaine Dependence Treatment Data: Methods for Measurement Error Problems With Predictors Derived From Stationary Stochastic Processes

In a cocaine dependence treatment study, we use linear and nonlinear regression models to model posttreatment cocaine craving scores and first cocaine relapse time. A subset of the covariates are summary statistics derived from baseline daily cocaine use trajectories, such as baseline cocaine use frequency and average daily use amount. These summary statistics are subject to estimation error and can therefore cause biased estimators for the regression coefficients. Unlike classical measurement error problems, the error we encounter here is heteroscedastic with an unknown distribution, and there are no replicates for the error-prone variables or instrumental variables. We propose two robust methods to correct for the bias: a computationally efficient method-of-moments-based method for linear regression models and a subsampling extrapolation method that is generally applicable to both linear and nonlinear regression models. Simulations and an application to the cocaine dependence treatment data are used to illustrate the efficacy of the proposed methods. Asymptotic theory and variance estimation for the proposed subsampling extrapolation method and some additional simulation results are described in the online supplementary material.     

Second‐order analysis of marked inhomogeneous spatiotemporal point processes: Applications to earthquake data

To analyze interactions in marked spatiotemporal point processes (MSTPPs), we introduce marked second‐order reduced moment measures and K‐functions for inhomogeneous second‐order intensity‐reweighted stationary MSTPPs. These summary statistics, which allow us to quantify dependence between different mark‐based classifications of points, depend on the specific mark space and mark reference measure chosen. Unbiased and consistent minus‐sampling estimators are derived for all statistics considered, and a test for random labeling is indicated. In addition, we treat Voronoi intensity estimators for MSTPPs. These new statistics are finally employed to analyze an Andaman Sea earthquake data set.     

SUBSAMPLING THE JOHANSEN TEST WITH STABLE INNOVATIONS

By taking into account the thick-tail property of the errors, cointegration analysis in vector error-correction models with infinite-variance stable errors is a natural generalization of cointegration analysis in error-correction models with normally distributed errors. We study the Johansen test for cointegrated systems under symmetric stable innovations with discrete spectral measures. The results show that the distributions of the Johansen test statistics under these innovations involve nuisance parameters. To overcome the problem of nuisance parameters, we implement a nonparametric subsampling procedure. We document some subsampling simulation results and demonstrate in an empirical example how the test can be used in practice.     

Subsampling Variance Estimation for Non-stationary Spatial Lattice Data

ABSTRACT.  Most proposed subsampling and resampling methods in the literature assume stationary data. In many empirical applications, however, the hypothesis of stationarity can easily be rejected. In this paper, we demonstrate that moment and variance estimators based on the subsampling methodology can also be employed for different types of non-stationarity data. Consistency of estimators are demonstrated under mild moment and mixing conditions. Rates of convergence are provided, giving guidance for the appropriate choice of subshape size. Results from a small simulation study on finite-sample properties are also reported.     

Subsampling,symmetrization, and robust interpolation

The recently developed subsampling methodology has been shown to be valid for the construction of large-sample confidence regions for a general unknown parameter 9 under very minimal conditions. Nevertheless, in some specific cases—e.g. in the case of the sample mean of i.i.d. data—it has been noted that the subsampling distribution estimator underperforms as compared to alternative estimators such as the bootstrap or the asymptotic normal distribution (with estimated variance). In the present report we introduce a (partially) symmetrized.     

Multivariate geometric anisotropic Cox processes

This paper introduces a new modeling and inference framework for multivariate and anisotropic point processes. Building on recent innovations in multivariate spatial statistics, we propose a new family of multivariate anisotropic random fields, and from them a family of anisotropic point processes. We give conditions that make the proposed models valid. We also propose a Palm likelihood-based inference method for this type of point process, circumventing issues of likelihood tractability. Finally we illustrate the utility of the proposed modeling framework by analyzing spatial ecological observations of plants and trees in the Barro Colorado Island data.     

Subsampling quantile estimators and uniformity criteria

Several asymptotically equivalent quantile estimators recently have been proposed as alternative to the conventional sample quantile. A variety of weight functions have been obtained either by subsampling considerations or by a kernel approach, analogous to density estimation techniques. Focusing on the former approach, a unified treatment of quantile estimators derived by subsampling is developed. Closely related to the generalized Harrell-Davis (HD) and Kaigh-Lachenbruch (KL) estimators, a new statistic performed well in Monte Carlo effiency comparisons presented here. Moreover, the new estimator shares certain desirable computational and finite-sample theeoretical properties with the KL estimator to yield convenient components representations for tests of uniformity and goodness-of-fit criteria. Similar analytic treatment for the HD statistics and kernel quantile estimators, however, is precluded by intractable eigenvalue problems.     

Indices of Dependence Between Types in Multivariate Point Patterns

We propose new summary statistics quantifying several forms of dependence between points of different types in a multi-type spatial point pattern. These statistics are the multivariate counterparts of the J -function for point processes of a single type, introduced by Lieshout & Baddeley (1996). They are based on comparing distances from a type i point to either the nearest type j point or to the nearest point in the pattern regardless of type to these distances seen from an arbitrary point in space. Information about the range of interaction can also be inferred. Our statistics can be computed explicitly for a range of well-known multivariate point process models. Some applications to bivariate and trivariate data sets are presented as well.     

Tests for an end-of-sample bubble in financial time series

ABSTRACT

In this paper, we examine the issue of detecting explosive behavior in economic and financial time series when an explosive episode is both ongoing at the end of the sample and of finite length. We propose a testing strategy based on a subsampling method in which a suitable test statistic is calculated on a finite number of end-of-sample observations, with a critical value obtained using subsample test statistics calculated on the remaining observations. This approach also has the practical advantage that, by virtue of how the critical values are obtained, it can deliver tests which are robust to, among other things, conditional heteroskedasticity and serial correlation in the driving shocks. We also explore modifications of the raw statistics to account for unconditional heteroskedasticity using studentization and a White-type correction. We evaluate the finite sample size and power properties of our proposed procedures and find that they offer promising levels of power, suggesting the possibility for earlier detection of end-of-sample bubble episodes compared to existing procedures.     

A family of test statistics for trend change in mean residual life with unknown turning point using censored data

In this paper, we develop a family of test statistics for testing exponentiality against increasing then decreasing mean residual life alternatives to accommodate the randomly censored data when neither the change point nor the proportion is known. We establish the asymptotic null distribution of test statistics. Monte Carlo simulations are conducted to investigate the speed of convergence of test statistics to the asymptotic null distribution and study the performance of test statistics by the power of tests.     

A model selection criterion for clustered survival analysis with informative cluster size

We propose a model selection criterion for correlated survival data when the cluster size is informative to the outcome. This approach, called Resampling Cluster Survival Information Criterion (RCSIC), uses the Cox proportional hazards model that is weighted with the inverse of the cluster size. The RCSIC based on the within-cluster resampling idea takes into account the possible variability of the within-cluster subsampling and the possible informativeness of cluster sizes. The RCSIC allows for easy execution for the within-cluster resampling idea without a large number of resamples of the data. In contrast with the traditional model selection method in survival analysis, the RCSIC has an additional penalization for the within-cluster subsampling variability. Our simulations show the satisfactory results where the RCSIC provides a more robust power for variable selection in terms of clustered survival analysis, regardless of whether informative cluster size exists or not. Applying the RCSIC method to a periodontal disease studies, we identify the tooth loss in patients associated with the risk factors, Age, Filled Tooth, Molar, Crown, Decayed Tooth, and Smoking Status, respectively.     

A Monte Carlo method for filtering a marked doubly stochastic Poisson process

The author provides an approximated solution for the filtering of a state-space model, where the hidden state process is a continuous-time pure jump Markov process and the observations come from marked point processes. Each state k corresponds to a different marked point process, defined by its conditional intensity function λ _k (t). When a state is visited by the hidden process, the corresponding marked point process is observed. The filtering equations are obtained by applying the innovation method and the integral representation theorem of a point process martingale. Since the filtering equations belong to the family of Kushner–Stratonovich equations, an iterative solution is calculated. The theoretical solution is approximated and a Monte Carlo integration technique employed to implement it. The sequential method has been tested on a simulated data set based on marked point processes widely used in the statistical analysis of seismic sequences: the Poisson model, the stress release model and the Etas model.     

政府统计生产及其行为研究

政府统计不仅是一种手段,更重要的是一种生产活动。该研究在分析了政府生产性过程的基础上,按照行为理论,以方法论上的个人主义和个人效用最大化作为立论的基点,比较全面地揭示出政府统计用户、统计数据生产者和统计调查对象的行为特点,对于人们总结政府统计规律,规范政府统计行为,提高政府统计数据质量具有重要的参考价值。     

A functional marked point process model for lupus data

Modelling for marked point processes is an important problem, but has received remarkably little attention in the statistical literature. The authors developed a marked point process model that incorporates the use of functional data analysis in a joint estimation of the frequency function of the point process and the intensity of the mark, with application to data from 22 lupus patients consisting of times of flares in symptom severity combined with a quantitative assessment of the severity. The data indicate that a rapid decrease in drug dose is significantly associated with a decrease in flare frequency. Experiments with simulated data designed to model the actual data further support this conclusion. The Canadian Journal of Statistics 40: 517–529; 2012 © 2012 Statistical Society of Canada     

FORECASTING VOLATILITY IN THE PRESENCE OF MODEL INSTABILITY

Recent advances in financial econometrics have allowed for the construction of efficient ex post measures of daily volatility. This paper investigates the importance of instability in models of realised volatility and their corresponding forecasts. Testing for model instability is conducted with a subsampling method. We show that removing structurally unstable data of a short duration has a negligible impact on the accuracy of conditional mean forecasts of volatility. In contrast, it does provide a substantial improvement in a model's forecast density of volatility. In addition, the forecasting performance improves, often dramatically, when we evaluate models on structurally stable data.     

Subsampling confidence intervals for parameters of atmospheric time series: block size choice and calibration

Problems of practical implementation of the computer intensive subsampling methodology are addressed by Monte Carlo simulations of a situation typical for atmospheric time series. The motivating data were collected under Lake-Effect Snow Studies Project in the winter of 1983–1984 over Lake Michigan. Certain enhancements of subsampling methodology are suggested specifically on the issue of optimal block size choice.     

基于高频数据的风险值度量IVar方法综述

IVar估计是基于高频金融数据的分笔数据构造风险测量的方法。IVar的估计依赖两个步骤：首先根据分笔数据信息以及假定的交易间持续期和交易收益的联合分布，得到10g—ACD模型和扩展的UHF—GARCH模型；根据这两个模型，利用蒙特卡罗方法构造测量风险的指标IVar的估计，并详细给出这两个步骤中的计算细节。     

Analysis of tidal data via the blockwise bootstrap

We analyze tidal data from Port Mansfield, TX, using Kunsch's blockwise bootstrap in the regression setting. In particular, we estimate the variability of parameter estimates in a harmonic analysis via block subsampling of residuals from a least-squares fit. We see that naive least-squares variance estimates can be either too large or too small, depending on the strength of correlation and the design matrix. We argue that the block bootstrap is a simple, omnibus method of accounting for correlation in a regression model with correlated errors.     

Monte Carlo integration with Markov chain

There are two conceptually distinct tasks in Markov chain Monte Carlo (MCMC): a sampler is designed for simulating a Markov chain and then an estimator is constructed on the Markov chain for computing integrals and expectations. In this article, we aim to address the second task by extending the likelihood approach of Kong et al. for Monte Carlo integration. We consider a general Markov chain scheme and use partial likelihood for estimation. Basically, the Markov chain scheme is treated as a random design and a stratified estimator is defined for the baseline measure. Further, we propose useful techniques including subsampling, regulation, and amplification for achieving overall computational efficiency. Finally, we introduce approximate variance estimators for the point estimators. The method can yield substantially improved accuracy compared with Chib's estimator and the crude Monte Carlo estimator, as illustrated with three examples.