On the dependent competing risks using Marshall–Olkin bivariate Weibull model: Parameter estimation with different methods

Competing risks models are of great importance in reliability and survival analysis. They are often assumed to have independent causes of failure in literature, which may be unreasonable. In this article, dependent causes of failure are considered by using the Marshall–Olkin bivariate Weibull distribution. After deriving some useful results for the model, we use ML, fiducial inference, and Bayesian methods to estimate the unknown model parameters with a parameter transformation. Simulation studies are carried out to assess the performances of the three methods. Compared with the maximum likelihood method, the fiducial and Bayesian methods could provide better parameter estimation.     

Maximum penalized likelihood estimation for skew-normal and skew-t distributions

The skew-normal and the skew-t distributions are parametric families which are currently under intense investigation since they provide a more flexible formulation compared to the classical normal and t distributions by introducing a parameter which regulates their skewness. While these families enjoy attractive formal properties from the probability viewpoint, a practical problem with their usage in applications is the possibility that the maximum likelihood estimate of the parameter which regulates skewness diverges. This situation has vanishing probability for increasing sample size, but for finite samples it occurs with non-negligible probability, and its occurrence has unpleasant effects on the inferential process. Methods for overcoming this problem have been put forward both in the classical and in the Bayesian formulation, but their applicability is restricted to simple situations. We formulate a proposal based on the idea of penalized likelihood, which has connections with some of the existing methods, but it applies more generally, including the multivariate case.     

Generalised data augmentation and posterior inferences

This paper explores the use of data augmentation in settings beyond the standard Bayesian one. In particular, we show that, after proposing an appropriate generalised data-augmentation principle, it is possible to extend the range of sampling situations in which fiducial methods can be applied by constructing Markov chains whose stationary distributions represent valid posterior inferences on model parameters. Some properties of these chains are presented and a number of open questions are discussed. We also use the approach to draw out connections between classical and Bayesian approaches in some standard settings.     

The fiducial argument and hacking's principle of irrelevance

In a searching analysis of the fiducial argument Hacking (1965) proposed the Principle of Irrelevance as a condition under which the argument is valid. His statement of the Principle was essentially non-mathematical and this paper presents a mathematical development of the Principle. The relationship with likelihood inference is explored and some of the proposed counter-examples to fiducial theory are considered. It is shown that even with the Principle of Irrelevance examples of non-uniqueness of fiducial distributions exist.     

Composite transformation models: A fiducial perspective

This paper discusses inferences for the parameters of a transformation model in the presence of a scalar nuisance parameter that describes the shape of the error distribution. The development is from the point of view of conditional inference and thus is an attempt to extend the classical fiducial (or structural inference) argument. For known shape parameter it is straightforward to derive a fiducial distribution of the transformation parameters from which confidence points can be obtained. For unknown shape parameter, the paper discusses a certain average of these fiducial distributions. The weights used in this averaging process are naturally induced by the action of the underlying group of transformations and correspond to a noninformative prior for the nuisance parameter. This results in a confidence distribution for the transformation parameters which in some cases has good frequentist properties. The method is illustrated by some examples.     

A fiducial approach to multiple comparisons

Comparing treatment means from populations that follow independent normal distributions is a common statistical problem. Many frequentist solutions exist to test for significant differences amongst the treatment means. A different approach would be to determine how likely it is that particular means are grouped as equal. We developed a fiducial framework for this situation. Our method provides fiducial probabilities that any number of means are equal based on the data and the assumed normal distributions. This methodology was developed when there is constant and non-constant variance across populations. Simulations suggest that our method selects the correct grouping of means at a relatively high rate for small sample sizes and asymptotic calculations demonstrate good properties. Additionally, we have demonstrated the flexibility in the methods ability to calculate the fiducial probability for any number of equal means. This was done by analyzing a simulated data set and a data set measuring the nitrogen levels of red clover plants that were inoculated with different treatments.     

Minimum distance lack-of-fit tests for fixed design

This paper discusses a class of tests of lack-of-fit of a parametric regression model when design is non-random and uniform on [0,1]. These tests are based on certain minimized distances between a nonparametric regression function estimator and the parametric model being fitted. We investigate asymptotic null distributions of the proposed tests, their consistency and asymptotic power against a large class of fixed and sequences of local nonparametric alternatives, respectively. The best fitted parameter estimate is seen to be n^1/2-consistent and asymptotically normal. A crucial result needed for proving these results is a central limit lemma for weighted degenerate U statistics where the weights are arrays of some non-random real numbers. This result is of an independent interest and an extension of a result of Hall for non-weighted degenerate U statistics.     

Inference for functions of parameters in discrete distributions based on fiducial approach: Binomial and Poisson cases

