Relative risk analysis of dengue cases using convolution extended into spatio-temporal model

Dengue Hemmorage Fever (DHF) cases have become a serious problem every year in tropical countries such as Indonesia. Understanding the dynamic spread of the disease is essential in order to find an effective strategy in controlling its spread. In this study, a convolution (Poisson-lognormal) model that integrates both uncorrelated and correlated random effects was developed. A spatial–temporal convolution model to accomodate both spatial and temporal variations of the disease spread dynamics was considered. The model was applied to the DHF cases in the city of Kendari, Indonesia. DHF data for 10 districts during the period 2007–2010 were collected from the health services. The data of rainfall and population density were obtained from the local offices in Kendari. The numerical experiments indicated that both the rainfall and the population density played an important role in the increasing DHF cases in the city of Kendari. The result suggested that DHF cases mostly occured in January, the wet session with high rainfall, and in Kadia, the densest district in the city. As people in the city have high mobility while dengue mosquitoes tend to stay localized in their area, the best intervention is in January and in the district of Kadia.     

Extended Hazard Regression Model for Reliability and Survival Analysis

We propose an extended hazard regression model which allows the spread parameter to be dependent on covariates. This allows a broad class of models which includes the most common hazard models, such as the proportional hazards model, the accelerated failure time model and a proportional hazards/accelerated failure time hybrid model with constant spread parameter. Simulations based on sub-classes of this model suggest that maximum likelihood performs well even when only small or moderate-size data sets are available and the censoring pattern is heavy. The methodology provides a broad framework for analysis of reliability and survival data. Two numerical examples illustrate the results.     

Acceptance control charts for non-normal data

Control charts are one of the most important methods in industrial process control. The acceptance control chart is generally applied in situations when an X¯ chart is used to control the fraction of conforming units produced by the process and where 6-sigma spread of the process is smaller than the spread in the specification limits. Traditionally, when designing control charts, one usually assumes that the data or measurements are normally distributed. However, this assumption may not be true in some processes. In this paper, we use the Burr distribution, which is employed to represent various non-normal distributions, to determine the appropriate control limits or sample size for the acceptance control chart under non-normality. Some numerical examples are given for illustration. From the presented examples, ignoring the effect of non-normality in the data leads to a higher type I or type II error probability.     

ANALYSIS OF DATA FROM A SINGLE EPIDEMIC1

A model with nonrandom latent and infectious periods is suggested for epidemics in a large community. This permits a relatively complete statistical analysis of data from the spread of a single epidemic. An attractive feature of such models is the possibility of exploring how the rate of spread of the disease depends on the number of susceptibles and infectives. An application to smallpox data is included.     

基于复杂网络的甲型H1N1流感病毒传播研究

基于复杂网络理论,研究了甲型HIN1流感病毒的传播机制。根据实证数据分析,建立了一个甲型HIN1流感病毒传播的小世界网络模型。通过对病毒传播中个体的接触程度、传染概率、潜伏期、增加个体免疫四个因素的研究,分析了当前应对甲型HIN1流感病毒防控措施的有效性。     

Perfect simulation for Reed-Frost epidemic models

The Reed-Frost epidemic model is a simple stochastic process with parameter q that describes the spread of an infectious disease among a closed population. Given data on the final outcome of an epidemic, it is possible to perform Bayesian inference for q using a simple Gibbs sampler algorithm. In this paper it is illustrated that by choosing latent variables appropriately, certain monotonicity properties hold which facilitate the use of a perfect simulation algorithm. The methods are applied to real data.     

Epidemics in heterogeneous communities: estimation of R0 and secure vaccination coverage

A stochastic multitype model for the spread of an infectious disease in a community of heterogeneous individuals is analysed. In particular, estimates of R ₀ (the basic reproduction number) and the critical vaccination coverage are derived, where estimation is based on final size data of an outbreak in the community. It is shown that these key parameters cannot be estimated consistently from data; only upper and lower bounds can be estimated. Confidence regions for the upper bounds are derived, thus giving conservative estimates of R ₀ and the fractions necessary to vaccinate.     

Statistical analysis of an endemic disease from a capture–recapture experiment

There are a number of statistical techniques for analysing epidemic outbreaks. However, many diseases are endemic within populations and the analysis of such diseases is complicated by changing population demography. Motivated by the spread of cowpox amongst rodent populations, a combined mathematical model for population and disease dynamics is introduced. A Markov chain Monte Carlo algorithm is then constructed to make statistical inference for the model based on data being obtained from a capture–recapture experiment. The statistical analysis is used to identify the key elements in the spread of the cowpox virus.     

