一种基于加权F-范数的半正定矩阵的逼近方法

正定性是许多金融预测模型的重要假设前提，然而从实际样本中得到的相关系数矩阵并不能保证其正定性。为此在介绍如何根据样本设定相关系数矩阵以及范数逼近原理的基础上，如何根据该原理找到与之最接近的相关系数矩阵，即最接近的单位对角半正定对称矩阵。通过实证，验证了其方法的有效性。     

On local influence for elliptical linear models

The local influence method plays an important role in regression diagnostics and sensitivity analysis. To implement it, we need the Delta matrix for the underlying scheme of perturbations, in addition to the observed information matrix under the postulated model. Galea, Paula and Bolfarine (1997) has recently given the observed information matrix and the Delta matrix for a scheme of scale perturbations and has assessed of local influence for elliptical linear regression models. In the present paper, we consider the same elliptical linear regression models. We study the schemes of scale, predictor and response perturbations, and obtain their corresponding Delta matrices, respectively. To illustrate the methodology for assessment of local influence for these schemes and the implementation of the obtained results, we give an example.     

Diagonalization matrix and its application in distribution theory

Some matrix representations of diverse diagonal arrays are studied in this work; the results allow new definitions of classes of elliptical distributions indexed by kernels mixing Hadamard and usual products. A number of applications are derived in the setting of prior densities from the Bayesian multivariate regression model and families of non-elliptical distributions, such as the matrix multivariate generalized Birnbaum–Saunders density. The philosophy of the research about matrix representations of quadratic and inverse quadratic forms can be extended as a methodology for exploring possible new applications in non-standard distributions, matrix transformations and inference.     

A NOTE ON THE INVERSE COVARIANCE MATRIX OF THE AUTOREGRESSIVE PROCESS1

The inverse covariance matrix of an autoregressive process of arbitrary order is found explicitly using the notion of the information matrix for the vector random variable, not the parameters. Any process for which a simple conditional representation exists, can be treated in the same way. The approach has merit in the teaching of statistics.     

Zonal Polynomials and Hypergeometric Functions of Quaternion Matrix Argument

We define zonal polynomials of quaternion matrix argument and deduce some impor-tant formulae of zonal polynomials and hypergeometric functions of quaternion matrix argument. As an application, we give the distributions of the largest and smallest eigenvalues of a quaternion central Wishart matrix W ～ ?W(n, Σ), respectively.     

On testing for an identity covariance matrix when the dimensionality equals or exceeds the sample size

This article explores the problem of testing the hypothesis that the covariance matrix is an identity matrix when the dimensionality is equal to the sample size or larger. Two new test statistics are proposed under comparable assumptions to those statistics in the literature. The asymptotic distribution of the proposed test statistics are found and are shown to be consistent in the general asymptotic framework. An extensive simulation study shows the newly proposed tests are comparable to, and in some cases more powerful than, the tests for an identity covariance matrix currently in the literature.     

Classification of difference matrices over cyclic groups

The existence of difference matrices over small cyclic groups is investigated in this computer-aided work. The maximum values of the parameters for which difference matrices exist as well as the number of inequivalent difference matrices in each case is determined up to the computational limit. Several new difference matrices have been found in this manner. The maximum number of rows is 9 for an r ×15 difference matrix over Z₃, 8 for an r ×15 difference matrix over Z₅, and 6 for an r ×12 difference matrix over Z₆; the number of inequivalent matrices with these parameters is 5, 2, and 7, respectively.     

Modified D-optimal design for logistic model

Asymptotic theory of using the Fisher information matrix may provide poor approximation to the exact variance matrix of maximum likelihood estimation in nonlinear models. This may be due to not obtaining an efficient D-optimal design. In this article, we propose a modified D-optimality criterion, using a more accurate information matrix, based on the Bhattacharyya matrix. The proposed information matrix and its properties are given for two parameters simple logistic model. It is shown that the resulted modified locally D-optimal design is more efficient than the previous one; particularly, for small sample size experiments.     

ESTIMATION OF REALIZED SIGNAL TO NOISE RATIO FOR A PAIR OF MULTIVARIATE SIGNALS

Consider the case of classifying an incoming message as one of two known p-dimension signals or as a pure noise. Let the noise co-variance matrix (assumed to be same in all the three cases) be unknown. We consider the problem of estimation of “realized signal to noise ratio matrix”, which is an index of discriminatory power, under various loss functions. Optimum estimators are obtained under these loss functions. Finally, an attempt is made to provide a lower confidence bound for the realized signal to noise ratio matrix. In the process, the probability distribution of the smaller eigenvalue of a 2 × 2 confluent hypergeometric random matrix is obtained.     

On the relationship between the matrix operators,vech and vecd

We introduce a matrix operator, which we call “vecd” operator. This operator stacks up “diagonals” of a symmetric matrix. This operator is more convenient for some statistical analyses than the commonly used “vech” operator. We show an explicit relationship between the vecd and vech operators. Using this relationship, various properties of the vecd operator are derived. As applications of the vecd operator, we derive concise and explicit expressions of the Wald and score tests for equal variances of a multivariate normal distribution and for the diagonality of variance coefficient matrices in a multivariate generalized autoregressive conditional heteroscedastic (GARCH) model, respectively.     

基于集值统计的区间综合评判及其应用

综合评判是根据评判对象的评判因素,采用综合评判模型对评判对象作出一个综合判断。为了克服传统综合评判的不足,研究了基本集值统计的区间综合评判,探讨了集值统计区间分析处理方法。运用了集值统计判断矩阵计算因素权重区间;集值统计和属性测度构建评判矩阵;综合评判区间构建可能度矩阵计算相对权重;得到了综合评判区间排序。该综合评判的优点在于充分体现了评价者的真实意图,充分融合了集值统计区间分析处理多种方法的优点。通过案例分析,基于集值统计的区间综合评判是可行的,具有较强的操作性。     

