Generalized Additive Modelling of Mixed Distribution Markov Models with Application to Melbourne's Rainfall

The paper considers the modelling of time series using a generalized additive model with first-order Markov structure and mixed transition density having a discrete component at zero and a continuous component with positive sample space. Such models have application, for example, in modelling daily occurrence and intensity of rainfall, and in modelling numbers and sizes of insurance claims. The paper shows how these methods extend the usual sinusoidal seasonal assumption in standard chain-dependent models by assuming a general smooth pattern of occurrence and intensity over time. These models can be fitted using standard statistical software. The methods of Grunwald & Jones (2000) can be used to combine these separate occurrence and intensity models into a single model for amount. The models are used to investigate the relationship between the Southern Oscillation Index and Melbourne's rainfall, illustrated with 36 years of rainfall data from Melbourne, Australia.     

The variability of rainfall acidity

Two years of rainfall acidity data for the eastern United States were analyzed. The data consist of rainfall-event pH measurements from a nine station monitoring network. A spatio-temporal stochastic model, including deterministic components for seasonal variation and rainfall washout, and stochastic components for spatial, temporal, and measurement variation, was fitted to the data. The fitted autocorrelation structure from this model was used, in the process known as Kriging, to obtain BLUE contour maps of seasonal and rainfall adjusted yearly average pH over the monitoring region.     

Venezuelan Rainfall Data Analysed by Using a Bayesian Space–time Model

We consider a set of data from 80 stations in the Venezuelan state of Guárico consisting of accumulated monthly rainfall in a time span of 16 years. The problem of modelling rainfall accumulated over fixed periods of time and recorded at meteorological stations at different sites is studied by using a model based on the assumption that the data follow a truncated and transformed multivariate normal distribution. The spatial correlation is modelled by using an exponentially decreasing correlation function and an interpolating surface for the means. Missing data and dry periods are handled within a Markov chain Monte Carlo framework using latent variables. We estimate the amount of rainfall as well as the probability of a dry period by using the predictive density of the data. We considered a model based on a full second-degree polynomial over the spatial co-ordinates as well as the first two Fourier harmonics to describe the variability during the year. Predictive inferences on the data show very realistic results, capturing the typical rainfall variability in time and space for that region. Important extensions of the model are also discussed.     

Mixed probability models for dry and wet spells

When identifying the best model for representing the behavior of rainfall distribution based on a sequence of dry (wet) days, focus is usually given on the fitted model with the least number of estimated parameters. If the model with lesser number of parameters is found not adequate for describing a particular data distribution, the model with a higher number of parameters is recommended. Based on several probability models developed by previous researchers in this field, we propose five types of mixed probability models as the alternative to describe the distribution of dry (wet) spells for daily rainfall events. The mixed probability models comprise of the combination of log series distribution with three other types of models, which are Poisson distribution (MLPD), truncated Poisson distribution (MLTPD), and geometric distribution (MLGD). In addition, the combination of the two log series distributions (MLSD) and the mixed geometric with the truncated Poisson distribution (MGTPD) are also introduced as the alternative models. Daily rainfall data from 14 selected rainfall stations in Peninsular Malaysia for the periods of 1975 to 2004 were used in this present study. When selecting the best probability model to describe the observed distribution of dry (wet) spells, the Akaike’s Information Criterion (AIC) was considered. The results revealed that MLGD was the best probability model to represent the distribution of dry spells over the Peninsular.     

Spatial modeling of extreme rainfall in northeast Thailand

It is well recognized that the generalized extreme value (GEV) distribution is widely used for any extreme events. This notion is based on the study of discrete choice behavior; however, there is a limit for predicting the distribution at ungauged sites. Hence, there have been studies on spatial dependence within extreme events in continuous space using recorded observations. We model the annual maximum daily rainfall data consisting of 25 locations for the period from 1982 to 2013. The spatial GEV model that is established under observations is assumed to be mutually independent because there is no spatial dependency between the stations. Furthermore, we divide the region into two regions for a better model fit and identify the best model for each region. We show that the regional spatial GEV model reflects the spatial pattern well compared with the spatial GEV model over the entire region as the local GEV distribution. The advantage of spatial extreme modeling is that more robust return levels and some indices of extreme rainfall can be obtained for observed stations as well as for locations without observed data. Thus, the model helps to determine the effects and assessment of vulnerability due to heavy rainfall in northeast Thailand.     

