基于倾斜-t Copula的半参数马尔可夫链

在不指定时间序列结构的情况下,我们的分布模型是基于多变量离散时间的相应马尔可夫族和相关变量一维的边际分布.这样的模型可以同时处理时间序列之间的相互依赖和每个时间序列沿时间方向的依赖.具体的参数copula被指定为倾斜-t. 倾斜-t Copla能够处理不对称,偏斜和粗尾的数据分布.三个股票指数日均收益的实证研究表明,倾斜-t copula的马尔可夫模型要比以下模型更好:倾斜正态Copula马可夫, t-copula马可夫, 倾斜-t copula但无马尔可夫特性.     

Classification of brain activation via spatial Bayesian variable selection in fMRI regression

Functional magnetic resonance imaging (fMRI) is the most popular technique in human brain mapping, with statistical parametric mapping (SPM) as a classical benchmark tool for detecting brain activity. Smith and Fahrmeir (J Am Stat Assoc 102(478):417–431, 2007) proposed a competing method based on a spatial Bayesian variable selection in voxelwise linear regressions, with an Ising prior for latent activation indicators. In this article, we alternatively link activation probabilities to two types of latent Gaussian Markov random fields (GMRFs) via a probit model. Statistical inference in resulting high-dimensional hierarchical models is based on Markov chain Monte Carlo approaches, providing posterior estimates of activation probabilities and enhancing formation of activation clusters. Three algorithms are proposed depending on GMRF type and update scheme. An application to an active acoustic oddball experiment and a simulation study show a substantial increase in sensitivity compared to existing fMRI activation detection methods like classical SPM and the Ising model.     

Nonparametric conditional autoregressive expectile model via neural network with applications to estimating financial risk

The parametric conditional autoregressive expectiles (CARE) models have been developed to estimate expectiles, which can be used to assess value at risk and expected shortfall. The challenge lies in parametric CARE modeling is the specification of a parametric form. To avoid any model misspecification, we propose a nonparametric CARE model via neural network. The nonparametric CARE model can be estimated by a classical gradient based nonlinear optimization algorithm, and the consistency of nonparametric conditional expectile estimators is established. We then apply the nonparametric CARE model to estimating value at risk and expected shortfall of six stock indices. Empirical results for the new model is competitive with those classical models and parametric CARE models. Copyright © 2016 John Wiley & Sons, Ltd.     

Bayesian Semiparametric Density Estimation and Model Verification Using a Logistic–Gaussian Process

This article proposes a semiparametric model, which consists of parametric and nonparametric components, for density estimation. The parametric component represents the researcher's a priori beliefs about a likely family of density functions. The nonparametric component, which is modeled by a logistic–Gaussian process, allows the predictive distribution to deviate from the parametric family if it is inadequate. Bayesian hypothesis testing is used to examine the adequacy of the parametric model relative to the flexible alternative provided by the semiparametric model. The article presents a Markov chain Monte Carlo algorithm that efficiently handles the large number of parameters.     

Nonparametric Bayesian Assessment of the Order of Dependence for Binary Sequences

This article discusses inference on the order of dependence in binary sequences. The proposed approach is based on the notion of partial exchangeability of order k. A partially exchangeable binary sequence of order k can be represented as a mixture of Markov chains. The mixture is with respect to the unknown transition probability matrix θ. We use this defining property to construct a semiparametric model for binary sequences by assuming a nonparametric prior on the transition matrix θ. This enables us to consider inference on the order of dependence without constraint to a particular parametric model. Implementing posterior simulation in the proposed model is complicated by the fact that the dimension of θ changes with the order of dependence k. We discuss appropriate posterior simulation schemes based on a pseudo prior approach. We extend the model to include covariates by considering an alternative parameterization as an autologistic regression which allows for a straightforward introduction of covariates. The regression on covariates raises the additional inference problem of variable selection. We discuss appropriate posterior simulation schemes, focusing on inference about the order of dependence. We discuss and develop the model with covariates only to the extent needed for such inference.     

Semiparametric Bayesian measurement error modeling

This work presents a Bayesian semiparametric approach for dealing with regression models where the covariate is measured with error. Given that (1) the error normality assumption is very restrictive, and (2) assuming a specific elliptical distribution for errors (Student-t for example), may be somewhat presumptuous; there is need for more flexible methods, in terms of assuming only symmetry of errors (admitting unknown kurtosis). In this sense, the main advantage of this extended Bayesian approach is the possibility of considering generalizations of the elliptical family of models by using Dirichlet process priors in dependent and independent situations. Conditional posterior distributions are implemented, allowing the use of Markov Chain Monte Carlo (MCMC), to generate the posterior distributions. An interesting result shown is that the Dirichlet process prior is not updated in the case of the dependent elliptical model. Furthermore, an analysis of a real data set is reported to illustrate the usefulness of our approach, in dealing with outliers. Finally, semiparametric proposed models and parametric normal model are compared, graphically with the posterior distribution density of the coefficients.     

