Using non-homogeneous Poisson models with multiple change-points to estimate the number of ozone exceedances in Mexico City

In this paper, we consider some non-homogeneous Poisson models to estimate the probability that an air quality standard is exceeded a given number of times in a time interval of interest. We assume that the number of exceedances occurs according to a non-homogeneous Poisson process (NHPP). This Poisson process has rate function lambda(t), t >= 0, which depends on some parameters that must be estimated. We take into account two cases of rate functions: the Weibull and the Goel-Okumoto. We consider models with and without change-points. When the presence of change-points is assumed, we may have the presence of either one, two or three change-points, depending of the data set. The parameters of the rate functions are estimated using a Gibbs sampling algorithm. Results are applied to ozone data provided by the Mexico City monitoring network. In a first instance, we assume that there are no change-points present. Depending on the adjustment of the model, we assume the presence of either one, two or three change-points. Copyright (C) 2009 John Wiley & Sons, Ltd.

Estimating the number of ozone peaks in Mexico City using a non-homogeneous Poisson model

In this paper, we consider the problem of estimating the number of times an air quality standard is exceeded in a given period of time. A non-homogeneous Poisson model is proposed to analyse this issue. The rate at which the Poisson events occur is given by a rate function lambda(t), t >= 0. This rate function also depends on some parameters that need to be estimated. Two forms of lambda(t), t >= 0 are considered. One of them is of the Weibull form and the other is of the exponentiated-Weibull form. The parameters estimation is made using a Bayesian formulation based on the Gibbs sampling algorithm. The assignation of the prior distributions for the parameters is made in two stages. In the first stage, non-informative prior distributions are considered. Using the information provided by the first stage, more informative prior distributions are used in the second one. The theoretical development is applied to data provided by the monitoring network of Mexico City. The rate function that best fit the data varies according to the region of the city and/or threshold that is considered. In some cases the best fit is the Weibull form and in other cases the best option is the exponentiated-Weibull. Copyright (C) 2007 John Wiley & Sons, Ltd.

Sequential sampling of Euschistus heros (Heteroptera: Pentatomidae) in soybean

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

A Model-Based Approach for Making Ecological Inference from Distance Sampling Data

We consider a fully model-based approach for the analysis of distance sampling data. Distance sampling has been widely used to estimate abundance (or density) of animals or plants in a spatially explicit study area. There is, however, no readily available method of making statistical inference on the relationships between abundance and environmental covariates. Spatial Poisson process likelihoods can be used to simultaneously estimate detection and intensity parameters by modeling distance sampling data as a thinned spatial point process. A model-based spatial approach to distance sampling data has three main benefits: it allows complex and opportunistic transect designs to be employed, it allows estimation of abundance in small subregions, and it provides a framework to assess the effects of habitat or experimental manipulation on density. We demonstrate the model-based methodology with a small simulation study and analysis of the Dubbo weed data set. In addition, a simple ad hoc method for handling overdispersion is also proposed. The simulation study showed that the model-based approach compared favorably to conventional distance sampling methods for abundance estimation. In addition, the overdispersion correction performed adequately when the number of transects was high. Analysis of the Dubbo data set indicated a transect effect on abundance via Akaike’s information criterion model selection. Further goodness-of-fit analysis, however, indicated some potential confounding of intensity with the detection function.

Proteins in silico-modeling and sampling

Proteins are linear chain molecules made out of amino acids. Only when they fold to their native states, they become functional. This dissertation aims to model the solvent (environment) effect and to develop & implement enhanced sampling methods that enable a reliable study of the protein folding problem in silico. We have developed an enhanced solvation model based on the solution to the Poisson-Boltzmann equation in order to describe the solvent effect. Following the quantum mechanical Polarizable Continuum Model (PCM), we decomposed net solvation free energy into three physical terms– Polarization, Dispersion and Cavitation. All the terms were implemented, analyzed and parametrized individually to obtain a high level of accuracy. In order to describe the thermodynamics of proteins, their conformational space needs to be sampled thoroughly. Simulations of proteins are hampered by slow relaxation due to their rugged free-energy landscape, with the barriers between minima being higher than the thermal energy at physiological temperatures. In order to overcome this problem a number of approaches have been proposed of which replica exchange method (REM) is the most popular. In this dissertation we describe a new variant of canonical replica exchange method in the context of molecular dynamic simulation. The advantage of this new method is the easily tunable high acceptance rate for the replica exchange. We call our method Microcanonical Replica Exchange Molecular Dynamic (MREMD). We have described the theoretical frame work, comment on its actual implementation, and its application to Trp-cage mini-protein in implicit solvent. We have been able to correctly predict the folding thermodynamics of this protein using our approach.

Importance sampling approximations to various probabilities of ruin of spectrally negative Lévy risk processes

