Influence diagnostics in Gaussian spatial linear models

Spatial linear models have been applied in numerous fields such as agriculture, geoscience and environmental sciences, among many others. Spatial dependence structure modelling, using a geostatistical approach, is an indispensable tool to estimate the parameters that define this structure. However, this estimation may be greatly affected by the presence of atypical observations in the sampled data. The purpose of this paper is to use diagnostic techniques to assess the sensitivity of the maximum-likelihood estimators, covariance functions and linear predictor to small perturbations in the data and/or the spatial linear model assumptions. The methodology is illustrated with two real data sets. The results allowed us to conclude that the presence of atypical values in the sample data have a strong influence on thematic maps, changing the spatial dependence structure.

A multiple time scale survival model with a cure fraction

Many recent survival studies propose modeling data with a cure fraction, i.e., data in which part of the population is not susceptible to the event of interest. This event may occur more than once for the same individual (recurrent event). We then have a scenario of recurrent event data in the presence of a cure fraction, which may appear in various areas such as oncology, finance, industries, among others. This paper proposes a multiple time scale survival model to analyze recurrent events using a cure fraction. The objective is analyzing the efficiency of certain interventions so that the studied event will not happen again in terms of covariates and censoring. All estimates were obtained using a sampling-based approach, which allows information to be input beforehand with lower computational effort. Simulations were done based on a clinical scenario in order to observe some frequentist properties of the estimation procedure in the presence of small and moderate sample sizes. An application of a well-known set of real mammary tumor data is provided.

A new long-term lifetime distribution induced by a latent complementary risk framework

In this paper, we proposed a new three-parameter long-term lifetime distribution induced by a latent complementary risk framework with decreasing, increasing and unimodal hazard function, the long-term complementary exponential geometric distribution. The new distribution arises from latent competing risk scenarios, where the lifetime associated scenario, with a particular risk, is not observable, rather we observe only the maximum lifetime value among all risks, and the presence of long-term survival. The properties of the proposed distribution are discussed, including its probability density function and explicit algebraic formulas for its reliability, hazard and quantile functions and order statistics. The parameter estimation is based on the usual maximum-likelihood approach. A simulation study assesses the performance of the estimation procedure. We compare the new distribution with its particular cases, as well as with the long-term Weibull distribution on three real data sets, observing its potential and competitiveness in comparison with some usual long-term lifetime distributions.

Integration of spatial and temporal models to predict extreme water table depths

The objective of this work was to evaluate extreme water table depths in a watershed, using methods for geographical spatial data analysis. Groundwater spatio-temporal dynamics was evaluated in an outcrop of the Guarani Aquifer System. Water table depths were estimated from monitoring of water levels in 23 piezometers and time series modeling available from April 2004 to April 2011. For generation of spatial scenarios, geostatistical techniques were used, which incorporated into the prediction ancillary information related to the geomorphological patterns of the watershed, using a digital elevation model. This procedure improved estimates, due to the high correlation between water levels and elevation, and aggregated physical sense to predictions. The scenarios showed differences regarding the extreme levels - too deep or too shallow ones - and can subsidize water planning, efficient water use, and sustainable water management in the watershed.

Spatial analysis of egg distribution and geographic changes in the spawning habitat of the Brazilian sardine Sardinella brasiliensis

This paper establishes the spawning habitat of the Brazilian sardine Sardinella brasiliensis and investigates the spatial variability of egg density and its relation with oceanographic conditions in the shelf of the south-east Brazil Bight (SBB). The spawning habitats of S. brasiliensis have been defined in terms of spatial models of egg density, temperature-salinity plots, quotient (Q) analysis and remote sensing data. Quotient curves (Q(C)) were constructed using the geographic distribution of egg density, temperature and salinity from samples collected during nine survey cruises between 1976 and 1993. The interannual sea surface temperature (SST) variability was determined using principal component analysis on the SST anomalies (SSTA) estimated from remote sensing data over the period between 1985 and 2007. The spatial pattern of egg occurrences in the SBB indicated that the largest concentration occurred between Paranagua and Sao Sebastiao. Spawning habitat expanded and contracted during the years, fluctuating around Paranagua. In January 1978 and January 1993, eggs were found nearly everywhere along the inner shelf of the SBB, while in January 1988 and 1991 spawning had contracted to their southernmost position. The SSTA maps for the spawning periods showed that in the case of habitat expansion (1993 only) anomalies over the SBB were zero or slightly negative, whereas for the contraction period anomalies were all positive. Sardinella brasiliensis is capable of exploring suitable spawning sites provided by the entrainment of the colder and less-saline South Atlantic Central Water onto the shelf by means of both coastal wind-driven (to the north-east of the SBB) and meander-induced (to the south-west of the SBB) upwelling.

