Spatio‐temporal modelling of disease mapping of rates

This paper studies generalized linear mixed models (GLMMs) for the analysis of geographic and temporal variability of disease rates. This class of models adopts spatially correlated random effects and random temporal components. Spatio‐temporal models that use conditional autoregressive smoothing across the spatial dimension and autoregressive smoothing over the temporal dimension are developed. The model also accommodates the interaction between space and time. However, the effect of seasonal factors has not been previously addressed and in some applications (e.g., health conditions), these effects may not be negligible. The authors incorporate the seasonal effects of month and possibly year as part of the proposed model and estimate model parameters through generalized estimating equations. The model provides smoothed maps of disease risk and eliminates the instability of estimates in low‐population areas while maintaining geographic resolution. They illustrate the approach using a monthly data set of the number of asthma presentations made by children to Emergency Departments (EDs) in the province of Alberta, Canada, during the period 2001–2004. The Canadian Journal of Statistics 38: 698–715; 2010 © 2010 Statistical Society of Canada     

Spatio-temporal modelling using B-spline for disease mapping: analysis of childhood cancer trends

To examine childhood cancer diagnoses in the province of Alberta, Canada during 1983–2004, we construct a generalized additive mixed model for the analysis of geographic and temporal variability of cancer ratios. In this model, spatially correlated random effects and temporal components are adopted. The interaction between space and time is also accommodated. Spatio-temporal models that use conditional autoregressive smoothing across the spatial dimension and B-spline over the temporal dimension are considered. We study the patterns of incidence ratios over time and identify areas with consistently high ratio estimates as areas for potential further investigation. We apply the method of penalized quasi-likelihood to estimate the model parameters. We illustrate this approach using a yearly data set of childhood cancer diagnoses in the province of Alberta, Canada during 1983–2004.     

Interlitter Response Variability in a Threshold Dose-Response Model

In this article, we extend a previously formulated threshold dose-response model with random litter effects that was applied to a data set from a developmental toxicity study. The dose-response pattern of the data indicates that a threshold dose level may exist. Additionally, there is noticeable variation between the responses across the dose levels. With threshold estimation being critical, the assumed variability structure should adequately model the variation while not taking away from the estimation of the threshold as well as the other parameters directly involved in the dose-response relationship. In the prior formulation, the random effect was modeled assuming identical variation in the interlitter response probabilities across all dose levels, that is, the model had a single parameter to account for the interlitter variability. In this new model, the random effect is modeled as having different response variability across dose levels, that is, multiple interlitter variability parameters. We performed the likelihood ratio test (LRT) to compare our extended model to the previous model. We conducted a simulation study to compare the bias of each model when fit to data generated with the underlying parametric structure of the opposing model. The extended threshold dose-response model with multiple response variation was less biased.     

Bayesian Spatial Analysis for the Evaluation of Breast Cancer Detection Methods

This study investigated the impact of spatial location on the effectiveness of population‐based breast screening in reducing breast cancer mortality compared to other detection methods among Queensland women. The analysis was based on linked population‐based datasets from BreastScreen Queensland and the Queensland Cancer Registry for the period of 1997–2008 for women aged less than 90 years at diagnosis. A Bayesian hierarchical regression modelling approach was adopted and posterior estimation was performed using Markov Chain Monte Carlo techniques. This approach accommodated sparse data resulting from rare outcomes in small geographic areas, while allowing for spatial correlation and demographic influences to be included. A relative survival model was chosen to evaluate the relative excess risk for each breast cancer related factor. Several models were fitted to examine the influence of demographic information, cancer stage, geographic information and detection method on women's relative survival. Overall, the study demonstrated that including the detection method and geographic information when assessing the relative survival of breast cancer patients helped capture unexplained and spatial variability. The study also found evidence of better survival among women with breast cancer diagnosed in a screening program than those detected otherwise, as well as lower risk for those residing in a more urban or socio‐economically advantaged region, even after adjusting for tumour stage, environmental factors and demographics. However, no evidence of dependency between method of detection and geographic location was found. This project provides a sophisticated approach to examining the benefit of a screening program while considering the influence of geographic factors.     