In this article, we propose a simple method of constructing confidence intervals for a function of binomial success probabilities and for a function of Poisson means. The method involves finding an approximate fiducial quantity (FQ) for the parameters of interest. A FQ for a function of several parameters can be obtained by substitution. For the binomial case, the fiducial approach is illustrated for constructing confidence intervals for the relative risk and the ratio of odds. Fiducial inferential procedures are also provided for estimating functions of several Poisson parameters. In particular, fiducial inferential approach is illustrated for interval estimating the ratio of two Poisson means and for a weighted sum of several Poisson means. Simple approximations to the distributions of the FQs are also given for some problems. The merits of the procedures are evaluated by comparing them with those of existing asymptotic methods with respect to coverage probabilities, and in some cases, expected widths. Comparison studies indicate that the fiducial confidence intervals are very satisfactory, and they are comparable or better than some available asymptotic methods. The fiducial method is easy to use and is applicable to find confidence intervals for many commonly used summary indices. Some examples are used to illustrate and compare the results of fiducial approach with those of other available asymptotic methods.     

A new generalized confidence interval for the among-group variance in the heteroscedastic one-way random effects model

Based on the generalized inference idea, a new kind of generalized confidence intervals is derived for the among-group variance component in the heteroscedastic one-way random effects model. We construct structure equations of all variance components in the model based on their minimal sufficient statistics; meanwhile, the fiducial generalized pivotal quantity (FGPQ) can be obtained through solving an implicit equation of the parameter of interest. Then, the confidence interval is derived naturally from the FGPQ. Simulation results demonstrate that the new procedure performs very well in terms of both empirical coverage probability and average interval length.     

Confidence Intervals of Performance Measures for an M/G/1 Queueing System

In this article, we discuss constructing confidence intervals (CIs) of performance measures for an M/G/1 queueing system. Fiducial empirical distribution is applied to estimate the service time distribution. We construct fiducial empirical quantities (FEQs) for the performance measures. The relationship between generalized pivotal quantity and fiducial empirical quantity is illustrated. We also present numerical examples to show that the FEQs can yield new CIs dominate the bootstrap CIs in relative coverage (defined as the ratio of coverage probability to average length of CI) for performance measures of an M/G/1 queueing system in most of the cases.     

Bias-corrected estimations in varying-coefficient partially nonlinear models with measurement error in the nonparametric part

In this paper, we consider the statistical inference for the varying-coefficient partially nonlinear model with additive measurement errors in the nonparametric part. The local bias-corrected profile nonlinear least-squares estimation procedure for parameter in nonlinear function and nonparametric function is proposed. Then, the asymptotic normality properties of the resulting estimators are established. With the empirical likelihood method, a local bias-corrected empirical log-likelihood ratio statistic for the unknown parameter, and a corrected and residual adjusted empirical log-likelihood ratio for the nonparametric component are constructed. It is shown that the resulting statistics are asymptotically chi-square distribution under some suitable conditions. Some simulations are conducted to evaluate the performance of the proposed methods. The results indicate that the empirical likelihood method is superior to the profile nonlinear least-squares method in terms of the confidence regions of parameter and point-wise confidence intervals of nonparametric function.     

Confidence intervals for variance components in unbalanced one-way random effects model using non-normal distributions

In scenarios where the variance of a response variable can be attributed to two sources of variation, a confidence interval for a ratio of variance components gives information about the relative importance of the two sources. For example, if measurements taken from different laboratories are nine times more variable than the measurements taken from within the laboratories, then 90% of the variance in the responses is due to the variability amongst the laboratories and 10% of the variance in the responses is due to the variability within the laboratories. Assuming normally distributed sources of variation, confidence intervals for variance components are readily available. In this paper, however, simulation studies are conducted to evaluate the performance of confidence intervals under non-normal distribution assumptions. Confidence intervals based on the pivotal quantity method, fiducial inference, and the large-sample properties of the restricted maximum likelihood (REML) estimator are considered. Simulation results and an empirical example suggest that the REML-based confidence interval is favored over the other two procedures in unbalanced one-way random effects model.     

Nonparametric estimation of the hazard function under dependence conditions

The estimation of the hazard rate has a great number of practical appli¬cations in dependence situations (seismicity analysis, reliability, economics), Based on kernel estimates of the density and the distribution function, we study the properties of the nonparametric estimator of the hazard function as-sociated with a strongly mixing time series. We prove consistency and asymp¬totic normality properties, and a cross-validation method for the smoothing parameter selection is studied. Some simulations and a practical application to real data are also shown.     

Fiducial distribution in a power series family

Abstract

