Fitting genetic models to twin data with binary and ordered categorical responses: A comparison of structural equation modelling and Bayesian hierarchical models

We compare Bayesian methodology utilizing free-ware BUGS (Bayesian Inference Using Gibbs Sampling) with the traditional structural equation modelling approach based on another free-ware package, Mx. Dichotomous and ordinal (three category) twin data were simulated according to different additive genetic and common environment models for phenotypic variation. Practical issues are discussed in using Gibbs sampling as implemented by BUGS to fit subject-specific Bayesian generalized linear models, where the components of variation may be estimated directly. The simulation study (based on 2000 twin pairs) indicated that there is a consistent advantage in using the Bayesian method to detect a correct model under certain specifications of additive genetics and common environmental effects. For binary data, both methods had difficulty in detecting the correct model when the additive genetic effect was low (between 10 and 20%) or of moderate range (between 20 and 40%). Furthermore, neither method could adequately detect a correct model that included a modest common environmental effect (20%) even when the additive genetic effect was large (50%). Power was significantly improved with ordinal data for most scenarios, except for the case of low heritability under a true ACE model. We illustrate and compare both methods using data from 1239 twin pairs over the age of 50 years, who were registered with the Australian National Health and Medical Research Council Twin Registry (ATR) and presented symptoms associated with osteoarthritis occurring in joints of the hand.

Bridging the gap between different statistical approaches: An integrated framework for modelling

This paper proposes a template for modelling complex datasets that integrates traditional statistical modelling approaches with more recent advances in statistics and modelling through an exploratory framework. Our approach builds on the well-known and long standing traditional idea of 'good practice in statistics' by establishing a comprehensive framework for modelling that focuses on exploration, prediction, interpretation and reliability assessment, a relatively new idea that allows individual assessment of predictions. The integrated framework we present comprises two stages. The first involves the use of exploratory methods to help visually understand the data and identify a parsimonious set of explanatory variables. The second encompasses a two step modelling process, where the use of non-parametric methods such as decision trees and generalized additive models are promoted to identify important variables and their modelling relationship with the response before a final predictive model is considered. We focus on fitting the predictive model using parametric, non-parametric and Bayesian approaches. This paper is motivated by a medical problem where interest focuses on developing a risk stratification system for morbidity of 1,710 cardiac patients given a suite of demographic, clinical and preoperative variables. Although the methods we use are applied specifically to this case study, these methods can be applied across any field, irrespective of the type of response.

Indirect evidence of density-dependent population regulation in Aponomma hydrosauri (Acari : Ixodidae), an ectoparasite of reptiles

The extent to which density-dependent processes regulate natural populations is the subject of an ongoing debate. We contribute evidence to this debate showing that density-dependent processes influence the population dynamics of the ectoparasite Aponomma hydrosauri (Acari: Ixodidae), a tick species that infests reptiles in Australia. The first piece of evidence comes from an unusually long-term dataset on the distribution of ticks among individual hosts. If density-dependent processes are influencing either host mortality or vital rates of the parasite population, and those distributions can be approximated with negative binomial distributions, then general host-parasite models predict that the aggregation coefficient of the parasite distribution will increase with the average intensity of infections. We fit negative binomial distributions to the frequency distributions of ticks on hosts, and find that the estimated aggregation coefficient k increases with increasing average tick density. This pattern indirectly implies that one or more vital rates of the tick population must be changing with increasing tick density, because mortality rates of the tick's main host, the sleepy lizard, Tiliqua rugosa, are unaffected by changes in tick burdens. Our second piece of evidence is a re-analysis of experimental data on the attachment success of individual ticks to lizard hosts using generalized linear modelling. The probability of successful engorgement decreases with increasing numbers of ticks attached to a host. This is direct evidence of a density-dependent process that could lead to an increase in the aggregation coefficient of tick distributions described earlier. The population-scale increase in the aggregation coefficient is indirect evidence of a density-dependent process or processes sufficiently strong to produce a population-wide pattern, and thus also likely to influence population regulation. The direct observation of a density-dependent process is evidence of at least part of the responsible mechanism.

Air pollution effects on myocardial infarction

