A stable manifold MCMC method for high dimensions

We combine two important recent advancements of MCMC algorithms: first, methods utilizing the intrinsic manifold structure of the parameter space; then, algorithms effective for targets in infinite-dimensions with the critical property that their mixing time is robust to mesh refinement.     

散度偏大计数数据回归模型的变量选择与模型比较

讨论了具有散度偏大特征计数数据的建模与拟合问题.针对导致数据散度偏大的原因和常用的几类候选模型的结构,分别给出了关于嵌套模型的模型与变量同时选择的Bayes方法和关于非嵌套模型的模型检验与比较方法,并在此基础上进一步完善,提出了较为系统完整的模型与变量选择方法.实际例子说明了方法的具体实现过程和有效性.     

Estimating Expected Information Gains for Experimental Designs With Application to the Random Fatigue-Limit Model

Expected gain in Shannon information is commonly suggested as a Bayesian design evaluation criterion. Because estimating expected information gains is computationally expensive, examples in which they have been successfully used in identifying Bayes optimal designs are both few and typically quite simplistic. This article discusses in general some properties of estimators of expected information gains based on Markov chain Monte Carlo (MCMC) and Laplacian approximations. We then investigate some issues that arise when applying these methods to the problem of experimental design in the (technically nontrivial) random fatigue-limit model of Pascual and Meeker. An example comparing follow-up designs for a laminate panel study is provided.     

MEXIT: Maximal un-coupling times for stochastic processes

Classical coupling constructions arrange for copies of the same Markov process started at two different initial states to become equal as soon as possible. In this paper, we consider an alternative coupling framework in which one seeks to arrange for two different Markov (or other stochastic) processes to remain equal for as long as possible, when started in the same state. We refer to this “un-coupling” or “maximal agreement” construction as MEXIT, standing for “maximal exit”. After highlighting the importance of un-coupling arguments in a few key statistical and probabilistic settings, we develop an explicit MEXIT construction for stochastic processes in discrete time with countable state-space. This construction is generalized to random processes on general state-space running in continuous time, and then exemplified by discussion of MEXIT for Brownian motions with two different constant drifts.     

基于SV-M模型的全球股市交易指数波动与收益率的多分辨分析

为了有效揭示收益率与波动的关系,本文采用极大重叠离散小波变换,将收益率分解在不同的交易周期上,建立各自的SV-M模型,考察各周期收益率与波动的关系及其他参数随交易周期的变化情况。实证表明,全球各大股票交易市场的收益率,收益率与波动的在较小的尺度上关系不显著,但在较大尺度上,存在明显的正相关关系;且随交易周期的增长,收益率的相关性及波动持续性逐渐增强。     

Pairwise clustering using a Monte Carlo Markov Chain

In this work an application of MCMC is proposed for unsupervised data classification, in conjunction with a novel pairwise objective function, which is shown to work well in situations where clusters to be identified have a strong overlap, and the centroid oriented methods (such as K-means) fail by construction. In particular, an exceptionally simple but difficult situation is addressed when cluster centroids coincide, and one can differentiate between the clusters only on the basis of their variance. Performance of the proposed approach is tested on synthetic and real datasets.     

矿难的贝叶斯分析——基于变点与事件分析结合的方法

本文利用贝叶斯方法与事件分析结合,建立了矿山死亡事故变点预测模型,并且文中将采用了变点与重大事件相关分析的方法,探究了矿难变化的规律,分析了事故发生原因,通过对每次矿难发生数的拐点处的事件分析,发现价格、产量、产权等的变化飙升固然会驱动矿难事故上升,但是通过加大媒体监督透明度与深度,加大官员问责力度,强化法律监督,设计合理的管理机制并严格执行可以有效地减少矿难的发生有效途径,文章采用我国安全生产局数据,利用蒙特卡洛马尔科夫方法进行了验证,同时也证明了该方法的预测突发事件的有效性.本文利用博弈思想对矿难的发生进行了经济学分析。     

Heteroscedastic normal-exponential mixture models: Bayesian and classical approaches

In this paper, we introduce a Bayesian analysis for mixture of distributions belonging to the exponential family. As a special case we consider a mixture of normal exponential distributions including joint modeling of the mean and variance. We also consider joint modeling of the mean and variance heterogeneity. Markov Chain Monte Carlo (MCMC) methods are used to obtain the posterior summaries of interest. We also introduce and apply an EM algorithm, where the maximization is obtained applying the Fisher scoring algorithm. Finally, we also include analysis of real data sets to illustrate the proposed methodology.     

Bayesian Inference for the RC(m) Association Model

Describing the structure in a two-way contingency table in terms of an RC(m) association model, we are concerned with the computation of posterior distributions of the model parameters using prior distributions which take into account the nonlinear restrictions of the model. We are further involved with the determination of the order of association m, based on Bayesian arguments. Using projection methods, a prior distribution over the parameters of the simpler RC(m) model is induced from a prior of the parameters of the saturated model. The fit of the assumed RC(m) model is evaluated using the posterior distribution of its distance from the full model. Our methods are illustrated with a popular dataset.