This article provides importance sampling algorithms for computing the probabilities of various types ruin of spectrally negative Lévy risk processes, which are ruin over the infinite time horizon, ruin within a finite time horizon and ruin past a finite time horizon. For the special case of the compound Poisson process perturbed by diffusion, algorithms for computing probabilities of ruins by creeping (i.e. induced by the diffusion term) and by jumping (i.e. by a claim amount) are provided. It is shown that these algorithms have either bounded relative error or logarithmic efficiency, as t,x→∞t,x→∞, where t>0t>0 is the time horizon and x>0x>0 is the starting point of the risk process, with y=t/xy=t/x held constant and assumed either below or above a certain constant.

Saddlepoint Approximations to the Probability of Ruin in Finite Time for the Compound Poisson Risk Process Perturbed by Diffusion

A large deviations type approximation to the probability of ruin within a finite time for the compound Poisson risk process perturbed by diffusion is derived. This approximation is based on the saddlepoint method and generalizes the approximation for the non-perturbed risk process by Barndorff-Nielsen and Schmidli (Scand Actuar J 1995(2):169–186, 1995). An importance sampling approximation to this probability of ruin is also provided. Numerical illustrations assess the accuracy of the saddlepoint approximation using importance sampling as a benchmark. The relative deviations between saddlepoint approximation and importance sampling are very small, even for extremely small probabilities of ruin. The saddlepoint approximation is however substantially faster to compute.

Sampling sites, sampled birds and plants of young larch plantations in Tokachi district, eastern Hokkaido, northern Japan, 2011

Models and data used to describe species-area relationships confound sampling with ecological process as they fail to acknowledge that estimates of species richness arise due to sampling. This compromises our ability to make ecological inferences from and about species-area relationships. We develop and illustrate hierarchical community models of abundance and frequency to estimate species richness. The models we propose separate sampling from ecological processes by explicitly accounting for the fact that sampled patches are seldom completely covered by sampling plots and that individuals present in the sampling plots are imperfectly detected. We propose a multispecies abundance model in which community assembly is treated as the summation of an ensemble of species-level Poisson processes and estimate patch-level species richness as a derived parameter. We use sampling process models appropriate for specific survey methods. We propose a multispecies frequency model that treats the number of plots in which a species occurs as a binomial process. We illustrate these models using data collected in surveys of early-successional bird species and plants in young forest plantation patches. Results indicate that only mature forest plant species deviated from the constant density hypothesis, but the null model suggested that the deviations were too small to alter the form of species-area relationships. Nevertheless, results from simulations clearly show that the aggregate pattern of individual species density-area relationships and occurrence probability-area relationships can alter the form of species-area relationships. The plant community model estimated that only half of the species present in the regional species pool were encountered during the survey. The modeling framework we propose explicitly accounts for sampling processes so that ecological processes can be examined free of sampling artefacts. Our modeling approach is extensible and could be applied to a variety of study designs and allows the inclusion of additional environmental covariates.

Finite Sampling Exponential Bounds

Thesis (Ph.D.)--University of Washington, 2016-08

Falls Among The Elderly: Risk Factors In A Population-based Study.

The aim of the present study was to identify factors associated with the occurrence of falls among elderly adults in a population-based study (ISACamp 2008). A population-based cross-sectional study was carried out with two-stage cluster sampling. The sample was composed of 1,520 elderly adults living in the urban area of the city of Campinas, São Paulo, Brazil. The occurrence of falls was analyzed based on reports of the main accident occurred in the previous 12 months. Data on socioeconomic/demographic factors and adverse health conditions were tested for possible associations with the outcome. Prevalence ratios (PR) were estimated and adjusted for gender and age using the Poisson multiple regression analysis. Falls were more frequent, after adjustment for gender and age, among female elderly participants (PR = 2.39; 95% confidence interval (95% CI) 1.47 - 3.87), elderly adults (80 years old and older) (PR = 2.50; 95% CI 1.61 - 3.88), widowed (PR = 1.74; 95% CI 1.04 - 2.89) and among elderly adults who had rheumatism/arthritis/arthrosis (PR = 1.58; 95% CI 1.00 - 2.48), osteoporosis (PR = 1.71; 95% CI 1.18 - 2.49), asthma/bronchitis/emphysema (PR = 1,73; 95% CI 1.09 - 2.74), headache (PR = 1.59; 95% CI 1.07 - 2.38), mental common disorder (PR = 1.72; 95% CI 1.12 - 2.64), dizziness (PR = 2.82; 95% CI 1.98 - 4.02), insomnia (PR = 1.75; 95% CI 1.16 - 2.65), use of multiple medications (five or more) (PR = 2.50; 95% CI 1.12 - 5.56) and use of cane/walker (PR = 2.16; 95% CI 1.19 - 3,93). The present study shows segments of the elderly population who are more prone to falls through the identification of factors associated with this outcome. The findings can contribute to the planning of public health policies and programs addressed to the prevention of falls.

Factors associated with the prevalence of hypertension and control practices among elderly residents of São Paulo city, Brazil

