A semiparametric Bayesian approach to extreme value estimation

This paper is concerned with extreme value density estimation. The generalized Pareto distribution (GPD) beyond a given threshold is combined with a nonparametric estimation approach below the threshold. This semiparametric setup is shown to generalize a few existing approaches and enables density estimation over the complete sample space. Estimation is performed via the Bayesian paradigm, which helps identify model components. Estimation of all model parameters, including the threshold and higher quantiles, and prediction for future observations is provided. Simulation studies suggest a few useful guidelines to evaluate the relevance of the proposed procedures. They also provide empirical evidence about the improvement of the proposed methodology over existing approaches. Models are then applied to environmental data sets. The paper is concluded with a few directions for future work.     

Semiparametric bivariate modelling with flexible extremal dependence

Statistics and Computing - Inference over multivariate tails often requires a number of assumptions which may affect the assessment of the extreme dependence structure. Models are usually...     

Bayesian dynamic models for survival data with a cure fraction

In this paper, we propose a new class of semi-parametric cure rate models. Specifically, we construct dynamic models for piecewise hazard functions over a finite partition of the time axis. Allowing the size of partition and the levels of baseline hazard to be random, our proposed models provide a great flexibility in controlling the degree of parametricity in the right tail of the survival distribution and the amount of correlations among the log-baseline hazard levels. Several properties of the proposed models are derived, and propriety of the implied posteriors with improper noninformative priors for regression coefficients based on the proposed models is established for the fixed partition of the time axis. In addition, an efficient reversible jump computational algorithm is developed for carrying out posterior computation. A real data set from a melanoma clinical trial is analyzed in detail to further demonstrate the proposed methodology.     

An adaptive resampling scheme for cycle estimation

Bayesian dynamic linear models (DLMs) are useful in time series modelling, because of the flexibility that they off er for obtaining a good forecast. They are based on a decomposition of the relevant factors which explain the behaviour of the series through a series of state parameters. Nevertheless, the DLM as developed by West and Harrison depend on additional quantities, such as the variance of the system disturbances, which, in practice, are unknown. These are referred to here as 'hyper-parameters' of the model. In this paper, DLMs with autoregressive components are used to describe time series that show cyclic behaviour. The marginal posterior distribution for state parameters can be obtained by weighting the conditional distribution of state parameters by the marginal distribution of hyper-parameters. In most cases, the joint distribution of the hyperparameters can be obtained analytically but the marginal distributions of the components cannot, so requiring numerical integration. We propose to obtain samples of the hyperparameters by a variant of the sampling importance resampling method. A few applications are shown with simulated and real data sets.     

Dynamic survival models with spatial frailty

In many survival studies, covariates effects are time-varying and there is presence of spatial effects. Dynamic models can be used to cope with the variations of the effects and spatial components are introduced to handle spatial variation. This paper proposes a methodology to simultaneously introduce these components into the model. A number of specifications for the spatial components are considered. Estimation is performed via a Bayesian approach through Markov chain Monte Carlo methods. Models are compared to assess relevance of their components. Analysis of a real data set is performed, showing the relevance of both time-varying covariate effects and spatial components. Extensions to the methodology are proposed along with concluding remarks.     

Sampling from the posterior distribution in generalized linear mixed models

Generalized linear mixed models provide a unified framework for treatment of exponential family regression models, overdispersed data and longitudinal studies. These problems typically involve the presence of random effects and this paper presents a new methodology for making Bayesian inference about them. The approach is simulation-based and involves the use of Markov chain Monte Carlo techniques. The usual iterative weighted least squares algorithm is extended to include a sampling step based on the Metropolis–Hastings algorithm thus providing a unified iterative scheme. Non-normal prior distributions for the regression coefficients and for the random effects distribution are considered. Random effect structures with nesting required by longitudinal studies are also considered. Particular interests concern the significance of regression coefficients and assessment of the form of the random effects. Extensions to unknown scale parameters, unknown link functions, survival and frailty models are outlined.