Bayesian analysis of multivariate mortality data with large families

This paper presents a Bayesian method for the analysis of toxicological multivariate mortality data when the discrete mortality rate for each family of subjects at a given time depends on familial random effects and the toxicity level experienced by the family. Our aim is to model and analyse one set of such multivariate mortality data with large family sizes: the potassium thiocyanate (KSCN) tainted fish tank data of O'Hara Hines. The model used is based on a discretized hazard with additional time-varying familial random effects. A similar previous study (using sodium thiocyanate (NaSCN)) is used to construct a prior for the parameters in the current study. A simulation-based approach is used to compute posterior estimates of the model parameters and mortality rates and several other quantities of interest. Recent tools in Bayesian model diagnostics and variable subset selection have been incorporated to verify important modelling assumptions regarding the effects of time and heterogeneity among the families on the mortality rate. Further, Bayesian methods using predictive distributions are used for comparing several plausible models.     

Estimation in the multipath change point problem for correlated data

In many experiments, several measurements on the same variable are taken over time, a geographic region, or some other index set. It is often of interest to know if there has been a change over the index set in the parameters of the distribution of the variable. Frequently, the data consist of a sequence of correlated random variables, and there may also be several experimental units under observation, each providing a sequence of data. A problem in ascertaining the boundaries between the layers in geological sedimentary beds is used to introduce the model and then to illustrate the proposed methodology. It is assumed that, conditional on the change point, the data from each sequence arise from an autoregressive process that undergoes a change in one or more of its parameters. Unconditionally, the model then becomes a mixture of nonstationary autoregressive processes. Maximum-likelihood methods are used, and results of simulations to evaluate the performance of these estimators under practical conditions are given.     

Frailty modeling for clustered survival data: an application to birth interval in Bangladesh

The present work demonstrates an application of random effects model for analyzing birth intervals that are clustered into geographical regions. Observations from the same cluster are assumed to be correlated because usually they share certain unobserved characteristics between them. Ignoring the correlations among the observations may lead to incorrect standard errors of the estimates of parameters of interest. Beside making the comparisons between Cox's proportional hazards model and random effects model for analyzing geographically clustered time-to-event data, important demographic and socioeconomic factors that may affect the length of birth intervals of Bangladeshi women are also reported in this paper.     

A new behavioural model for fertility schedules

Modelling age-specific fertility rates is of great importance in demography because of their influence on population growth. Although we have a variety of fertility models in the demographic literature, most of them do not have any demographic interpretation for their parameters. It is generally expected that models with behavioural interpretation are more universal than those without any interpretation. Even though the famous Gompertz model has some behavioural interpretation it suffers from other drawbacks. In the present work, we propose a new fertility model, which has its genesis in the generalization of logistic law. The proposed model has good behavioural interpretation, alongside having nice parameter interpretations.     

Stochastic Differential Mixed-Effects Models

Abstract.  Stochastic differential equations have been shown useful in describing random continuous time processes. Biomedical experiments often imply repeated measurements on a series of experimental units and differences between units can be represented by incorporating random effects into the model. When both system noise and random effects are considered, stochastic differential mixed-effects models ensue. This class of models enables the simultaneous representation of randomness in the dynamics of the phenomena being considered and variability between experimental units, thus providing a powerful modelling tool with immediate applications in biomedicine and pharmacokinetic/pharmacodynamic studies. In most cases the likelihood function is not available, and thus maximum likelihood estimation of the unknown parameters is not possible. Here we propose a computationally fast approximated maximum likelihood procedure for the estimation of the non-random parameters and the random effects. The method is evaluated on simulations from some famous diffusion processes and on real data sets.     

Modelling Trends and Inequality in Small Area Mortality

