�������ֲ������¶�Ԫ����ģ���к����������

In Bayesian analysis, the Markov Chain Monte Carlo (MCMC) algorithm is an efficient and simple method to compute posteriors. However, the chain may appear to converge while the posterior is improper, which will leads to incorrect statistical inferences. In this paper, we focus on the necessary and sufficient conditions for which improper hierarchical priors can yield proper posteriors in a multivariate linear model. In addition, we carry out a simulation study to illustrate the theoretical results, in which the Gibbs sampling and Metropolis-Hasting sampling are employed to generate the posteriors.     

ij������ֵ����ʱCopula���Ž��Ⱥ�ṹ�Ϳ�խ��

本文证明了copula , 生存copula , 对偶copula 和伴随copula 关于copula的复合运算构成一个四元群, 给出了当某个单点值给定时它们的最优上下界. 计算了, 时copula最优上下界的宽窄度, 并与时的宽窄度进行了比较.     

������ݽ����̶�ЧӦģ�͵�Э����������

??In this paper, we focus on the tests for covariance matrices in panel data model with interactive fixed effects. For the problem of testing identity and sphericity of covariance matrices, we first propose test statistics based on the estimators of the trace of covariance matrices. Under both the null hypothesis and the alternatives, we establish the asymptotic distributions of the proposed test statistics under some regularity conditions, and we further show that the proposed tests are distribution free. Subsequently simulation studies suggest that the proposed tests perform well under the high dimensional panel data.     

һ�������ΥԼ�ʲ�����Ornstein-Uhlenbeck���̻̿��Ĺ�Ʊ�������ٱ��պ�Ͷ������

In this paper, the insurer is allowed to buy reinsurance and allocate his money among three financial securities: a defaultable corporate zero-coupon bond, a default-free bank account, and a stock, while the instantaneous rate of the stock is described by an Ornstein-Uhlenbeck process. The objective is to maximize the exponential utility of the terminal wealth. We decompose the original optimization problem into two subproblems: a pre-default case and a post-default case. Using dynamic programming principle, and then solving the corresponding HJB equations, we derive the closed-form solutions for the optimal reinsurance and investment strategies and the corresponding value functions     

�лָ������ɷ�֧���̵�Feller����

首先, 当$Q$是一个拟单调的q矩阵的时候, 我们找出最小的$Q$函数是一个Feller的转移函数的准则. 然后我们把这个结论应用于生成分支q矩阵并得到相应的生成分支过程的Feller准则. 特别地, 设$\theta$是分支q矩阵中的非线性数, 总是存在一个分点$\theta_0$满足$1\leq\theta_0\leq2$或$\theta_0<+\infty$使得 生成分支过程是否是Feller的要依据$\theta<\theta_0$或者$\theta>\theta_0$.     

��λ���ع�ı�Ҷ˹����ѡ�񷽷�

??When the data has heavy tail feature or contains outliers, conventional variable selection methods based on penalized least squares or likelihood functions perform poorly. Based on Bayesian inference method, we study the Bayesian variable selection problem for median linear models. The Bayesian estimation method is proposed by using Bayesian model selection theory and Bayesian estimation method through selecting the Spike and Slab prior for regression coefficients, and the effective posterior Gibbs sampling procedure is also given. Extensive numerical simulations and Boston house price data analysis are used to illustrate the effectiveness of the proposed method.     

��������֦��������ǿ�������ɺ�Shannon-McMillan����

In this paper, we first study the strong law of large numbers for the frequencies of occurrence of random ordered couples of states for nonhomogeneous bifurcating Markov chains indexed by a binary tree. Then the strong law of large numbers are studied for functions of the nonhomogeneous bifurcating Markov chains indexed by a binary tree. As a corollary, we obtain the shannon-McMillan theorem for these Markov chains with finite state space.     

�����췽��ģ����̫���ܵ���˻����ɿ��Է����е�Ӧ��

Solar cell is the basic component of satellite photovoltaic panels with complicated redundant system structure. Its reliability plays an important role in the system, and its performance shows a degradation trend over time. In this paper, study is conducted for the solar cell degradation modeling and reliability analysis based on practical testing results. Specifically, we illustrate an accelerated test for the attenuation ratio character test under different accumulative irradiation levels, focusing on the heteroscedasticity of the collected testing data. A heteroscedastic linear model is proposed, and the life distribution of the photovoltaic panel is obtained by using Fiducial method. A numerical example is shown for the purpose of illustration.     

����O-U���̵�һ��ע��

