Objective Bayesian Inference for the Half-Normal and Half-t Distributions

In this article, Bayesian inference for the half-normal and half-t distributions using uninformative priors is considered. It is shown that exact Bayesian inference can be undertaken for the half-normal distribution without the need for Gibbs sampling. Simulation is then used to compare the sampling properties of Bayesian point and interval estimators with those of their maximum likelihood based counterparts. Inference for the half-t distribution based on the use of Gibbs sampling is outlined, and an approach to model comparison based on the use of Bayes factors is discussed. The fitting of the half-normal and half-t models is illustrated using real data on the body fat measurements of elite athletes.     

A Generalized Poisson Distribution

Cumulative probabilities of a Poisson distribution can be written in terms of incomplete gamma function where the parameter of the gamma function is an integer. From this definition a new generalization of the Poisson distribution is obtained with two parameters. Asymptotic behavior of this distribution is shown to be normal. Some order properties of this distribution are also studied.     

Matrix variate Macdonald distribution

Abstract

In this article, we generalize the univariate Macdonald distribution to the matrix case and give its derivation using matrix variate gamma distribution. We study several properties such as cumulative distribution function, marginal distribution of submatrix, triangular factorization, moment generating function, and expected values of several functions of the Macdonald matrix. Some of these results are expressed in terms of special functions of matrix arguments and zonal polynomials.     

Empirical Bayes models of Poisson clinical trials and sample size determination

Bayesian methods are often used to reduce the sample sizes and/or increase the power of clinical trials. The right choice of the prior distribution is a critical step in Bayesian modeling. If the prior not completely specified, historical data may be used to estimate it. In the empirical Bayesian analysis, the resulting prior can be used to produce the posterior distribution. In this paper, we describe a Bayesian Poisson model with a conjugate Gamma prior. The parameters of Gamma distribution are estimated in the empirical Bayesian framework under two estimation schemes. The straightforward numerical search for the maximum likelihood (ML) solution using the marginal negative binomial distribution is unfeasible occasionally. We propose a simplification to the maximization procedure. The Markov Chain Monte Carlo method is used to create a set of Poisson parameters from the historical count data. These Poisson parameters are used to uniquely define the Gamma likelihood function. Easily computable approximation formulae may be used to find the ML estimations for the parameters of gamma distribution. For the sample size calculations, the ML solution is replaced by its upper confidence limit to reflect an incomplete exchangeability of historical trials as opposed to current studies. The exchangeability is measured by the confidence interval for the historical rate of the events. With this prior, the formula for the sample size calculation is completely defined. Published in 2009 by John Wiley & Sons, Ltd.     

运用生存分析与变点理论对上证指数的研究

本文运用生存分析与变点理论对上证指数进行了研究,发现连涨和连跌的股指服从伽玛分布,表明股市是有记忆的并且不服从随机游动假说。当股市发生变化时,伽玛分布的参数也会发生变化,而且伽玛分布的两个参数很好地反映了牛市、熊市的变化以及股市波动的变化。我们运用变点理论对这一变化进行了实证分析,得出了一些有意义的结论。     

On k-th record times,record values and their moments

We give new formula for moments of k-th record values in terms of Stirling numbers of the first kind. In particular, the formulae allow to derive the explicit formulae for moments of k-th lower record values from exponential distribution which have not been known yet. Moreover, some interesting identities involving harmonic numbers are also obtained as corollaries to presented results.     

The Additive Genetic Gamma Frailty Model

In this paper the additive genetic gamma frailty model is defined. Individual frailties are correlated as a result of an additive genetic model. An algorithm to construct additive genetic gamma frailties for any pedigree is given so that the variance–covariance structure among individual frailties equals the numerator relationship matrix times a variance. The EM algorithm can be used to estimate the parameters in the model. Calculations are similar using the EM algorithm in the shared frailty model, however the E step is not correspondingly simple. This is illustrated re-analysing data, analysed by the shared frailty model in Nielsen et al . (1992), from the Danish adoptive register. Goodness of fit of the additive genetic gamma frailty model can be tested after analysing data with the correlated frailty model. Doing so, a "defect" in the often used and otherwise well behaving likelihood was found     

On improved estimation of a gamma shape parameter

This paper deals with improved estimation of a gamma shape parameter from a decision-theoretic point of view. First we study the second-order properties of three estimators – (i) the maximum-likelihood estimator (MLE), (ii) a bias corrected version of the MLE, and (iii) an improved version (in terms of mean squared error) of the MLE. It is shown that all the three estimators mentioned above are second-order inadmissible. Next, we obtain superior estimators which are second order better than the above three estimators. Simulation results are provided to study the relative risk improvement of each improved estimator over the MLE.     

Modified likelihood ratio tests for unit gamma regressions

Regression analyses are commonly performed with doubly limited continuous dependent variables; for instance, when modeling the behavior of rates, proportions and income concentration indices. Several models are available in the literature for use with such variables, one of them being the unit gamma regression model. In all such models, parameter estimation is typically performed using the maximum likelihood method and testing inferences on the model''s parameters are usually based on the likelihood ratio test. Such a test can, however, deliver quite imprecise inferences when the sample size is small. In this paper, we propose two modified likelihood ratio test statistics for use with the unit gamma regressions that deliver much more accurate inferences when the number of data points in small. Numerical (i.e. simulation) evidence is presented for both fixed dispersion and varying dispersion models, and also for tests that involve nonnested models. We also present and discuss two empirical applications.     

Some classes of distributions and their application in queueing

New classes of distributions are defined in analogy to the properties NBU, NBUE known from reliability. They are applied to obtain bounds on certain parameters of the GI/G/1 queue, such as the mean and the variance of the stationary waiting time, the probability of waiting, and the covariances of waiting times.