Algebraic Markov Bases and MCMC for Two-Way Contingency Tables

ABSTRACT.  The Diaconis–Sturmfels algorithm is a method for sampling from conditional distributions, based on the algebraic theory of toric ideals. This algorithm is applied to categorical data analysis through the notion of Markov basis. An application of this algorithm is a non-parametric Monte Carlo approach to the goodness of fit tests for contingency tables. In this paper, we characterize or compute the Markov bases for some log-linear models for two-way contingency tables using techniques from Computational Commutative Algebra, namely Gröbner bases. This applies to a large set of cases including independence, quasi-independence, symmetry, quasi-symmetry. Three examples of quasi-symmetry and quasi-independence from Fingleton ( Models of category counts , Cambridge University Press, Cambridge, 1984) and Agresti ( An Introduction to categorical data analysis , Wiley, New York, 1996) illustrate the practical applicability and the relevance of this algebraic methodology.     

非对称随机波动模型在中国股市的应用

为了探测随机波动模型的非对称特征,修改传统的随机波动模型建立非对称的随机波动模型,采用基于马尔可夫链蒙特卡洛(MCMC)模拟的贝叶斯分析对模型进行参数估计。对中国深圳、上海股市波动进行实证研究发现,非对称随机波动模型能较好地探测波动存在的非对称波动。与GJR-GARCH模型相比,非对称随机波动模型预测效果更好。     

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.     

A distance-based diagnostic for trans-dimensional Markov chains

Over the last decade the use of trans-dimensional sampling algorithms has become endemic in the statistical literature. In spite of their application however, there are few reliable methods to assess whether the underlying Markov chains have reached their stationary distribution. In this article we present a distance-based method for the comparison of trans-dimensional Markov chain sample output for a broad class of models. This diagnostic will simultaneously assess deviations between and within chains. Illustration of the analysis of Markov chain sample-paths is presented in simulated examples and in two common modelling situations: a finite mixture analysis and a change-point problem.     

Bayesian automatic polynomial wavelet regression

This paper considers the problem of Bayesian automatic polynomial wavelet regression (PWR). We propose three different Bayesian methods based on integrated likelihood, conditional empirical Bayes, and reversible jump Markov chain Monte Carlo (MCMC). From the simulation results, we find that the proposed methods are similar to or superior to the existing ones.     

Transform Analysis and Asset Pricing for Affine Jump‐diffusions

In the setting of ‘affine’ jump‐diffusion state processes, this paper provides an analytical treatment of a class of transforms, including various Laplace and Fourier transforms as special cases, that allow an analytical treatment of a range of valuation and econometric problems. Example applications include fixed‐income pricing models, with a role for intensity‐based models of default, as well as a wide range of option‐pricing applications. An illustrative example examines the implications of stochastic volatility and jumps for option valuation. This example highlights the impact on option ‘smirks’ of the joint distribution of jumps in volatility and jumps in the underlying asset price, through both jump amplitude as well as jump timing.     

Bayesian quantile regression and variable selection for partial linear single-index model: Using free knot spline

Partial linear single-index model (PLSIM) has both the flexibility of nonparametric treatment and interpretability of linear term, yet existing literatures about it mainly focused on mean regression, and quantile regression analysis is scarce. Based on free knot spline approximation, we apply asymmetric Laplace distribution to implement Bayesian quantile regression, and perform variable selection in linear term and index vector via binary indicators. Our approach is exempt from regularity conditions in frequentist method, and could execute variable selection and quantile regression under mutual posterior correction, which is also the first work to implement them jointly for PLSIM in fully Bayesian framework. The numerical simulation manifests the superiority of our approach to previous methods, which embodied in better efficiency of variable selection, index vector estimates and link function approximation with different error distributions. For illustration of its application, we build a power consumption model of A₂/O process in wastewater treatment and emphatically analyze the impact of water quality factors.     

污染减排、外部性与环境技术创新:来自省级环境专利的证据

基于污染外部性机制,推导出含有污染减排外部性的空间Tobit设定。利用中国31个省级环境专利数据,在控制技术外部性的情况下,分离出污染减排外部性对环境技术创新的影响。贝叶斯MCMC方法估计结果显示,污染减排对环境技术创新有显著的正的直接效应,而且还存在显著的负外部效应,从而污染减排的总效应为不显著的正影响,表明忽略污染减排外部性对环境技术创新有重要影响。进而探讨避免污染减排外部性的政策含义。     

Dynamic generalized extreme value modeling via particle filters

This article addresses how particle filters compare to MCMC methods for posterior density approximations of a model that allows for a dynamic state with fixed parameters and where the observation equation is nonlinear. This is a problem that was not been well studied in the specialized literature. We prove that these state and parameter estimations can be achieved via particle filter methods without the need of more expensive Forward Filtering Backward Sampling (FFBS) simulation. Estimation of a time-varying extreme value model via the generalized extreme value distribution is considered using these particle filter methods and compared to a MCMC algorithm that involves a variety of Metropolis-Hastings steps. We illustrate and compare the different methodologies with simulated data and some minimum daily stock returns occurring monthly from January 4, 1990 to December 28, 2007 using the Tokyo Stock Price Index (TOPIX).     

Simultaneous Bayesian modeling of longitudinal and survival data in breast cancer patients

Abstract

Using simultaneous Bayesian modeling, an attempt is made to analyze data on the size of lymphedema occurring in the arms of breast cancer patients after breast cancer surgery (as the longitudinal data) and the time interval for disease progression (as the time-to-event occurrence). A model based on a multivariate skew t distribution is shown to provide the best fit. This outcome was confirmed by simulation studies too.