Prediction and improved estimation in linear models /

Includes index.

How well can models predict changes in species distributions? A 13-year-old otter model revisited

Species distribution and ecological niche models are increasingly used in biodiversity management and conservation. However, one thing that is important but rarely done is to follow up on the predictive performance of these models over time, to check if their predictions are fulfilled and maintain accuracy, or if they apply only to the set in which they were produced. In 2003, a distribution model of the Eurasian otter (Lutra lutra) in Spain was published, based on the results of a country-wide otter survey published in 1998. This model was built with logistic regression of otter presence-absence in UTM 10 km2 cells on a diverse set of environmental, human and spatial variables, selected according to statistical criteria. Here we evaluate this model against the results of the most recent otter survey, carried out a decade later and after a significant expansion of the otter distribution area in this country. Despite the time elapsed and the evident changes in this species’ distribution, the model maintained a good predictive capacity, considering both discrimination and calibration measures. Otter distribution did not expand randomly or simply towards vicinity areas,m but specifically towards the areas predicted as most favourable by the model based on data from 10 years before. This corroborates the utility of predictive distribution models, at least in the medium term and when they are made with robust methods and relevant predictor variables.

Age-related differences in exercise and quality of life among breast cancer survivors

Purpose: Physical activity has become a focus of cancer recovery research as it has the potential to reduce treatment-related burden and optimize health-related quality of life (HRQoL). However, the potential for physical activity to influence recovery may be age-dependent. This paper describes physical activity levels and HRQoL among younger and older women after surgery for breast cancer and explores the correlates of physical inactivity. Methods: A population-based sample of breast cancer patients diagnosed in South-East Queensland, Australia, (n=287) were assessed once every three months, from 6 to 18 months post-surgery. The Functional Assessment of Cancer Therapy-Breast questionnaire (FACTB+4) and items from the Behavioral Risk Factor Surveillance System (BRFSS) questionnaire were used to measure HRQoL and physical activity, respectively. Physical activity was assigned metabolic equivalent task (MET) values, and categorized as < 3, 3 to 17.9 and 18+ MET-hours/weeks. Descriptive statistics, generalized linear models with age stratification (<50 years versus 50+ years), and logistic regression were used for analyses (p=0.05, two-tailed). Results: Younger women who engaged in 3 or more MET-hours/week of physical activity reported a higher HRQoL at 18 months compared to their more sedentary counterparts (p<0.05). Older women reported similar HRQoL irrespective of activity level and consistently reported clinically higher HRQoL than younger women. Increasing age, being overweight or obese, and restricting use of the treated side at six months post-surgery increased the likelihood of sedentary behavior (OR>3, p<0.05). Conclusions: Age influences the potential to observe HRQoL benefits related to physical activity participation. These results also provide relevant information for the design of exercise interventions for breast cancer survivors and highlights that some groups of women are at greater risk of long-term sedentary behavior.

Analysing Seasonal Health Data

Seasonal patterns have been found in a remarkable range of health conditions, including birth defects, respiratory infections and cardiovascular disease. Accurately estimating the size and timing of seasonal peaks in disease incidence is an aid to understanding the causes and possibly to developing interventions. With global warming increasing the intensity of seasonal weather patterns around the world, a review of the methods for estimating seasonal effects on health is timely. This is the first book on statistical methods for seasonal data written for a health audience. It describes methods for a range of outcomes (including continuous, count and binomial data) and demonstrates appropriate techniques for summarising and modelling these data. It has a practical focus and uses interesting examples to motivate and illustrate the methods. The statistical procedures and example data sets are available in an R package called ‘season’. Adrian Barnett is a senior research fellow at Queensland University of Technology, Australia. Annette Dobson is a Professor of Biostatistics at The University of Queensland, Australia. Both are experienced medical statisticians with a commitment to statistical education and have previously collaborated in research in the methodological developments and applications of biostatistics, especially to time series data. Among other projects, they worked together on revising the well-known textbook "An Introduction to Generalized Linear Models," third edition, Chapman Hall/CRC, 2008. In their new book they share their knowledge of statistical methods for examining seasonal patterns in health.

Design of experiments for bivariate binary responses modelled by Copula functions

Optimal design for generalized linear models has primarily focused on univariate data. Often experiments are performed that have multiple dependent responses described by regression type models, and it is of interest and of value to design the experiment for all these responses. This requires a multivariate distribution underlying a pre-chosen model for the data. Here, we consider the design of experiments for bivariate binary data which are dependent. We explore Copula functions which provide a rich and flexible class of structures to derive joint distributions for bivariate binary data. We present methods for deriving optimal experimental designs for dependent bivariate binary data using Copulas, and demonstrate that, by including the dependence between responses in the design process, more efficient parameter estimates are obtained than by the usual practice of simply designing for a single variable only. Further, we investigate the robustness of designs with respect to initial parameter estimates and Copula function, and also show the performance of compound criteria within this bivariate binary setting.

