Comparison of Two Methods for Calculating the Partition Functions of Various Spatial Statistical Models

Likelihood computation in spatial statistics requires accurate and efficient calculation of the normalizing constant (i.e. partition function) of the Gibbs distribution of the model. Two available methods to calculate the normalizing constant by Markov chain Monte Carlo methods are compared by simulation experiments for an Ising model, a Gaussian Markov field model and a pairwise interaction point field model.     

Algorithmic Inference of Two-Parameter Gamma Distribution

We provide an estimation procedure of the two-parameter Gamma distribution based on the Algorithmic Inference approach. As a key feature of this approach, we compute the joint probability distribution of these parameters without assuming any prior. To this end, we propose a numerical algorithm which is often beneficial of a highly efficient speed up based on an approximate analytical expression of the probability distribution. We contrast our interval and point estimates with those recently obtained in Son and Oh (2006 Son , Y. S. , Oh , M. ( 2006 ). Bayesian estimation of the two-parameter gamma distribution . Communications in Statistics – Simulation and Computation 35 : 285 – 293 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) for the same problem. From this benchmark we realize that our estimates are both unbiased and more accurate, albeit more dispersed, in some cases, than the competitors' methods, where the dispersion drawback is notably mitigated w.r.t. Bayesian methods by a greater estimate decorrelation. We also briefly discuss the theoretical novelty of the adopted inference paradigm which actually represents a brush up on a Fisher perspective dating to almost a century, made feasible today by the available computational tools.     

Nonlinear mixed‐effect models with nonignorably missing covariates

Nonlinear mixed‐effect models are often used in the analysis of longitudinal data. However, it sometimes happens that missing values for some of the model covariates are not purely random. Motivated by an application to HTV viral dynamics, where this situation occurs, the author considers likelihood inference for this type of problem. His approach involves a Monte Carlo EM algorithm, along with a Gibbs sampler and rejection/importance sampling methods. A concrete application is provided.     

Bayesian Estimation of the Two-Parameter Gamma Distribution

In this article, we present and discuss an original price index being a special case of a general formula for price indices. We show that the discussed formula satisfies most postulates coming from the axiomatic price index theory. We compare this index to other known and popular price indices in a simulation study.     

Binary probit,tobit and hazard rate

This short note drows attention to an obvious although so far only partly noted relation between an inequality involving the hazard rate of the normal distribution and maximum likelihood estimation in the binary probit and in the Tobit model. Global concavity of the likelihood function for both the binary probit and the Tobit model can be proved by means of an inequality concerning the hazard rate of the standard normal distribution. As not yet noted in the literature this inequality may also be used to show that the hazard rate is monotonically increasing. Results of this paper are related to research which is financially supported by Thyssen-Stiftung. Comments on an earlier version by an anonymous referee are gratefully acknowledged.     

Likelihood-Based Inference for the Shape Parameter of a Two Parameter Weibull Distribution

Based on a Type 2 censored sample, we use the likelihood-based approach to draw likelihood inference on the shape parameter gamma of a two-parameter Weibull distribution. In particular, we derive the profile, conditional and marginal likelihoods of gamma. Numerical results along with some concluding remarks regarding the use of likelihood-based methods for inference are provided.     

Information on sexual behaviour when some data are missing

When tables contain missing values, statistical models that allow the non-response probability to be a function of the intended response have been proposed by several researchers. We investigate the properties of these methods in the context of a survey of the sexual behaviour of university students. Profile likelihoods can be computed, even when models are not identified and saturated profile likelihoods (making no assumptions about the non-response mechanism) are derived. Bayesian approaches are investigated and it is shown that their results may be highly sensitive to the prior specification. The proportion of responders answering 'yes' to the question 'have you ever had sexual intercourse?' was 73%. However, different assumptions about the non-responders gave proportions as low as 46% or as high as 83%. Our preferred estimate, derived from the response-saturated profile likelihood, is 67% with a 95% confidence interval of 58–74%. This is in line with other studies on response bias in the reports of young people's sexual behaviour that suggest that the respondents overrepresent the sexually active.     

Estimation of conditional mode with truncated,censored, and dependent data

In this paper, we studied the uniform convergence with rates for the kernel estimator of the conditional mode function for a left truncated and right censored model. It is assumed that the lifetime observations with multivariate covariates form a stationary α-mixing sequence. Also, the asymptotic normality of the estimator is established.     

A Bayesian analysis for the Block and Basu bivariate exponential distribution in the presence of covariates and censored data