Model Selection For A Family Of Discrete Failure Time Distributions

In recent years, a large number of new discrete distributions have appeared in the literature. However, flexible discrete models which, at the same time, allow for easy statistical inference, are still an exception. This paper makes a detailed analysis of a family of discrete failure time distributions which meets both requirements. It examines the maximum likelihood estimation of the unknown parameters and presents a goodness-of-fit test for this model. The test is used for the selection of an appropriate model for datasets of frequencies of the duration of atmospheric circulation patterns.     

A note comparing component-slope,Scheffé and Cox parameterizations of the linear mixture experiment model

A mixture experiment involves combining two or more components in various proportions and collecting data on one or more responses. A linear mixture model may adequately represent the relationship between a response and mixture component proportions and be useful in screening the mixture components. The Scheffé and Cox parameterizations of the linear mixture model are commonly used for analyzing mixture experiment data. With the Scheffé parameterization, the fitted coefficient for a component is the predicted response at that pure component (i.e. single-component mixture). With the Cox parameterization, the fitted coefficient for a mixture component is the predicted difference in response at that pure component and at a pre-specified reference composition. This article presents a new component-slope parameterization, in which the fitted coefficient for a mixture component is the predicted slope of the linear response surface along the direction determined by that pure component and at a pre-specified reference composition. The component-slope, Scheffé, and Cox parameterizations of the linear mixture model are compared and their advantages and disadvantages are discussed.     

Functional regression concurrent model with spatially correlated errors: application to rainfall ground validation

In this paper, we give an extension of the functional regression concurrent model to the case of spatially correlated errors. We propose estimating the spatial correlation structure by using functional geostatistics. The estimation of the regression parameters is carried out by feasible generalized least squares. This modeling approach is motivated by the problem of validating rainfall data retrieved from satellite sensors. In this sense, we use the methodology to study the relationship between satellite and ground rainfall time series recorded in 82 weather stations from Department of Valle del Cauca, Colombia. The model obtained allows predicting pentadal rainfall curves in many sites of the region of interest by using as input the satellite information. A residual analysis shows a good performance of the methodology proposed.     

Parametric relative survival regression using generalized linear models with application to Hodgkin's lymphoma

Changes in survival rates during 1940–1992 for patients with Hodgkin's disease are studied by using population-based data. The aim of the analysis is to identify when the breakthrough in clinical trials of chemotherapy treatments started to increase population survival rates, and to find how long it took for the increase to level off, indicating that the full population effect of the breakthrough had been realized. A Weibull relative survival model is used because the model parameters are easily interpretable when assessing the effect of advances in clinical trials. However, the methods apply to any relative survival model that falls within the generalized linear models framework. The model is fitted by using modifications of existing software (SAS, GLIM) and profile likelihood methods. The results are similar to those from a cause-specific analysis of the data by Feuer and co-workers. Survival started to improve around the time that a major chemotherapy breakthrough (nitrogen mustard, Oncovin, prednisone and procarbazine) was publicized in the mid 1960s but did not level off for 11 years. For the analysis of data where the cause of death is obtained from death certificates, the relative survival approach has the advantage of providing the necessary adjustment for expected mortality from causes other than the disease without requiring information on the causes of death.     

A latent Gaussian Markov random-field model for spatiotemporal rainfall disaggregation

Summary. Rainfall data are often collected at coarser spatial scales than required for input into hydrology and agricultural models. We therefore describe a spatiotemporal model which allows multiple imputation of rainfall at fine spatial resolutions, with a realistic dependence structure in both space and time and with the total rainfall at the coarse scale consistent with that observed. The method involves the transformation of the fine scale rainfall to a thresholded Gaussian process which we model as a Gaussian Markov random field. Gibbs sampling is then used to generate realizations of rainfall efficiently at the fine scale. Results compare favourably with previous, less elegant methods.     

Factors influencing Soay sheep survival