This paper considers the modelling of mortality rates classified by age, time, and small area with a view to developing life table parameters relevant to assessing trends in inequalities in life chances. In particular, using a fully Bayes perspective, one may assess the stochastic variation in small area life table parameters, such as life expectancies, and also formally assess whether trends in indices of inequality in mortality are significant. Modelling questions include choice between random walk priors for age and time effects as against non-linear regression functions, questions of identifiability when several random effects are present in the death rates model, and the choice of model when both within and out-of-sample performance may be important. A case study application involves 44 small areas in North East London and mortality in five sub-periods (1986-88, 1989-91, 1992-94, 1995-97, 1998-2000) between 1986 and 2000, with the final period used for assessing out-of-sample performance.     

Impact of misspecified residual correlation structure on the parameter estimates in a shared spatial frailty model

In practice, survival data are often collected over geographical regions. Shared spatial frailty models have been used to model spatial variation in survival times, which are often implemented using the Bayesian Markov chain Monte Carlo method. However, this method comes at the price of slow mixing rates and heavy computational cost, which may render it impractical for data-intensive application. Alternatively, a frailty model assuming an independent and identically distributed (iid) random effect can be easily and efficiently implemented. Therefore, we used simulations to assess the bias and efficiency loss in the estimated parameters, if residual spatial correlation is present but using an iid random effect. Our simulations indicate that a shared frailty model with an iid random effect can estimate the regression coefficients reasonably well, even with residual spatial correlation present, when the percentage of censoring is not too high and the number of clusters and cluster size are not too low. Therefore, if the primary goal is to assess the covariate effects, one may choose the frailty model with an iid random effect; whereas if the goal is to predict the hazard, additional care needs to be given due to the efficiency loss in the parameter(s) for the baseline hazard.     

Multivariate time series projections of parameterized age-specific fertility rates

1 solution to the dimensionality problem raised by projection of individual age-specific fertility rates is the use of parametric curves to approximate the annual age-specific rates and a multivariate time series model to forecast the curve parameters. Such a method reduces the number of time series to be modeled for women 14-45 years of age from 32 to 40 (the number of curve parameters). In addition, the curves force even longterm fertility projections to exhibit the same smooth distribution across age as historical data. The data base used to illustrate this approach was age-specific fertility rates for US white women in 1921-84. An important advantage of this model is that it permits investigation of the interactions among the total fertility rate, the mean age of childbearing, and the standard deviation of age at childbearing. In the analysis of this particular data base, the contemporaneous relationship between the mean and standard deviation of age at childbearing was the only significant relationship. The addition of bias forecasts to the forecast gamma curve improves forecast accuracy, especially 1-2 years ahead. The most recent US Census Bureau projections have combined a time series model with longterm projections based on demographic judgment. These official projections yielded a slightly higher ultimate mean age and slightly lower standard deviation than those resulting from the model described in this paper.     

A Comparison of Methods for Analysing a Nested Frailty Model to Child Survival in Malawi

Demographic and Health Surveys collect child survival times that are clustered at the family and community levels. It is assumed that each cluster has a specific, unobservable, random frailty that induces an association in the survival times within the cluster. The Cox proportional hazards model, with family and community random frailties acting multiplicatively on the hazard rate, is presented. The estimation of the fixed effect and the association parameters of the modified model is then examined using the Gibbs sampler and the expectation–maximization (EM) algorithm. The methods are compared using child survival data collected in the 1992 Demographic and Health Survey of Malawi. The two methods lead to very similar estimates of fixed effect parameters. However, the estimates of random effect variances from the EM algorithm are smaller than those of the Gibbs sampler. Both estimation methods reveal considerable family variation in the survival of children, and very little variability over the communities.     

Accurate estimation for extra-Poisson variability assuming random effect models

In this study, the components of extra-Poisson variability are estimated assuming random effect models under a Bayesian approach. A standard existing methodology to estimate extra-Poisson variability assumes a negative binomial distribution. The obtained results show that using the proposed random effect model it is possible to get more accurate estimates for the extra-Poisson variability components when compared to the use of a negative binomial distribution where it is possible to estimate only one component of extra-Poisson variability. Some illustrative examples are introduced considering real data sets.     

Estimating the Prediction Mean Squared Error in Gaussian Stochastic Processes with Exponential Correlation Structure

