Bayesian estimation and prediction based on generalized Type-I hybrid censored sample

In this paper, we consider an exponential form for the underlying distributionand a conjugate prior, and develop a procedure for deriving the maximum likelihood and Bayesian estimators based on an observed generalized Type-I hybrid censored sample. The problems of predicting the future order statistics from the same sample and that from a future sample are also discussed from a Bayesian viewpoint. For the illustration of the developed results, the exponential and Pareto distributions are used as examples. Finally, two numerical examples are presented for illustrating all the inferential procedures developed here.     

Bayesian prediction of order statistics with fixed and random sample sizes based on k-record values from Pareto distribution

In this paper, the two-parameter Pareto distribution is considered and the problem of prediction of order statistics from a future sample and that of its geometric mean are discussed. The Bayesian approach is applied to construct predictors based on observed k-record values for the cases when the future sample size is fixed and when it is random. Several Bayesian prediction intervals are derived. Finally, the results of a simulation study and a numerical example are presented for illustrating all the inferential procedures developed here.     

Bayesian prediction of future number of failures based on finite mixture of general class of distributions

This paper is concerned with the Bayesian prediction problem of the number of components which will fail in a future time interval. The failure times are distributed according to a finite mixture of a general class of distributions. Type-I censored sample from this nonhomogeneous population and a general class of prior density functions are used. A one-sample scheme is used to predict the number of failures in a future time interval. An example of a finite mixture of k exponential components is given to illustrate our results.     

Robust Bayesian analysis for parallel system with masked data under inverse Weibull lifetime distribution

Abstract

This paper considers the statistical analysis of masked data in a parallel system with inverse Weibull distributed components under type II censoring. Based on Gamma conjugate prior, the Bayesian estimation as well as the hierarchical Bayesian estimation for the parameters and the reliability function of system are obtained by using the Bayesian theory and the hierarchical Bayesian method. Finally, Monte Carlo simulations are provided to compare the performances of the estimates under different masking probabilities and effective sample sizes.     

Prediction bounds for the burr model

Given a type 2 censored sample from the Burr life time distribution, Bayesian prediction bounds are derived for future observations. An approximate Bayesian method has been used to simplify the computation of the prediction bounds. Numerical examples are used to illustrate the procedures.     

Bayesian inference based on a jointly type-II censored sample from two exponential populations

In this paper, based on a jointly type-II censored sample from two exponential populations, the Bayesian inference for the two unknown parameters are developed with the use of squared-error, linear-exponential and general entropy loss functions. The problem of predicting the future failure times, both point and interval prediction, based on the observed joint type-II censored data, is also addressed from a Bayesian viewpoint. A Monte Carlo simulation study is conducted to compare the Bayesian estimators with the maximum likelihood estimator developed by Balakrishnan and Rasouli [Exact likelihood inference for two exponential populations under joint type-II censoring. Comput Stat Data Anal. 2008;52:2725–2738]. Finally, a numerical example is utilized for the purpose of illustration.     

A unified approach to marginal equivalence in the general framework of group invariance

ABSTRACT

Two Bayesian models with different sampling densities are said to be marginally equivalent if the joint distribution of observables and the parameter of interest is the same for both models. We discuss marginal equivalence in the general framework of group invariance. We introduce a class of sampling models and derive marginal equivalence when the prior for the nuisance parameter is relatively invariant. We also obtain some robustness properties of invariant statistics under our sampling models. Besides the prototypical example of v-spherical distributions, we apply our general results to two examples—analysis of affine shapes and principal component analysis.     

Bayesian Causality

Abstract

Although no universally accepted definition of causality exists, in practice one is often faced with the question of statistically assessing causal relationships in different settings. We present a uniform general approach to causality problems derived from the axiomatic foundations of the Bayesian statistical framework. In this approach, causality statements are viewed as hypotheses, or models, about the world and the fundamental object to be computed is the posterior distribution of the causal hypotheses, given the data and the background knowledge. Computation of the posterior, illustrated here in simple examples, may involve complex probabilistic modeling but this is no different than in any other Bayesian modeling situation. The main advantage of the approach is its connection to the axiomatic foundations of the Bayesian framework, and the general uniformity with which it can be applied to a variety of causality settings, ranging from specific to general cases, or from causes of effects to effects of causes.     

Laplace approximations and Bayesian information criteria in possibly misspecified models

We provide general conditions to ensure the valid Laplace approximations to the marginal likelihoods under model misspecification, and derive the Bayesian information criteria including all terms of order O_p(1). Under conditions in theorem 1 of Lv and Liu [J. R. Statist. Soc. B, 76, (2014), 141–167] and a continuity condition for prior densities, asymptotic expansions with error terms of order o_p(1) are derived for the log-marginal likelihoods of possibly misspecified generalized linear models. We present some numerical examples to illustrate the finite sample performance of the proposed information criteria in misspecified models.     

Simple and flexible Bayesian inferences for standardized regression coefficients

ABSTRACT

In statistical practice, inferences on standardized regression coefficients are often required, but complicated by the fact that they are nonlinear functions of the parameters, and thus standard textbook results are simply wrong. Within the frequentist domain, asymptotic delta methods can be used to construct confidence intervals of the standardized coefficients with proper coverage probabilities. Alternatively, Bayesian methods solve similar and other inferential problems by simulating data from the posterior distribution of the coefficients. In this paper, we present Bayesian procedures that provide comprehensive solutions for inferences on the standardized coefficients. Simple computing algorithms are developed to generate posterior samples with no autocorrelation and based on both noninformative improper and informative proper prior distributions. Simulation studies show that Bayesian credible intervals constructed by our approaches have comparable and even better statistical properties than their frequentist counterparts, particularly in the presence of collinearity. In addition, our approaches solve some meaningful inferential problems that are difficult if not impossible from the frequentist standpoint, including identifying joint rankings of multiple standardized coefficients and making optimal decisions concerning their sizes and comparisons. We illustrate applications of our approaches through examples and make sample R functions available for implementing our proposed methods.     

