一般偏差右删失数据下剩余寿命分位数回归

剩余寿命是刻画个体预期寿命的一个重要度量,对剩余寿命的早期研究主要集中在剩余均值上.然而当总体生存函数偏态或厚尾时剩余均值函数可能不存在,因此统计学者建议用剩余寿命分位数来刻画预期寿命.在完全数据和右删失数据下,剩余寿命分位数的建模和理论已经很完善.但是,在实际的调查研究中经常会遇到偏差抽样数据.例如,临床医学中的左截断数据,流行病学中的病例队列抽样数据,医学大型队列研究中的长度偏差抽样数据等等.忽略抽样偏差会导致参数估计有偏和不合理的推断结果.本文考虑一般偏差右删失数据下剩余寿命分位数回归的统计推断问题.首先,我们提出了一个一般偏差右删失数据下的剩余寿命分位数回归模型,并利用一般估计方程方法对模型中的参数进行了估计.针对已有文献常用的删失变量与协变量独立性假设,本文重点考虑了删失变量依赖于协变量场合.其次,由于估计量的渐近方差中涉及非参密度函数,在估计渐近方差时,本文采用Bootstrap方法.最后,数值模拟显示本文提出的方法有限样本性质表现很好.     

长度偏差右删失数据剩余寿命的分位数回归

本文研究长度偏差数据下剩余寿命分位数模型的估计方法,充分考虑有偏抽样机制对模型估计的影响.如果忽略这种有偏性会导致估计产生严重偏差甚至错误的结果.本文首先针对长度偏差右删失数据的剩余寿命分位数提出了对数形式的线性回归模型,对删失变量与协变量独立和不独立的两种情况利用估计方程给出了模型参数的估计.其次,通过经验过程和弱收敛理论给出了参数估计的相合性和渐近正态性.最后,本文对提出的估计方法进行了数值模拟并用该方法对奥斯卡奖数据进行分析.     

对数稳定分布加速寿命试验的贝叶斯分析

对数稳定分布与对数正态分布相比具有厚尾等良好的特性,由于没有显式的密度函数,其参数估计有一定的困难.运用Gibbs抽样方法给出了对数稳定分布的Bayes参数估计,与对数正态分布参数估计进行了比较,并应用到恒加寿命试验数据的分析中.模拟研究说明了该方法的有效性,对数稳定分布作为新的寿命分布是合适的.     

风险度量ES半参数模型估计及其应用

预期不足或称期损(Expected Shortfall,ES)是近几年发展起来的重要风险度量工具,对其进行建模和估计是统计学和金融计量经济学研究的前沿问题之一.本文基于平均剩余寿命模型提出一种ES估计的半参数模型,并使用广义估计方程(GEE)的方法估计参数.同时建立了严平稳α混合相依序列下参数估计的大样本理论.本文模型的意义在于可以研究资产组合的风险来源以及各风险因素对ES大小的影响程度.最后,将本文的模型应用到金融股票市场的风险评估中,结果表明此模型可以对某些金融市场现象作出合理的解释,是一个灵活且合理的金融计量统计模型.     

左截断右删失数据下离散型寿命失效率变点模型的Bayes参数估计

利用EM算法和MCMC方法得到了左截断右删失数据下离散型寿命失效率变点模型的参数估计.利用筛选法对缺失数据进行填充,对各参数进行Gibbs抽样.随机模拟证实方法可行且参数估计的精度较高.     

长度偏差右删失数据下比例均值剩余寿命模型的估计方法

本文考虑了长度偏差右删失数据下均值剩余寿命模型的统计推断.当截断变量满足平稳性假设时,长度偏差右删失数据比左截断右删失数据具有更多的信息.为了提高参数估计的效率,我们在估计方程构造中添加了额外信息,通过组合方法获得了新的估计.模拟研究的结果也表明,组合估计方程的方法比仅考虑左截断右删失数据的方法更有效,结果表现更好.     

长度偏差右删失数据下[1mm]比例均值剩余寿命模型的估计方法

本文考虑了长度偏差右删失数据下均值剩余寿命模型的统计推断.当截断变量满足平稳性假设时，长度偏差右删失数据比左截断右删失数据具有更多的信息.为了提高参数估计的效率，我们在估计方程构造中添加了额外信息，通过组合方法获得了新的估计.模拟研究的结果也表明，组合估计方程的方法比仅考虑左截断右删失数据的方法更有效，结果表现更好.     

