Bayesian Estimation for Poisson-exponential Model under Progressive Type-II Censoring Data with Binomial Removal and Its Application to Ovarian Cancer Data

In this article, we propose Maximum likelihood estimators (MLEs) and Bayes estimators of parameters of Poisson-exponential distribution (PED) under General entropy loss function (GELF) and Squared error loss function (SELF) for Progressive type-II censored data with binomial removals (PT-II CBRs). The MLEs and corresponding Bayes estimators are compared in terms of their risks based on simulated samples from PED. The proposed methodology is illustrated on a real dataset of ovarian cancer.     

An accurate,tractable, and analytically integrable polynomial expansion of the skewed Student’s t-distribution

An integrable polynomial expansion is derived for the skewed Student’s t-distribution. Analytic integration of the density function results in a polynomial expression for the cumulative distribution function too, which is easily manageable from a numerical and mathematical point of view and enables computing quantiles with minimum numerical effort.     

Non-parametric predictive inference for future order statistics

This article presents non-parametric predictive inference for future order statistics. Given the data consisting of n real-valued observations, m future observations are considered and predictive probabilities are presented for the rth-ordered future observation. In addition, joint and conditional probabilities for events involving multiple future order statistics are presented. The article further presents the use of such predictive probabilities for order statistics in statistical inference, in particular considering pairwise and multiple comparisons based on two or more independent groups of data.     

Making Low Probabilities Useful

This paper explores how people process information on low probability-high consequence negative events and what it will take to get individuals to be sensitive to the likelihood of these types of accidents or disasters. In a set of experiments, information is presented to individuals on the likelihood of serious accidents from a chemical facility. Comparisons are made with other risks, such as fatalities from automobile accidents, to see whether laypersons can determine the relative safety of different plants. We conclude that fairly rich context information must be available for people to be able to judge differences between low probabilities. In particular, it appears that one needs to present comparison scenarios that are located on the probability scale to evoke people's own feelings of risk. The concept of evaluability recently introduced by Hsee and his colleagues provides a useful explanation of these findings.     

Global and local tests to assess stationarity of Markov transition models

We present global and local likelihood-based tests to evaluate stationarity in transition models. Three motivational studies are considered. A simulation study was carried out to assess the performance of the proposed tests. The results showed that they present good performance with the control of the type-I error, especially for ordinal responses, and control of the type-II error, especially for the nominal case. Asymptotically they are close to the classical test performance. They can be executed in a single framework without the need to estimate the transition probabilities, incorporating both categorical and continuous covariates, and used to identify sources of non-stationarity.     

An alternative method for evaluating stationarity in transition models

Transition models are an important framework that can be used to model longitudinal categorical data. A relevant issue in applying these models is the condition of stationarity, or homogeneity of transition probabilities over time. We propose two tests to assess stationarity in transition models: Wald and likelihood-ratio tests, which do not make use of transition probabilities, using only the estimated parameters of the models in contrast to the classical test available in the literature. In this paper, we present two motivating studies, with ordinal longitudinal data, to which proportional odds transition models are fitted and the two proposed tests are applied as well as the classical test. Additionally, their performances are assessed through simulation studies. The results show that the proposed tests have good performance, being better for control of type-I error and they present equivalent power functions asymptotically. Also, the correlations between the Wald, likelihood-ratio and the classical test statistics are positive and large, an indicator of general concordance. Additionally, both of the proposed tests are more flexible and can be applied in studies with qualitative and quantitative covariates.     

Goodness of fit of generalized linear models to sparse data

We derive approximations to the first three moments of the conditional distribution of the deviance statistic, for testing the goodness of fit of generalized linear models with non-canonical links, by using an estimating equations approach, for data that are extensive but sparse. A supplementary estimating equation is proposed from which the modified deviance statistic is obtained. An application of a modified deviance statistic is shown to binomial and Poisson data. We also conduct a performance study of the modified Pearson statistic derived by Farrington and the modified deviance statistic derived in this paper, in terms of size and power, through a small scale simulation experiment. Both statistics are shown to perform well in terms of size. The deviance statistic, however, shows an advantage of power. Two examples are given.     

Elicitation of Subjective Probabilities when the Initial Endowment is Unobservable

In the framework of subjective expected utility theory we develop two distinct procedures for the elicitation of a person's subjective probabilities when the initial endowment is random and unobservable. Procedures of the first kind rely on the boundedness of the utility function to elicit the person's subjective probabilities directly. Procedures of the second kind infer the person's initial endowment from his choice behavior prior to the application of probability elicitation.     

The Effect of Non-Stationarity on Extreme Sea-Level Estimation

The sea-level is the composition of astronomical tidal and meteorological surge processes. It exhibits temporal non-stationarity due to a combination of long-term trend in the mean level, the deterministic tidal component, surge seasonality and interactions between the tide and surge. We assess the effect of these non-stationarities on the estimation of the distribution of extreme sea-levels. This is important for coastal flood assessment as the traditional method of analysis assumes that, once the trend has been removed, extreme sea-levels are from a stationary sequence. We compare the traditional approach with a recently proposed alternative that incorporates the knowledge of the tidal component and its associated interactions, by applying them to 22 UK data sites and through a simulation study. Our main finding is that if the tidal non-stationarity is ignored then a substantial underestimation of extreme sea-levels results for most sites. In contrast, if surge seasonality and the tide–surge interaction are not modelled the traditional approach produces little additional bias. The alternative method is found to perform well but requires substantially more statistical modelling and better data quality.     

On Bayesian analyses and finite mixtures for proportions

When the results of biological experiments are tested for a possible difference between treatment and control groups, the inference is only valid if based upon a model that fits the experimental results satisfactorily. In dominant-lethal testing, foetal death has previously been assumed to follow a variety of models, including a Poisson, Binomial, Beta-binomial and various mixture models. However, discriminating between models has always been a particularly difficult problem. In this paper, we consider the data from 6 separate dominant-lethal assay experiments and discriminate between the competing models which could be used to describe them. We adopt a Bayesian approach and illustrate how a variety of different models may be considered, using Markov chain Monte Carlo (MCMC) simulation techniques and comparing the results with the corresponding maximum likelihood analyses. We present an auxiliary variable method for determining the probability that any particular data cell is assigned to a given component in a mixture and we illustrate the value of this approach. Finally, we show how the Bayesian approach provides a natural and unique perspective on the model selection problem via reversible jump MCMC and illustrate how probabilities associated with each of the different models may be calculated for each data set. In terms of estimation we show how, by averaging over the different models, we obtain reliable and robust inference for any statistic of interest.