Cation Distribution and Local Order in Mixed Sodium Metaphosphate Glasses

Two structural properties in mixed alkali metal phosphate glasses that seem to be crucial to the development of the mixed ion effect in dc conductivity were systematically analyzed in Na mixed metaphosphates: the local order around the mobile species, and their distribution and mixing in the glass network. The set of glasses considered here, Na1-xMxPO3 with M = Li, Ag, K, Rb, and Cs and 0 <= x <= 1, encompass a broad degree of size mismatch between the mixed cation species. A comprehensive solid-state nuclear magnetic resonance study was carried out using P-31 MAS, Na-23 triple quantum MAS, Rb-87 QCPMG, P-31-Na-23 REDOR, Na-23-Li-7 and Li-7-Li-6 SEDOR, and Na-23 spin echo decay. It was observed that the arrangement of P atoms around Na in the mixed glasses was indistinguishable from that observed in the NaPO3 glass. However, systematic distortions in the local structure of the 0 environments around Na were observed, related to the presence of the second cation. The average Na-O distances show an expansion/compression When Na+ ions are replaced by cations with respectively smaller/bigger radii. The behavior of the nuclear electric quadrupole coupling. constants indicates that this expansion reduces the local symmetry, while the compression produces the opposite effect These effects become marginally small when the site mismatch between the cations is small, as in Na-Ag mixed glasses. The present study confirms the intimate mixing of cation species at the atomic scale, but clear deviations from random mixing were detected in systems with larger alkali metal ions (Cs-Na, K-Na, Rb-Na). In contrast, no deviations from the statistical ion mixture were found in the systems Ag-Na and Li-Na, where mixed cations are either of radii comparable to (Ag+) or smaller than (Li+) Na+. The set of results supports two fundamental structural features of the models proposed to explain the mixed ion effect: the. structural specificity of the sites occupied by each cation species and their mixing at the atomic scale.

Stochastic simulation of time-series models combined with geostatistics to predict water-table scenarios in a Guarani Aquifer System outcrop area, Brazil

Stochastic methods based on time-series modeling combined with geostatistics can be useful tools to describe the variability of water-table levels in time and space and to account for uncertainty. Monitoring water-level networks can give information about the dynamic of the aquifer domain in both dimensions. Time-series modeling is an elegant way to treat monitoring data without the complexity of physical mechanistic models. Time-series model predictions can be interpolated spatially, with the spatial differences in water-table dynamics determined by the spatial variation in the system properties and the temporal variation driven by the dynamics of the inputs into the system. An integration of stochastic methods is presented, based on time-series modeling and geostatistics as a framework to predict water levels for decision making in groundwater management and land-use planning. The methodology is applied in a case study in a Guarani Aquifer System (GAS) outcrop area located in the southeastern part of Brazil. Communication of results in a clear and understandable form, via simulated scenarios, is discussed as an alternative, when translating scientific knowledge into applications of stochastic hydrogeology in large aquifers with limited monitoring network coverage like the GAS.

Higher spatial derivative field theories

In this work, we employ renormalization group methods to study the general behavior of field theories possessing anisotropic scaling in the spacetime variables. The Lorentz symmetry breaking that accompanies these models are either soft, if no higher spatial derivative is present, or it may have a more complex structure if higher spatial derivatives are also included. Both situations are discussed in models with only scalar fields and also in models with fermions as a Yukawa-like model.

A general class of zero-or-one inflated beta regression models

This paper proposes a general class of regression models for continuous proportions when the data contain zeros or ones. The proposed class of models assumes that the response variable has a mixed continuous-discrete distribution with probability mass at zero or one. The beta distribution is used to describe the continuous component of the model, since its density has a wide range of different shapes depending on the values of the two parameters that index the distribution. We use a suitable parameterization of the beta law in terms of its mean and a precision parameter. The parameters of the mixture distribution are modeled as functions of regression parameters. We provide inference, diagnostic, and model selection tools for this class of models. A practical application that employs real data is presented. (C) 2011 Elsevier B.V. All rights reserved.

The importance of plant life form on spatial associations along a subtropical coastal dune gradient