We present a survival analysis of Soay sheep mark recapture and recovery data. Unlike previous conditional analyses, it is not necessary to assume equality of recovery and recapture probabilities; instead these are estimated by maximum likelihood. Male and female sheep are treated separately, with the higher numbers and survival probabilities of the females resulting in a more complex model than that used for the males. In both cases, however, age and time aspects need to be included and there is a strong indication of a reduction in survival for sheep aged 7 years or more. Time variation in survival is related to the size of the population and selected weather variables, by using logistic regression. The size of the population significantly affects the survival probabilities of male and female lambs, and of female sheep aged 7 or more years. March rainfall and a measure of the North Atlantic oscillation are found to influence survival significantly for all age groups considered, for both males and females. Either of these weather variables can be used in a model. Several phenotypic and genotypic individual covariates are also fitted. The only covariate which is found to influence survival significantly is the type of horn of first-year female sheep. There is a substantial variation in the recovery probabilities over time, reflecting in part the increased effort when a population crash was expected. The goodness of fit of the model is checked by using graphical procedures.     

A note on diagnostic checking of the double autoregressive model

The double autoregressive model finds its use in the modelling of conditional heteroscedasticity of time series data. In view of its growing popularity, the goodness-of-fit of the model is examined. The asymptotic distributions of the residual and squared residual autocorrelations are derived. Two test statistics are then constructed which can be used to measure the adequacy of the conditional mean and conditional variance components of a fitted model. Our goodness-of-fit tests out-perform other benchmark tests such as the Ljung–Box test in simulation studies. To illustrate the testing procedure, the model is fitted to the weekly log-return series of the Hang Seng Index.     

Multivariate extreme-value distributions with applications to environmental data

Some parametric families of multivariate extreme-value distributions have been proposed in recent years; several additional parametric families are derived here. The parametric models are fitted, using numerical maximum likelihood, to some environmental multivariate extreme data sets consisting of extreme concentrations of a pollutant at several monitoring stations in a region. Some multivariate nonnormal data analysis techniques are proposed to aid in the likelihood analysis. The new models, together with previous models, appear to be adequate for inferences in that they cover a wide range of possible dependence patterns.     

Generalized linear spatial models in epidemiology: A case study of zoonotic cutaneous leishmaniasis in Tunisia

Generalized linear spatial models (GLSM) are used here to study spatial characters of zoonotic cutaneous leishmaniasis (ZCL) in Tunisia. The response variable stands for the number of affected by district during the period 2001–2002. The model covariates are: climates (temperature and rainfall), humidity and surrounding vegetation status. As the environmental and weather data are not available for all the studied districts, Kriging based on linear interpolation was used to estimate the missing data. To account for unexplained spatial variation in the model, we include a stationary Gaussian process S with a powered exponential spatial correlation function. Moran coefficient, DIC criterion and residuals variograms are used to show the high goodness-of-fit of the GLSM. When compared with the statistical tools used in the previous ZCL studies, the optimal GLSM found here yields a better assessment of the impact of the risk factors, a better prediction of ZCL evolution and a better comprehension of the disease transmission. The statistical results show the progressive increase in the number of affected in zones with high temperature, low rainfall and high surrounding vegetation index. Relative humidity does not seem to affect the distribution of the disease in Tunisia. The results of the statistical analyses stress the important risk of misleading epidemiological conclusions when non-spatial models are used to analyse spatially structured data.     

A Bayesian kriged Kalman model for short-term forecasting of air pollution levels

Summary.  Short-term forecasts of air pollution levels in big cities are now reported in news-papers and other media outlets. Studies indicate that even short-term exposure to high levels of an air pollutant called atmospheric particulate matter can lead to long-term health effects. Data are typically observed at fixed monitoring stations throughout a study region of interest at different time points. Statistical spatiotemporal models are appropriate for modelling these data. We consider short-term forecasting of these spatiotemporal processes by using a Bayesian kriged Kalman filtering model. The spatial prediction surface of the model is built by using the well-known method of kriging for optimum spatial prediction and the temporal effects are analysed by using the models underlying the Kalman filtering method. The full Bayesian model is implemented by using Markov chain Monte Carlo techniques which enable us to obtain the optimal Bayesian forecasts in time and space. A new cross-validation method based on the Mahalanobis distance between the forecasts and observed data is also developed to assess the forecasting performance of the model implemented.     

中国省域人口总数修订

