Open Access
Improving the quality of satellite imagery based on ground-truth data from rain gauge stations
(MDPI, 2018) Militino, Ana F.; Ugarte Martínez, María Dolores; Pérez Goya, Unai; Estatistika eta Ikerketa Operatiboa; Institute for Advanced Materials and Mathematics - INAMAT2; Estadística e Investigación Operativa; Gobierno de Navarra / Nafarroako Gobernua
Multitemporal imagery is by and large geometrically and radiometrically accurate, but the residual noise arising from removal clouds and other atmospheric and electronic effects can produce outliers that must be mitigated to properly exploit the remote sensing information. In this study, we show how ground-truth data from rain gauge stations can improve the quality of satellite imagery. To this end, a simulation study is conducted wherein different sizes of outlier outbreaks are spread and randomly introduced in the normalized difference vegetation index (NDVI) and the day and night land surface temperature (LST) of composite images from Navarre (Spain) between 2011 and 2015. To remove outliers, a new method called thin-plate splines with covariates (TpsWc) is proposed. This method consists of smoothing the median anomalies with a thin-plate spline model, whereby transformed ground-truth data are the external covariates of the model. The performance of the proposed method is measured with the square root of the mean square error (RMSE), calculated as the root of the pixel-by-pixel mean square differences between the original data and the predicted data with the TpsWc model and with a state-space model with and without covariates. The study shows that the use of ground-truth data reduces the RMSE in both the TpsWc model and the state-space model used for comparison purposes. The new method successfully removes the abnormal data while preserving the phenology of the raw data. The RMSE reduction percentage varies according to the derived variables (NDVI or LST), but reductions of up to 20% are achieved with the new proposal.
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
Interpolation of the mean anomalies for cloud filling in land surface temperature and normalized difference vegetation index
(IEEE, 2019) Militino, Ana F.; Ugarte Martínez, María Dolores; Pérez Goya, Unai; Genton, Marc G.; Estatistika, Informatika eta Matematika; Institute for Advanced Materials and Mathematics - INAMAT2; Estadística, Informática y Matemáticas
When monitoring time series of remote sensing data, it is advisable to fill gaps, i.e., missing or distorted data, caused by atmospheric effects or technical failures. In this paper, a new method for filling these gaps called interpolation of the mean anomalies (IMA) is proposed and compared with some competitors. The method consists of: 1) defining a neighborhood for the target image from previous and subsequent images across previous and subsequent years; 2) computing the mean target image of the neighborhood; 3) estimating the anomalies in the target image by subtracting the mean image from the target image; 4) filtering the anomalies; 5) averaging the anomalies over a predefined window; 6) interpolating the averaged anomalies; and 7) adding the interpolated anomalies to the mean image. To assess the performance of the IMA method, both a real example and a simulation study are conducted with a time series of Moderate Resolution Imaging Spectroradiometer (MODIS) TERRA and MODIS AQUA images captured over the region of Navarre (Spain) from 2011 to 2013. We analyze the land surface temperature (LST) day and night, and the normalized difference vegetation index (NDVI). In the simulation study, seven sizes of artificial clouds are randomly introduced to each image in the studied time series. The square root of the mean-squared prediction error (RMSE) between the original and the filled data is chosen as an indicator of the goodness of fit. The results show that the IMA method outperforms Timesat, Hants, and Gapfill (GF) in filling small, moderate, and big cloud gaps in both the day and night LST and NDVI data, reaching RMSE reductions of up to 23%.
Open Access
A scalable approach for short-term disease forecasting in high spatial resolution areal data
(Wiley-VCH, 2023) Orozco Acosta, Erick; Riebler, Andrea; Adin Urtasun, Aritz; Ugarte Martínez, María Dolores; Estadística, Informática y Matemáticas; Estatistika, Informatika eta Matematika; Institute for Advanced Materials and Mathematics - INAMAT2; Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa
Short-term disease forecasting at specific discrete spatial resolutions has become a high-impact decision-support tool in health planning. However, when the number of areas is very large obtaining predictions can be computationally intensive or even unfeasible using standard spatiotemporal models. The purpose of this paper is to provide a method for short-term predictions in high-dimensional areal data based on a newly proposed ¿divide-and-conquer¿ approach. We assess the predictive performance of this method and other classical spatiotemporal models in a validation study that uses cancer mortality data for the 7907 municipalities of continental Spain. The new proposal outperforms traditional models in terms of mean absolute error, root mean square error, and interval score when forecasting cancer mortality 1, 2, and 3 years ahead. Models are implemented in a fully Bayesian framework using the well-known integrated nested Laplace estimation technique.
Open Access
Estimación del desempleo por comarcas en Navarra
(Gobierno de Navarra, Departamento de Economía y Hacienda, 2005) Ugarte Martínez, María Dolores; Militino, Ana F.; González Ramajo, Begoña; Goicoa Mangado, Tomás; Sagaseta López, M.; Estadística e Investigación Operativa; Estatistika eta Ikerketa Operatiboa
El conocimiento del desempleo en una región es un indicador potente del ritmo de crecimiento de una economía, ya que de forma indirecta mide su capacidad para generar empleo. El Instituto de Estadística de Navarra está apostando por proporcionar en un futuro cercano estimaciones del desempleo a un nivel cada vez más desagregado. La heterogeneidad de las comarcas navarras y el interés mostrado por administraciones locales y sindicatos, hace necesario tener un conocimiento de la situación de desempleo a nivel comarcal, evitando así descansar únicamente en el resultado global para toda Navarra tal y como lo proporciona la Encuesta de Población Activa (EPA). La tarea es compleja, pero está incardinada además en uno de los objetivos prioritarios del proyecto europeo EURAREA, del cual ha formado parte el Instituto Nacional de Estadística (INE), y por ende, el Instituto de Estadística de Navarra. Es decir, hay un interés real en Europa por proporcionar estimaciones a nivel comarcal. En Navarra esta tarea ya ha comenzado y en este congreso presentamos algunos de los resultados obtenidos. En particular se ilustran las estimaciones preliminares derivadas de la aplicación de diversos estimadores basados en el diseño para obtener la proporción de parados por sexo en las siete comarcas de Navarra. Se compara además el comportamiento de diversos estimadores en términos del sesgo relativo y del error cuadrático medio relativo. Los estimadores ofrecidos permiten calcular además la estimación del número de ocupados e inactivos, así como de sus correspondientes tasas.
Open Access
Alleviating confounding in spatio-temporal areal models with an application on crimes against women in India
(SAGE Publications, 2021) Adin Urtasun, Aritz; Goicoa Mangado, Tomás; Hodges, James S.; Schnell, Patrick M.; Ugarte Martínez, María Dolores; Estatistika, Informatika eta Matematika; Institute for Advanced Materials and Mathematics - INAMAT2; Estadística, Informática y Matemáticas
Assessing associations between a response of interest and a set of covariates in spatial areal models is the leitmotiv of ecological regression. However, the presence of spatially correlated random effects can mask or even bias estimates of such associations due to confounding effects if they are not carefully handled. Though potentially harmful, confounding issues have often been ignored in practice leading to wrong conclusions about the underlying associations between the response and the covariates. In spatio-temporal areal models, the temporal dimension may emerge as a new source of confounding, and the problem may be even worse. In this work, we propose two approaches to deal with confounding of fixed effects by spatial and temporal random effects, while obtaining good model predictions. In particular, restricted regression and an apparently—though in fact not—equivalent procedure using constraints are proposed within both fully Bayes and empirical Bayes approaches. The methods are compared in terms of fixed-effect estimates and model selection criteria. The techniques are used to assess the association between dowry deaths and certain socio-demographic covariates in the districts of Uttar Pradesh, India.
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
Two-level resolution of relative risk of dengue disease in a hyperendemic city of Colombia
(Public Library of Science, 2018) Adin Urtasun, Aritz; Martínez Bello, Daniel Adyro; López Quílez, Antonio; Ugarte Martínez, María Dolores; Estatistika, Informatika eta Matematika; Institute for Advanced Materials and Mathematics - INAMAT2; Estadística, Informática y Matemáticas
Risk maps of dengue disease offer to the public health officers a tool to model disease risk in space and time. We analyzed the geographical distribution of relative incidence risk of dengue disease in a high incidence city from Colombia, and its evolution in time during the period January 2009—December 2015, identifying regional effects at different levels of spatial aggregations. Cases of dengue disease were geocoded and spatially allocated to census sectors, and temporally aggregated by epidemiological periods. The census sectors are nested in administrative divisions defined as communes, configuring two levels of spatial aggregation for the dengue cases. Spatio-temporal models including census sector and commune-level spatially structured random effects were fitted to estimate dengue incidence relative risks using the integrated nested Laplace approximation (INLA) technique. The final selected model included two-level spatial random effects, a global structured temporal random effect, and a census sector-level interaction term. Risk maps by epidemiological period and risk profiles by census sector were generated from the modeling process, showing the transmission dynamics of the disease. All the census sectors in the city displayed high risk at some epidemiological period in the outbreak periods. Relative risk estimation of dengue disease using INLA offered a quick and powerful method for parameter estimation and inference.
Open Access
Flexible Bayesian P-splines for smoothing age-specific spatio-temporal mortality patterns
(SAGE, 2019) Goicoa Mangado, Tomás; Adin Urtasun, Aritz; Etxeberria Andueza, Jaione; Militino, Ana F.; Ugarte Martínez, María Dolores; Estadística, Informática y Matemáticas; Estatistika, Informatika eta Matematika; Institute for Advanced Materials and Mathematics - INAMAT2
In this paper age-space-time models based on one and two-dimensional P-splines with B-spline bases are proposed for smoothing mortality rates, where both xed relative scale and scale invariant two-dimensional penalties are examined. Model tting and inference are carried out using integrated nested Laplace approximations (INLA), a recent Bayesian technique that speeds up computations compared to McMC methods. The models will be illustrated with Spanish breast cancer mortality data during the period 1985-2010, where a general decline in breast cancer mortality has been observed in Spanish provinces in the last decades. The results reveal that mortality rates for the oldest age groups do not decrease in all provinces.
Open Access
Detecting change-points in the time series of surfaces occupied by pre-defined NDVI categories in continental Spain from 1981 to 2015
(Springer, 2018) Militino, Ana F.; Ugarte Martínez, María Dolores; Pérez Goya, Unai; Estadística, Informática y Matemáticas; Estatistika, Informatika eta Matematika; Institute for Advanced Materials and Mathematics - INAMAT2
The free access to satellite images since more than 40 years ago has provoked a rapid increase of multitemporal derived information of remote sensing data that should be summarized and analyzed for future inferences. In particular, the study of trends and trend changes is of crucial interest in many studies of phenology, climatology, agriculture, hydrology, geology or many other environmental disciplines. Overall, the normalized dierence vegetation index (NDVI), as a satellite derived variable, plays a crucial role because of its usefulness for vegetation and landscape characterization, land use and land cover mapping, environmental monitoring, climate change or crop prediction models. Since the eighties, it can be retrieved all over the world from dierent satellites. In this work we propose to analyze its temporal evolution, looking for breakpoints or change-points in trends of the surfaces occupied by four NDVI classications made in Spain from 1981 to 2015. The results show a decrease of bare soils and semi-bare soils starting in the middle nineties or before, and a slight increase of middle-vegetation and high-vegetation soils starting in 1990 and 2000 respectively.