A Benchmarking Approach to Forecast Combination

This article is concerned with the development of a statistical model-based approach to optimally combine forecasts derived from an extrapolative model, such as an autoregressive integrated moving average (ARIMA) time series model, with forecasts of a particular characteristic of the same series obtained from independent sources. The methods derived combine the strengths of all forecasting approaches considered in the combination scheme. The implications of the general theory are investigated in the context of some commonly encountered seasonal ARIMA models. An empirical example to illustrate the method is included.     

Estimability,testability, and connectedness in the cell means model

This article discusses the conocepts of estimation of parameters and tests of linear hyp hypotheses in the fixed effects (ANOVA) model with some zero cell frequencies. The concept of connected experiments is generalized ot he cel1-means model and relatedtotheest imation of and test of hy pothe sesa bout the cell means. A simple test for connectedness is developed as a n atur alre sult of thea nalysis of the constrai nedcell means model. The a nalysis of unconnected experimentsisbriefly discussed. Examples are presented toillustrate the ideas.     

Statistical Reversion toward the Mean: More Universal than Regression toward the Mean

Schmittlein discussed the lack of universality of regression toward the mean. The present note emphasizes the universality of a similar effect, dubbed “reversion” toward the mean, defined as the shift in conditional expectation of the upper or lower portion of a distribution. Reversion toward the mean is a useful concept for statistical reasoning in applications and is more self-evidently plausible than regression toward the mean.     

Comparing alternating logistic regressions to other approaches to modelling correlated binary data

Alternating logistic regressions (ALRs) seem to offer some of the advantages of marginal models estimated via generalized estimating equations (GEE) and generalized linear mixed models (GLMMs). Via simulation study we compared ALRs to marginal models estimated via GEE and subject-specific models estimated via GLMMs, with a focus on estimation of the correlation structure in three-level data sets (e.g. students in classes in schools). Data set size and structure, and amount of correlation in the data sets were varied. For simple correlation structures, ALRs performed well. For three-level correlation structures, all approaches, but especially ALRs, had difficulty assigning the correlation to the correct level, though sample sizes used were small. In addition, ALRs and GEEs had trouble attaching correct inference to the mean effects, though this improved as overall sample size improved. ALRs are a valuable addition to the data analyst's toolkit, though care should be taken when modelling data with three-level structures.     

Testing in semiparametric models with interaction, with applications to gene–environment interactions

Summary.  Motivated from the problem of testing for genetic effects on complex traits in the presence of gene–environment interaction, we develop score tests in general semiparametric regression problems that involves Tukey style 1 degree-of-freedom form of interaction between parametrically and non-parametrically modelled covariates. We find that the score test in this type of model, as recently developed by Chatterjee and co-workers in the fully parametric setting, is biased and requires undersmoothing to be valid in the presence of non-parametric components. Moreover, in the presence of repeated outcomes, the asymptotic distribution of the score test depends on the estimation of functions which are defined as solutions of integral equations, making implementation difficult and computationally taxing. We develop profiled score statistics which are unbiased and asymptotically efficient and can be performed by using standard bandwidth selection methods. In addition, to overcome the difficulty of solving functional equations, we give easy interpretations of the target functions, which in turn allow us to develop estimation procedures that can be easily implemented by using standard computational methods. We present simulation studies to evaluate type I error and power of the method proposed compared with a naive test that does not consider interaction. Finally, we illustrate our methodology by analysing data from a case–control study of colorectal adenoma that was designed to investigate the association between colorectal adenoma and the candidate gene NAT2 in relation to smoking history.     

Inferring an Augmented Bayesian Network to Confront a Complex Quantitative Microbial Risk Assessment Model with Durability Studies: Application to Bacillus Cereus on a Courgette Purée Production Chain

The Monte Carlo (MC) simulation approach is traditionally used in food safety risk assessment to study quantitative microbial risk assessment (QMRA) models. When experimental data are available, performing Bayesian inference is a good alternative approach that allows backward calculation in a stochastic QMRA model to update the experts’ knowledge about the microbial dynamics of a given food‐borne pathogen. In this article, we propose a complex example where Bayesian inference is applied to a high‐dimensional second‐order QMRA model. The case study is a farm‐to‐fork QMRA model considering genetic diversity of Bacillus cereus in a cooked, pasteurized, and chilled courgette purée. Experimental data are Bacillus cereus concentrations measured in packages of courgette purées stored at different time‐temperature profiles after pasteurization. To perform a Bayesian inference, we first built an augmented Bayesian network by linking a second‐order QMRA model to the available contamination data. We then ran a Markov chain Monte Carlo (MCMC) algorithm to update all the unknown concentrations and unknown quantities of the augmented model. About 25% of the prior beliefs are strongly updated, leading to a reduction in uncertainty. Some updates interestingly question the QMRA model.     

A Class of Markov Chain Models for Average Run Length Computations for Autocorrelated Processes

The average run length (ARL) of conventional control charts is typically computed assuming temporal independence. However, this assumption is frequently violated in practical applications. Alternative ARL computations have often been conducted via time consuming and yet not necessarily very accurate simulations. In this article, we develop a class of Markov chain models for evaluating the run length performance of traditional control charts for autocorrelated processes. We show extensions from the univariate AR(1) model to the general multivariate VARMA(p, q) time series. The results of the proposed method are highly comparable to those of simulations and with significantly less computational overhead.     

基于三种模型的江苏沿海人口比较测算

本文运用马尔萨斯人口模型、Logistic 增长模型和线性函数三种理论模型，结合江苏省第六次全国人口普查研究统计数据对江苏沿海2012—2020年的人口发展规模进行预测，研究江苏沿海人口演变特征，预测未来人口增长率及规模，结果显示：基于三种模型运行的五种结果差异较大，根据江苏沿海的实际情况和前人研究结果，确认江苏沿海到2020年的总人口为2125.47万人。研究结果以期能为江苏沿海的科学发展提供参考依据。     

Smoothing Spline Estimation for Partially Linear Single-index Models

In this article, the partially linear single-index models are discussed based on smoothing spline and average derivative estimation method. This proposed technique consists of two stages: one is to estimate the vector parameter in the linear part using the smoothing cubic spline method, simultaneously, obtaining the estimator of unknown single-index function; the other is to estimate the single-index coefficients in the single-index part by the using average derivative estimator procedure. Some simulated and real examples are presented to illustrate the performance of this method.     

Optimal designs to improve the power of multiresponse hypothesis tests

In this paper we develop a design criterion to improve the power of the multiresponse hypothesis tests, and give a sequential procedure for generating designs which satisfy the criterion.