Bayesian updating of atmospheric dispersion after a nuclear accident

Summary.  We consider a Bayesian forecasting system to predict the dispersal of contamination on a large scale grid in the event of an accidental release of radioactivity. The statistical model is built on a physical model for atmospheric dispersion and transport called MATCH. Our spatiotemporal model is a dynamic linear model where the state parameters are the (essentially, deterministic) predictions of MATCH; the distributions of these are updated sequentially in the light of monitoring data. One of the distinguishing features of the model is that the number of these parameters is very large (typically several hundreds of thousands) and we discuss practical issues arising in its implementation as a realtime model. Our procedures have been checked against a variational approach which is used widely in the atmospheric sciences. The results of the model are applied to test data from a tracer experiment.     

Temporal pattern in number of staff on sick leave: the effect of downsizing

Summary.  The pattern of absenteeism in the downsizing process of companies is a topic in focus in economics and social science. A general question is whether employees who are frequently absent are more likely to be selected to be laid off or in contrast whether employees to be dismissed are more likely to be absent for the remaining time of their working contract. We pursue an empirical and microeconomic investigation of these theses. We analyse longitudinal data that were collected in a German company over several years. We fit a semiparametric transition model based on a mixture Poisson distribution for the days of absenteeism per month. Prediction intervals are considered and the primary focus is on the period of downsizing. The data reveal clear evidence for the hypothesis that employees who are to be laid off are more frequently absent before leaving the company. Interestingly, though, no clear evidence is seen that employees being selected to leave the company are those with a bad absenteeism profile.     

Causal inference with generalized structural mean models

Summary.  We estimate cause–effect relationships in empirical research where exposures are not completely controlled, as in observational studies or with patient non-compliance and self-selected treatment switches in randomized clinical trials. Additive and multiplicative structural mean models have proved useful for this but suffer from the classical limitations of linear and log-linear models when accommodating binary data. We propose the generalized structural mean model to overcome these limitations. This is a semiparametric two-stage model which extends the structural mean model to handle non-linear average exposure effects. The first-stage structural model describes the causal effect of received exposure by contrasting the means of observed and potential exposure-free outcomes in exposed subsets of the population. For identification of the structural parameters, a second stage 'nuisance' model is introduced. This takes the form of a classical association model for expected outcomes given observed exposure. Under the model, we derive estimating equations which yield consistent, asymptotically normal and efficient estimators of the structural effects. We examine their robustness to model misspecification and construct robust estimators in the absence of any exposure effect. The double-logistic structural mean model is developed in more detail to estimate the effect of observed exposure on the success of treatment in a randomized controlled blood pressure reduction trial with self-selected non-compliance.     

言语生成及理解的心理模型

言语生成 (speechproduction)和言语理解 (speechunderstanding)是语言交际中十分复杂的心理认知过程 ,也是心理语言学研究中的一个重要内容。本文拟对有关言语生成和理解的几种心理模型予以讨论 ,探讨言语生成及理解的过程和实质。     

Decision analysis for practical negotiation application

The family of decision analysis techniques can be applied effectively to support practical negotiators in international settings. These techniques are most appropriate in support of the prenegotiation phase, when parties are diagnosing the situation, assessing their own plans and strategies, and evaluating likely reactions and outcomes. The paper identifies how these approaches have and can be used to assist negotiation practitioners, offers a rationale for the application of decision analytic approaches in terms of the particular analytical requirements of the prenegotiation period, suggests how these process-oriented tools can be integrated with substantive tools, and discusses ways in which these tools can be presented and delivered to practitioners in a practical and confidence-building manner.     

密度分布非均匀的反馈控制Logistic单种群模型平衡态的稳定性

本文对种群密度在非均匀分布情形下,考虑了具反馈控制的滞后 Logistic 生态模型平衡位置的稳定性;分别给出了在常时滞和弱连续时滞以及强连续时滞情况下的稳定性条件;其结果是对 Gopalsamy 在密度均匀分布情形下相应结果的推广.     