Questions Does the spatial association between isolated adult trees and understorey plants change along a gradient of sand dunes? Does this association depend on the life form of the understorey plant? Location Coastal sand dunes, southeast Brazil. Methods We recorded the occurrence of understorey plant species in 100 paired 0.25 m2 plots under adult trees and in adjacent treeless sites along an environmental gradient from beach to inland. Occurrence probabilities were modelled as a function of the fixed variables of the presence of a neighbour, distance from the seashore and life form, and a random variable, the block (i.e. the pair of plots). Generalized linear mixed models (GLMM) were fitted in a backward step-wise procedure using Akaike's information criterion (AIC) for model selection. Results The occurrence of understorey plants was affected by the presence of an adult tree neighbour, but the effect varied with the life form of the understorey species. Positive spatial association was found between isolated adult neighbour and young trees, whereas a negative association was found for shrubs. Moreover, a neutral association was found for lianas, whereas for herbs the effect of the presence of an adult neighbour ranged from neutral to negative, depended on the subgroup considered. The strength of the negative association with forbs increased with distance from the seashore. However, for the other life forms, the associational pattern with adult trees did not change along the gradient. Conclusions For most of the understorey life forms there is no evidence that the spatial association between isolated adult trees and understorey plants changes with the distance from the seashore, as predicted by the stress gradient hypothesis, a common hypothesis in the literature about facilitation in plant communities. Furthermore, the positive spatial association between isolated adult trees and young trees identified along the entire gradient studied indicates a positive feedback that explains the transition from open vegetation to forest in subtropical coastal dune environments.

Spatial-temporal variations in phytoplankton size and colored detrital matter absorption at global and regional scales, as derived from twelve years of SeaWiFS data (1998-2009)

A procedure has been proposed by Ciotti and Bricaud (2006) to retrieve spectral absorption coefficients of phytoplankton and colored detrital matter (CDM) from satellite radiance measurements. This was also the first procedure to estimate a size factor for phytoplankton, based on the shape of the retrieved algal absorption spectrum, and the spectral slope of CDM absorption. Applying this method to the global ocean color data set acquired by SeaWiFS over twelve years (1998-2009), allowed for a comparison of the spatial variations of chlorophyll concentration ([Chl]), algal size factor (S-f), CDM absorption coefficient (a(cdm)) at 443 nm, and spectral slope of CDM absorption (S-cdm). As expected, correlations between the derived parameters were characterized by a large scatter at the global scale. We compared temporal variability of the spatially averaged parameters over the twelve-year period for three oceanic areas of biogeochemical importance: the Eastern Equatorial Pacific, the North Atlantic and the Mediterranean Sea. In all areas, both S-f and a(cdm)(443) showed large seasonal and interannual variations, generally correlated to those of algal biomass. The CDM maxima appeared in some occasions to last longer than those of [Chl]. The spectral slope of CDM absorption showed very large seasonal cycles consistent with photobleaching, challenging the assumption of a constant slope commonly used in bio-optical models. In the Equatorial Pacific, the seasonal cycles of [Chl], S-f, a(cdm)(443) and S-cdm, as well as the relationships between these parameters, were strongly affected by the 1997-98 El Ni o/La Ni a event.

Spatial variability of leaf nutrient contents in a drip irrigated citrus orchard

This study aimed to evaluate the spatial variability of leaf content of macro and micronutrients. The citrus plants orchard with 5 years of age, planted at regular intervals of 8 x 7 m, was managed under drip irrigation. Leaf samples were collected from each plant to be analyzed in the laboratory. Data were analyzed using the software R, version 2.5.1 Copyright (C) 2007, along with geostatistics package GeoR. All contents of macro and micronutrients studied were adjusted to normal distribution and showed spatial dependence.The best-fit models, based on the likelihood, for the macro and micronutrients were the spherical and matern. It is suggest for the macronutrients nitrogen, phosphorus, potassium, calcium, magnesium and sulfur the minimum distances between samples of 37; 58; 29; 63; 46 and 15 m respectively, while for the micronutrients boron, copper, iron, manganese and zinc, the distances suggests are 29; 9; 113; 35 and 14 m, respectively.

The role of asymmetries in rock-paper-scissors biodiversity models.

The maintenance of biodiversity is a long standing puzzle in ecology. It is a classical result that if the interactions of the species in an ecosystem are chosen in a random way, then complex ecosystems can't sustain themselves, meaning that the structure of the interactions between the species must be a central component on the preservation of biodiversity and on the stability of ecosystems. The rock-paper-scissors model is one of the paradigmatic models that study how biodiversity is maintained. In this model 3 species dominate each other in a cyclic way (mimicking a trophic cycle), that is, rock dominates scissors, that dominates paper, that dominates rock. In the original version of this model, this dominance obeys a 'Z IND 3' symmetry, in the sense that the strength of dominance is always the same. In this work, we break this symmetry, studying the effects of the addition of an asymmetry parameter. In the usual model, in a two dimensional lattice, the species distribute themselves according to spiral patterns, that can be explained by the complex Landau-Guinzburg equation. With the addition of asymmetry, new spatial patterns appear during the transient and the system either ends in a state with spirals, similar to the ones of the original model, or in a state where unstable spatial patterns dominate or in a state where only one species survives (and biodiversity is lost).