Approximate bayes estimation of reliability of two parameter inverse gaussian distribution

This paper extends the result of Padgett (1981) and gives a Bayes estimate of the reliability function of two-parameter inverse Gaussian distribution using Jeffreys' non-informative joint prior and a squared error loss fun ction . A numerical example is given. Based on a Monte Carlo simulation, Bayes estimator of reliability is compared with its maximum likelihood counterpart.     

On multiple change-point estimation for Poisson process

This work is devoted to the problem of change-point parameter estimation in the case of the presence of multiple changes in the intensity function of the Poisson process. It is supposed that the observations are independent inhomogeneous Poisson processes with the same intensity function and this intensity function has two jumps separated by a known quantity. The asymptotic behavior of the maximum-likelihood and Bayesian estimators are described. It is shown that these estimators are consistent, have different limit distributions, the moments converge and that the Bayesian estimators are asymptotically efficient. The numerical simulations illustrate the obtained results.     

Sequential estimation for a family of counting processes in the nuisance parameter case

The paper deals with the problem of sequential estimation for stochastic processes in the presence of a nuisance parameter. Using the approach to estimation through estimating equations, optimum estimating functions based on a random observation time are investigated in some models for processes appearing in reliability systems theory.     

Software reliability assessment using accelerated testing methods

The use of operational profiles and usage-based testing has received considerable attention recently in the software engineering literature. Testing under the actual operational profile can, however, be expensive, time consuming or even infeasible in situations where the performance of a system is dominated by infrequent but highly critical events. We consider a real application that deals with telecommunications network restoration after network failure caused by cuts in fibre optic cables. We use this application to demonstrate the usefulness of traditional accelerated testing methods to test and estimate software reliability. These methods, which have been extensively used in hardware reliability, have an important role to play in software reliability assessment as well.     

Pairwise likelihood estimation for multivariate mixed Poisson models generated by Gamma intensities

Estimating the parameters of multivariate mixed Poisson models is an important problem in image processing applications, especially for active imaging or astronomy. The classical maximum likelihood approach cannot be used for these models since the corresponding masses cannot be expressed in a simple closed form. This paper studies a maximum pairwise likelihood approach to estimate the parameters of multivariate mixed Poisson models when the mixing distribution is a multivariate Gamma distribution. The consistency and asymptotic normality of this estimator are derived. Simulations conducted on synthetic data illustrate these results and show that the proposed estimator outperforms classical estimators based on the method of moments. An application to change detection in low-flux images is also investigated.     

Robust Poisson likelihood estimation for frailty Cox models: A simulation study

Frailty models can be fit as mixed-effects Poisson models after transforming time-to-event data to the Poisson model framework. We assess, through simulations, the robustness of Poisson likelihood estimation for Cox proportional hazards models with log-normal frailties under misspecified frailty distribution. The log-gamma and Laplace distributions were used as true distributions for frailties on a natural log scale. Factors such as the magnitude of heterogeneity, censoring rate, number and sizes of groups were explored. In the simulations, the Poisson modeling approach that assumes log-normally distributed frailties provided accurate estimates of within- and between-group fixed effects even under a misspecified frailty distribution. Non-robust estimation of variance components was observed in the situations of substantial heterogeneity, large event rates, or high data dimensions.     

A New Goodness-of-Fit Test Based on the Laplace Statistic for a Large Class of NHPP Models

ABSTRACT

Nonhomogeneous Poisson processes (NHPP) provide many models for hardware and software reliability analysis. In order to get an appropriate NHPP model, goodness-of-Fit (GOF for short) tests have to be carried out. For the power-law processes, lots of GOF tests have been developed. For other NHPP models, only the Conditional Probability Integral Transformation (CPIT) test has been proposed. However, the CPIT test is less powerful and cannot be applied to some NHPP models. This article proposes a general GOF test based on the Laplace statistic for a large class of NHPP models with intensity functions of the form αλ(t, β). The simulation results show that this test is more powerful than CPIT test.     

Bayes estimation for exponential distributions with common location parameter and applications to multi-state reliability models