Statement of Withdrawal: Statistical analysis of an endemic disease from a capture–recapture experiment

There are a number of statistical techniques for analysing epidemic outbreaks. However, many diseases are endemic within populations and the analysis of such diseases are complicated by changing population demography. Motivated by the spread of cowpox among rodent populations, a combined mathematical model for population and disease dynamics is introduced. An MCMC algorithm is then constructed to make statistical inference for the model based on data being obtained from a capture–recapture experiment. The statistical analysis is used to identify the key elements in the spread of the cowpox virus.     

Predicting recessions using trends in the yield spread

The yield spread, measured as the difference between long- and short-term interest rates, is widely regarded as one of the strongest predictors of economic recessions. In this paper, we propose an enhanced recession prediction model that incorporates trends in the value of the yield spread. We expect our model to generate stronger recession signals because a steadily declining value of the yield spread typically indicates growing pessimism associated with the reduced future business activity. We capture trends in the yield spread by considering both the level of the yield spread at a lag of 12 months as well as its value at each of the previous two quarters leading up to the forecast origin, and we evaluate its predictive abilities using both logit and artificial neural network models. Our results indicate that models incorporating information from the time series of the yield spread correctly predict future recession periods much better than models only considering the spread value as of the forecast origin. Furthermore, the results are strongest for our artificial neural network model and logistic regression model that includes interaction terms, which we confirm using both a blocked cross-validation technique as well as an expanding estimation window approach.     

收益率曲面预测及其在信用债投资组合管理中的应用

收益率曲线是信用债投资者的重要参考。在市场套利活动的作用下,跨期限和信用等级的债券收益率具有很强的内生联系,其变动规律具有整体性和连续性。以往研究将信用债收益率曲线拆分成无风险利率曲线和信用利差两个部分,前者关注期限的影响,后者关注信用等级的影响,但很少同时考虑两者的共同影响。本文在收益率曲线的基础上增加信用等级维度,将AAA+级到AA级收益率曲线视为一个相互关联的整体,定义为信用债的收益率曲面(Yield Surface)。相比收益率曲线,收益率曲面包含了跨等级的系统性预测信息,通过预测收益率曲面能够构建具有较高收益的投资管理策略。本文以中期票据市场为例,采用卷积神经网络模型预测1周后的收益率曲面,在此基础上计算债券的预测理论价格,发现该价格对未来交易价格的预测有显著作用。基于预测价差,本文提出了信用债投资管 理策略,应用该策略的投资组合能够获得显著的正收益。业绩归因分析发现,该策略取得的收益同时来 自投资组合在信用风险和久期风险上的暴露,预测价差可以刻画债券市场的“风险前沿”。本文采用了中债估值价格进行稳健性检验,主要结论均保持一致,具有较强的稳健性。     

An overview of model-robust regression

The proper combination of parametric and nonparametric regression procedures can improve upon the shortcomings of each when used individually. Considered is the situation where the researcher has an idea of which parametric model should explain the behavior of the data, but this model is not adequate throughout the entire range of the data. An extension of partial linear regression and two other methods of model-robust regression are developed and compared in this context. The model-robust procedures each involve the proportional mixing of a parametric fit to the data and a nonparametric fit to either the data or residuals. The emphasis of this work is on fitting in the small-sample situation, where nonparametric regression alone has well-known inadequacies. Performance is based on bias and variance considerations, and theoretical mean squared error formulas are developed for each procedure. An example is given that uses generated data from an underlying model with defined misspecification to provide graphical comparisons of the fits and to show the theoretical benefits of the model-robust procedures. Simulation results are presented which establish the accuracy of the theoretical formulas and illustrate the potential benefits of the model-robust procedures. Simulations are also used to illustrate the advantageous properties of a data-driven selector developed in this work for choosing the smoothing and mixing parameters. It is seen that the model-robust procedures (the final proposed method, in particular) give much improved fits over the individual parametric and nonparametric fits.     

A permutation test for the spread of three-dimensional rotation data

The permutation test is a nonparametric test that can be used to compare measures of spread for two data sets, but is yet to be explored in the context of three-dimensional rotation data. A permutation test for such data is developed and the statistical power of this test is considered under various conditions. The test is then used in a brief application comparing movement around the calcaneocuboid joint for a human, chimpanzee, and baboon.     

