Transformed generalized linear models

The estimation of data transformation is very useful to yield response variables satisfying closely a normal linear model. Generalized linear models enable the fitting of models to a wide range of data types. These models are based on exponential dispersion models. We propose a new class of transformed generalized linear models to extend the Box and Cox models and the generalized linear models. We use the generalized linear model framework to fit these models and discuss maximum likelihood estimation and inference. We give a simple formula to estimate the parameter that index the transformation of the response variable for a subclass of models. We also give a simple formula to estimate the r$r$th moment of the original dependent variable. We explore the possibility of using these models to time series data to extend the generalized autoregressive moving average models discussed by Benjamin et al. [Generalized autoregressive moving average models. J. Amer. Statist. Assoc. 98, 214–223]. The usefulness of these models is illustrated in a simulation study and in applications to three real data sets.     

应用同伦正则化算法反演二维溶质运移模型中的弥散系数

利用同伦正则化算法探讨了二维对流弥散方程的依赖空间变量的弥散系数反演问题.讨论了初始迭代值、数值微分步长、以及收敛精度对算法实现的影响.数值模拟表明,同伦正则化算法对于此类参数反演问题是一种有效的方法.     

On the computation of the noncentral <Emphasis Type="Italic">F</Emphasis> and noncentral beta distribution

Unfortunately many of the numerous algorithms for computing the comulative distribution function (cdf) and noncentrality parameter of the noncentral F and beta distributions can produce completely incorrect results as demonstrated in the paper by examples. Existing algorithms are scrutinized and those parts that involve numerical difficulties are identified. As a result, a pseudo code is presented in which all the known numerical problems are resolved. This pseudo code can be easily implemented in programming language C or FORTRAN without understanding the complicated mathematical background. Symbolic evaluation of a finite and closed formula is proposed to compute exact cdf values. This approach makes it possible to check quickly and reliably the values returned by professional statistical packages over an extraordinarily wide parameter range without any programming knowledge. This research was motivated by the fact that a very useful table for calculating the size of detectable effects for ANOVA tables contains suspect values in the region of large noncentrality parameter values compared to the values obtained by Patnaik’s 2-moment central-F approximation. The cause is identified and the corrected form of the table for ANOVA purposes is given. The accuracy of the approximations to the noncentral-F distribution is also discussed. The authors wish to thank Mr. Richárd Király for his preliminary work. The authors are grateful to the Editor and Associate Editor of STCO and the unknown reviewers for their helpful suggestions.     

The Bayes rule of the variance parameter of the hierarchical normal and inverse gamma model under Stein's loss

For the variance parameter of the hierarchical normal and inverse gamma model, we analytically calculate the Bayes rule (estimator) with respect to a prior distribution IG (alpha, beta) under Stein's loss function. This estimator minimizes the posterior expected Stein's loss (PESL). We also analytically calculate the Bayes rule and the PESL under the squared error loss. Finally, the numerical simulations exemplify that the PESLs depend only on alpha and the number of observations. The Bayes rules and PESLs under Stein's loss are unanimously smaller than those under the squared error loss.     

Extended scaled prediction variance optimality for modified central composite design

Robust parameter designs (RPDs) enable the experimenter to discover how to modify the design of the product to minimize the effect due to variation from noise sources. The aim of this article is to show how this amount of work can be reduced under modified central composite design (MCCD). We propose a measure of extended scaled prediction variance (ESPV) for evaluation of RPDs on MCCD. Using these measures, we show that we can check the error or bias associated with estimating the model parameters and suggest the values of α recommended for MCCS under minimum ESPV.     

A new confidence interval in measurement error model with the reliability ratio known

For the slope parameter of the measurement error model with the reliability ratio known, this article constructs a fiducial generalized confidence interval (FGCI) which is proved to have correct asymptotic coverage. Simulation results demonstrate that the FGCI often outperforms the existing intervals in terms of empirical coverage probability, average interval length, and false parameter coverage rate. Two examples are also provided to illustrate our approach.     

一种简单有效的二项分布比例参数似然比置信区间的求法

对二项分布比例参数p的似然比置信区间,提出一种简便求解方法。在平均覆盖率、平均区间长度及区间长度的95%置信区间准则下与WScore、Plus4、Jeffreys置信区间进行模拟比较。试验表明,在二项分布b(n,p)的参数n≥20且p∈(0.1,0.9)时,该方法获取的似然比置信区间性能优良。当点估计p值不是接近于0或1且n≥20时,推荐使用本方法获取p的置信区间。     

The discrepancy of P-values and posterior probability: Two-sided hypothesis of variance of normal distribution with unknown mean

This article is concerned with the comparison of P-value and Bayesian measure in point null hypothesis for the variance of Normal distribution with unknown mean. First, using fixed prior for test parameter, the posterior probability is obtained and compared with the P-value when an appropriate prior is used for the mean parameter. In the second, lower bounds of the posterior probability of H₀ under a reasonable class of prior are compared with the P-value. It has been shown that even in the presence of nuisance parameters, these two approaches can lead to different results in the statistical inference.     

Assessing reporting patterns of child sexual abuse within the Catholic Church using discontinuities in model parameter timeseries

Prior to 2002, little was known about sexual abuse within the Catholic Church. After the Boston Globe broke the story about John Geoghan - a priest in the Boston Archdiocese who was accused of abusing numerous children, convicted of one count of indecent assault, and eventually murdered in prison - the Church had many questions to answer. To this end, the United States Conference of Catholic Bishops (USCCB) commissioned John Jay College of Criminal Justice to research the nature and scope, as well as the causes and context of child sexual abuse within the Catholic Church.This research analyzes the data from the John Jay studies using a new quantitative technique, capable of adjusting for distortions introduced by delays in abuse reporting. By isolating discontinuities in model parameter timeseries, we determine changes in reporting patterns occurred during the period 1982-1988. A posteriori to the analysis, we provide some possible explanations for the changes in abuse reporting associated with the change-point. While the scope of this paper is limited to presenting a new methodological approach within the frame of a particular case study, the techniques are more broadly applicable in settings where reporting lag is manifested.     

TESTS FOR EXPONENTIALITY AGAINST GAMMA ALTERNATIVES USING NORMALIZED WAITING TIMES

In this article, we present some tests for Exponentiality against Gamma alternatives by using normalized waiting times. The test is constructed by using a quadratic form. The asymptotic distribution of the proposed test is derived. The power of the test is computed through Monte Carlo simulation and is compared with Linhart (1965) test, Bain and Engelhardt (1975) test and Keating et al. (1990) test.