Given one or more realizations from the finite dimensional marginal distribution of a stochastic process, we consider the problem of estimating the squared prediction error when predicting the process at unobserved locations. An approximation taking into account the additional variability due to estimating parameters involved in the correlation structure was developed by Kackar & Harville (1984) and was revisited by Harville & Jeske (1992) as well as Zimmerman & Cressie (1992). The present paper discusses an extension of these methods. The approaches will be compared via an extensive simulation study for models with and without random error term. Effects due to the designs used for prediction and for model fitting as well as due to the strength of the correlation between neighbouring observations of the stochastic process are investigated. The results show that considering the additional variability in the predictor due to estimating the covariance structure is of great importance and should not be neglected in practical applications.     

A multivariate two-factor skew model

It is well known that many data, such as the financial or demographic data, exhibit asymmetric distributions. In recent years, researchers have concentrated their efforts to model this asymmetry. Skew normal model is one of such models that are skew and yet possess many properties of the normal model. In this paper, a new multivariate skew model is proposed, along with its statistical properties. It includes the multivariate normal distribution and multivariate skew normal distribution as special cases. The quadratic form of this random vector follows a χ² distribution. The roles of the parameters in the model are investigated using contour plots of bivariate densities.     

Estimation of the additive hazards model with interval-censored data and missing covariates

The additive hazards model is one of the most commonly used regression models in the analysis of failure time data and many methods have been developed for its inference in various situations. However, no established estimation procedure exists when there are covariates with missing values and the observed responses are interval-censored; both types of complications arise in various settings including demographic, epidemiological, financial, medical and sociological studies. To address this deficiency, we propose several inverse probability weight-based and reweighting-based estimation procedures for the situation where covariate values are missing at random. The resulting estimators of regression model parameters are shown to be consistent and asymptotically normal. The numerical results that we report from a simulation study suggest that the proposed methods work well in practical situations. An application to a childhood cancer survival study is provided. The Canadian Journal of Statistics 48: 499–517; 2020 © 2020 Statistical Society of Canada     

Bayesian Spatial Modeling of Housing Prices Subject to a Localized Externality

This work proposes a non stationary random field model to describe the spatial variability of housing prices that are affected by a localized externality. The model allows for the effect of the localized externality on house prices to be represented in the mean function and/or the covariance function of the random field. The correlation function of the proposed model is a mixture of an isotropic correlation function and a correlation function that depends on the distances between home sales and the localized externality. The model is fit using a Bayesian approach via a Markov chain Monte Carlo algorithm. A dataset of 437 single family home sales during 2001 in the city of Cedar Falls, Iowa, is used to illustrate the model.     

The Influence of Random Effects on Univariate and Bivariate Discrete Proportional Hazards Models

The random effects survival model has been widely used in the recent literature as a generalization of the continuous proportional hazards model. When a random effect is present, it is known that the hazard rates are generally underestimated in the context of continuous proportional hazards models. This article establishes theorems for the influence of random effects on both univariate and bivariate discrete proportional hazards models.     

A monte carlo study of variable selection and inferences in a two stage random coefficient linear regression model with unbalanced data

A growth curve analysis is often applied to estimate patterns of changes in a given characteristic of different individuals. It is also used to find out if the variations in the growth rates among individuals are due to effects of certain covariates. In this paper, a random coefficient linear regression model, as a special case of the growth curve analysis, is generalized to accommodate the situation where the set of influential covariates is not known a priori. Two different approaches for seleaing influential covariates (a weighted stepwise selection procedure and a modified version of Rao and Wu’s selection criterion) for the random slope coefficient of a linear regression model with unbalanced data are proposed. Performances of these methods are evaluated by means of Monte-Carlo simulation. In addition, several methods (Maximum Likelihood, Restricted Maximum Likelihood, Pseudo Maximum Likelihood and Method of Moments) for estimating the parameters of the selected model are compared Proposed variable selection schemes and estimators are appliedtotheactualindustrial problem which motivated this investigation.