OBJECTIVE: Myocardial infarction is an acute and severe cardiovascular disease that generally leads to patient admissions to intensive care units and few cases are initially admitted to infirmaries. The objective of the study was to assess whether estimates of air pollution effects on myocardial infarction morbidity are modified by the source of health information. METHODS: The study was carried out in hospitals of the Brazilian Health System in the city of São Paulo, Southern Brazil. A time series study (1998-1999) was performed using two outcomes: infarction admissions to infirmaries and to intensive care units, both for people older than 64 years of age. Generalized linear models controlling for seasonality (long and short-term trends) and weather were used. The eight-day cumulative effects of air pollutants were assessed using third degree polynomial distributed lag models. RESULTS: Almost 70% of daily hospital admissions due to myocardial infarction were to infirmaries. Despite that, the effects of air pollutants on infarction were higher for intensive care units admissions. All pollutants were positively associated with the study outcomes but SO2 presented the strongest statistically significant association. An interquartile range increase on SO2 concentration was associated with increases of 13% (95% CI: 6-19) and 8% (95% CI: 2-13) of intensive care units and infirmary infarction admissions, respectively. CONCLUSIONS: It may be assumed there is a misclassification of myocardial infarction admissions to infirmaries leading to overestimation. Also, despite the absolute number of events, admissions to intensive care units data provides a more adequate estimate of the magnitude of air pollution effects on infarction admissions.

Habitat predictive modelling of demersal fish species in the Azores

Dissertação de Mestrado, Estudos Integrados dos Oceanos, 25 de Março de 2013, Universidade dos Açores.

O impacto do excesso de gordura corporal sobre a remodelação óssea de adolescentes

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Factors associated with vaccination coverage in children < 5 years in Angola

OBJECTIVE To analyze vaccination coverage and factors associated with a complete immunization scheme in children < 5 years old. METHODS This cross-sectional household census survey evaluated 1,209 children < 5 years old living in Bom Jesus, Angola, in 2010. Data were obtained from interviews, questionnaires, child immunization histories, and maternal health histories. The statistical analysis used generalized linear models, in which the dependent variable followed a binary distribution (vaccinated, unvaccinated) and the association function was logarithmic and had the children’s individual, familial, and socioeconomic factors as independent variables. RESULTS Vaccination coverage was 37.0%, higher in children < 1 year (55.0%) and heterogeneous across neighborhoods; 52.0% of children of both sexes had no immunization records. The prevalence rate of vaccination significantly varied according to child age, mother’s level of education, family size, ownership of household appliances, and destination of domestic waste. CONCLUSIONS Vulnerable groups with vaccination coverage below recommended levels continue to be present. Some factors indicate inequalities that represent barriers to full immunization, indicating the need to implement more equitable policies. The knowledge of these factors contributes to planning immunization promotion measures that focus on the most vulnerable groups.

Cervical cancer screening coverage in a high-incidence region

OBJECTIVE To analyze the coverage of a cervical cancer screening program in a city with a high incidence of the disease in addition to the factors associated with non-adherence to the current preventive program.METHODS A cross-sectional study based on household surveys was conducted. The sample was composed of women between 25 and 59 years of age of the city of Boa Vista, RR, Northern Brazil who were covered by the cervical cancer screening program. The cluster sampling method was used. The dependent variable was participation in a women’s health program, defined as undergoing at least one Pap smear in the 36 months prior to the interview; the explanatory variables were extracted from individual data. A generalized linear model was used.RESULTS 603 women were analyzed, with an mean age of 38.2 years (SD = 10.2). Five hundred and seventeen women underwent the screening test, and the prevalence of adherence in the last three years was up to 85.7% (95%CI 82.5;88.5). A high per capita household income and recent medical consultation were associated with the lower rate of not being tested in multivariate analysis. Disease ignorance, causes, and prevention methods were correlated with chances of non-adherence to the screening system; 20.0% of the women were reported to have undergone opportunistic and non-routine screening.CONCLUSIONS The informed level of coverage is high, exceeding the level recommended for the control of cervical cancer. The preventive program appears to be opportunistic in nature, particularly for the most vulnerable women (with low income and little information on the disease). Studies on the diagnostic quality of cervicovaginal cytology and therapeutic schedules for positive cases are necessary for understanding the barriers to the control of cervical cancer.

Occurrence of severe dengue in Rio de Janeiro: an ecological study

Introduction This study aimed to analyze the relationship between the incidence of severe dengue during the 2008 epidemic in Rio de Janeiro, Brazil, and socioeconomic indicators, as well as indicators of health service availability and previous circulation of the dengue virus serotype-3 (DENV-3). Methods In this ecological study, the units of analysis were the districts of Rio de Janeiro. The data were incorporated into generalized linear models, and the incidence of severe dengue in each district was the outcome variable. Results The districts with more cases of dengue fever in the 2001 epidemic and a higher percentage of residents who declared their skin color or race as black had higher incidence rates of severe dengue in the 2008 epidemic [incidence rate ratio (IRR)= 1.21; 95% confidence interval (95%CI)= 1.05-1.40 and IRR= 1.34; 95%CI= 1.16-1.54, respectively]. In contrast, the districts with Family Health Strategy (FHS) clinics were more likely to have lower incidence rates of severe dengue in the 2008 epidemic (IRR= 0.81; 95%CI= 0.70-0.93). Conclusions At the ecological level, our findings suggest the persistence of health inequalities in this region of Brazil that are possibly due to greater social vulnerability among the self-declared black population. Additionally, the protective effect of FHS clinics may be due to the ease of access to other levels of care in the health system or to a reduced vulnerability to dengue transmission that is afforded by local practices to promote health.