The aim of this study was to analyze the prevalence of hypertension and control practices among the elderly. The survey analyzed data from 872 elderly people in São Paulo, Brazil, through a cluster sampling, stratified according to education and income. A Poisson multiple regression model checked for the existence of factors associated with hypertension. The prevalence of self-reported hypertension among the elderly was 46.9%. Variables associated with hypertension were self-rated health, alcohol consumption, gender, and hospitalization in the last year, regardless of age. The three most common measures taken to control hypertension, but only rarely, are oral medication, routine salt-free diet and physical activity. Lifestyle and socioeconomic status did not affect the practice of control, but knowledge about the importance of physical activity was higher among those older people with higher education and greater income. The research suggests that health policies that focus on primary care to encourage lifestyle changes among the elderly are necessary.

Spatial pattern detection modeling of thrips (Thrips tabaci) on onion fields

Onion (Allium cepa) is one of the most cultivated and consumed vegetables in Brazil and its importance is due to the large laborforce involved. One of the main pests that affect this crop is the Onion Thrips (Thrips tabaci), but the spatial distribution of this insect, although important, has not been considered in crop management recommendations, experimental planning or sampling procedures. Our purpose here is to consider statistical tools to detect and model spatial patterns of the occurrence of the onion thrips. In order to characterize the spatial distribution pattern of the Onion Thrips a survey was carried out to record the number of insects in each development phase on onion plant leaves, on different dates and sample locations, in four rural properties with neighboring farms under different infestation levels and planting methods. The Mantel randomization test proved to be a useful tool to test for spatial correlation which, when detected, was described by a mixed spatial Poisson model with a geostatistical random component and parameters allowing for a characterization of the spatial pattern, as well as the production of prediction maps of susceptibility to levels of infestation throughout the area.

Temporal variation of the phytoplankton community at short sampling intervals in the Mundaú reservoir, Northeastern Brazil

The aim of this study was to determine how abiotic factors drive the phytoplankton community in a water supply reservoir within short sampling intervals. Samples were collected at the subsurface (0.1 m) and bottom of limnetic (8 m) and littoral (2 m) zones in both the dry and rainy seasons. The following abiotic variables were analyzed: water temperature, dissolved oxygen, electrical conductivity, total dissolved solids, turbidity, pH, total nitrogen, nitrite, nitrate, total phosphorus, total dissolved phosphorus and orthophosphate. Phytoplankton biomass was determined from biovolume values. The role abiotic variables play in the dynamics of phytoplankton species was determined by means of Canonical Correspondence Analysis. Algae biomass ranged from 1.17×10(4) to 9.21×10(4) µg.L-1; cyanobacteria had biomass values ranging from 1.07×10(4) to 8.21×10(4) µg.L-1. High availability of phosphorous, nitrogen limitation, alkaline pH and thermal stability all favored cyanobacteria blooms, particularly during the dry season. Temperature, pH, total phosphorous and turbidity were key factors in characterizing the phytoplankton community between sampling times and stations. Of the species studied, Cylindrospermopsis raciborskii populations were dominant in the phytoplankton in both the dry and rainy seasons. We conclude that the phytoplankton was strongly influenced by abiotic variables, particularly in relation to seasonal distribution patterns.

Autopercepção da saúde bucal em idosos e fatores associados em Campinas, SP, 2008-2009

OBJETIVO: Descrever a autopercepção de saúde bucal em idosos e analisar fatores sociodemográficos e clínicos associados. MÉTODOS: Estudo transversal com 876 participantes em amostra representativa de idosos (65 anos ou mais) de Campinas, SP, em 2008-2009. Os exames odontológicos seguiram critérios padronizados pela Organização Mundial da Saúde para levantamentos epidemiológicos de saúde bucal. A autopercepção da saúde bucal foi avaliada pelo índice Geriatric Oral Health Assessment Index (GOHAI). Os indivíduos foram classificados segundo características sociodemográficas, odontológicas e prevalência de fragilidade biológica. O estudo de associações utilizou análise de regressão de Poisson; a análise considerou os pesos amostrais e a estrutura complexa da amostra por conglomerados. RESULTADOS: A média de idade dos indivíduos foi de 72,8 anos; 70,1% eram mulheres. A proporção de indivíduos com mais de 20 dentes presentes foi 17,2%; 38,2% usavam prótese dentária total em ambos os arcos; 8,5% necessitavam desse recurso em ao menos um arco dentário. Em média, o índice GOHAI foi elevado: 33,9 (máximo possível 36,0). Manter 20 dentes ou mais, usar prótese total nos dois arcos, não necessitar desse tratamento, não apresentar alterações de mucosa oral e não apresentar fragilidade biológica foram os fatores significantemente associados com melhor autopercepção de saúde bucal (p < 0,05). CONCLUSÕES: A avaliação de autopercepção em saúde bucal permitiu identificar os principais fatores associados a esse desfecho. Esse instrumento pode contribuir para o planejamento de serviços odontológicos, orientando estratégias de promoção em saúde voltadas à melhora da qualidade de vida das pessoas desse grupo etário.