Chronological age and age-related cognitive deficits are associated with an increase in multiple types of driving errors in late life

Objective: Older driver research has mostly focused on identifying that small proportion of older drivers who are unsafe. Little is known about how normal cognitive changes in aging affect driving in the wider population of adults who drive regularly. We evaluated the association of cognitive function and age, with driving errors. Method: A sample of 266 drivers aged 70 to 88 years were assessed on abilities that decline in normal aging (visual attention, processing speed, inhibition, reaction time, task switching) and the UFOV® which is a validated screening instrument for older drivers. Participants completed an on-road driving test. Generalized linear models were used to estimate the associations of cognitive factor with specific driving errors and number of errors in self-directed and instructor navigated conditions. Results: All errors types increased with chronological age. Reaction time was not associated with driving errors in multivariate analyses. A cognitive factor measuring Speeded Selective Attention and Switching was uniquely associated with the most errors types. The UFOV predicted blindspot errors and errors on dual carriageways. After adjusting for age, education and gender the cognitive factors explained 7% of variance in the total number of errors in the instructor navigated condition and 4% of variance in the self-navigated condition. Conclusion: We conclude that among older drivers errors increase with age and are associated with speeded selective attention particularly when that requires attending to the stimuli in the periphery of the visual field, task switching, errors inhibiting responses and visual discrimination. These abilities should be the target of cognitive training.

Co-morbidity of depression and anxiety in common age-related eye diseases : a population-based study of 662 adults

This study examined the prevalence of co-morbid age-related eye disease and symptoms of depression and anxiety in late life, and the relative roles of visual function and disease in explaining symptoms of depression and anxiety. A community-based sample of 662 individuals aged over 70 years was recruited through the electoral roll. Vision was measured using a battery of tests including high and low contrast visual acuity, contrast sensitivity, motion sensitivity, stereoacuity, Useful Field of View, and visual fields. Depression and anxiety symptoms were measured using the Goldberg scales. The prevalence of self-reported eye disease [cataract, glaucoma, or age-related macular degeneration (AMD)] in the sample was 43.4%, with 7.7% reporting more than one form of ocular pathology. Of those with no eye disease, 3.7% had clinically significant depressive symptoms. This rate was 6.7% among cataract patients, 4.3% among those with glaucoma, and 10.5% for AMD. Generalized linear models adjusting for demographics, general health, treatment, and disability examined self-reported eye disease and visual function as correlates of depression and anxiety. Depressive symptoms were associated with cataract only, AMD, comorbid eye diseases and reduced low contrast visual acuity. Anxiety was significantly associated with self-reported cataract, and reduced low contrast visual acuity, motion sensitivity and contrast sensitivity. We found no evidence for elevated rates of depressive or anxiety symptoms associated with self-reported glaucoma. The results support previous findings of high rates of depression and anxiety in cataract and AMD, and in addition show that mood and anxiety are associated with objective measures of visual function independently of self-reported eye disease. The findings have implications for the assessment and treatment of mental health in the context of late-life visual impairment...

Quality of life of women with lower limb swelling or lymphedema 3-5 years following endometrial cancer

Objective To quantitatively assess and compare the quality of life (QoL) of women with a self-reported diagnosis of lower limb lymphedema (LLL), to women with lower limb swelling (LLS), and to women without LLL or LLS following treatment for endometrial cancer. Methods 1399 participants in the Australian National Endometrial Cancer Study were sent a follow-up questionnaire 3–5 years after diagnosis. Women were asked if they had experienced swelling in the lower limbs and, if so, whether they had received a diagnosis of lymphedema by a health professional. The 639 women who responded were categorised as: Women with LLL (n = 68), women with LLS (n = 177) and women without LLL or LLS (n = 394). Multivariable-adjusted generalized linear models were used to compare women’s physical and mental QoL by LLL status. Results On average, women were 65 years of age and 4 years after diagnosis. Women with LLL had clinically lower physical QoL (M= 41.8, SE= 1.4) than women without LLL or LLS (M= 45.1, SE= 0.8, p = .07), however, their mental QoL was within the normative range (M= 49.6; SE= 1.1 p = 1.0). Women with LLS had significantly lower physical (M= 41.0, SE= 1.0, p = .003) and mental QoL (M= 46.8; SE= 0.8, p < .0001) than women without LLL or LLS (Mental QoL: M= 50.6, SE= 0.8). Conclusion Although LLL was associated with reductions in physical QoL, LLS was related to reductions in both physical and mental QoL 3-5 years after cancer treatment. Early referral to evidence-based lymphedema programs may prevent long-term impairments to women’s QoL.

Violence in and around entertainment districts : a longitudinal analysis of the impact of late-night lockout legislation

