Predicting integrals of diffusion processes

Consider predicting the integral of a diffusion process Z in a bounded interval A, based on the observations Z(t_1n),…,Z(t_nn), where t_1n,…,t_nn is a dense triangular array of points (the step of discretization tends to zero as n increases) in the bounded interval. The best linear predictor is generally not asymptotically optimal. Instead, we predict using the conditional expectation of the integral of the diffusion process, the optimal predictor in terms of minimizing the mean squared error, given the observed values of the process. We obtain that, conditioning on the observed values, the order of convergence in probability to zero of the mean squared prediction error is O_p(n⁻²). We prove that the standardized conditional prediction error is approximately Gaussian with mean zero and unit variance, even though the underlying diffusion is generally non-Gaussian. Because the optimal predictor is hard to calculate exactly for most diffusions, we present an easily computed approximation that is asymptotically optimal. This approximation is a function of the diffusion coefficient.     

Testing equality of proportions with incomplete correlated data

Let (ψ_i,φ_i) be independent, identically distributed pairs of zero-one random variables with (possible) dependence of ψ_i and φ_i within the pair. For n pairs, both variables are observed, but for m₁ additional pairs only ψ_i is observed and for m₂ others φ_i is observed. If π_1· = P{ψ_i = 1} and π_·1=P{φ_i, the problem is to test π_1·=π_·1. Maximum likelihood estimates of π_1· and π_·1 are obtained via the EM algorithm. A test statistic is developed whose null distribution is asymptotically chi-square with one degree of freedom (as n and either m₁ or m₂ tend to infinity). If m₁ = m₂ = 0 the statistic reduces to that of McNemar's test; if n = 0, it is equivalent to the statistic for testing equality of two independent proportions. This test is compared with other tests by means of Pitman efficiency. Examples are presented.     

The flower intersection problem for Kirkman triple systems

The flower at a point x in a Steiner triple system is the set of all triples containing x. Denote by I_R^*[r] the set of all integers k such that there exists a pair of KTS(2r+1) having k+r triples in common, r of them being the triples of a common flower. In this article we determine the set I_R^*[r] for any positive integer r≡1 (mod 3) (only nine cases are left undecided for r=7,13,16,19), and establish that I_R^*[r]=J[r] for r≡1 (mod 3) and r22 where J[r]={0,1,…,2r(r−1)/3−6,2r(r−1)/3−4,2r(r−1)/3}.     

Asymptotic normality of a weighted integrated squared error of kernel regression estimates with data-dependent bandwidth

Let (X, Y) be a bivariate random vector and let be the regression function of Y on X that has to be estimated from a sample of i.i.d. random vectors (X₁, Y₁),…,(X_n, Y_n) having the same distribution as (X, Y). In the present paper it is shown that the normalized integrated squared error of a kernel estimator with data-driven bandwidth is asymptotically normally distributed.     

Bayesian Inference in Cointegrated I (2) Systems: A Generalization of the Triangular Model

This paper generalizes the cointegrating model of Phillips (1991) to allow for I (0), I (1) and I (2) processes. The model has a simple form that permits a wider range of I (2) processes than are usually considered, including a more flexible form of polynomial cointegration. Further, the specification relaxes restrictions identified by Phillips (1991) on the I (1) and I (2) cointegrating vectors and restrictions on how the stochastic trends enter the system. To date there has been little work on Bayesian I (2) analysis and so this paper attempts to address this gap in the literature. A method of Bayesian inference in potentially I (2) processes is presented with application to Australian money demand using a Jeffreys prior and a shrinkage prior.     

A critical look at Lo's modified R/S statistic

We report on an empirical investigation of the modified rescaled adjusted range or R/S statistic that was proposed by Lo, 1991. Econometrica 59, 1279–1313, as a test for long-range dependence with good robustness properties under ‘extra’ short-range dependence. In contrast to the classical R/S statistic that uses the standard deviation S to normalize the rescaled range R, Lo's modified R/S-statistic V_q is normalized by a modified standard deviation S_q which takes into account the covariances of the first q lags, so as to discount the influence of the short-range dependence structure that might be present in the data. Depending on the value of the resulting test-statistic V_q, the null hypothesis of no long-range dependence is either rejected or accepted. By performing Monte-Carlo simulations with ‘truly’ long-range- and short-range dependent time series, we study the behavior of V_q, as a function of q, and uncover a number of serious drawbacks to using Lo's method in practice. For example, we show that as the truncation lag q increases, the test statistic V_q has a strong bias toward accepting the null hypothesis (i.e., no long-range dependence), even in ideal situations of ‘purely’ long-range dependent data.     

Linearity of regression for non-adjacent record values

Let X₁,X₂,… be a sequence of iid random variables having a continuous distribution; by R₁,R₂,… denote the corresponding record values. All the distributions allowing linearity of regressions either E(R_m+k|R_m) or E(R_m|R_m+k) are identified.     

图解Bootstrap方法用于空间相关性Moran检验的有效性

因为区域间经济收敛、外商直接投资和知识溢出等领域的空间经济计量研究依赖于空间关系的存在,所以进行空间相关性Moran’s I检验是关键。然而,已有空间相关性Moran’s I检验理论受到众多假设条件限制。利用＂名义水平—实际水平＂图和＂名义水平—功效＂图,解析非对称Wild Bootstrap方法用于空间相关性Moran’s I检验的有限样本性质,发现即使模型不满足经典的分布假设条件,与渐近检验相比,Bootstrap方法也能够有效地检验研究对象间的空间相关性。     

A-optimal diallel crosses for test versus control comparisons

A-optimality of block designs for control versus test comparisons in diallel crosses is investigated. A sufficient condition for designs to be A-optimal is derived. Type S₀ designs are defined and A-optimal type S₀ designs are characterized. A lower bound to the A-efficiency of type S₀ designs is also given. Using the lower bound to A-efficiency, type S₀ designs are shown to yield efficient designs for test versus control comparisons.     

Adaptive wavelet empirical Bayes estimation of a location or a scale parameter

Assume that in independent two-dimensional random vectors (X₁,θ₁),…,(X_n,θ_n), each θ_i is distributed according to some unknown prior density function g. Also, given θ_i=θ, X_i has the conditional density function q(x−θ), x,θ(−∞,∞) (a location parameter case), or θ⁻¹q(x/θ), x,θ(0,∞) (a scale parameter case). In each pair the first component is observable, but the second is not. After the (n+1)th pair (X_n+1,θ_n+1) is obtained, the objective is to construct an empirical Bayes (EB) estimator of θ. In this paper we derive the EB estimators of θ based on a wavelet approximation with Meyer-type wavelets. We show that these estimators provide adaptation not only in the case when g belongs to the Sobolev space H with an unknown , but also when g is supersmooth.     

Likelihood ratio ordering of order statistics

This paper provides an improvement on the work of Bapat and Kochar (1994, Linear Algebra Appl., 199, 281–291) and strengthens the literature on the likelihood ratio ordering of order statistics. For independent (but possibly nonidentically distributed) absolutely continuous random variables X₁,…,X_n, it is shown under some weak conditions that

X_1:n^lrlrX_n:n,

where ^lr stands for the likelihood ratio ordering and X_k:n represents the kth-order statistic.     

Some new partial orderings on Sn and Sn×Sn

A new class of partial orderings on S_n, the set of permutations of {1,…,n}, is given. Each of these partial orderings is shown to be a subordering of a recently described partial ordering on S_n (Block, Chhetry, Fang and Sampson (1990)) which is related to Schriever's (1987) more associated ordering on bivariate distributions. Also given is an extension of three known partial orders on S_n to partial orders on S_n×S_n. These extensions facilitate the study of functions from S_n×S_n into , which preser these partial orderings, thereby, providing a methodology for extending the notion of arrangement increasing functions.     

Nonparametric estimation of the lower tail dependence λL in bivariate copulas

The lower tail dependence λ_L is a measure that characterizes the tendency of extreme co-movements in the lower tails of a bivariate distribution. It is invariant with respect to strictly increasing transformations of the marginal distribution and is therefore a function of the copula of the bivariate distribution. λ_L plays an important role in modelling aggregate financial risk with copulas. This paper introduces three non-parametric estimators for λ_L. They are weakly consistent under mild regularity conditions on the copula and under the assumption that the number k = k(n) of observations in the lower tail, used for estimation, is asymptotically k ≈ √n. The finite sample properties of the estimators are investigated using a Monte Carlo simulation in special cases. It turns out that these estimators are biased, where amount and sign of the bias depend on the underlying copula, on the sample size n, on k, and on the true value of λ_L.     

The effect of fat-tailed error terms on the properties of systemwise RESET test

The small sample properties of the systemwise RESET (Regression Specification Error Test) test for functional misspecification are investigated using normal and non-normal error terms. When using normally distributed or less heavy tailed error terms, we find the Rao's multivariate F-test to be best among all other alternative test methods (i.e. Wald, Lagrange Multiplier and Likelihood Ratio). Using the bootstrap critical values, however, all test methods perform satisfactorily in almost all situations. However, the test methods perform extremely badly (even the RAO test) when the error terms are very heavy tailed.     

On the non-existence of generalised Hadamard matrices

The determinant of a generalized Hadamard matrix over its group ring factored out by the relation Σ_gεG G = 0 is shown to have certain number theoretic properties. These are exploited to prove the non-existence of many generalised Hadamard matrices for groups whose orders are divisible by 3, 5 or 7. For example the GH(15, C₁₅), GH(15, C₃) and GH(15, C₅) do not exist. Also for certain n and G we find the set of determinants of the GH(n, G) matrices.     

空间回归模型选择的反思

空间计量经济学存在两种最基本的模型:空间滞后模型和空间误差模型,这里旨在重新思考和探讨这两种空间回归模型的选择,结论为:Moran’s I指数可以用来判断回归模型后的残差是否存在空间依赖性;在实证分析中,采用拉格朗日乘子检验判断两种模型优劣是最常见的做法。然而,该检验仅仅是基于统计推断而忽略了理论基础,因此,可能导致选择错误的模型;在实证分析中,空间误差模型经常被选择性遗忘,而该模型的适用性较空间滞后模型更为广泛;实证分析大多缺乏空间回归模型设定的探讨,Anselin提出三个统计量,并且,如果模型设定正确,应该遵从Wald统计量>Log likelihood统计量>LM统计量的排列顺序。     

Dynamic Graphics for Exploring Spatial Data with Application to Locating Global and Local Anomalies

We explore the application of dynamic graphics to the exploratory analysis of spatial data. We introduce a number of new tools and illustrate their use with prototype software, developed at Trinity College, Dublin. These tools are used to examine local variability—anomalies—through plots of the data that display its marginal and multivariate distributions, through interactive smoothers, and through plots motivated by the spatial auto-covariance ideas implicit in the variogram. We regard these as alternative and linked views of the data. We conclude that the most important single view of the data is the Map View: All other views must be cross-referred to this, and the software must encourage this. The view can be enriched by overlaying on other pertinent spatial information. We draw attention to the possibilities of one-many linking, and to the use of line-objects to link pairs of data points. We draw attention to the parallels with work on Geographical Information Systems.     

OPTIMAL DESIGN FOR EARLY‐GENERATION PLANT‐BREEDING TRIALS WITH UNREPLICATED OR PARTIALLY REPLICATED TEST LINES

Some subtle difficulties in optimal design are highlighted by the example of unreplicated field trials laid out on plots with spatial errors defined by uniformity trials. There is a dual problem of the arrangement of control plots and maximizing the number of test‐line entries. A simulation study is conducted by randomizing the allocation of genotypes to the plots of four uniformity trials in accordance with the rules defining a number of competing designs. Results are summarized in terms of the ‘SE ratio’, which reflects the improvement in precision of a given design relative to a completely random design on the same plots. The definition of the SE ratio overcomes problems induced by differential shrinkage and consequent precision of test and control lines. A general result applying to all designs shows a curvilinear improvement in SE ratio with increasing error degrees of freedom of the design. The actual arrangement of check plots is of less importance than their increasing number, which contributes to increasing error degrees of freedom. Overall measures, including expected genetic gain, are used to illustrate the choice of a balance between the total number of test‐line entries and the error degrees of freedom.     

The use of spatial and randomisation-based methods for analysis of trials with treatments randomised into rows and columns

The effect of restricted randomisations on the validity and efficiency of using spatial model as well as more common analysis of variance methods for analysing field trials was examined by simulating yields in agricultural fields with known spatial variation and analysing those using eight different statistical models. Some of the models took into account the restrictions on the randomisation and some matched the true generated spatial variation. Two different types of spatial variations were generated. The results showed that the type I error was controlled if the correct spatial model was applied or if the model took into account the restrictions on the randomisation. If the model did not include the correct spatial model nor reflected the randomisation, the type I error could be much higher than the nominal value. The power of the tests could be increased if the model paid attention to the spatial variation. However, it may not always be possible to find the correct and true model. Therefore it is recommended that the random part of the model should reflect both the randomisation process and include terms that take into account spatial variation in the data in order to increase the power.