For the Ornstein-Uhlenbeck process, we prove that if the classical Kalman's condition doesn't hold, then it can be divided into a deterministic process plus a stochastic process after a linear and invertible transformation.     

�������������ݰ����ģ�͵Ĺ���

We propose a new estimation method for the parameters of a partial functional linear model when the parameter curve is subject to monotone constraint. The proposed estimators are implemented under the nonlinear mixed effects model framework. The small sample properties are illustrated through a simulation experiment.     

Reference Bayesian analysis of inverse Gaussian degradation process

In this paper, objective Bayesian method is applied to analyze degradation model based on the inverse Gaussian process. Noninformative priors (Jefferys prior and two reference priors) for model parameters are obtained and their properties are discussed. Moreover, we propose a class of modified reference priors to remedy weaknesses of the usual reference priors and show that the modified reference priors not only have proper posterior distributions but also have probability matching properties for model parameters. Gibbs sampling algorithms for Bayesian inference based on the Jefferys prior and the modified reference priors are studied. Simulations are conducted to compare the objective Bayesian estimates with the maximum likelihood estimates and subjective Bayesian estimates and shows better performance of the objective method than the other two estimates especially for the case of small sample size. Finally, two real data examples are analyzed for illustration.     

Objective Bayesian analysis for competing risks model with Wiener degradation phenomena and catastrophic failures

In this paper, the objective Bayesian method is applied to investigate the competing risks model involving both catastrophic and degradation failures. By modeling soft failure as the Wiener degradation process, and hard failures as a Weibull distribution, we obtain the noninformative priors (Jefferys prior and two reference priors) for the parameters. Moreover, we show that their posterior distributions have good properties and we propose Gibbs sampling algorithms for the Bayesian inference based on the Jefferys prior and two reference priors. Some simulation studies are conducted to illustrate the superiority of objective Bayesian method. Finally, we apply our methods to two real data examples and compare the objective Bayesian estimates with the other estimates.     

Objective Bayesian analysis for the accelerated degradation model based on the inverse Gaussian process

The inverse Gaussian process is an attractive stochastic process to model monotone degradation paths in degradation analysis. In this paper, we propose an objective Bayesian method to analyze the accelerated degradation model based on the inverse Gaussian process. Noninformative priors including the Jeffreys prior and reference priors are derived, and the propriety of the posteriors under each prior is validated. A simulation study is carried out to investigate the performance of the approach compared with the maximum likelihood method and the Bootstrap method. Numerical results show that the proposed method has better performance in terms of the mean squared error and the frequentist coverage probability. The reference prior ${\pi }_{R_2}$ is recommended to use in practice. Finally, the Bayesian approach is applied to a real data.     

Priors for Bayesian adaptive spline smoothing

Adaptive smoothing has been proposed for curve-fitting problems where the underlying function is spatially inhomogeneous. Two Bayesian adaptive smoothing models, Bayesian adaptive smoothing splines on a lattice and Bayesian adaptive P-splines, are studied in this paper. Estimation is fully Bayesian and carried out by efficient Gibbs sampling. Choice of prior is critical in any Bayesian non-parametric regression method. We use objective priors on the first level parameters where feasible, specifically independent Jeffreys priors (right Haar priors) on the implied base linear model and error variance, and we derive sufficient conditions on higher level components to ensure that the posterior is proper. Through simulation, we demonstrate that the common practice of approximating improper priors by proper but diffuse priors may lead to invalid inference, and we show how appropriate choices of proper but only weakly informative priors yields satisfactory inference.     

Bayesian analysis of conditional autoregressive models

Conditional autoregressive (CAR) models have been extensively used for the analysis of spatial data in diverse areas, such as demography, economy, epidemiology and geography, as models for both latent and observed variables. In the latter case, the most common inferential method has been maximum likelihood, and the Bayesian approach has not been used much. This work proposes default (automatic) Bayesian analyses of CAR models. Two versions of Jeffreys prior, the independence Jeffreys and Jeffreys-rule priors, are derived for the parameters of CAR models and properties of the priors and resulting posterior distributions are obtained. The two priors and their respective posteriors are compared based on simulated data. Also, frequentist properties of inferences based on maximum likelihood are compared with those based on the Jeffreys priors and the uniform prior. Finally, the proposed Bayesian analysis is illustrated by fitting a CAR model to a phosphate dataset from an archaeological region.     

Computation of reference Bayesian inference for variance components in longitudinal studies