人口统计数据是制定宏观经济政策和规划人口发展的基础。中国非普查年度人口总数常常与普查年度人口总数差异较大,而关于中国非普查年度人口统计数据修订的研究较少。建立省域人口增长率的固定效应模型,描述中国各省域人口总数的变化趋势,基于拟合模型和普查数据,提出中国省域人口总数修订方法,并对2001年至2009年中国省域人口总数进行修订,结果表明:2001年至2009年中国人口总数修订值与《中国统计年鉴》公布的人口总数之间存在差异,差异率低于0.5%;省域人口总数修订值变化相对平稳,在年度之间的变化幅度相对小,能够更合理反映各省域人口总数的发展趋势。     

Constrained clustering of irregularly sampled spatial data

Many spatial data such as those in climatology or environmental monitoring are collected over irregular geographical locations. Furthermore, it is common to have multivariate observations at each location. We propose a method of segmentation of a region of interest based on such data that can be carried out in two steps: (1) clustering or classification of irregularly sample points and (2) segmentation of the region based on the classified points.

We develop a spatially-constrained clustering algorithm for segmentation of the sample points by incorporating a geographical-constraint into the standard clustering methods. Both hierarchical and nonhierarchical methods are considered. The latter is a modification of the seeded region growing method known in image analysis. Both algorithms work on a suitable neighbourhood structure, which can for example be defined by the Delaunay triangulation of the sample points. The number of clusters is estimated by testing the significance of successive change in the within-cluster sum-of-squares relative to a null permutation distribution. The methodology is validated on simulated data and used in construction of a climatology map of Ireland based on meteorological data of daily rainfall records from 1294 stations over the period of 37 years.     

GENERALIZED EXTREME VALUE ADDITIVE MODEL ANALYSIS VIA MEAN FIELD VARIATIONAL BAYES

We develop Mean Field Variational Bayes methodology for fast approximate inference in Bayesian Generalized Extreme Value additive model analysis. Such models are useful for flexibly assessing the impact of continuous predictor variables on sample extremes. The new methodology allows large Bayesian models to be fitted and assessed without the significant computing costs of Markov Chain Monte Carlo methods. We illustrate our new methodology with maximum rainfall data from the Sydney, Australia, hinterland. Comparisons are made between the Mean Field Variational Bayes and Markov Chain Monte Carlo approaches.     

A Bayesian compartmental model for the evaluation of 1,3-butadiene metabolism

Summary. We propose a Bayesian model for physiologically based pharmacokinetics of 1,3-butadiene (BD). BD is classified as a suspected human carcinogen and exposure to it is common, especially through cigarette smoke as well as in urban settings. The main aim of the methodology and analysis that are presented here is to quantify variability in the rates of BD metabolism by human subjects. A three-compartmental model is described, together with informative prior distributions for the population parameters, all of which represent real physiological variables. The model is described in detail along with the meanings and interpretations of the associated parameters. A four-compartment model is also given for comparison. Markov chain Monte Carlo methods are described for fitting the model proposed. The model is fitted to toxicokinetic data obtained from 133 healthy subjects (males and females) from the four major racial groups in the USA, with ages ranging from 19 to 62 years. Subjects were exposed to 2 parts per million of BD for 20 min through a face mask by using a computer-controlled exposure and respiratory monitoring system. Stratification by ethnic group results in major changes in the physiological parameters. Sex and age were also tested but not found to have a significant effect.     

Stochastic modelling of regional annual rainfall anomalies in East Africa

ARIMA (p, d, q) models were fitted to areal annual rainfall of two homogeneous regions in East Africa with rainfall records extending between the period 1922–80. The areal estimates of the regional rainfall were derived from the time series of the first eigenvector, which was significantly dominant at each of the two regions. The first eigenvector accounted for about 80% of the total rainfall variance in each region.

The class of ARIMA (p, d, q) models which best fitted the areal indices of relative wetness/dryness were the A R M A (3, 1) models. Tests of forecasting skill however indicated low skill in the forecasts given by these models. In all cases the models accounted for less than 50% of the total variance.

Spectral analysis of the indices time series indicated dominant quasi-periodic fluctuations around 2.2–2.8 years, 3–3.7 years, 5–6 years and 10–13 years. These spectral bands however accounted for very low proportion of the total rainfall variance.