块Jacobi迭代阵的收敛性

研究大型线性方程组迭代解法中分块ＪＡＣＯＢＩ迭代阵的收敛性。采用块矩阵分析方法和谱半径降维估计法得到块Ｊａｃｏｂｉ迭代阵收敛的实用充分条件。     

Generalized Leverage and its Applications

The generalized leverage of an estimator is defined in regression models as a measure of the importance of individual observations. We derive a simple but powerful result, developing an explicit expression for leverage in a general M -estimation problem, of which the maximum likelihood problems are special cases. A variety of applications are considered, most notably to the exponential family non-linear models. The relationship between leverage and local influence is also discussed. Numerical examples are given to illustrate our results     

A Unified Approach to the Approximation of Multivariate Densities

Approximation of a density by another density is considered in the case of different dimensionalities of the distributions. The results have been derived by inverting expansions of characteristic functions with the help of matrix techniques. The approximations obtained are all functions of cumulant differences and derivatives of the approximating density. The multivariate Edgeworth expansion follows from the results as a special case. Furthermore, the density functions of the trace and eigenvalues of the sample covariance matrix are approximated by the multivariate normal density and a numerical example is given     

Standard errors for EM estimation

The EM algorithm is a popular method for computing maximum likelihood estimates. One of its drawbacks is that it does not produce standard errors as a by-product. We consider obtaining standard errors by numerical differentiation. Two approaches are considered. The first differentiates the Fisher score vector to yield the Hessian of the log-likelihood. The second differentiates the EM operator and uses an identity that relates its derivative to the Hessian of the log-likelihood. The well-known SEM algorithm uses the second approach. We consider three additional algorithms: one that uses the first approach and two that use the second. We evaluate the complexity and precision of these three and the SEM in algorithm seven examples. The first is a single-parameter example used to give insight. The others are three examples in each of two areas of EM application: Poisson mixture models and the estimation of covariance from incomplete data. The examples show that there are algorithms that are much simpler and more accurate than the SEM algorithm. Hopefully their simplicity will increase the availability of standard error estimates in EM applications. It is shown that, as previously conjectured, a symmetry diagnostic can accurately estimate errors arising from numerical differentiation. Some issues related to the speed of the EM algorithm and algorithms that differentiate the EM operator are identified.     

Toeplitz矩阵的性质及其应用

运用矩阵方法,研究了特殊的Toeplitz矩阵A和n阶幂零矩阵在相似变换下的广义逆,得到了Toeplitz矩阵的一些性质,并用这些性质刻画了幂零矩阵A在相似变换下的广义逆的表达形式.     

Analysis of some ruin-related quantities in a Markov-modulated risk model

In this paper, we study the joint Laplace transform and probability generating function of some random quantities that occur in each environment state by the time of ruin in a Markov-modulated risk process. These quantities include the duration spent in each state, the number of claims and the aggregate amount of claims that occurred in each state by the time of ruin. Explicit formulae for the joint transforms, given the initial surplus, and the initial and terminal environment states, are expressed in terms of a matrix version of the scale function. Moments and covariances of these ruin-related quantities are obtained and numerical illustrations are presented. The joint transform of the duration spent in each state, the number of claims, and the aggregate amount of claims that occurred in each state by the time the surplus attains a certain level are also investigated.     

Multivariate process dispersion monitoring without subgrouping

The memory-type adaptive and non-adaptive control charts are among the best control charts for detecting small-to-moderate changes in the process parameter(s). In this paper, we propose the Crosier CUSUM (CCUSUM), EWMA, adaptive CCUSUM (ACCUSUM) and adaptive EWMA (AEWMA) charts for efficiently monitoring the changes in the covariance matrix of a multivariate normal process without subgrouping. Using extensive Monte Carlo simulations, the length characteristics of these control charts are computed. It turns out that the ACCUSUM and AEWMA charts perform uniformly and substantially better than the CCUSUM and EWMA charts when detecting a range of shift sizes in the covariance matrix. Moreover, the AEWMA chart outperforms the ACCUSUM chart. A real dataset is used to explain the implementation of the proposed control charts.     

Accurate computation of single and product moments of order statistics under progressive censoring

An accurate procedure is proposed to calculate approximate moments of progressive order statistics in the context of statistical inference for lifetime models. The study analyses the performance of power series expansion to approximate the moments for location and scale distributions with high precision and smaller deviations with respect to the exact values. A comparative analysis between exact and approximate methods is shown using some tables and figures. The different approximations are applied in two situations. First, we consider the problem of computing the large sample variance–covariance matrix of maximum likelihood estimators. We also use the approximations to obtain progressively censored sampling plans for log-normal distributed data. These problems illustrate that the presented procedure is highly useful to compute the moments with precision for numerous censoring patterns and, in many cases, is the only valid method because the exact calculation may not be applicable.