Modelling Dependence Uncertainty in the Extremes of Markov Chains

General theory on the extremes of stationary processes leads only to a limited representation for extreme-state behaviour, usually summarised by the extremal index. In practice this means that other quantities such as the duration of extreme episodes or aggregate of threshold exceedances within a cluster require stronger model assumptions. In this paper we propose a model based on a Markov assumption for the underlying process, with high-level transitions determined by an asymptotically motivated distribution. This idea is not new: Smith et al. (1997) first developed the statistical basis for such a procedure, which was subsequently extended by Bortot and Tawn (1998) to better handle the case of weak extremal temporal dependence for which the extremal index is unity. We adopt similar procedures to each of these earlier works, but suggest a different model for the Markov transitions. The model we use was developed by Coles and Pauli (2002) to enable a Bayesian inference of multivariate extremes that provides a posterior distribution on the status of asymptotic independence. By adopting this model in the Markov framework, we show here that the model has all the flexibility of the model developed by Bortot and Tawn (1998), but with the additional advantage of providing a posterior probability on the extremal index and inferences that take full account of the uncertainty in the extremal index. We demonstrate the methodology on both simulated data and a time series of daily rainfall that exhibit weak temporal dependence at extreme levels.     

Bayesian estimation of system reliability in Brownian stress-strength models

A stress-strength system fails as soon as the applied stress,X, is at least as much as the strength,Y, of the system. Stress and strength are time-varying in many real-life systems but typical statistical models for stress-strength systems are static. In this article, the stress and strength processes are dynamically modeled as Brownian motions. The resulting stress-strength system is then governed by a time-homogeneous Markov process with an absorption barrier at O. Conjugate as well as non-informative priors are developed for the model parameters and Markov chain sampling methods are used for posterior inference of the reliability of the stress-strength system. A generalization of this model is described next where the different stress-strength systems are assumed to be exchangeable. The proposed Bayesian analyses are illustrated in two examples where we obtain posterior estimates as well as perform model checking by cross-validation.     

A general autoregressive model with Markov switching: Estimation and consistency

In this paper, a general autoregressive model with Markov switching is considered, where the autoregression may be of an infinite order. The consistency of the maximum likelihood estimators for this model is obtained under regularity assumptions. Examples of finite and infinite order autoregressive models with Markov switching are discussed. Simulation studies with these examples illustrate the consistency and asymptotic normality of the estimators.      

Semiparametric Bayesian optimal replacement policies: application to railroad tracks

We present a Bayesian decision theoretic approach for developing replacement strategies. In so doing, we consider a semiparametric model to describe the failure characteristics of systems by specifying a nonparametric form for cumulative intensity function and by taking into account effect of covariates by a parametric form. Use of a gamma process prior for the cumulative intensity function complicates the Bayesian analysis when the updating is based on failure count data. We develop a Bayesian analysis of the model using Markov chain Monte Carlo methods and determine replacement strategies. Adoption of Markov chain Monte Carlo methods involves a data augmentation algorithm. We show the implementation of our approach using actual data from railroad tracks. Copyright © 2016 John Wiley & Sons, Ltd.     

Assessing Convergence and Mixing of MCMC Implementations via Stratification

Some posterior distributions lead to Markov chain Monte Carlo (MCMC) chains that are naturally viewed as collections of subchains. Examples include mixture models, regime-switching models, and hidden Markov models. We obtain MCMC-based estimators of posterior expectations by combining different subgroup (subchain) estimators using stratification and poststratification methods. Variance estimates of the limiting distributions of such estimators are developed. Based on these variance estimates, we propose a test statistic to aid in the assessment of convergence and mixing of chains. We compare our diagnostic with other commonly used methods. The approach is illustrated in two examples: a latent variable model for arsenic concentration in public water systems in Arizona and a Bayesian hierarchical model for Pacific sea surface temperatures. Supplementary materials, which include MATLAB codes for the proposed method, are available online.     

Estimation of the nuisance parameter for a seminartingale regression model

A nuisance parameter is introduced to the semimartingale regression model proposed by Aalen (1980), and we construct two estimators for this nuisance parameter based on the results of parametric estimation which were given by Mckeague (1986) using the method of sieves. The consistency of the estimators is also provided.     

Sampling Schemes for Bayesian Variable Selection in Generalized Linear Models

