Pupillary behavior in relation to wavelength and age

Pupil light reflex can be used as a non-invasive ocular predictor of cephalic autonomic nervous system integrity. Spectral sensitivity of the pupil's response to light has, for some time, been an interesting issue. It has generally, however, only been investigated with the use of white light and studies with monochromatic wavelengths are scarce. This study investigates the effects of wavelength and age within three parameters of the pupil light reflex (amplitude of response, latency, and velocity of constriction) in a large sample of younger and older adults (N = 97), in mesopic conditions. Subjects were exposed to a single light stimulus at four different wavelengths: white (5600° K), blue (450 nm), green (510 nm), and red (600 nm). Data was analyzed appropriately, and, when applicable, using the General Linear Model (GLM), Randomized Complete Block Design (RCBD), Student's t-test and/or ANCOVA. Across all subjects, pupillary response to light had the greatest amplitude and shortest latency in white and green light conditions. In regards to age, older subjects (46-78 years) showed an increased latency in white light and decreased velocity of constriction in green light compared to younger subjects (18-45 years old). This study provides data patterns on parameters of wavelength-dependent pupil reflexes to light in adults and it contributes to the large body of pupillometric research. It is hoped that this study will add to the overall evaluation of cephalic autonomic nervous system integrity.

Realized Matrix-Exponential Stochastic Volatility with Asymmetry, Long Memory and Spillovers

The paper develops a novel realized matrix-exponential stochastic volatility model of multivariate returns and realized covariances that incorporates asymmetry and long memory (hereafter the RMESV-ALM model). The matrix exponential transformation guarantees the positivedefiniteness of the dynamic covariance matrix. The contribution of the paper ties in with Robert Basmann’s seminal work in terms of the estimation of highly non-linear model specifications (“Causality tests and observationally equivalent representations of econometric models”, Journal of Econometrics, 1988, 39(1-2), 69–104), especially for developing tests for leverage and spillover effects in the covariance dynamics. Efficient importance sampling is used to maximize the likelihood function of RMESV-ALM, and the finite sample properties of the quasi-maximum likelihood estimator of the parameters are analysed. Using high frequency data for three US financial assets, the new model is estimated and evaluated. The forecasting performance of the new model is compared with a novel dynamic realized matrix-exponential conditional covariance model. The volatility and co-volatility spillovers are examined via the news impact curves and the impulse response functions from returns to volatility and co-volatility.

Avaliação comparativa da remodelação óssea periimplantar em implantes com junção cone-morse e hexágono externo: estudo clínico randomizado controlado, duplo cego, boca dividida

