Bayesian System Reliability Prediction

In system reliability studies a common problem is the coherent assessment of system reliability, based on generic database information from components and on failure data from a system and its components in a common, but unknown, environment. A solution to this problem is given.     

大厚高压油藏物理模型饱和度分布测量方法

当前针对物理模型饱和度分布的研究主要针对于一维小岩样和二维小薄模型,针对大厚高压油藏物理模型的方法无人涉及。针对大厚高压油藏物理模型的特点,选取 射线吸收法和lambert 定律进行单个测试点的饱和度测量和计算,然后利用二维数控移动定位系统完成对整个物理模型所有测试点饱和度的测量和计算,最终利用数学插分方法得到油藏物理模型的饱和度分布等值线图,实现对大厚高压油藏二维饱和度分布的预测。通过单点精确度实验及简单验证实验研究表明,利用 射线吸收法预测大厚高压油藏二维饱和度分布可行性强。     

A survey on multivariate chi-square distributions and their applications in testing multiple hypotheses

We are concerned with three different types of multivariate chi-square distributions. Their members play important roles as limiting distributions of vectors of test statistics in several applications of multiple hypotheses testing. We explain these applications and consider the computation of multiplicity-adjusted p-values under the respective global hypothesis. By means of numerical examples, we demonstrate how much gain in level exhaustion or, equivalently, power can be achieved with corresponding multivariate multiple tests compared with approaches which are only based on univariate marginal distributions and do not take the dependence structure among the test statistics into account. As a further contribution of independent value, we provide an overview of essentially all analytic formulas for computing multivariate chi-square probabilities of the considered types which are available up to present. These formulas were scattered in the previous literature and are presented here in a unified manner.     

Estimating the Prevalence of Illegal Drug Use

A generic globally applicable integral equation allows the estimation of hidden drug use. Formulas for estimating the prevalence of drug use are provided for both the assumed and known latency period. Estimates of prevalence of opiate use in Ireland ranging from 7,955 [5,789, 10,122] to 11,986 [8,792, 15,181] with prevalence rates ranging from 2.7 [2.0, 3.5] per 1,000 to 4.1 [3.0, 5.2] were obtained. The model can be implemented by service providers and planners if numbers of first treatments for drug use and the latency period are known or assumed.     

Diagnostic tools in generalized Weibull linear regression models

We propose some statistical tools for diagnosing the class of generalized Weibull linear regression models [A.A. Prudente and G.M. Cordeiro, Generalized Weibull linear models, Comm. Statist. Theory Methods 39 (2010), pp. 3739–3755]. This class of models is an alternative means of analysing positive, continuous and skewed data and, due to its statistical properties, is very competitive with gamma regression models. First, we show that the Weibull model induces ma-ximum likelihood estimators asymptotically more efficient than the gamma model. Standardized residuals are defined, and their statistical properties are examined empirically. Some measures are derived based on the case-deletion model, including the generalized Cook's distance and measures for identifying influential observations on partial F-tests. The results of a simulation study conducted to assess behaviour of the global influence approach are also presented. Further, we perform a local influence analysis under the case-weights, response and explanatory variables perturbation schemes. The Weibull, gamma and other Weibull-type regression models are fitted into three data sets to illustrate the proposed diagnostic tools. Statistical analyses indicate that the Weibull model fitted into these data yields better fits than other common alternative models.     

The behavior of estimated measures of association in small and moderate sample sizes for 2×3 tables

We consider the distributions of Goodman and Kruskal's G, Kendall's tau-b, and correlation coefficients rho and rho-s for sample sizes 10‘10’40 from 2×3 tables. The results are compared with asymptotic theory. It is found that the convergence of G to its asymptotic normal distribution is much slower than the convergence of the other measures to theirs, and that G is more likely to be significantly biased. However, the variances and biases of all four measures come close to their asymptotic values for quite moderate sample sizes.     

A simple monotone process with application to radiocarbon-dated depth chronologies

Summary.  We propose a new and simple continuous Markov monotone stochastic process and use it to make inference on a partially observed monotone stochastic process. The process is piecewise linear, based on additive independent gamma increments arriving in a Poisson fashion. An independent increments variation allows very simple conditional simulation of sample paths given known values of the process. We take advantage of a reparameterization involving the Tweedie distribution to provide efficient computation. The motivating problem is the establishment of a chronology for samples taken from lake sediment cores, i.e. the attribution of a set of dates to samples of the core given their depths, knowing that the age–depth relationship is monotone. The chronological information arises from radiocarbon (¹⁴C) dating at a subset of depths. We use the process to model the stochastically varying rate of sedimentation.     

Testing Homogeneity in Gamma Mixture Models

This paper characterizes the asymptotic behaviour of the likelihood ratio test statistic (LRTS) for testing homogeneity (i.e. no mixture) against gamma mixture alternatives. Under the null hypothesis, the LRTS is shown to be asymptotically equivalent to the square of Davies's Gaussian process test statistic and diverges at a log n rate to infinity in probability. Based on the asymptotic analysis, we propose and demonstrate a computationally efficient method to simulate the null distributions of the LRTS for small to moderate sample sizes.     

Bayesian estimation in Kibble's bivariate gamma distribution

The authors describe Bayesian estimation for the parameters of the bivariate gamma distribution due to Kibble (1941). The density of this distribution can be written as a mixture, which allows for a simple data augmentation scheme. The authors propose a Markov chain Monte Carlo algorithm to facilitate estimation. They show that the resulting chain is geometrically ergodic, and thus a regenerative sampling procedure is applicable, which allows for estimation of the standard errors of the ergodic means. They develop Bayesian hypothesis testing procedures to test both the dependence hypothesis of the two variables and the hypothesis of equal means. They also propose a reversible jump Markov chain Monte Carlo algorithm to carry out the model selection problem. Finally, they use sets of real and simulated data to illustrate their methodology.     

伽马分布的形状参数的一种估计

本文利用服从伽马分布的随机变量 X 与其逆 X~(-1)的协方差只与形状参数有关这一性质,给出伽马分布形状参数的所谓自逆协方差估计,进而构造了相应的无偏估计,并证明了这类估计的大样本性质:强相合性以及渐近正态性。