Bayesian approaches to prediction and the assessment of predictive uncertainty in generalized linear models are often based on averaging predictions over different models, and this requires methods for accounting for model uncertainty. When there are linear dependencies among potential predictor variables in a generalized linear model, existing Markov chain Monte Carlo algorithms for sampling from the posterior distribution on the model and parameter space in Bayesian variable selection problems may not work well. This article describes a sampling algorithm based on the Swendsen-Wang algorithm for the Ising model, and which works well when the predictors are far from orthogonality. In problems of variable selection for generalized linear models we can index different models by a binary parameter vector, where each binary variable indicates whether or not a given predictor variable is included in the model. The posterior distribution on the model is a distribution on this collection of binary strings, and by thinking of this posterior distribution as a binary spatial field we apply a sampling scheme inspired by the Swendsen-Wang algorithm for the Ising model in order to sample from the model posterior distribution. The algorithm we describe extends a similar algorithm for variable selection problems in linear models. The benefits of the algorithm are demonstrated for both real and simulated data.     

Variable selection for generalized varying coefficient partially linear models with diverging number of parameters

Semiparametric models with diverging number of predictors arise in many contemporary scientific areas.Variable selection for these models consists of two components:model selection for non-parametric components and selection of significant variables for the parametric portion.In this paper,we consider a variable selection procedure by combining basis function approximation with SCAD penalty.The proposed procedure simultaneously selects significant variables in the parametric components and the nonparametric components.With appropriate selection of tuning parameters,we establish the consistency and sparseness of this procedure.     

A comparison of some criteria for states selection in the latent Markov model for longitudinal data

We compare different selection criteria to choose the number of latent states of a multivariate latent Markov model for longitudinal data. This model is based on an underlying Markov chain to represent the evolution of a latent characteristic of a group of individuals over time. Then, the response variables observed at different occasions are assumed to be conditionally independent given this chain. Maximum likelihood estimation of the model is carried out through an Expectation–Maximization algorithm based on forward–backward recursions which are well known in the hidden Markov literature for time series. The selection criteria we consider are based on penalized versions of the maximum log-likelihood or on the posterior probabilities of belonging to each latent state, that is, the conditional probability of the latent state given the observed data. Among the latter criteria, we propose an appropriate entropy measure tailored for the latent Markov models. We show the results of a Monte Carlo simulation study aimed at comparing the performance of the above states selection criteria on the basis of a wide set of model specifications.     

Local consistency of Markov chain Monte Carlo methods

In this paper, we introduce the notion of efficiency (consistency) and examine some asymptotic properties of Markov chain Monte Carlo methods. We apply these results to the data augmentation (DA) procedure for independent and identically distributed observations. More precisely, we show that if both the sample size and the running time of the DA procedure tend to infinity, the empirical distribution of the DA procedure tends to the posterior distribution. This is a local property of the DA procedure, which may be, in some cases, more helpful than the global properties to describe its behavior. The advantages of using the local properties are the simplicity and the generality of the results. The local properties provide useful insight into the problem of how to construct efficient algorithms.     

基于马氏决策向量过程模型的有限阶段期望总报酬准则及其最优方程

在马氏决策向量过程模型的理论基础上,结合决策向量和相合度等新定义,进一步提出有限阶段期望总报酬准则和最优方程,并证明最优方程的解的存在性.     

随机截尾下污染数据部分线性模型的强相合估计

考虑一类新的污染数据部分线性模型,当受污染后的因变量被随机右截断时,就截断分布已知的情形,利用所获得截断观测数据构造了模型中的参数分量,非参数分量及污染系数的估计量,并在适当的条件下,证明了这些估计量的强相合性.     

随机删失下半参数估计的强相合性

主要讨论了随机删失下的部分线性模型,利用基于分布函数的核估计和最小二乘法,给出了删失情况下参数和非参数部分的估计,并证明了它们的强相合性.     

Relevance of functional flexibility for heterogeneous sales response models: A comparison of parametric and semi-nonparametric models

So far studies estimating sales response functions on the basis of store-specific data either consider heterogeneity or functional flexibility. That is why in this contribution a model is developed possessing both these features. It is a multilayer perceptron with store-specific coefficients which follows a hierarchical Bayesian framework. An appropriate Markov Chain Monte Carlo estimation technique is introduced capable to satisfy theoretical constraints (e.g. sign constraints on elasticities). The empirical study refers to a data base consisting of weekly observations of sales and prices for nine leading brands of a packaged consumer good category. The data were acquired in 81 stores over a time span of at least 61 weeks. The multilayer perceptron is compared to a strict parametric multiplicative model and turns out to be clearly superior in terms of posterior model probability. This result indicates the benefits of using a flexible model even if heterogeneity is dealt with. Estimated sales curves and elasticities demonstrate that both models differ with regard to implications on price response.