Background: Several theories, such as the biological width formation, the inflammatory reactions due to the implant-abutment microgap contamination, and the periimplant stress/strain concentration causing bone microdamage accumulation, have been suggested to explain early periimplant bone loss. However, it is yet not well understood to which extent the implant-abutment connection type may influence the remodeling process around dental implants. Aim: to evaluate clinical, bacteriological, and biomechanical parameters related to periimplant bone loss at the crestal region, comparing external hexagon (EH) and Morse-taper (MT) connections. Materials and methods: Twelve patients with totally edentulous mandibles received four custom made Ø 3.8 x 13 mm implants in the interforaminal region of the mandible, with the same design, but different prosthetic connections (two of them EH or MT, randomly placed based on a split-mouth design), and a immediate implant- supported prosthesis. Clinical parameters (periimplant probing pocket depth, modified gingival index and mucosal thickness) were evaluated at 6 sites around the implants, at a 12 month follow-up. The distance from the top of the implant to the first bone-to-implant contact – IT-FBIC was evaluated on standardized digital peri-apical radiographs acquired at 1, 3, 6 and 12 months follow-up. Samples of the subgingival microbiota were collected 1, 3 and 6 months after implant loading. DNA were extracted and used for the quantification of Tanerella forsythia, Porphyromonas gingivalis, Aggragatibacter actinomycetemcomitans, Prevotella intermedia and Fusobacterium nucleatum. Comparison among multiple periods of observation were performed using repeated-measures Analysis of Variance (ANOVA), followed by a Tukey post-hoc test, while two-period based comparisons were made using paired t- test. Further, 36 computer-tomographic based finite element (FE) models were accomplished, simulating each patient in 3 loading conditions. The results for the peak EQV strain in periimplant bone were interpreted by means of a general linear model (ANOVA). Results: The variation in periimplant bone loss assessed by means of radiographs was significantly different between the connection types (P<0.001). Mean IT-FBIC was 1.17±0.44 mm for EH, and 0.17±0.54 mm for MT, considering all evaluated time periods. All clinical parameters presented not significant differences. No significant microbiological differences could be observed between both connection types. Most of the collected samples had very few pathogens, meaning that these regions were healthy from a microbiological point of view. In FE analysis, a significantly higher peak of EQV strain (P=0.005) was found for EH (mean 3438.65 µ∑) compared to MT (mean 840.98 µ∑) connection. Conclusions: Varying implant-abutment connection type will result in diverse periimplant bone remodeling, regardless of clinical and microbiological conditions. This fact is more likely attributed to the singular loading transmission through different implant-abutment connections to the periimplant bone. The present findings suggest that Morse-taper connection is more efficient to prevent periimplant bone loss, compared to an external hexagon connection.

Facial affect processing inbipolar disorder:a magnetoencephalography study

Background: Identifying biological markers to aid diagnosis of bipolar disorder (BD) is critically important. To be considered a possible biological marker, neural patterns in BD should be discriminant from those in healthy individuals (HI). We examined patterns of neuromagnetic responses revealed by magnetoencephalography (MEG) during implicit emotion-processing using emotional (happy, fearful, sad) and neutral facial expressions, in sixteen BD and sixteen age- and gender-matched healthy individuals. Methods: Neuromagnetic data were recorded using a 306-channel whole-head MEG ELEKTA Neuromag System, and preprocessed using Signal Space Separation as implemented in MaxFilter (ELEKTA). Custom Matlab programs removed EOG and ECG signals from filtered MEG data, and computed means of epoched data (0-250ms, 250-500ms, 500-750ms). A generalized linear model with three factors (individual, emotion intensity and time) compared BD and HI. A principal component analysis of normalized mean channel data in selected brain regions identified principal components that explained 95% of data variation. These components were used in a quadratic support vector machine (SVM) pattern classifier. SVM classifier performance was assessed using the leave-one-out approach. Results: BD and HI showed significantly different patterns of activation for 0-250ms within both left occipital and temporal regions, specifically for neutral facial expressions. PCA analysis revealed significant differences between BD and HI for mild fearful, happy, and sad facial expressions within 250-500ms. SVM quadratic classifier showed greatest accuracy (84%) and sensitivity (92%) for neutral faces, in left occipital regions within 500-750ms. Conclusions: MEG responses may be used in the search for disease specific neural markers.

Bayesian Inference in Large-scale Problems

Many modern applications fall into the category of "large-scale" statistical problems, in which both the number of observations n and the number of features or parameters p may be large. Many existing methods focus on point estimation, despite the continued relevance of uncertainty quantification in the sciences, where the number of parameters to estimate often exceeds the sample size, despite huge increases in the value of n typically seen in many fields. Thus, the tendency in some areas of industry to dispense with traditional statistical analysis on the basis that "n=all" is of little relevance outside of certain narrow applications. The main result of the Big Data revolution in most fields has instead been to make computation much harder without reducing the importance of uncertainty quantification. Bayesian methods excel at uncertainty quantification, but often scale poorly relative to alternatives. This conflict between the statistical advantages of Bayesian procedures and their substantial computational disadvantages is perhaps the greatest challenge facing modern Bayesian statistics, and is the primary motivation for the work presented here.