This paper considers the estimation of the stress–strength reliability of a multi-state component or of a multi-state system where its states depend on the ratio of the strength and stress variables through a kernel function. The article presents a Bayesian approach assuming the stress and strength as exponentially distributed with a common location parameter but different scale parameters. We show that the limits of the Bayes estimators of both location and scale parameters under suitable priors are the maximum likelihood estimators as given by Ghosh and Razmpour [15 M. Ghosh and A. Razmpour, Estimation of the common location parameter of several exponentials, Sankhyā, Ser. A 46 (1984), pp. 383–394. [Google Scholar]]. We use the Bayes estimators to determine the multi-state stress–strength reliability of a system having states between 0 and 1. We derive the uniformly minimum variance unbiased estimators of the reliability function. Interval estimation using the bootstrap method is also considered. Under the squared error loss function and linex loss function, risk comparison of the reliability estimators is carried out using extensive simulations.     

Isotonic estimation of the intensity of a nonhomogeneous Poisson process: The multiple realization setup

In the 1950s Brunk and Van Eeden each obtained maximum-likelihood estimators of a finite product of probability density functions under partial or complete ordering of their parameters. Their results play an important role in the general theory of inference under order restrictions and lead to an isotonic estimator of the intensity of a nonhomogeneous Poisson process. Here an elementary derivation of the maximum likelihood estimator (m.l.e.) for the intensity of a nonhomogeneous Poisson process is given when several (possibly censored) realizations are available. Boswell obtained the m.l.e. based on a single realization as well as a conditional m.l.e. under the same conditions. An example is given to show that in the multirealization setup a conditional m.l.e. may not exist; the proofs are, we believe, new and elementary. An illustrative application is given.     

Improved estimation for Poisson INAR(1) models

We consider the first-order Poisson autoregressive model proposed by McKenzie [Some simple models for discrete variate time series. Water Resour Bull. 1985;21:645–650] and Al-Osh and Alzaid [First-order integer valued autoregressive (INAR(1)) process. J Time Ser Anal. 1987;8:261–275], which may be suitable in situations where the time series data are non-negative and integer valued. We derive the second-order bias of the squared difference estimator [Weiß. Process capability analysis for serially dependent processes of Poisson counts. J Stat Comput Simul. 2012;82:383–404] for one of the parameters and show that this bias can be used to define a bias-reduced estimator. The behaviour of a modified conditional least-squares estimator is also studied. Furthermore, we access the asymptotic properties of the estimators here discussed. We present numerical evidence, based upon Monte Carlo simulation studies, showing that the here proposed bias-adjusted estimator outperforms the other estimators in small samples. We also present an application to a real data set.     

Bias correction for parameter estimates of spatial point process models

When a spatial point process model is fitted to spatial point pattern data using standard software, the parameter estimates are typically biased. Contrary to folklore, the bias does not reflect weaknesses of the underlying mathematical methods, but is mainly due to the effects of discretization of the spatial domain. We investigate two approaches to correcting the bias: a Newton–Raphson-type correction and Richardson extrapolation. In simulation experiments, Richardson extrapolation performs best.     

Uniqueness of estimation and identifiability in mixture models

Further properties of the nonparametric maximum-likelihood estimator of a mixing distribution are obtained by exploiting the properties of totally positive kernels. Sufficient conditions for uniqueness of the estimator are given. This result is more general, and the proof is substantially simpler, than given previously. When the component density has support on N points, it is shown that all identifiable mixing distributions have support on no more than N/2 points. Identifiable mixtures are shown to lie on the boundary of the mixture model space. The maximum-likelihood estimate is shown to be unique if the vector of observations lies outside this space.     

Non-homogeneous Markov models in the analysis of survival after breast cancer

A cohort of 300 women with breast cancer who were submitted for surgery is analysed by using a non-homogeneous Markov process. Three states are onsidered: no relapse, relapse and death. As relapse times change over time, we have extended previous approaches for a time homogeneous model to a non omogeneous multistate process. The trends of the hazard rate functions of transitions between states increase and then decrease, showing that a changepoint can be considered. Piecewise Weibull distributions are introduced as transition intensity functions. Covariates corresponding to treatments are incorporated in the model multiplicatively via these functions. The likelihood function is built for a general model with k changepoints and applied to the data set, the parameters are estimated and life-table and transition probabilities for treatments in different periods of time are given. The survival probability functions for different treatments are plotted and compared with the corresponding function for the homogeneous model. The survival functions for the various cohorts submitted for treatment are fitted to the mpirical survival functions.     