资金流量矩阵表预测中的DRAS法研究

在投入产出分析中RAS法已被广泛应用,主要用于对预测期的投入产出消耗系数矩阵进行修订和预测。随着研究的深入,RAS法有了很多改进,但如果需要两个矩阵共同平衡预测的时候,常规的RAS法或改进的RAS法都难以完成,而对资金流量矩阵延长表的预测,正是需要同时预测两个平衡矩阵的情况。鉴于此,讨论了一种双矩阵RAS(Double Matrix RAS,DRAS)的平衡方法,对这种方法进行了数学表述,并就其在资金流量矩阵表预测中的应用给出了说明。     

A counterexample to Beder's conjectures about Hadamard matrices

In this note we provide a counterexample which resolves conjectures about Hadamard matrices made in this journal. Beder [1998. Conjectures about Hadamard matrices. Journal of Statistical Planning and Inference 72, 7–14] conjectured that if H$\mathit{H}$ is a maximal m×n$m\times n$ row-Hadamard matrix then m is a multiple of 4; and that if n   is a power of 2 then every row-Hadamard matrix can be extended to a Hadamard matrix. Using binary integer programming we obtain a maximal 13×32$13\times 32$ row-Hadamard matrix, which disproves both conjectures. Additionally for n being a multiple of 4 up to 64, we tabulate values of m   for which we have found a maximal row-Hadamard matrix. Based on the tabulated results we conjecture that a m×n$m\times n$ row-Hadamard matrix with m?n-7$m?n-7$ can be extended to a Hadamard matrix.     

Maximum likelihood estimation of the linearly structured correlation matrix by a Jacobi-type iterative scheme

Mukherjee and Maiti [Q-procedure for solving likelihood equations in the analysis of covariance structures, Comput. Statist. Quart. 2 (1988), pp. 105–128] proposed an iterative scheme to derive the maximum likelihood estimates of the parameters involved in the population covariance matrix when it is linearly structured. The present investigation provides a Jacobi-type of iterative scheme, MSIII, when the underlying correlation matrix is linearly structured. Such scheme is shown to be quite competent and efficient compared to the prevalent Fisher-scoring (FS) and the Newton–Raphson iterative scheme (NR). An illustrative example is provided for a numerical comparison of the iterates of MSIII, FS and NR choosing the Toeplitz matrix as the population correlation matrix. Numerical behaviour of such schemes is studied in the context of ‘bad’ initial try-out vectors. Additionally a simulation experiment is performed to judge the superiority of MSIII over FS.     

Rank covariance matrix estimation of a partially known covariance matrix

Classical multivariate methods are often based on the sample covariance matrix, which is very sensitive to outlying observations. One alternative to the covariance matrix is the affine equivariant rank covariance matrix (RCM) that has been studied in Visuri et al. [2003. Affine equivariant multivariate rank methods. J. Statist. Plann. Inference 114, 161–185]. In this article we assume that the covariance matrix is partially known and study how to estimate the corresponding RCM. We use the properties that the RCM is affine equivariant and that the RCM is proportional to the inverse of the regular covariance matrix, and hence reduce the problem of estimating the original RCM to estimating marginal rank covariance matrices. This is a great computational advantage when the dimension of the original data vector is large.     

Unbiased Estimator for a Covariance Matrix Under Two-Step Monotone Incomplete Sample

In this article, we consider an inference for a covariance matrix under two-step monotone incomplete sample. The maximum likelihood estimator of the mean vector is unbiased but that of the covariance matrix is biased. We derive an unbiased estimator for the covariance matrix using some fundamental properties of the Wishart matrix. The properties of the estimators are investigated and the accuracies are checked by a numerical simulation.     

Profile information matrix for nonlinear transformation models

For semiparametric models, interval estimation and hypothesis testing based on the information matrix for the full model is a challenge because of potentially unlimited dimension. Use of the profile information matrix for a small set of parameters of interest is an appealing alternative. Existing approaches for the estimation of the profile information matrix are either subject to the curse of dimensionality, or are ad-hoc and approximate and can be unstable and numerically inefficient. We propose a numerically stable and efficient algorithm that delivers an exact observed profile information matrix for regression coefficients for the class of Nonlinear Transformation Models [A. Tsodikov (2003) J R Statist Soc Ser B 65:759-774]. The algorithm deals with the curse of dimensionality and requires neither large matrix inverses nor explicit expressions for the profile surface.     

On the independence between two generalized second degree polynomial statistics and their covariance

The condition of independence between two generalized second degree polynomial statistics is worked out in this paper under the assumption that the covariance matrix of the underlying normal distribution may be singular. The covariance between two such forms is also worked out without the assumption of normality. The results are indicated in compact matrix form.     

On the use of the selection matrix in the maximum likelihood estimation of normal distribution models with missing data

In this article, by using the constant and random selection matrices, several properties of the maximum likelihood (ML) estimates and the ML estimator of a normal distribution with missing data are derived. The constant selection matrix allows us to obtain an explicit form of the ML estimates and the exact relationship between the EM algorithm and the score function. The random selection matrix allows us to clarify how the missing-data mechanism works in the proof of the consistency of the ML estimator, to derive the asymptotic properties of the sequence by the EM algorithm, and to derive the information matrix.     

The Density Functions of the Singular Quaternion Normal Matrix and the Singular Quaternion Wishart Matrix

In this article, we give the density functions of the singular quaternion normal matrix and the singular quaternion Wishart matrix. Furthermore, we also give the density functions of certain singular quaternion β-matrix and the singular quaternion F-matrix in terms of the density function of the singular quaternion Wishart matrix and hypergeometric functions of quaternion matrix argument.