Weibull 分布基于恒加试验尺度参数估计

对 Weibull 分布场合恒加寿命试验的参数估计进行了讨论,在现有参数估计的基础上作了进一步改进,从而使改进后的参数估计更优.     

左截断右删失数据下百分剩余寿命函数的非参数估计

在左截断右删失数据的模型中，文章讨论３了可靠性中一类重要的α－百分剩余寿命函数的非参数估计，证明了该估计的强一致相合性并获得了该仗垢弱收敛性结果。     

Weibull分布基于恒加试验的参数估计

讨论了Weibull分布场合恒加寿命试验的参数估计,在文[1]对参数估计改进的基础上作了进一步改进,从而使改进后的参数估计更优.     

Rolling element bearing diagnosis based on probability box theory

Feature extraction leads to the loss of statistical information of raw data and ignores the sampling uncertainty and the fluctuations in the signal over time in mechanical fault diagnosis. In this paper, novel modeling methods for mechanical signals based on probability box theory were proposed to solve the above problem. First, the type of random distribution of the bearing signals were analyzed. Then, a Dempster-Shafer structure was obtained to establish a probability box model. To address the identification difficulty of the type of random distribution for the bearing signals, a second probability box model was established based on a vector consisting of features from the bearing signals. If the data are not found to follow a random distribution, a third modeling method based on the definition of probability boxes was proposed. The effectiveness and applicability of the three proposed models were compared with experimental data from rolling element bearings. The combination of probability box theory and mechanical fault diagnosis theory can open up a new research direction for mechanical fault diagnosis.     

A semiparametric panel approach to mortality modeling

During the past twenty years, there has been a rapid growth in life expectancy and an increased attention on funding for old age. Attempts to forecast improving life expectancy have been boosted by the development of stochastic mortality modeling, for example the Cairns–Blake–Dowd (CBD) 2006 model. The most common optimization method for these models is maximum likelihood estimation (MLE) which relies on the assumption that the number of deaths follows a Poisson distribution. However, several recent studies have found that the true underlying distribution of death data is overdispersed in nature (see Cairns et al. 2009 and Dowd et al. 2010). Semiparametric models have been applied to many areas in economics but there are very few applications of such models in mortality modeling. In this paper we propose a local linear panel fitting methodology to the CBD model which would free the Poisson assumption on number of deaths. The parameters in the CBD model will be considered as smooth functions of time instead of being treated as a bivariate random walk with drift process in the current literature. Using the mortality data of several developed countries, we find that the proposed estimation methods provide comparable fitting results with the MLE method but without the need of additional assumptions on number of deaths. Further, the 5-year-ahead forecasting results show that our method significantly improves the accuracy of the forecast.     

Statistical Analysis Technique of Two-parameter Generalized Exponential Sum Distribution

A new life distribution is proposed, known as ``two-parameter generalized exponential sum distribution". We study the density function and failure rate function, the average failure rate function, the image features and the numerical characteristics of the mean residual life of the distribution. Several methods of calculating point estimation of parameters are discussed. Through the Monte-Carlo simulation, we compare the precision of the point estimations. In our opinion, the best linear unbiased estimation is the most optimal solution of these methods. At the same time, several methods of calculating parameters of interval estimations are given. We also discuss the precision of interval estimations by Monte-Carlo simulation and use the best linear unbiased estimation and the best linear invariant estimation to construct interval estimations which are better than other estimation method. Finally, several simulation examples and a case of maintaining tanks is used to illustrate the application of the methods presented in this paper.     

A two-stage prognosis model in condition based maintenance

We often observe in practice that the life of a piece of production equipment can be divided into two stages. The first stage is referred to as the normal working stage where no significant deviation from the normal operating state is observed. The second stage is called the failure delay period, since a defect may be initiated, and progressively develop into an actual failure, i.e., the equipment is in a defective stage but still working during this stage. With the help of condition monitoring, hidden defects already present in the equipment may be detected, but for maintenance planning purposes, the prediction of the initiation point of the second stage, and more importantly, the residual life thereafter is important. This paper reports on the development of a probability model to predict the initiation point of the second stage and the remaining life based on available condition monitoring information. The method for model parameters estimation is discussed and applied to real data.     