Two general strategies for scaling Bayesian inference are considered. The first is the development of methods that lend themselves to faster computation, and the second is design and characterization of computational algorithms that scale better in n or p. In the first instance, the focus is on joint inference outside of the standard problem of multivariate continuous data that has been a major focus of previous theoretical work in this area. In the second area, we pursue strategies for improving the speed of Markov chain Monte Carlo algorithms, and characterizing their performance in large-scale settings. Throughout, the focus is on rigorous theoretical evaluation combined with empirical demonstrations of performance and concordance with the theory.

One topic we consider is modeling the joint distribution of multivariate categorical data, often summarized in a contingency table. Contingency table analysis routinely relies on log-linear models, with latent structure analysis providing a common alternative. Latent structure models lead to a reduced rank tensor factorization of the probability mass function for multivariate categorical data, while log-linear models achieve dimensionality reduction through sparsity. Little is known about the relationship between these notions of dimensionality reduction in the two paradigms. In Chapter 2, we derive several results relating the support of a log-linear model to nonnegative ranks of the associated probability tensor. Motivated by these findings, we propose a new collapsed Tucker class of tensor decompositions, which bridge existing PARAFAC and Tucker decompositions, providing a more flexible framework for parsimoniously characterizing multivariate categorical data. Taking a Bayesian approach to inference, we illustrate empirical advantages of the new decompositions.

Latent class models for the joint distribution of multivariate categorical, such as the PARAFAC decomposition, data play an important role in the analysis of population structure. In this context, the number of latent classes is interpreted as the number of genetically distinct subpopulations of an organism, an important factor in the analysis of evolutionary processes and conservation status. Existing methods focus on point estimates of the number of subpopulations, and lack robust uncertainty quantification. Moreover, whether the number of latent classes in these models is even an identified parameter is an open question. In Chapter 3, we show that when the model is properly specified, the correct number of subpopulations can be recovered almost surely. We then propose an alternative method for estimating the number of latent subpopulations that provides good quantification of uncertainty, and provide a simple procedure for verifying that the proposed method is consistent for the number of subpopulations. The performance of the model in estimating the number of subpopulations and other common population structure inference problems is assessed in simulations and a real data application.

In contingency table analysis, sparse data is frequently encountered for even modest numbers of variables, resulting in non-existence of maximum likelihood estimates. A common solution is to obtain regularized estimates of the parameters of a log-linear model. Bayesian methods provide a coherent approach to regularization, but are often computationally intensive. Conjugate priors ease computational demands, but the conjugate Diaconis--Ylvisaker priors for the parameters of log-linear models do not give rise to closed form credible regions, complicating posterior inference. In Chapter 4 we derive the optimal Gaussian approximation to the posterior for log-linear models with Diaconis--Ylvisaker priors, and provide convergence rate and finite-sample bounds for the Kullback-Leibler divergence between the exact posterior and the optimal Gaussian approximation. We demonstrate empirically in simulations and a real data application that the approximation is highly accurate, even in relatively small samples. The proposed approximation provides a computationally scalable and principled approach to regularized estimation and approximate Bayesian inference for log-linear models.

Another challenging and somewhat non-standard joint modeling problem is inference on tail dependence in stochastic processes. In applications where extreme dependence is of interest, data are almost always time-indexed. Existing methods for inference and modeling in this setting often cluster extreme events or choose window sizes with the goal of preserving temporal information. In Chapter 5, we propose an alternative paradigm for inference on tail dependence in stochastic processes with arbitrary temporal dependence structure in the extremes, based on the idea that the information on strength of tail dependence and the temporal structure in this dependence are both encoded in waiting times between exceedances of high thresholds. We construct a class of time-indexed stochastic processes with tail dependence obtained by endowing the support points in de Haan's spectral representation of max-stable processes with velocities and lifetimes. We extend Smith's model to these max-stable velocity processes and obtain the distribution of waiting times between extreme events at multiple locations. Motivated by this result, a new definition of tail dependence is proposed that is a function of the distribution of waiting times between threshold exceedances, and an inferential framework is constructed for estimating the strength of extremal dependence and quantifying uncertainty in this paradigm. The method is applied to climatological, financial, and electrophysiology data.

