Open Access
Estimating unemployment in very small areas
(Institut d'Estadística de Catalunya-IDESCAT, 2009) Ugarte Martínez, María Dolores; Goicoa Mangado, Tomás; Militino, Ana F.; Sagaseta López, M.; Estadística e Investigación Operativa; Estatistika eta Ikerketa Operatiboa
In the last few years, European countries have shown a deep interest in applying small area techniques to produce reliable estimates at county level. However, the specificity of every European country and the heterogeneity of the available auxiliary information, make the use of a common methodology a very difficult task. In this study, the performance of several design-based, model-assisted, and model-based estimators using different auxiliary information for estimating unemployment at small area level is analyzed. The results are illustrated with data from Navarre, an autonomous region located at the north of Spain and divided into seven small areas. After discussing pros and cons of the different alternatives, a composite estimator is chosen, because of its good trade-off between bias and variance. Several methods for estimating the prediction error of the proposed estimator are also provided.
Open Access
Temporal evolution of brain cancer incidence in the municipalities of Navarre and the Basque Country, Spain
(BioMed Central, 2015) Ugarte Martínez, María Dolores; Adin Urtasun, Aritz; Goicoa Mangado, Tomás; Casado, Itziar; Ardanaz, Eva; Larrañaga, Nerea; Estatistika eta Ikerketa Operatiboa; Institute for Advanced Materials and Mathematics - INAMAT2; Estadística e Investigación Operativa; Gobierno de Navarra / Nafarroako Gobernua: proyecto 113 Res. 2186/2014
Background: Brain cancer incidence rates in Spain are below the European’s average. However, there are two regions in the north of the country, Navarre and the Basque Country, ranked among the European regions with the highest incidence rates for both males and females. Our objective here was two-fold. Firstly, to describe the temporal evolution of the geographical pattern of brain cancer incidence in Navarre and the Basque Country, and secondly, to look for specific high risk areas (municipalities) within these two regions in the study period (1986–2008). Methods: A mixed Poisson model with two levels of spatial effects is used. The model also included two levels of spatial effects (municipalities and local health areas). Model fitting was carried out using penalized quasi-likelihood. High risk regions were detected using upper one-sided confidence intervals. Results: Results revealed a group of high risk areas surrounding Pamplona, the capital city of Navarre, and a few municipalities with significant high risks in the northern part of the region, specifically in the border between Navarre and the Basque Country (Gipuzkoa). The global temporal trend was found to be increasing. Differences were also observed among specific risk evolutions in certain municipalities. Conclusions: Brain cancer incidence in Navarre and the Basque Country (Spain) is still increasing with time. The number of high risk areas within those two regions is also increasing. Our study highlights the need of continuous surveillance of this cancer in the areas of high risk. However, due to the low percentage of cases explained by the known risk factors, primary prevention should be applied as a general recommendation in these populations.
Open Access
Bayesian inference in multivariate spatio-temporal areal models using INLA: analysis of gender-based violence in small areas
(Springer, 2020) Vicente Fuenzalida, Gonzalo; Goicoa Mangado, Tomás; Ugarte Martínez, María Dolores; Estatistika, Informatika eta Matematika; Institute for Advanced Materials and Mathematics - INAMAT2; Estadística, Informática y Matemáticas
Multivariate models for spatial count data are currently receiving attention in disease mapping to model two or more diseases jointly. They have been thoroughly studied from a theoretical point of view, but their use in practice is still limited because they are computationally expensive and, in general, they are not implemented in standard software to be used routinely. Here, a new multivariate proposal, based on the recently derived M models for spatial data, is developed for spatio-temporal areal data. The model takes account of the correlation between the spatial and temporal patterns of the phenomena being studied, and it also includes spatio-temporal interactions. Though multivariate models have been traditionally fitted using Markov chain Monte Carlo techniques, here we propose to adopt integrated nested Laplace approximations to speed up computations as results obtained using both fitting techniques were nearly identical. The techniques are used to analyse two forms of crimes against women in India. In particular, we focus on the joint analysis of rapes and dowry deaths in Uttar Pradesh, the most populated Indian state, during the years 2001-2014.
Open Access
Online relative risks/rates estimation in spatial and spatio-temporal disease mapping
(Elsevier, 2019) Adin Urtasun, Aritz; Goicoa Mangado, Tomás; Ugarte Martínez, María Dolores; Estatistika, Informatika eta Matematika; Institute for Advanced Materials and Mathematics - INAMAT2; Estadística, Informática y Matemáticas
Background and objective: Spatial and spatio-temporal analyses of count data are crucial in epidemiology and other fields to unveil spatial and spatio-temporal patterns of incidence and/or mortality risks. However, fitting spatial and spatio-temporal models is not easy for non-expert users. The objective of this paper is to present an interactive and user-friendly web application (named SSTCDapp) for the analysis of spatial and spatio-temporal mortality or incidence data. Although SSTCDapp is simple to use, the underlying statistical theory is well founded and all key issues such as model identifiability, model selection, and several spatial priors and hyperpriors for sensitivity analyses are properly addressed. Methods: The web application is designed to fit an extensive range of fairly complex spatio-temporal models to smooth the very often extremely variable standardized incidence/mortality risks or crude rates. The application is built with the R package shiny and relies on the well founded integrated nested Laplace approximation technique for model fitting and inference. Results: The use of the web application is shown through the analysis of Spanish spatio-temporal breast cancer data. Different possibilities for the analysis regarding the type of model, model selection criteria, and a range of graphical as well as numerical outputs are provided. Conclusions: Unlike other software used in disease mapping, SSTCDapp facilitates the fit of complex statistical models to non-experts users without the need of installing any software in their own computers, since all the analyses and computations are made in a powerful remote server. In addition, a desktop version is also available to run the application locally in those cases in which data confidentiality is a serious issue.
Open Access
A two-stage approach to estimate spatial and spatio-temporal disease risks in the presence of local discontinuities and clusters
(SAGE, 2018-04-13) Adin Urtasun, Aritz; Lee, Duncan; Goicoa Mangado, Tomás; Ugarte Martínez, María Dolores; Estadística, Informática y Matemáticas; Estatistika, Informatika eta Matematika; Institute for Advanced Materials and Mathematics - INAMAT2
Disease risk maps for areal unit data are often estimated from Poisson mixed models with local spatial smoothing, for example by incorporating random effects with a conditional autoregressive prior distribution. However, one of the limitations is that local discontinuities in the spatial pattern are not usually modelled, leading to over-smoothing of the risk maps and a masking of clusters of hot/coldspot areas. In this paper, we propose a novel two-stage approach to estimate and map disease risk in the presence of such local discontinuities and clusters. We propose approaches in both spatial and spatio-temporal domains, where for the latter the clusters can either be fixed or allowed to vary over time. In the first stage, we apply an agglomerative hierarchical clustering algorithm to training data to provide sets of potential clusters, and in the second stage, a two-level spatial or spatio-temporal model is applied to each potential cluster configuration. The superiority of the proposed approach with regard to a previous proposal is shown by simulation, and the methodology is applied to two important public health problems in Spain, namely stomach cancer mortality across Spain and brain cancer incidence in the Navarre and Basque Country regions of Spain.
Open Access
A simplified spatial+ approach to mitigate spatial confounding in multivariate spatial areal models
(Elsevier, 2024) Urdangarin Iztueta, Arantxa; Goicoa Mangado, Tomás; Kneib, Thomas; Ugarte Martínez, María Dolores; Estadística, Informática y Matemáticas; Estatistika, Informatika eta Matematika; Institute for Advanced Materials and Mathematics - INAMAT2
Spatial areal models encounter the well-known and challenging problem of spatial confounding. This issue makes it arduous to distinguish between the impacts of observed covariates and spatial random effects. Despite previous research and various proposed methods to tackle this problem, finding a definitive solution remains elusive. In this paper, we propose a simplified version of the spatial+ approach that involves dividing the covariate into two components. One component captures large-scale spatial dependence, while the other accounts for short-scale dependence. This approach eliminates the need to separately fit spatial models for the covariates. We apply this method to analyse two forms of crimes against women, namely rapes and dowry deaths, in Uttar Pradesh, India, exploring their relationship with socio-demographic covariates. To evaluate the performance of the new approach, we conduct extensive simulation studies under different spatial confounding scenarios. The results demonstrate that the proposed method provides reliable estimates of fixed effects and posterior correlations between different responses.
Open Access
In spatio-temporal disease mapping models, identifiability constraints affect PQL and INLA results
(Springer, 2018) Goicoa Mangado, Tomás; Adin Urtasun, Aritz; Ugarte Martínez, María Dolores; Hodges, James S.; Institute for Advanced Materials and Mathematics - INAMAT2
Disease mapping studies the distribution of relative risks or rates in space and time, and typically relies on generalized linear mixed models (GLMMs) including fixed effects and spatial, temporal, and spatio-temporal random effects. These GLMMs are typically not identifiable and constraints are required to achieve sensible results. However, automatic specification of constraints can sometimes lead to misleading results. In particular, the penalized quasi-likelihood fitting technique automatically centers the random effects even when this is not necessary. In the Bayesian approach, the recently-introduced integrated nested Laplace approximations computing technique can also produce wrong results if constraints are not wellspecified. In this paper the spatial, temporal, and spatiotemporal interaction random effects are reparameterized using the spectral decompositions of their precision matrices to establish the appropriate identifiability constraints. Breast cancer mortality data from Spain is used to illustrate the ideas.
Open Access
Using mortality to predict incidence for rare and lethal cancers in very small areas
(VCH Publishers, 2022) Etxeberria Andueza, Jaione; Goicoa Mangado, Tomás; Ugarte Martínez, María Dolores; Estadística, Informática y Matemáticas; Estatistika, Informatika eta Matematika; Institute for Advanced Materials and Mathematics - INAMAT2
Incidence and mortality figures are needed to get a comprehensive overview of cancer burden. In many countries, cancer mortality figures are routinely recorded by statistical offices, whereas incidence depends on regional cancer registries. However, due to the complexity of updating cancer registries, incidence numbers become available 3 or 4 years later than mortality figures. It is, therefore, necessary to develop reliable procedures to predict cancer incidence at least until the period when mortality data are available. Most of the methods proposed in the literature are designed to predict total cancer (except nonmelanoma skin cancer) or major cancer sites. However, less frequent lethal cancers, such as brain cancer, are generally excluded from predictions because the scarce number of cases makes it difficult to use univariate models. Our proposal comes to fill this gap and consists of modeling jointly incidence and mortality data using spatio-temporal models with spatial and age shared components. This approach allows for predicting lethal cancers improving the performance of individual models when data are scarce by taking advantage of the high correlation between incidence and mortality. A fully Bayesian approach based on integrated nested Laplace approximations is considered for model fitting and inference. A validation process is also conducted to assess the performance of alternative models. We use the new proposals to predict brain cancer incidence rates by gender and age groups in the health units of Navarre and Basque Country (Spain) during the period 2005-2008.
Open Access
Space-time interactions in bayesian disease mapping with recent tools: making things easier for practitioners
(Edward Arnold, 2022) Urdangarin Iztueta, Arantxa; Ugarte Martínez, María Dolores; Goicoa Mangado, Tomás; Estatistika, Informatika eta Matematika; Institute for Advanced Materials and Mathematics - INAMAT2; Estadística, Informática y Matemáticas
Spatio-temporal disease mapping studies the distribution of mortality or incidence risks in space and its evolution in time, and it usually relies on fitting hierarchical Poisson mixed models. These models are complex for practitioners as they generally require adding constraints to correctly identify and interpret the different model terms. However, including constraints may not be straightforward in some recent software packages. This paper focuses on NIMBLE, a library of algorithms that contains among others a configurable system for Markov chain Monte Carlo (MCMC) algorithms. In particular, we show how to fit different spatio-temporal disease mapping models with NIMBLE making emphasis on how to include sum-to-zero constraints to solve identifiability issues when including spatio-temporal interactions. Breast cancer mortality data in Spain during the period 1990-2010 is used for illustration purposes. A simulation study is also conducted to compare NIMBLE with R-INLA in terms of parameter estimates and relative risk estimation. The results are very similar but differences are observed in terms of computing time.
Open Access
Effect of thermal stress on tissue ultrastructure and metabolite profiles during initiation of Radiata pine somatic embryogenesis
(Frontiers Media, 2019) Castander Olarieta, Ander; Montalbán, Itziar A.; Oliveira, Eliana de Medeiros; Dell'Aversana, Emilia; D'Amelia, Luisa; Carillo, Petronia; Steiner, Neusa; Fraga, Hugo Pacheco de Freitas; Guerra, Miguel Pedro; Goicoa Mangado, Tomás; Ugarte Martínez, María Dolores; Estadística, Informática y Matemáticas; Estatistika, Informatika eta Matematika
Climate change will inevitably lead to environmental variations, thus plant drought tolerance will be a determinant factor in the success of plantations and natural forestry recovery. Some metabolites, such as soluble carbohydrates and amino acids, have been described as being the key to both embryogenesis efficiency and abiotic stress response, contributing to phenotypic plasticity and the adaptive capacity of plants. For this reason, our main objectives were to evaluate if the temperature during embryonal mass initiation in radiata pine was critical to the success of somatic embryogenesis, to alter the morphological and ultrastructural organization of embryonal masses at cellular level and to modify the carbohydrate, protein, or amino acid contents. The first SE initiation experiments were carried out at moderate and high temperatures for periods of different durations prior to transfer to the control temperature of 23 degrees C. Cultures initiated at moderate temperatures (30 degrees C, 4 weeks and 40 degrees C, 4 days) showed significantly lower initiation and proliferation rates than those at the control temperature or pulse treatment at high temperatures (50 degrees C, 5 min). No significant differences were observed either for the percentage of embryogenic cell lines that produced somatic embryos, or for the number of somatic embryos per gram of embryonal mass. Based on the results from the first experiments, initiation was carried out at 40 degrees C 4 h; 50 degrees C, 30 min; and a pulse treatment of 60 degrees C, 5 min. No significant differences were found for the initiation or number of established lines or for the maturation of somatic embryos. However, large morphological differences were observed in the mature somatic embryos. At the same time, changes observed at cellular level suggested that strong heat shock treatments may trigger the programmed cell death of embryogenic cells, leading to an early loss of embryogenic potential, and the formation of supernumerary suspensor cells. Finally, among all the differences observed in the metabolic profile, it is worth highlighting the accumulation of tyrosine and isoleucine, both amino acids involved in the synthesis of abiotic stress response-related secondary metabolites.