Spatial measurement error in infectious disease models

Individual-level models (ILMs) for infectious disease can be used to model disease spread between individuals while taking into account important covariates. One important covariate in determining the risk of infection transfer can be spatial location. At the same time, measurement error is a concern in many areas of statistical analysis, and infectious disease modelling is no exception. In this paper, we are concerned with the issue of measurement error in the recorded location of individuals when using a simple spatial ILM to model the spread of disease within a population. An ILM that incorporates spatial location random effects is introduced within a hierarchical Bayesian framework. This model is tested upon both simulated data and data from the UK 2001 foot-and-mouth disease epidemic. The ability of the model to successfully identify both the spatial infection kernel and the basic reproduction number (R ₀) of the disease is tested.     

The Design of and R Control Charts for Skew Normal Distributed Data

The Shewhart-type control chart is traditionally developed under the normality assumption. In practice, however, this assumption may not hold. Because the skew normal distribution represents a broad distribution class and is more flexible than is the normal distribution, we propose two new control charts to monitor process mean and spread for skew normal distributed data. Moreover, to facilitate practical implementation, tables of charting constants are provided. We conducted simulation studies to compare the false alarm rates, and the results show that new proposed charts perform better than others as skewness increases. Finally, an illustrative example is provided.     

Markov Chain Monte Carlo Methods for Fitting Spatiotemporal Stochastic Models in Plant Epidemiology

Strategies for controlling plant epidemics are investigated by fitting continuous time spatiotemporal stochastic models to data consisting of maps of disease incidence observed at discrete times. Markov chain Monte Carlo methods are used for fitting two such models to data describing the spread of citrus tristeza virus (CTV) in an orchard. The approach overcomes some of the difficulties encountered when fitting stochastic models to infrequent observations of a continuous process. The results of the analysis cast doubt on the effectiveness of a strategy identified from a previous spatial analysis of the CTV data. Extensions of the approaches to more general models and other problems are also considered.     

An influence method for outlier detection applied to time series traffic data

The work of Chernick et al. (1982) is extended to form a quantitative outlier detection statistic for use with time series data. The statistic is formed from the squared elements of the influence function matrix, where each element of the matrix gives the influence on the theoretical autocorrelation function at lag k (pk) of a pair of obser vations at time lag k. The approximate first four moments for the statistic are derived and, by fitting Johnson curves to these theoretical moments, critical points are also produced. The statistic is also used to form the basis of an adjustment procedure to treat outliers or estimate missing values in the time series. The nuclear power data of Chernick et al. and the traffic count data of the Department of Transport are used for practical illustration.     

Predicting infectious disease outbreak risk via migratory waterfowl vectors

The spread of an emerging infectious disease is a major public health threat. Given the uncertainties associated with vector-borne diseases, in terms of vector dynamics and disease transmission, it is critical to develop statistical models that address how and when such an infectious disease could spread throughout a region such as the USA. This paper considers a spatio-temporal statistical model for how an infectious disease could be carried into the USA by migratory waterfowl vectors during their seasonal migration and, ultimately, the risk of transmission of such a disease to domestic fowl. Modeling spatio-temporal data of this type is inherently difficult given the uncertainty associated with observations, complexity of the dynamics, high dimensionality of the underlying process, and the presence of excessive zeros. In particular, the spatio-temporal dynamics of the waterfowl migration are developed by way of a two-tiered functional temporal and spatial dimension reduction procedure that captures spatial and seasonal trends, as well as regional dynamics. Furthermore, the model relates the migration to a population of poultry farms that are known to be susceptible to such diseases, and is one of the possible avenues toward transmission to domestic poultry and humans. The result is a predictive distribution of those counties containing poultry farms that are at the greatest risk of having the infectious disease infiltrate their flocks assuming that the migratory population was infected. The model naturally fits into the hierarchical Bayesian framework.     

A model for the spread and control of citrus tristeza virus disease

Citrus Tristeza Virus (CTV) causes a most destructive citrus disease in many parts of the world. The strategy for controlling the disease in Israel is based on the eradication of virus-infected trees, detected by an immunological method. A model for the spread and control of CTV is described. The disease spread depends on an infection rate parameter and on a detection coefficient. Both parameters vary with time. Simulation of the model under different conditions assists us to make decisions regarding the control of CTV disease.