The remainder of this thesis focuses on posterior computation by Markov chain Monte Carlo. The Markov Chain Monte Carlo method is the dominant paradigm for posterior computation in Bayesian analysis. It has long been common to control computation time by making approximations to the Markov transition kernel. Comparatively little attention has been paid to convergence and estimation error in these approximating Markov Chains. In Chapter 6, we propose a framework for assessing when to use approximations in MCMC algorithms, and how much error in the transition kernel should be tolerated to obtain optimal estimation performance with respect to a specified loss function and computational budget. The results require only ergodicity of the exact kernel and control of the kernel approximation accuracy. The theoretical framework is applied to approximations based on random subsets of data, low-rank approximations of Gaussian processes, and a novel approximating Markov chain for discrete mixture models.

Data augmentation Gibbs samplers are arguably the most popular class of algorithm for approximately sampling from the posterior distribution for the parameters of generalized linear models. The truncated Normal and Polya-Gamma data augmentation samplers are standard examples for probit and logit links, respectively. Motivated by an important problem in quantitative advertising, in Chapter 7 we consider the application of these algorithms to modeling rare events. We show that when the sample size is large but the observed number of successes is small, these data augmentation samplers mix very slowly, with a spectral gap that converges to zero at a rate at least proportional to the reciprocal of the square root of the sample size up to a log factor. In simulation studies, moderate sample sizes result in high autocorrelations and small effective sample sizes. Similar empirical results are observed for related data augmentation samplers for multinomial logit and probit models. When applied to a real quantitative advertising dataset, the data augmentation samplers mix very poorly. Conversely, Hamiltonian Monte Carlo and a type of independence chain Metropolis algorithm show good mixing on the same dataset.

Advances in Dynamic Modeling and Computation for Count Data

A class of multi-process models is developed for collections of time indexed count data. Autocorrelation in counts is achieved with dynamic models for the natural parameter of the binomial distribution. In addition to modeling binomial time series, the framework includes dynamic models for multinomial and Poisson time series. Markov chain Monte Carlo (MCMC) and Po ́lya-Gamma data augmentation (Polson et al., 2013) are critical for fitting multi-process models of counts. To facilitate computation when the counts are high, a Gaussian approximation to the P ́olya- Gamma random variable is developed.

Three applied analyses are presented to explore the utility and versatility of the framework. The first analysis develops a model for complex dynamic behavior of themes in collections of text documents. Documents are modeled as a “bag of words”, and the multinomial distribution is used to characterize uncertainty in the vocabulary terms appearing in each document. State-space models for the natural parameters of the multinomial distribution induce autocorrelation in themes and their proportional representation in the corpus over time.

The second analysis develops a dynamic mixed membership model for Poisson counts. The model is applied to a collection of time series which record neuron level firing patterns in rhesus monkeys. The monkey is exposed to two sounds simultaneously, and Gaussian processes are used to smoothly model the time-varying rate at which the neuron’s firing pattern fluctuates between features associated with each sound in isolation.

The third analysis presents a switching dynamic generalized linear model for the time-varying home run totals of professional baseball players. The model endows each player with an age specific latent natural ability class and a performance enhancing drug (PED) use indicator. As players age, they randomly transition through a sequence of ability classes in a manner consistent with traditional aging patterns. When the performance of the player significantly deviates from the expected aging pattern, he is identified as a player whose performance is consistent with PED use.

All three models provide a mechanism for sharing information across related series locally in time. The models are fit with variations on the P ́olya-Gamma Gibbs sampler, MCMC convergence diagnostics are developed, and reproducible inference is emphasized throughout the dissertation.

Work-to-family spillover: Gender differences in Hungary

