Reliability estimation in stress–strength models: an MCMC approach

In this paper, we consider the estimation of the stress–strength parameter R=P(Y<X) when X and Y are independent and both are modified Weibull distributions with the common two shape parameters but different scale parameters. The Markov Chain Monte Carlo sampling method is used for posterior inference of the reliability of the stress–strength model. The maximum-likelihood estimator of R and its asymptotic distribution are obtained. Based on the asymptotic distribution, the confidence interval of R can be obtained using the delta method. We also propose a bootstrap confidence interval of R. The Bayesian estimators with balanced loss function, using informative and non-informative priors, are derived. Different methods and the corresponding confidence intervals are compared using Monte Carlo simulations.     

Non-ergodic martingale estimating functions and related asymptotics

Well-known estimation methods such as conditional least squares, quasilikelihood and maximum likelihood (ML) can be unified via a single framework of martingale estimating functions (MEFs). Asymptotic distributions of estimates for ergodic processes use constant norm (e.g. square root of the sample size) for asymptotic normality. For certain non-ergodic-type applications, however, such as explosive autoregression and super-critical branching processes, one needs a random norm in order to get normal limit distributions. In this paper, we are concerned with non-ergodic processes and investigate limit distributions for a broad class of MEFs. Asymptotic optimality (within a certain class of non-ergodic MEFs) of the ML estimate is deduced via establishing a convolution theorem using a random norm. Applications to non-ergodic autoregressive processes, generalized autoregressive conditional heteroscedastic-type processes, and super-critical branching processes are discussed. Asymptotic optimality in terms of the maximum random limiting power regarding large sample tests is briefly discussed.     

Fitting of self-exciting threshold autoregressive moving average nonlinear time-series model through genetic algorithm and development of out-of-sample forecasts

Self-exciting threshold autoregressive moving average (SETARMA) nonlinear time-series model is considered here. Sufficient conditions for invertibility and stationarity are derived. Parameter estimation algorithm is developed by employing real-coded genetic algorithm stochastic optimization procedure. A significant feature of the work done is that optimal out-of-sample forecasts up to three-step ahead and their forecast error variances are derived analytically. Relevant computer programs are written in statistical analysis system (SAS) and C. As an illustration, annual mackerel catch time-series data are considered. Forecast performance of the fitted model for hold-out data is evaluated by using Naive and Monte Carlo approaches. It is found that optimal out-of-sample forecast values are quite close to actual values and estimated variances are quite close to theoretical values. Superiority of the SETARMA model over the SETAR model for equal predictive ability through Diebold–Mariano test is also established.     

Generalized varying-coefficient single-index model

In this paper, the generalized varying-coefficient single-index model is discussed based on penalized likelihood. All the unknown functions are fitted by penalized spline. The estimates of the unknown parameters and the unknown coefficient functions are obtained and the estimation approach is rapid and computationally stable. Under some mild conditions, the consistency and the asymptotic normality of these resulting estimators are given. Two simulation studies are carried out to illustrate the performance of the estimates. An application of the model to the Hong Kong environmental data further demonstrates the potential of the proposed modelling procedures.     

Model averaging for M-estimation

M-estimation is a widely used technique for robust statistical inference. In this paper, we study model selection and model averaging for M-estimation to simultaneously improve the coverage probability of confidence intervals of the parameters of interest and reduce the impact of heavy-tailed errors or outliers in the response. Under general conditions, we develop robust versions of the focused information criterion and a frequentist model average estimator for M-estimation, and we examine their theoretical properties. In addition, we carry out extensive simulation studies as well as two real examples to assess the performance of our new procedure, and find that the proposed method produces satisfactory results.     

Semiparametric estimation of plane similarities: application to fast computation of aeronautic loads

In the big data era, it is often needed to resolve the problem of parsimonious data representation. In this paper, the data under study are curves and the sparse representation is based on a semiparametric model. Indeed, we propose an original registration model for noisy curves. The model is built transforming an unknown function by plane similarities. We develop a statistical method that allows to estimate the parameters characterizing the plane similarities. The properties of the statistical procedure are studied. We show the convergence and the asymptotic normality of the estimators. Numerical simulations and a real-life aeronautic example illustrate and demonstrate the strength of our methodology.     

Bayes designs for multiple linear regression on the unit sphere

We deal with experimental designs minimizing the mean square error of the linear BAYES estimator for the parameter vector of a multiple linear regression model where the experimental region is the k-dimensional unit sphere. After computing the uniquely determined optimum information matrix, we construct, separately for the homogeneous and the inhomogeneous model, both approximate and exact designs having such an information matrix.     

An Exponential Model for Damage Accumulation

In this article, we present a general model for predicting the fatigue behavior for any stress level and amplitude using the exponential model. Based on the Wöhler field for fixed stress level, a compatibility functional equation enables us to derive the general model with eight parameters. The problem of parameter estimation is then discussed and some methods are described. Some examples are finally presented to illustrate the derived model and the proposed methods of estimation.     

Dirichlet Processes in Nonlinear Mixed Effects Models

In this article, we use two efficient approaches to deal with the difficulty in computing the intractable integrals when implementing Gibbs sampling in the nonlinear mixed effects model (NLMM) based on Dirichlet processes (DP). In the first approach, we compute the Laplace's approximation to the integral for its high accuracy, low cost, and ease of implementation. The second approach uses the no-gaps algorithm of MacEachern and Müller (1998 MacEachern , S. , Müller , P. ( 1998 ). Estimating mixtures of Dirichlet process models . Journal of Computational and Graphical Statistics 7 : 223 – 238 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) to perform Gibbs sampling without evaluating the difficult integral. We apply both approaches to real problems and simulations. Results show that both approaches perform well in density estimation and prediction and are superior to the parametric analysis in that they can detect important model features, such as skewness, long tails, and multimodality, whereas the parametric analysis cannot.     

A comparison of several analysis of variance programs with unequal cell size

Various programs in statistical packages for analysis of variance with unequal cell size give different results to the same data because of nonorthogonality of the main effects and interactions. This paper explains how these programs treat the problem of analysis of variance of unbalanced data.