Maximum likelihood estimation of ARFIMA models with a Box-Cox transformation

In this paper we study the interaction between the estimation of the fractional differencing parameter d of ARFIMA models and the common practice of instantaneous transformation of the observed time series. At this aim, we first discuss the effect of a nonlinear transformation of the data on the identification of the process and on the estimate of d. Thus, we propose a joint estimation of the Box-Cox parameter and d by means of a modified normalized version of the Whittle likelihood. Then, the variance and covariance matrix of the parameters estimates is obtained. Finally, a Monte Carlo study is performed in order to check the behaviour of the proposed estimators in finite samples.The paper is the result of a joint research of the two authors. As far as it concerns this version of the work, A. DElia wrote Sects. 2, 3, 4, while D. Piccolo wrote Sects. 1, 5, 6.     

Performance of adaptive estimators in slowly varying parameter models

This paper analyzes the MSE of the exponentially weighted least squares (EWLS) estimator in dynamic regression models with time-varying parameters. Under the assumption of differentiable parameter functions, it is derived an asymptotic expression which is the sum of a stationary and of an evolutionary component. The validity of the analytical expression is illustrated with simulation experiments, and its usefulness in designing the exponential discounting factor is illustrated on a real case-study. The practical finding is similar to the plug-in bandwidth selection in nonparametric smoothers.     

General orthogonal designs for parameter estimation of ANOVA models under weighted sum-to-zero constraints

ABSTRACT

Traditional studies on optimal designs for ANOVA parameter estimation are based on the framework of equal probabilities of appearance for each factor's levels. However, this premise does not hold in a variety of experimental problems, and it is of theoretical and practical interest to investigate optimal designs for parameters with unequal appearing odds. In this paper, we propose a general orthogonal design via matrix image, in which all columns’ matrix images are orthogonal with each other. Our main results show that such designs have A- and E-optimalities on the estimation of ANOVA parameters which have unequal appearing odds. In addition, we develop two simple methods to construct the proposed designs. The optimality of the design is also validated by a simulation study.     

Some remarks on classical and bayesian reliability estimation of binomial and poisson distributions

The problems of estimating the reliability function and Pr{X₁+...+X_k ≤ Y} are considered. The random variables X’s and Y are assumed to follow binomial and Poisson distributions. Classical estimators available in the literature are discussed and Bayes estimators are derived. In order to obtain the estimators of these parametric functions, the basic role is played by the estimators of factorial moments of the two distributions.     

Doubly periodic non-homogeneous Poisson models for hurricane data

Non-homogeneous Poisson processes with periodic claim intensity rate have been proposed as claim counts in risk theory. Here a doubly periodic Poisson model with short- and long-term trends is studied. Beta-type intensity functions are presented as illustrations. The likelihood function and the maximum likelihood estimates of the model parameters are derived.Doubly periodic Poisson models are appropriate when the seasonality does not repeat exactly the same short-term pattern every year, but has a peak intensity that varies over a longer period. This reflects periodic environments like those forming hurricanes, in alternating El Niño/La Niña years. An application of the model to the data set of Atlantic hurricanes affecting the United States (1899–2000) is discussed in detail.     

On estimation of conditional density models with two-phase sampling

Suppose that the conditional density of a response variable given a vector of explanatory variables is parametrically modelled, and that data are collected by a two-phase sampling design. First, a simple random sample is drawn from the population. The stratum membership in a finite number of strata of the response and explanatory variables is recorded for each unit. Second, a subsample is drawn from the phase-one sample such that the selection probability is determined by the stratum membership. The response and explanatory variables are fully measured at this phase. We synthesize existing results on nonparametric likelihood estimation and present a streamlined approach for the computation and the large sample theory of profile likelihood in four different situations. The amount of information in terms of data and assumptions varies depending on whether the phase-one data are retained, the selection probabilities are known, and/or the stratum probabilities are known. We establish and illustrate numerically the order of efficiency among the maximum likelihood estimators, according to the amount of information utilized, in the four situations.     

Sample size choice for reliability verification in strength-stress models

This paper discusses the choice of sample size for experiments concerned with inference on R = P(Y < X), where X and Y are normal variates, in an acceptance-sampling-theory framework. A conservative approach is derived, and the properties of this solution examined by simulation.