It is crucial to understand the role that labor market positions might play in creating gender differences in work–life balance. One theoretical approach to understanding this relationship is the spillover theory. The spillover theory argues that an individual’s life domains are integrated; meaning that well-being can be transmitted between life domains. Based on data collected in Hungary in 2014, this paper shows that work-to-family spillover does not affect both genders the same way. The effect of work on family life tends to be more negative for women than for men. Two explanations have been formulated in order to understand this gender inequality. According to the findings of the analysis, gender is conditionally independent of spillover if financial status and flexibility of work are also incorporated into the analysis. This means that the relative disadvantage for women in terms of spillover can be attributed to their lower financial status and their relatively low access to flexible jobs. In other words, the gender inequalities in work-to-family spillover are deeply affected by individual labor market positions. The observation of the labor market’s effect on work–life balance is especially important in Hungary since Hungary has one of the least flexible labor arrangements in Europe. A marginal log-linear model, which is a method for categorical multivariate analysis, has been applied in this analysis.

Modelagem dinâmica e simulação computacional de poços de petróleo verticais e direcionais com elevação por bombeio mecânico

LINS, Filipe C. A. et al. Modelagem dinâmica e simulação computacional de poços de petróleo verticais e direcionais com elevação por bombeio mecânico. In: CONGRESSO BRASILEIRO DE PESQUISA E DESENVOLVIMENTO EM PETRÓLEO E GÁS, 5. 2009, Fortaleza, CE. Anais... Fortaleza: CBPDPetro, 2009.

Learning environment, preparedness and satisfaction in osteopathy in Europe: The PreSS study

Objective: 1) to assess the preparedness to practice and satisfaction in learning environment amongst new graduates from European osteopathic institutions; 2) to compare the results of preparedness to practice and satisfaction in learning environment between and within countries where osteopathy is regulated and where regulation is still to be achieved; 3) to identify possible correlations between learning environment and preparedness to practice. Method: Osteopathic education providers of full-time education located in Europe were enrolled, and their final year students were contacted to complete a survey. Measures used were: Dundee Ready Educational Environment Measure (DREEM), the Association of American Medical Colleges (AAMC) and a demographic questionnaire. Scores were compared across institutions using one-way ANOVA and generalised linear model. Results: Nine European osteopathic education institutions participated in the study (4 located in Italy, 2 in the UK, 1 in France, 1 in Belgium and 1 in the Netherlands) and 243 (77%) of their final-year students completed the survey. The DREEM total score mean was 121.4 (SEM: 1.66) whilst the AAMC was 17.58 (SEM:0.35). A generalised linear model found a significant association between not-regulated countries and total score as well as subscales DREEM scores (p<0.001). Learning environment and preparedness to practice were significantly positively correlated (r=0.76; p<0.01). Discussion: A perceived higher level of preparedness and satisfaction was found amongst students from osteopathic institutions located in countries without regulation compared to those located in countries where osteopathy is regulated; however, all institutions obtained a 'more positive than negative' result. Moreover, in general, cohorts with fewer than 20 students scored significantly higher compared to larger student cohorts. Finally, an overall positive correlation between students' preparedness and satisfaction were found across all institutions recruited.

Analyse des tendances des séries temporelles de la concentration en azote minéral dissous dans les trois principaux estuaires de la façade Manche-Atlantique

Estuaries are areas which, from their structure, their fonctioning, and their localisation, are subject to significant contribution of nutrients. One of the objectif of the RNO, the French network for coastal water quality monitoring, is to assess the levels and trends of nutrient concentrations in estuaries. A linear model was used in order to describe and to explain the total dissolved nitrogen concentration evolution in the three most important estuaries on the Chanel-Atlantic front (Seine, Loire and Gironde). As a first step, the selection of a reliable data set was performed. Then total dissolved nitrogen evolution schemes in estuary environment were graphically studied, and allowed a resonable choice of covariables. The salinity played a major role in explaining nitrogen concentration variability in estuary, and dilution lines were proved to be a useful tool to detect outlying observations and to model the nitrogenlsalinity relation. Increasing trends were detected by the model, with a high magnitude in Seine, intermediate in Loire, and lower in Gironde. The non linear trend estimated in Loire and Seine estuaries could be due to important interannual variations as suggest in graphics. In the objective of the QUADRIGE database valorisation, a discussion on the statistical model, and on the RNO hydrological data sampling strategy, allowed to formulate suggestions towards a better exploitation of nutrient data.

