Bayesian analysis of mixture modelling using the multivariate t distribution

A finite mixture model using the multivariate t distribution has been shown as a robust extension of normal mixtures. In this paper, we present a Bayesian approach for inference about parameters of t-mixture models. The specifications of prior distributions are weakly informative to avoid causing nonintegrable posterior distributions. We present two efficient EM-type algorithms for computing the joint posterior mode with the observed data and an incomplete future vector as the sample. Markov chain Monte Carlo sampling schemes are also developed to obtain the target posterior distribution of parameters. The advantages of Bayesian approach over the maximum likelihood method are demonstrated via a set of real data.     

Estimation of Hazard, Density and Survivor Functions for Randomly Censored Data

Maximum likelihood estimation and goodness-of-fit techniques are used within a competing risks framework to obtain maximum likelihood estimates of hazard, density, and survivor functions for randomly right-censored variables. Goodness-of- fit techniques are used to fit distributions to the crude lifetimes, which are used to obtain an estimate of the hazard function, which, in turn, is used to construct the survivor and density functions of the net lifetime of the variable of interest. If only one of the crude lifetimes can be adequately characterized by a parametric model, then semi-parametric estimates may be obtained using a maximum likelihood estimate of one crude lifetime and the empirical distribution function of the other. Simulation studies show that the survivor function estimates from crude lifetimes compare favourably with those given by the product-limit estimator when crude lifetimes are chosen correctly. Other advantages are discussed.     

On Optimality of Bayesian Wavelet Estimators

Abstract.  We investigate the asymptotic optimality of several Bayesian wavelet estimators, namely, posterior mean, posterior median and Bayes Factor, where the prior imposed on wavelet coefficients is a mixture of a mass function at zero and a Gaussian density. We show that in terms of the mean squared error, for the properly chosen hyperparameters of the prior, all the three resulting Bayesian wavelet estimators achieve optimal minimax rates within any prescribed Besov space     for p  ≥ 2. For 1 ≤  p  < 2, the Bayes Factor is still optimal for (2 s +2)/(2 s +1) ≤  p  < 2 and always outperforms the posterior mean and the posterior median that can achieve only the best possible rates for linear estimators in this case.     

Robust testing with generalized partial linear models for longitudinal data

By approximating the nonparametric component using a regression spline in generalized partial linear models (GPLM), robust generalized estimating equations (GEE), involving bounded score function and leverage-based weighting function, can be used to estimate the regression parameters in GPLM robustly for longitudinal data or clustered data. In this paper, score test statistics are proposed for testing the regression parameters with robustness, and their asymptotic distributions under the null hypothesis and a class of local alternative hypotheses are studied. The proposed score tests reply on the estimation of a smaller model without the testing parameters involved, and perform well in the simulation studies and real data analysis conducted in this paper.     

低浓度三分子反应模型双曲型反应－扩散方程的稳定性与化学振荡

考虑浓度扩散和扩散流随时间的变化，得到低浓度三分子反应模型的双曲型反应－扩散方程；讨论了方程的稳定性及存在行波解的条件，得到了化学振荡的色散关系和相速度．     

将分段函数表示为初等函数的方法

本文给出了一种化具有多个分段点，且在各区间段上的表达式为多项式函数的分段函数为初等函数的具体方法．     

二阶非线性微分方程解的振动性和渐近性

讨论了二阶非线性常微分方程(α(t)(?)(x)x′)′+q(t)f(x)=r(t)的解的振动性和渐近性.获得了有关该方程的五个新的定理.     

A study of simulated annealing and a revised cascade algorithm for image reconstruction

We describe an image reconstruction problem and the computational difficulties arising in determining the maximum a posteriori (MAP) estimate. Two algorithms for tackling the problem, iterated conditional modes (ICM) and simulated annealing, are usually applied pixel by pixel. The performance of this strategy can be poor, particularly for heavily degraded images, and as a potential improvement Jubb and Jennison (1991) suggest the cascade algorithm in which ICM is initially applied to coarser images formed by blocking squares of pixels. In this paper we attempt to resolve certain criticisms of cascade and present a version of the algorithm extended in definition and implementation. As an illustration we apply our new method to a synthetic aperture radar (SAR) image. We also carry out a study of simulated annealing, with and without cascade, applied to a more tractable minimization problem from which we gain insight into the properties of cascade algorithms.     

Estimation of noncentrality parameters

We propose some estimators of noncentrality parameters which improve upon usual unbiased estimators under quadratic loss. The distributions we consider are the noncentral chi-square and the noncentral F. However, we give more general results for the family of elliptically contoured distributions and propose a robust dominating estimator.     

Some remarks on the physical models concerning two different approaches to inference in statistical process control

Summary For technological applications it can be useful to identify some simple physical mechanisms, which, on the basis of the available knowledge of the production process, may suggest the most appropriate approach to statistical control of the random quantities of interest. For this purpose the notion of rupture point is introduced firstly. A rupture point is characterized bym randomly arising out of control states, assumed to be mutually exclusive and stochastically independent. Shewhart's control charts seem to represent the natural statistical tool for controlling a rupture point; however it is shown that they are fully justified only when the hazard rates attached to the causes of failure are constant. Otherwise, typically in the presence of time increasing hazard rates, Shewhart's control charts should be completed by a preventive intervention rule (preventive maintenance). In the second place, the notion of dynamic instability point is introduced, which is specifically characterized by assuming that the random quantity of interest is ruled by a stochastic differential equation with constant coefficients. By discretization, developed according to a possibly new approach, it is shown that the former model reduces to an equation error model, which is among the simplest used in adaptive control, and thus particularly easy to deal with in regard to parameter estimation and the definition of the optimum control rule.