Resolução da ambivalência em terapia narrativa

Dissertação de mestrado integrado em Psicologia

Momentos de inovação: emergência e evolução das protonarrativas em terapia narrativa

Dissertação de mestrado em Psicologia Aplicada

Artificial nest predation rates vary depending on visibility in the eastern Brazilian Amazon

Observational and experimental studies have shown that increased concealment of bird nests reduces nest predation rates. The objective of the present study was to evaluate differences in predation rates between two experimental manipulations of artificial ground nests (i.e., clearing an area around the artificial nest or leaving it as natural as possible), and test whether environmental variables also affected nest predation in an undisturbed area of Amazonian forest in eastern Brazil. A generalized linear model was used to examine the influence of five variables (manipulation type, perpendicular distance from the main trail, total basal area of trees surrounding the nest site, understorey density, and liana quantity) on nest predation rates. Model results, showed that manipulation type was the only variable that significantly affected nest predation rates. Thus, to avoid systematic biases, the influence of nest site manipulation must be taken into consideration when conducting experiments with artificial nests.

Avaliação do impacto na circulação e risco de acidente em estradas em obras

Dissertação de mestrado integrado em Engenharia Civil

The ghost of past species occurrence: improving species distribution models for presence-only data

1. Model-based approaches have been used increasingly in conservation biology over recent years. Species presence data used for predictive species distribution modelling are abundant in natural history collections, whereas reliable absence data are sparse, most notably for vagrant species such as butterflies and snakes. As predictive methods such as generalized linear models (GLM) require absence data, various strategies have been proposed to select pseudo-absence data. However, only a few studies exist that compare different approaches to generating these pseudo-absence data. 2. Natural history collection data are usually available for long periods of time (decades or even centuries), thus allowing historical considerations. However, this historical dimension has rarely been assessed in studies of species distribution, although there is great potential for understanding current patterns, i.e. the past is the key to the present. 3. We used GLM to model the distributions of three 'target' butterfly species, Melitaea didyma, Coenonympha tullia and Maculinea teleius, in Switzerland. We developed and compared four strategies for defining pools of pseudo-absence data and applied them to natural history collection data from the last 10, 30 and 100 years. Pools included: (i) sites without target species records; (ii) sites where butterfly species other than the target species were present; (iii) sites without butterfly species but with habitat characteristics similar to those required by the target species; and (iv) a combination of the second and third strategies. Models were evaluated and compared by the total deviance explained, the maximized Kappa and the area under the curve (AUC). 4. Among the four strategies, model performance was best for strategy 3. Contrary to expectations, strategy 2 resulted in even lower model performance compared with models with pseudo-absence data simulated totally at random (strategy 1). 5. Independent of the strategy model, performance was enhanced when sites with historical species presence data were not considered as pseudo-absence data. Therefore, the combination of strategy 3 with species records from the last 100 years achieved the highest model performance. 6. Synthesis and applications. The protection of suitable habitat for species survival or reintroduction in rapidly changing landscapes is a high priority among conservationists. Model-based approaches offer planning authorities the possibility of delimiting priority areas for species detection or habitat protection. The performance of these models can be enhanced by fitting them with pseudo-absence data relying on large archives of natural history collection species presence data rather than using randomly sampled pseudo-absence data.

Modelling the distribution of bats in relation to landscape structure in a temperate mountain environment

1. Landscape modification is often considered the principal cause of population decline in many bat species. Thus, schemes for bat conservation rely heavily on knowledge about species-landscape relationships. So far, however, few studies have quantified the possible influence of landscape structure on large-scale spatial patterns in bat communities. 2. This study presents quantitative models that use landscape structure to predict (i) spatial patterns in overall community composition and (ii) individual species' distributions through canonical correspondence analysis and generalized linear models, respectively. A geographical information system (GIS) was then used to draw up maps of (i) overall community patterns and (ii) distribution of potential species' habitats. These models relied on field data from the Swiss Jura mountains. 3. Fight descriptors of landscape structure accounted for 30% of the variation in bat community composition. For some species, more than 60% of the variance in distribution could be explained by landscape structure. Elevation, forest or woodland cover, lakes and suburbs, were the most frequent predictors. 4. This study shows that community composition in bats is related to landscape structure through species-specific relationships to resources. Due to their nocturnal activities and the difficulties of remote identification, a comprehensive bat census is rarely possible, and we suggest that predictive modelling of the type described here provides an indispensable conservation tool.