Assessing the influence of abiotic factors and leaf-level properties on the stability of growing-season canopy greenness in a deciduous forest

Maps depicting spatial pattern in the stability of summer greenness could advance understanding of how forest ecosystems will respond to global changes such as a longer growing season. Declining summer greenness, or “greendown”, is spectrally related to declining near-infrared reflectance and is observed in most remote sensing time series to begin shortly after peak greenness at the end of spring and extend until the beginning of leaf coloration in autumn,. Understanding spatial patterns in the strength of greendown has recently become possible with the advancement of Landsat phenology products, which show that greendown patterns vary at scales appropriate for linking these patterns to proposed environmental forcing factors. This study tested two non-mutually exclusive hypotheses for how leaf measurements and environmental factors correlate with greendown and decreasing NIR reflectance across sites. At the landscape scale, we used linear regression to test the effects of maximum greenness, elevation, slope, aspect, solar irradiance and canopy rugosity on greendown. Secondly, we used leaf chemical traits and reflectance observations to test the effect of nitrogen availability and intrinsic water use efficiency on leaf-level greendown, and landscape-level greendown measured from Landsat. The study was conducted using Quercus alba canopies across 21 sites of an eastern deciduous forest in North America between June and August 2014. Our linear model explained greendown variance with an R2=0.47 with maximum greenness as the greatest model effect. Subsequent models excluding one model effect revealed elevation and aspect were the two topographic factors that explained the greatest amount of greendown variance. Regression results also demonstrated important interactions between all three variables, with the greatest interaction showing that aspect had greater influence on greendown at sites with steeper slopes. Leaf-level reflectance was correlated with foliar δ13C (proxy for intrinsic water use efficiency), but foliar δ13C did not translate into correlations with landscape-level variation in greendown from Landsat. Therefore, we conclude that Landsat greendown is primarily indicative of landscape position, with a small effect of canopy structure, and no measureable effect of leaf reflectance. With this understanding of Landsat greendown we can better explain the effects of landscape factors on vegetation reflectance and perhaps on phenology, which would be very useful for studying phenology in the context of global climate change

Modelagem dinâmica e simulação computacional de poços de petróleo verticais e direcionais com elevação por bombeio mecânico

LINS, Filipe C. A. et al. Modelagem dinâmica e simulação computacional de poços de petróleo verticais e direcionais com elevação por bombeio mecânico. In: CONGRESSO BRASILEIRO DE PESQUISA E DESENVOLVIMENTO EM PETRÓLEO E GÁS, 5. 2009, Fortaleza, CE. Anais... Fortaleza: CBPDPetro, 2009.

Efeitos dos fatores ambientais na comunidade de macrófitas em pequenos cursos lóticos perenes subtropicais