Violence in entertainment districts is a major problem across urban landscapes throughout the world. Research shows that licensed premises are the third most common location for homicides and serious assaults, accounting for one in ten fatal and nonfatal assaults. One class of interventions that aims to reduce violence in entertainment districts involves the use of civil remedies: a group of strategies that use civil or regulatory measures as legal “levers” to reduce problem behavior. One specific civil remedy used to reduce problematic behavior in entertainment districts involves manipulation of licensed premise trading hours. This article uses generalized linear models to analyze the impact of lockout legislation on recorded violent offences in two entertainment districts in the Australian state of Queensland. Our research shows that 3 a.m. lockout legislation led to a direct and significant reduction in the number of violent incidents inside licensed premises. Indeed, the lockouts cut the level of violent crime inside licensed premises by half. Despite these impressive results for the control of violence inside licensed premises, we found no evidence that the lockout had any impact on violence on streets and footpaths outside licensed premises that were the site for more than 80 percent of entertainment district violence. Overall, however, our analysis suggests that lockouts are an important mechanism that helps to control the level of violence inside licensed premises but that finely grained contextual responses to alcohol-related problems are needed rather than one-size-fits-all solutions.

Developments of the total entropy utility function for the dual purpose of model discrimination and parameter estimation in Bayesian design

The total entropy utility function is considered for the dual purpose of Bayesian design for model discrimination and parameter estimation. A sequential design setting is proposed where it is shown how to efficiently estimate the total entropy utility for a wide variety of data types. Utility estimation relies on forming particle approximations to a number of intractable integrals which is afforded by the use of the sequential Monte Carlo algorithm for Bayesian inference. A number of motivating examples are considered for demonstrating the performance of total entropy in comparison to utilities for model discrimination and parameter estimation. The results suggest that the total entropy utility selects designs which are efficient under both experimental goals with little compromise in achieving either goal. As such, the total entropy utility is advocated as a general utility for Bayesian design in the presence of model uncertainty.

Spectrum of ankylosing spondylitis in Portugal. Development of BASDAI, BASFI, BASMI and mSASSS reference centile charts

The availability of population-specific normative data regarding disease severity measures is essential for patient assessment. The goals of the current study were to characterize the pattern of ankylosing spondylitis (AS) in Portuguese patients and to develop reference centile charts for BASDAI, BASFI, BASMI and mSASSS, the most widely used assessment tools in AS. AS cases were recruited from hospital outpatient clinics, with AS defined according to the modified New York criteria. Demographic and clinical data were recorded. All radiographs were evaluated by two independent experienced readers. Centile charts for BASDAI, BASFI, BASMI and mSASSS were constructed for both genders, using generalized linear models and regression models with duration of disease as independent variable. A total of 369 patients (62.3% male, mean ± (SD) age 45.4 ± 13.2 years, mean ± (SD) disease duration 11.4 ± 10.5 years, 70.7% B27-positive) were included. Family history of AS in a first-degree relative was reported in 17.6% of the cases. Regarding clinical disease pattern, at the time of assessment 42.3% had axial disease, 2.4% peripheral disease, 40.9% mixed disease and 7.1% isolated enthesopatic disease. Anterior uveitis (33.6%) was the most common extra-articular manifestation. The centile charts suggest that females reported greater disease activity and more functional impairment than males but had lower BASMI and mSASSS scores. Data collected through this study provided a demographic and clinical profile of patients with AS in Portugal. The development of centile charts constitutes a useful tool to assess the change of disease pattern over time and in response to therapeutic interventions.

Generalizing the use of geographical weights in biodiversity modelling

Aim Determining how ecological processes vary across space is a major focus in ecology. Current methods that investigate such effects remain constrained by important limiting assumptions. Here we provide an extension to geographically weighted regression in which local regression and spatial weighting are used in combination. This method can be used to investigate non-stationarity and spatial-scale effects using any regression technique that can accommodate uneven weighting of observations, including machine learning. Innovation We extend the use of spatial weights to generalized linear models and boosted regression trees by using simulated data for which the results are known, and compare these local approaches with existing alternatives such as geographically weighted regression (GWR). The spatial weighting procedure (1) explained up to 80% deviance in simulated species richness, (2) optimized the normal distribution of model residuals when applied to generalized linear models versus GWR, and (3) detected nonlinear relationships and interactions between response variables and their predictors when applied to boosted regression trees. Predictor ranking changed with spatial scale, highlighting the scales at which different species–environment relationships need to be considered. Main conclusions GWR is useful for investigating spatially varying species–environment relationships. However, the use of local weights implemented in alternative modelling techniques can help detect nonlinear relationships and high-order interactions that were previously unassessed. Therefore, this method not only informs us how location and scale influence our perception of patterns and processes, it also offers a way to deal with different ecological interpretations that can emerge as different areas of spatial influence are considered during model fitting.