WEIGHTED SUMS OF NEGATIVELY ASSOCIATED RANDOM VARIABLES

In this paper, we establish strong laws for weighted sums of negatively associated (NA) random variables which have a higher‐order moment condition. Some results of Bai Z.D. & Cheng P.E. (2000) [Marcinkiewicz strong laws for linear statistics. Statist. and Probab. Lett. 43, 105–112,] and Sung S.K. (2001) [Strong laws for weighted sums of i.i.d. random variables, Statist. and Probab. Lett. 52, 413–419] are sharpened and extended from the independent identically distributed case to the NA setting. Also, one of the results of Li D.L. et al. (1995) [Complete convergence and almost sure convergence of weighted sums of random variables. J. Theoret. Probab. 8, 49–76,] is complemented and extended.     

Residual analysis for spatial point processes (with discussion)

Summary.  We define residuals for point process models fitted to spatial point pattern data, and we propose diagnostic plots based on them. The residuals apply to any point process model that has a conditional intensity; the model may exhibit spatial heterogeneity, interpoint interaction and dependence on spatial covariates. Some existing ad hoc methods for model checking (quadrat counts, scan statistic, kernel smoothed intensity and Berman's diagnostic) are recovered as special cases. Diagnostic tools are developed systematically, by using an analogy between our spatial residuals and the usual residuals for (non-spatial) generalized linear models. The conditional intensity λ plays the role of the mean response. This makes it possible to adapt existing knowledge about model validation for generalized linear models to the spatial point process context, giving recommendations for diagnostic plots. A plot of smoothed residuals against spatial location, or against a spatial covariate, is effective in diagnosing spatial trend or co-variate effects. Q – Q -plots of the residuals are effective in diagnosing interpoint interaction.     

An individual measure of relative survival

Summary.  Relative survival techniques are used to compare survival experience in a study cohort with that expected if background population rates apply. The techniques are especially useful when cause-specific death information is not accurate or not available as they provide a measure of excess mortality in a group of patients with a certain disease. Whereas these methods are based on group comparisons, we present here a transformation approach which instead gives for each individual an outcome measure relative to the appropriate background population. The new outcome measure is easily interpreted and can be analysed by using standard survival analysis techniques. It provides additional information on relative survival and gives new options in regression analysis. For example, one can estimate the proportion of patients who survived longer than a given percentile of the respective general population or compare survival experience of individuals while accounting for the population differences. The regression models for the new outcome measure are different from existing models, thus providing new possibilities in analysing relative survival data. One distinctive feature of our approach is that we adjust for expected survival before modelling. The paper is motivated by a study into the survival of patients after acute myocardial infarction.     

On Block Updating in Markov Random Field Models for Disease Mapping

Gaussian Markov random field (GMRF) models are commonly used to model spatial correlation in disease mapping applications. For Bayesian inference by MCMC, so far mainly single-site updating algorithms have been considered. However, convergence and mixing properties of such algorithms can be extremely poor due to strong dependencies of parameters in the posterior distribution. In this paper, we propose various block sampling algorithms in order to improve the MCMC performance. The methodology is rather general, allows for non-standard full conditionals, and can be applied in a modular fashion in a large number of different scenarios. For illustration we consider three different applications: two formulations for spatial modelling of a single disease (with and without additional unstructured parameters respectively), and one formulation for the joint analysis of two diseases. The results indicate that the largest benefits are obtained if parameters and the corresponding hyperparameter are updated jointly in one large block. Implementation of such block algorithms is relatively easy using methods for fast sampling of Gaussian Markov random fields ( Rue, 2001 ). By comparison, Monte Carlo estimates based on single-site updating can be rather misleading, even for very long runs. Our results may have wider relevance for efficient MCMC simulation in hierarchical models with Markov random field components.