Neste estudo foi investigado como a distribuição das espécies e a produção de biomassa de macrófitas aquáticas são influenciadas pelas condições físico-químicas do ambiente. Também foi avaliado como uma espécie com maior potencial competitivo pode interferir na diversidade de espécies da comunidade macrofítica. Para tanto, em cada um dos três arroios, foram dispostos seis transecções, perpendiculares à margem. Em cada transecção foram demarcadas três unidades amostrais de 1m², nas quais foram registrados os parâmetros fitossociológicos cobertura e frequência relativas e valor de importância. A diversidade de espécies foi estimada pelo índice de Shannon, utilizando os valores de cobertura de espécies. Para determinar a biomassa das macrófitas aquáticas foram usados quadrats de 0,25m², alocados dentro da unidade amostral de 1m² usadas para quantificar os dados fitossociológicos, nos mesmos pontos onde foi feito o levantamento de cobertura da vegetação. Utilizamos como variáveis preditoras a velocidade da corrente, radiação solar incidente, coeficiente de sombreamento, vegetação ripária arbórea adjacente, nitrogênio orgânico dissolvido, carbono orgânico dissolvido e condutividade elétrica. Foram registradas 32 espécies de macrófitas aquáticas, distribuídas em 19 famílias e 28 gêneros. Conforme Análise de Correspondência Canônica (CCA), as espécies com maiores valores de biomassa foram relacionadas a unidades amostrais com alta incidência luminosa. As unidades amostrais com dominância de Pistia stratiotes apresentaram menor diversidade de espécies indicando que esta espécie, quando encontra condições que permitam sua proliferação, pode excluir espécies de menor potencial competitivo. De acordo com GLM (Generalized Linear Model), a ausência de vegetação ripária ou presente em apenas uma das margens e baixas velocidades de corrente configura-se em condições favoráveis ao estabelecimento e desenvolvimento de macrófitas aquáticas, possibilitando produção maiores valores de biomassa.

A generic method of engagement to elicit regional coastal management options

Stakeholder engagement is important for successful management of natural resources, both to make effective decisions and to obtain support. However, in the context of coastal management, questions remain unanswered on how to effectively link decisions made at the catchment level with objectives for marine biodiversity and fisheries productivity. Moreover, there is much uncertainty on how to best elicit community input in a rigorous manner that supports management decisions. A decision support process is described that uses the adaptive management loop as its basis to elicit management objectives, priorities and management options using two case studies in the Great Barrier Reef, Australia. The approach described is then generalised for international interest. A hierarchical engagement model of local stakeholders, regional and senior managers is used. The result is a semi-quantitative generic elicitation framework that ultimately provides a prioritised list of management options in the context of clearly articulated management objectives that has widespread application for coastal communities worldwide. The case studies show that demand for local input and regional management is high, but local influences affect the relative success of both engagement processes and uptake by managers. Differences between case study outcomes highlight the importance of discussing objectives prior to suggesting management actions, and avoiding or minimising conflicts at the early stages of the process. Strong contributors to success are a) the provision of local information to the community group, and b) the early inclusion of senior managers and influencers in the group to ensure the intellectual and time investment is not compromised at the final stages of the process. The project has uncovered a conundrum in the significant gap between the way managers perceive their management actions and outcomes, and community's perception of the effectiveness (and wisdom) of these same management actions.

Interaction between mercury (Hg), arsenic (As) and selenium (Se) affects the activity of glutathione S-transferase in breast milk: possible relationship with fish and shellfish intake

Breast milk is regarded as an ideal source of nutrients for the growth and development of neonates, but it can also be a potential source of pollutants. Mothers can be exposed to different contaminants as a result of their lifestyle and environmental pollution. Mercury (Hg) and arsenic (As) could adversely affect the development of fetal and neonatal nervous system. Some fish and shellfish are rich in selenium (Se), an essential trace element that forms part of several enzymes related to the detoxification process, including glutathione S-transferase (GST). The goal of this study was to determine the interaction between Hg, As and Se and analyze its effect on the activity of GST in breast milk. Milk samples were collected from women between day 7 and 10 postpartum. The GST activity was determined spectrophotometrically; total Hg, As and Se concentrations were measured by atomic absorption spectrometry. To explain the possible association of Hg, As and Se concentrations with GST activity in breast milk, generalized linear models were constructed. The model explained 44% of the GST activity measured in breast milk. The GLM suggests that GST activity was positively correlated with Hg, As and Se concentrations. The activity of the enzyme was also explained by the frequency of consumption of marine fish and shellfish in the diet of the breastfeeding women.