������ɷ��Ԫ��̬�ֲ������ļ�����Ȼ����

??Hidden Markov model is widely used in statistical modeling of time, space and state transition data. The definition of hidden Markov multivariate normal distribution is given. The principle of using cluster analysis to determine the hidden state of observed variables is introduced. The maximum likelihood estimator of the unknown parameters in the model is derived. The simulated observation data set is used to test the estimation effect and stability of the method. The characteristic is simple classical statistical inference such as cluster analysis and maximum likelihood estimation. The method solves the parameter estimation problem of complex statistical models.     

Maximum Likelihood Estimation of Hidden Markov Multivariate Normal Distribution Parameters

Hidden Markov model is widely used in statistical modeling of time, space and state transition data. The definition of hidden Markov multivariate normal distribution is given. The principle of using cluster analysis to determine the hidden state of observed variables is introduced. The maximum likelihood estimator of the unknown parameters in the model is derived. The simulated observation data set is used to test the estimation effect and stability of the method. The characteristic is simple classical statistical inference such as cluster analysis and maximum likelihood estimation. The method solves the parameter estimation problem of complex statistical models.     

截断回归模型参数的非参数最小二乘估计

本文考虑截断回归模型,给出了基于截断数据估计回归参数的一种新方法,此处并不设定残差分布.我们使用早先的关于误差分布非参数估计的结果,在某些正则条件下建立了估计量的相合性.并给出实例说明我们的结果是Heckman(1979)-项工作的本质改进.     

分位数变系数模型基于核光滑的变量选择

分位数变系数模型是一种稳健的非参数建模方法.使用变系数模型分析数据时,一个自然的问题是如何同时选择重要变量和从重要变量中识别常数效应变量.本文基于分位数方法研究具有稳健和有效性的估计和变量选择程序.利用局部光滑和自适应组变量选择方法,并对分位数损失函数施加双惩罚,我们获得了惩罚估计.通过BIC准则合适地选择调节参数,提出的变量选择方法具有oracle理论性质,并通过模拟研究和脂肪实例数据分析来说明新方法的有用性.数值结果表明,在不需要知道关于变量和误差分布的任何信息前提下,本文提出的方法能够识别不重要变量同时能区分出常数效应变量.     

Optimum design of compressive residual stress field caused by ultrasonic surface rolling with a mathematical model

Compressive residual stress fields induced by ultrasonic surface rolling process play a key role in determining the fatigue performance of machine parts. The present work is an analytical approach to conducting optimum design of this field to obtain an optimum fatigue resistance. Thus, a mathematical model was presented to predict residual stresses based on circular and elliptical Hertz contact areas. Moreover, to validate the proposed analytical approach, experimental verification was carried out on 18CrNiMo7-6 steel. Analytic solutions were derived from the mathematical model and optimum characteristic parameters were obtained by investigating the characteristic parameters in this field, such as surface residual stress, maximum residual stress and its depths. Results showed that increasing the total force, Hertz contact area and ratio of the radius of tool tip to that of target body could significantly enhance the peak of compressive residual stress and its depth.     

Efficient and robust estimation for autoregressive regression models using shape mixtures of skewt normal distribution

Multiple linear regression model based on normally distributed and uncorrelated errors is a popular statistical tool with application in various fields. But these assumptions of normality and no serial correlation are hardly met in real life. Hence, this study considers the linear regression time series model for series with outliers and autocorrelated errors. These autocorrelated errors are represented by a covariance-stationary autoregressive process where the independent innovations are driven by shape mixture of skew-t normal distribution. The shape mixture of skew-t normal distribution is a flexible extension of the skew-t normal with an additional shape parameter that controls skewness and kurtosis. With this error model, stochastic modeling of multiple outliers is possible with an adaptive robust maximum likelihood estimation of all the parameters. An Expectation Conditional Maximization Either algorithm is developed to carryout the maximum likelihood estimation. We derive asymptotic standard errors of the estimators through an information-based approximation. The performance of the estimation procedure developed is evaluated through Monte Carlo simulations and real life data analysis.