Bayesian prediction with a random sample size for the burr lifetime distribution

This paper is concerned with the problem of deriving Bayesian prediction bounds for the Burr distribution when the sample size is a random variable. Prediction bounds for both the future observations (the case of two-sample prediction) and the remaining observations in the same sample (the case of one-sample prediction) will be derived. The analysis will depend mainly on assuming that the size of the sample is a random variable having the Poisson distribution. Finally, numerical examples are given to illustrate the results.     

Bayesian prediction of the total time on test using doubly censored rayleigh data

In a Bayesian setting, and on the basis of a doubly censored random sample of failure times drawn from a Rayleigh distribution, Fernandez (2000, Statist. Probab. Lett. , 48 , 393-399) considered the problem of predicting an independent future sample from the same distribution. In this article, we extend his work to include the estimation of the predictive distribution of the total time on test up to a certain failure in a future sample, as well as that of the remaining testing time time until all the items in the original sample have failed. Two examples are used to illustrate the prediction procedure.     

A flexible procedure for formulating probability distributions on the unit interval with applications

Abstract

In this paper, we present a flexible mechanism for constructing probability distributions on a bounded intervals which is based on the composition of the baseline cumulative probability function and the quantile transformation from another cumulative probability distribution. In particular, we are interested in the (0, 1) intervals. The composite quantile family of probability distributions contains many models that have been proposed in the recent literature and new probability distributions are introduced on the unit interval. The proposed methodology is illustrated with two examples to analyze a poverty dataset in Peru from the Bayesian paradigm and Likelihood points of view.     

Blending Bayesian and Classical Tools to Define Optimal Sample-Size-Dependent Significance Levels

ABSTRACT

This article argues that researchers do not need to completely abandon the p-value, the best-known significance index, but should instead stop using significance levels that do not depend on sample sizes. A testing procedure is developed using a mixture of frequentist and Bayesian tools, with a significance level that is a function of sample size, obtained from a generalized form of the Neyman–Pearson Lemma that minimizes a linear combination of α, the probability of rejecting a true null hypothesis, and β, the probability of failing to reject a false null, instead of fixing α and minimizing β. The resulting hypothesis tests do not violate the Likelihood Principle and do not require any constraints on the dimensionalities of the sample space and parameter space. The procedure includes an ordering of the entire sample space and uses predictive probability (density) functions, allowing for testing of both simple and compound hypotheses. Accessible examples are presented to highlight specific characteristics of the new tests.     

Bayesian two-stage design for drug screening trials with switching hypothesis tests based on continuous endpoints

Abstract

In this paper, we propose a Bayesian two-stage design with changing hypothesis test by bridging a single-arm study and a double-arm randomized trial in one phase II clinical trial based on continuous endpoints rather than binary endpoints. We have also calibrated with respect to frequentist and Bayesian error rates. The proposed design minimizes the Bayesian expected sample size if the new candidate has low or high efficacy activity subject to the constraint upon error rates in both frequentist and Bayesian perspectives. Tables of designs for various combinations of design parameters are also provided.     

Equal-tailed and shortest Bayesian tolerance intervals based on exponential k-records

The problem of constructing the equal-tailed and shortest Bayesian tolerance intervals that control percentages in both tails of the exponential distribution based on k-record values is considered. Equal-tailed and shortest Bayesian tolerance factors are derived. Practical examples using real and simulated k-record values are given to illustrate the proposed results.     

Exact Likelihood Inference for k Exponential Populations Under Joint Type-II Censoring

In this paper, when a jointly Type-II censored sample arising from k independent exponential populations is available, the conditional MLEs of the k exponential mean parameters are derived. The moment generating functions and the exact densities of these MLEs are obtained using which exact confidence intervals are developed for the parameters. Moreover, approximate confidence intervals based on the asymptotic normality of the MLEs and credible confidence regions from a Bayesian viewpoint are also discussed. An empirical comparison of the exact, approximate, bootstrap, and Bayesian intervals is also made in terms of coverage probabilities. Finally, an example is presented in order to illustrate all the methods of inference developed here.     

Bayesian Elastic Net Tobit Quantile Regression

In this article, the problem of parameter estimation and variable selection in the Tobit quantile regression model is considered. A Tobit quantile regression with the elastic net penalty from a Bayesian perspective is proposed. Independent gamma priors are put on the l₁ norm penalty parameters. A novel aspect of the Bayesian elastic net Tobit quantile regression is to treat the hyperparameters of the gamma priors as unknowns and let the data estimate them along with other parameters. A Bayesian Tobit quantile regression with the adaptive elastic net penalty is also proposed. The Gibbs sampling computational technique is adapted to simulate the parameters from the posterior distributions. The proposed methods are demonstrated by both simulated and real data examples.     

Bayesian estimation for non zero inflated modified power series distribution under linex and generalized entropy loss functions

Abstract

In this paper, we derive Bayesian estimators of the parameters of modified power series distributions inflated at any of a support point under linex and general entropy loss function. We assume that the prior information can be summarized by a uniform, Beta, two-sided power, Gamma or generalized Pareto distributions. The obtained results are demonstrated on the generalized Poisson and the generalized negative binomial distribution inflated at a given point.