In this paper, we introduce a Bayesian Analysis for the Block and Basu bivariate exponential distribution using Markov Chain Monte Carlo (MCMC) methods and considering lifetimes in presence of covariates and censored data. Posterior summaries of interest are obtained using the popular WinBUGS software. Numerical illustrations are introduced considering a medical data set related to the recurrence times of infection for kidney patients and a medical data set related to bone marrow transplantation for leukemia.     

Estimating fault hitting rates by recapture sampling

For the recapture debugging design introduced by Nayak (1988) we consider the problem of estimating the hitting rates of the faults remaining in a system. In the context of a conditional likelihood, moment estimators are derived and are shown to be asymptotically normal and fully efficient. Fixed sample properties of the moment estimators are compared, through simulation, with those of the conditional maximum likelihood estimators. Also considered is a procedure for testing the assumption that faults have identical hitting rates; this provides a test of fit of the Jelinski-Moranda (1972) model. It is assumed that the residual hitting rates follow a log linear rate model and that the testing process is truncated when the gaps between the detection of new errors exceed a fixed amount of time.     

Asymptotic normality of conditional density estimation under truncated,censored and dependent data

Abstract

In this paper, we focus on the left-truncated and right-censored model, and construct the local linear and Nadaraya-Watson type estimators of the conditional density. Under suitable conditions, we establish the asymptotic normality of the proposed estimators when the observations are assumed to be a stationary α-mixing sequence. Finite sample behavior of the estimators is investigated via simulations too.     

The Two-interval Line-segment Problem

In this paper we define and study the non-parametric maximum likelihood estimator (NPMLE) in the one-dimensional line-segment problem, where we observe line-segments on the real line through an interval with a gap which is smaller than the two remaining intervals. We define the self-consistency equations for the NPMLE and provide a quick algorithm for solving them. We prove supremum norm weak convergence to a Gaussian process and efficiency of the NPMLE. The problem has a geological application in the study of the lifespan of species     

Markov chain Monte Carlo methods for parameter estimation of the modified Weibull distribution

In this paper, the Markov chain Monte Carlo (MCMC) method is used to estimate the parameters of a modified Weibull distribution based on a complete sample. While maximum-likelihood estimation (MLE) is the most used method for parameter estimation, MCMC has recently emerged as a good alternative. When applied to parameter estimation, MCMC methods have been shown to be easy to implement computationally, the estimates always exist and are statistically consistent, and their probability intervals are convenient to construct. Details of applying MCMC to parameter estimation for the modified Weibull model are elaborated and a numerical example is presented to illustrate the methods of inference discussed in this paper. To compare MCMC with MLE, a simulation study is provided, and the differences between the estimates obtained by the two algorithms are examined.     

Bayesian prediction of unobserved values for Type-II censored data

In this paper, we consider posterior predictive distributions of Type-II censored data for an inverse Weibull distribution. These functions are given by using conditional density functions and conditional survival functions. Although the conditional survival functions were expressed by integral forms in previous studies, we derive the conditional survival functions in closed forms and thereby reduce the computation cost. In addition, we calculate the predictive confidence intervals of unobserved values and coverage probabilities of unobserved values by using the posterior predictive survival functions.     

Random variate transformations in the Gibbs sampler: issues of efficiency and convergence

In the non-conjugate Gibbs sampler, the required sampling from the full conditional densities needs the adoption of black-box sampling methods. Recent suggestions include rejection sampling, adaptive rejection sampling, generalized ratio of uniforms, and the Griddy-Gibbs sampler. This paper describes a general idea based on variate transformations which can be tailored in all the above methods and increase the Gibbs sampler efficiency. Moreover, a simple technique to assess convergence is suggested and illustrative examples are presented.     

Computations via Auxiliary Random Functions for Survival Models

Abstract.  A new simulation method, auxiliary random functions is introduced. When used within a Gibbs sampler, this method enables a unified treatment of exact, right-censored, left-censored, left-truncated and interval censored data, with and without covariates in survival models. The models and methods are exemplified via illustrative analysis.     

Estimation from Censored Data with Incomplete Information

Phillips and Sweeting [J. R. Statist. Soc. B 58 (1996) 775–783.] considered estimation of the parameter of the exponential distribution with censored failure time data when there is incomplete knowledge of the censoring times. It was shown that, under particular models for the censoring mechanism and censoring errors, it will usually be safe to ignore such errors provided they are not expected to be too large. A flexible model is introduced which includes the extreme cases of no censoring errors and no information on the censoring values. The effect of alternative assumptions about knowledge of the censoring values on the estimation of failure rate is investigated.