Generalized linear mixed models (GLMMs) have been applied widely in the analysis of longitudinal data. This model confers two important advantages, namely, the flexibility to include random effects and the ability to make inference about complex covariances. In practice, however, the inference of variance components can be a difficult task due to the complexity of the model itself and the dimensionality of the covariance matrix of random effects. Here we first discuss for GLMMs the relation between Bayesian posterior estimates and penalized quasi-likelihood (PQL) estimates, based on the generalization of Harville’s result for general linear models. Next, we perform fully Bayesian analyses for the random covariance matrix using three different reference priors, two with Jeffreys’ priors derived from approximate likelihoods and one with the approximate uniform shrinkage prior. Computations are carried out via the combination of asymptotic approximations and Markov chain Monte Carlo methods. Under the criterion of the squared Euclidean norm, we compare the performances of Bayesian estimates of variance components with that of PQL estimates when the responses are non-normal, and with that of the restricted maximum likelihood (REML) estimates when data are assumed normal. Three applications and simulations of binary, normal, and count responses with multiple random effects and of small sample sizes are illustrated. The analyses examine the differences in estimation performance when the covariance structure is complex, and demonstrate the equivalence between PQL and the posterior modes when the former can be derived. The results also show that the Bayesian approach, particularly under the approximate Jeffreys’ priors, outperforms other procedures.     

Bayesian reference analysis for Gaussian Markov random fields

Gaussian Markov random fields (GMRF) are important families of distributions for the modeling of spatial data and have been extensively used in different areas of spatial statistics such as disease mapping, image analysis and remote sensing. GMRFs have been used for the modeling of spatial data, both as models for the sampling distribution of the observed data and as models for the prior of latent processes/random effects; we consider mainly the former use of GMRFs. We study a large class of GMRF models that includes several models previously proposed in the literature. An objective Bayesian analysis is presented for the parameters of the above class of GMRFs, where explicit expressions for the Jeffreys (two versions) and reference priors are derived, and for each of these priors results on posterior propriety of the model parameters are established. We describe a simple MCMC algorithm for sampling from the posterior distribution of the model parameters, and study frequentist properties of the Bayesian inferences resulting from the use of these automatic priors. Finally, we illustrate the use of the proposed GMRF model and reference prior for studying the spatial variability of lip cancer cases in the districts of Scotland over the period 1975-1980.     

Noninformative priors for the ratio of the shape parameters of two Weibull distributions

The Weibull distribution is one of the most widely used lifetime distributions in reliability engineering. Here, the noninformative priors for the ratio of the shape parameters of two Weibull models are introduced. The first criterion used is the asymptotic matching of the coverage probabilities of Bayesian credible intervals with the corresponding frequentist coverage probabilities. We develop the probability matching priors for the ratio of the shape parameters using the following matching criteria: quantile matching, matching of the distribution function, highest posterior density matching, and matching via inversion of the test statistics. We obtain one particular prior that meets all the matching criteria. Next, we derive the reference priors for different groups of ordering. Our findings show that some of the reference priors satisfy a first-order matching criterion and the one-at-a-time reference prior is a second-order matching prior. Lastly, we perform a simulation study and provide a real-world example.     

Bayesian shrinkage prediction for the regression problem

We consider Bayesian shrinkage predictions for the Normal regression problem under the frequentist Kullback-Leibler risk function.Firstly, we consider the multivariate Normal model with an unknown mean and a known covariance. While the unknown mean is fixed, the covariance of future samples can be different from that of training samples. We show that the Bayesian predictive distribution based on the uniform prior is dominated by that based on a class of priors if the prior distributions for the covariance and future covariance matrices are rotation invariant.Then, we consider a class of priors for the mean parameters depending on the future covariance matrix. With such a prior, we can construct a Bayesian predictive distribution dominating that based on the uniform prior.Lastly, applying this result to the prediction of response variables in the Normal linear regression model, we show that there exists a Bayesian predictive distribution dominating that based on the uniform prior. Minimaxity of these Bayesian predictions follows from these results.     

Bayesian analysis of Birnbaum–Saunders distribution via the generalized ratio-of-uniforms method

This paper deals with the Bayesian inference for the parameters of the Birnbaum–Saunders distribution. We adopt the inverse-gamma priors for the shape and scale parameters because the continuous conjugate joint prior distribution does not exist and the reference prior (or independent Jeffreys’ prior) results in an improper posterior distribution. We propose an efficient sampling algorithm via the generalized ratio-of-uniforms method to compute the Bayesian estimates and the credible intervals. One appealing advantage of the proposed procedure over other sampling techniques is that it efficiently generates independent samples from the required posterior distribution. Simulation studies are conducted to investigate the behavior of the proposed method, and two real-data applications are analyzed for illustrative purposes.