In this article the interest is on finding the fiducial distribution of the parameter, when the probability distribution belongs to the power series family, as in Johnson et al. (1992 Johnson, N. L., S. Kotz, and A. W. Kemp. 1992. Univariate discrete distributions. New York, NY: John Wiley and Sons. [Google Scholar]). Recently in Nájera and O’Reilly (2017 Nájera, E., and F. O’Reilly. 2017. On fiducial generators. Communications Statistics - Theory Methods 46 (5):2232–2248. doi:10.1080/03610926.2015.1040505.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) an argument is given to obtain a unique fiducial in the Bernoulli case. An attempt is made here to define some sort of invariance in a power series distribution so that, as was done in the Bernoulli case, one may find a unique invariant fiducial for the parameter. The Bernoulli case is reviewed in detail and the Poisson and negative binomial cases are addressed.     

Limits of experiments associated with sampling plans

Starting from Milbrodt (1985), the asymptotic behaviour of experiments associated with Poisson sampling, Rejective sampling and its Sampford-Durbin modification is investigated. As superpopulation models so-called L^r-generated regression parameter families (1⩽r⩽2) are considered, allowing also the presence of nuisance parameters. Under some assumptions on the first order probabilities of inclusion it can be shown that the sampling experiments converge weakly if the underlying shift parameter families do so. In case of convergence the limit of the sampling experiments is characterized in terms of its Hellinger transforms and its Lévy-Khintchine representation, leading to criteria for the limit to be a pure Gaussian or a pure Poisson experiment respectively. These results are then applied to the situation of sampling in the presence of random non-response, and to establish local asymptotic normality (LAN) under more restrictive conditions. Applications also include asymptotic optimality properties of tests based on Horvitz-Thompson-type statistics, and LAM bounds and criteria for adaptivity, when testing or estimating a continuous linear functional in LAN situations. They especially cover the case of sampling from an unknown symmetric distribution, which has been subject to detailed investigations in the i.i.d. case.     

Fiducial Inferences for the Truncated Exponential Distribution

Lindqvist and Taraldsen (2005 Lindqvist , B. H. , Taraldsen , G. ( 2005 ). Monte Carlo conditioning on a sufficient statistic . Biometrika 92 : 451 – 464 .[Crossref], [Web of Science ®] , [Google Scholar]) introduced an interesting parametric family of distributions in the unit interval. In this note, inference procedures are given, both from the classical and the Bayesian view point. It is shown numerically through various examples that the posterior distribution for the parameter and the induced fiducial distribution are almost equivalent. The parametric family under study is a regular member of the Natural Exponential Family and so use of this fact permits induction of a unique fiducial in terms of the minimal sufficient statistic.     

Fiducial prediction intervals

This paper presents an approach for constructing prediction intervals for any given distribution. The approach is based on the principle of fiducial inference. We use several examples, including the normal, binomial, exponential, gamma, and Weibull distributions, to illustrate the proposed procedure.     

Nonparametric Phase-II monitoring for detecting monotone trend based on inverse sampling

Recently, Mukherjee and Bandyopadhyay (J Stat Plan Inference, 2011, doi:10.1016/j.jspi.2011.02.017) introduced some partially sequential tests for detecting liner trend among the incoming series of observations when a training sample is available a-priori. Their work is very useful in econometric or environmental monitoring under certain situations. The present work is intended for generalization of their tests for any monotone trend. We develop two nonparametric tests for the identity of some unknown univariate continuous distribution functions against monotone or unidirectional trend in location. One of these two tests is based on usual ranks and the other is based on sequential ranks. These are typical nonparametric tests for monitoring structural changes. Performance of the two tests are compared using asymptotic studies as well as through some numerical results based on Monte-Carlo simulations. An illustration is offered using a real data on monthly production of certain beverage.     

Nonparametric model checking for k-out-of-n systems

It is an important problem in reliability analysis to decide whether for a given k-out-of-n system the static or the sequential k-out-of-n model is appropriate. Often components are redundantly added to a system to protect against failure of the system. If the failure of any component of the system induces a higher rate of failure of the remaining components due to increased load, the sequential k-out-of-n model is appropriate. The increase of the failure rate of the remaining components after a failure of some component implies that the effects of the component redundancy are diminished. On the other hand, if all the components have the same failure distribution and whenever a failure occurs, the remaining components are not affected, the static k-out-of-n model is adequate. In this paper, we consider nonparametric hypothesis tests to make a decision between these two models. We analyze test statistics based on the profile score process as well as test statistics based on a multivariate intensity ratio and derive their asymptotic distribution. Finally, we compare the different test statistics.     

Prediction intervals for hypergeometric distributions

Abstract

The problem of constructing prediction intervals (PIs) for a future sample from a hypergeometric distribution is addressed. Simple closed-form approximate PIs based on the Wald approach, the joint sampling approach, and a fiducial approach are proposed and compared in terms of coverage probability and precision. Construction of the proposed PIs are illustrated using an example.