On Wald Residuals in Generalized Linear Models

A rigorous asymptotic theory for Wald residuals in generalized linear models is not yet available. The authors provide matrix formulae of order O(n(-1)), where n is the sample size, for the first two moments of these residuals. The formulae can be applied to many regression models widely used in practice. The authors suggest adjusted Wald residuals to these models with approximately zero mean and unit variance. The expressions were used to analyze a real dataset. Some simulation results indicate that the adjusted Wald residuals are better approximated by the standard normal distribution than the Wald residuals.

Awareness and knowledge of HPV, cervical cancer, and vaccines in young women after first delivery in São Paulo, Brazil - a cross-sectional study

Abstract Background The success of HPV vaccination programs will require awareness regarding HPV associated diseases and the benefits of HPV vaccination for the general population. The aim of this study was to assess the level of awareness and knowledge of human papillomavirus (HPV) infection, cervical cancer prevention, vaccines, and factors associated with HPV awareness among young women after birth of the first child. Methods This analysis is part of a cross-sectional study carried out at Hospital Maternidade Leonor Mendes de Barros, a large public maternity hospital in Sao Paulo. Primiparous women (15-24 years) who gave birth in that maternity hospital were included. A questionnaire that included questions concerning knowledge of HPV, cervical cancer, and vaccines was applied. To estimate the association of HPV awareness with selected factors, prevalence ratios (PR) were estimated using a generalized linear model (GLM). Results Three hundred and one primiparous women were included; 37% of them reported that they "had ever heard about HPV", but only 19% and 7%, respectively, knew that HPV is a sexually transmitted infection (STI) and that it can cause cervical cancer. Seventy-four percent of interviewees mentioned the preventive character of vaccines and all participants affirmed that they would accept HPV vaccination after delivery. In the multivariate analysis, only increasing age (P for trend = 0.021) and previous STI (P < 0.001) were factors independently associated with HPV awareness ("had ever heard about HPV"). Conclusions This survey indicated that knowledge about the association between HPV and cervical cancer among primiparous young women is low. Therefore, these young low-income primiparous women could benefit greatly from educational interventions to encourage primary and secondary cervical cancer prevention programs.

Catch rates of artisanal fisheries targeting Dasyatis guttata (Bloch & Schneider, 1801) in Todos os Santos Bay, northeastern Brazil

Dasyatis guttata has been target of artisanal fisheries in the coast of Bahia (Northeast Brazil) mainly by “arraieira” (gillnet) and “grozeira” (bottom long-line), but until now there is no stock assessment study. One of the important data for this knowledge is reliable indices of abundance. The aims of the present work are to: (1) estimate the best predictor for relative abundance (catch-per-unit-of-effort, CPUE), examining whether catch (production – kg) was related to: soak time of the gear, size of the gillnet or number of hooks, applying generalized linear model (GLM); (2) estimate the annual CPUE (kg/hooks and kg/m) averaged by gear; and (3) assess the temporal CPUE variance. Based on monthly sampling between January 2012 and January 2013, 222 landings by grozeira and 76 by arraiaiera were recorded in the two landing sites in Todos os Santos Bay, Bahia. A total of 14,550 kg (average = 44 kg/month) of D. guttata was captured. Models for both gears were highly significant (P < 0.0001). The analysis indicated that the most appropriate variable for CPUE analysis was the size of the gillnet (P < 0.001) and the number of hooks (P < 0.0001). Soak time of the gear was not significant for both gears (P = 0.4). High residual deviance expresses the complexity of the relations between ecosystem factors and other fisheries factors affecting relative abundance, which were not considered in this study. The average CPUE by grozeira was 6.39 kg/100 hooks ± 8.89 and by arraieira, 1.47 kg/100 m ± 1.66 over the year. Kruskal-Wallis test showed effect of the month on the mean grozeira CPUE (P = <0.001), but no effect (P = 0.096) on the mean arraieira CPUE. Grozeira CPUE values were highest in December and March, and lowest between May to August

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.

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.

On the relationship between squid and the environment: artisanal jigging for Loligo plei at Sao Sebastiao Island (24 degrees S), southeastern Brazil

The squid Loligo plei concentrates in the southeastern Brazil Bight, where it has traditionally supported small-scale fisheries around Sao Sebastiao Island (SSI). Sea surface temperature (SST), chlorophyll-a (Chl a), windspeed, wave height, rainfall, and lunar phase are related to fishing records and to the results of a survey of local fishers to investigate how they believe environmental variables might affect catches of L. plei. Daily fishery-dependent data over the years 2005-2009 were obtained from a fishing cooperative and were matched with satellite and meteorological forecast data. Generalized linear models were used to explore the significance of environmental variables in relation to variability in catch and catch per unit effort (cpue). Squid are fished with jigs in water shallower than 20 m, generally where SST is warmer and Chl a and windspeed are lower. Cpue and monthly catches decreased from 2005 to 2008, followed by a slight increase in 2009. The correlations between fishery and environmental data relate well to fishers` oceanological knowledge, underscoring the potential of incorporating such knowledge into evaluations of the fishery.

Tank experiments with borehole self-potential data

Despite their importance in the evaluation of petroleum and gas reservoirs, measurements of self-potential data under borehole conditions (well-logging) have found only minor applications in aquifer and waste-site characterization. This can be attributed to lower signals from the diffusion fronts in near-surface environments because measurements are made long after the drilling of the well, when concentration fronts are already disappearing. Proportionally higher signals arise from streaming potentials that prevent using simple interpretation models that assume signals from diffusion only. Our laboratory experiments found that dual-source self-potential signals can be described by a simple linear model, and that contributions (from diffusion and streaming potentials) can be isolated by slightly perturbing the borehole conditions. Perturbations are applied either by changing the concentration of the borehole-filling solution or its column height. Parameters useful for formation evaluation can be estimated from data measured during perturbations, namely, pore water resistivity, pressure drop across the borehole wall, and electrokinetic coupling parameter. These are important parameters to assess, respectively, water quality, aquifer lateral continuity, and interfacial properties of permeable formations.

Atmospheric pollution: influence on hospital admissions in paediatric rheumatic diseases

Objective: To investigate the lag structure effects from exposure to atmospheric pollution in acute outbursts in hospital admissions of paediatric rheumatic diseases (PRDs). Methods: Morbidity data were obtained from the Brazilian Hospital Information System in seven consecutive years, including admissions due to seven PRDs (juvenile idiopathic arthritis, systemic lupus erythematosus, dermatomyositis, Henoch-Schonlein purpura, polyarteritis nodosa, systemic sclerosis and ankylosing spondylitis). Cases with secondary diagnosis of respiratory diseases were excluded. Daily concentrations of inhaled particulate matter (PM10), sulphur dioxide (SO2) nitrogen dioxide (NO2), ozone (O-3) and carbon monoxide (CO) were evaluated. Generalized linear Poisson regression models controlling for short-term trend, seasonality, holidays, temperature and humidity were used. Lag structures and magnitude of air pollutants' effects were adopted to estimate restricted polynomial distributed lag models. Results: The total number of admissions due to acute outbursts PRD was 1,821. The SO2 interquartile range (7.79 mu g/m(3)) was associated with an increase of 1.98% (confidence interval 0.25-3.69) in the number of hospital admissions due to outcome studied after 14 days of exposure. This effect was maintained until day 17. Of note, the other pollutants, with the exception of O-3, showed an increase in the number of hospital admissions from the second week. Conclusion: This study is the first to demonstrate a delayed association between SO2 and PRD outburst, suggesting that oxidative stress reaction could trigger the inflammation of these diseases. Lupus (2012) 21, 526-533.

Bartlett corrections in Birnbaum-Saunders nonlinear regression models

Lemonte and Cordeiro [Birnbaum-Saunders nonlinear regression models, Comput. Stat. Data Anal. 53 (2009), pp. 4441-4452] introduced a class of Birnbaum-Saunders (BS) nonlinear regression models potentially useful in lifetime data analysis. We give a general matrix Bartlett correction formula to improve the likelihood ratio (LR) tests in these models. The formula is simple enough to be used analytically to obtain several closed-form expressions in special cases. Our results generalize those in Lemonte et al. [Improved likelihood inference in Birnbaum-Saunders regressions, Comput. Stat. DataAnal. 54 (2010), pp. 1307-1316], which hold only for the BS linear regression models. We consider Monte Carlo simulations to show that the corrected tests work better than the usual LR tests.

Temporal Trends of Asthma Mortality Rates in Brazil from 1980 to 2010

Objective. Mortality from asthma has varied among countries during the last several decades. This study aimed to identify temporal trends of asthma mortality in Brazil from 1980 to 2010. Method. We analyzed 6840 deaths of patients aged 5-34 years that occurred in Brazil with the underlying cause of asthma. We applied a log-linear model using Poisson regression to verify peaks and trends. We also calculated the point estimation and 95% confidence interval (CI 95%) of the annual percent change (APC) of the mortality rates, and the average annual percent change (AAPC) for 2001-2010. Results. A decline was observed from 1980 to 1992 [APC = -3.4 (-5.0 to -1.8)], followed by a nonsignificant rise until 1996 [APC = 6.8 (-1.4 to 15.6)], and a new downward trend from 1997 to 2010 [APC = -2.7 (-3.9 to -1.6)]. The APCs varied according to age strata: 5-14 years from 1980 to 2010 [-0.3 (-1.1 to 0.5)]; 15-24 years from 1980 to 1991 [-2.1 (-5.0 to 0.9)], from 1992 to 1996 [6.8 (-6.7 to 22.2)], and from 1997 to 2010 [-3.9 (-5.7 to -2.0)]; 24-25 years from 1980 to 1992 [-2.5 (-4.6 to -0.3)], from 1993 to 1995 [12.0 (-21.1 to 59.1)], and from 1996-2010 [-1.7 (-3.0 to -0.4)]. AAPC from 2001 to 2010 was -1.7 (-3.0 to -0.4); the decline for this period was significant for patients over 15 years old, women, and those living in the Southeast region. Conclusion. Asthma mortality rates in Brazil have been declining since the late 1990s.

Competing regression models for longitudinal data

The choice of an appropriate family of linear models for the analysis of longitudinal data is often a matter of concern for practitioners. To attenuate such difficulties, we discuss some issues that emerge when analyzing this type of data via a practical example involving pretestposttest longitudinal data. In particular, we consider log-normal linear mixed models (LNLMM), generalized linear mixed models (GLMM), and models based on generalized estimating equations (GEE). We show how some special features of the data, like a nonconstant coefficient of variation, may be handled in the three approaches and evaluate their performance with respect to the magnitude of standard errors of interpretable and comparable parameters. We also show how different diagnostic tools may be employed to identify outliers and comment on available software. We conclude by noting that the results are similar, but that GEE-based models may be preferable when the goal is to compare the marginal expected responses.

The impact of asymptomatic vertebral fractures on quality of life in older community-dwelling women: the Sao Paulo Ageing & Health Study

OBJECTIVES: The aim of this study was to investigate the impact of asymptomatic vertebral fractures on the quality of life in older women as part of the Sao Paulo Ageing & Health Study. METHODS: This study was a cross-sectional study with a random sample of 180 women 65 years of age or older with or without vertebral fractures. The Quality of Life Questionnaire of the European Foundation for Osteoporosis was administered to all subjects. Anthropometric data were obtained by physical examination, and the body mass index was calculated. Lateral thoracic and lumbar spine X-ray scans were obtained to identify asymptomatic vertebral fractures using a semi-quantitative method. RESULTS: Women with asymptomatic vertebral fractures had lower total scores [61.4(15.3) vs. 67.1(14.2), p = 0.03] and worse physical function domain scores [69.5(20.1) vs. 77.3(17.1), p = 0.02] for the Quality of Life Questionnaire of the European Foundation for Osteoporosis compared with women without fractures. The total score of this questionnaire was also worse in women classified as obese than in women classified as overweight or normal. High physical activity was related to a better total score for this questionnaire (p = 0.01). Likewise, lower physical function scores were observed in women with higher body mass index values (p < 0.05) and lower physical activity levels (p < 0.05). Generalized linear models with gamma distributions and logarithmic link functions, adjusted for age, showed that lower total scores and physical function domain scores for the Quality of Life Questionnaire of the European Foundation for Osteoporosis were related to a high body mass index, lower physical activity, and the presence of vertebral fractures (p < 0.05). CONCLUSION: Vertebral fractures are associated with decreased quality of life mainly physical functioning in older community-dwelling women regardless of age, body mass index, and physical activity. Therefore, the results highlight the importance of preventing and controlling asymptomatic vertebral fractures to reduce their impact on quality of life among older women.

Structure-function Relationships Using the Cirrus Spectral Domain Optical Coherence Tomograph and Standard Automated Perimetry

Purpose: To evaluate the relationship between glaucomatous structural damage assessed by the Cirrus Spectral Domain OCT (SDOCT) and functional loss as measured by standard automated perimetry (SAP). Methods: Four hundred twenty-two eyes (78 healthy, 210 suspects, 134 glaucomatous) of 250 patients were recruited from the longitudinal Diagnostic Innovations in Glaucoma Study and from the African Descent and Glaucoma Evaluation Study. All eyes underwent testing with the Cirrus SDOCT and SAP within a 6-month period. The relationship between parapapillary retinal nerve fiber layer thickness (RNFL) sectors and corresponding topographic SAP locations was evaluated using locally weighted scatterplot smoothing and regression analysis. SAP sensitivity values were evaluated using both linear as well as logarithmic scales. We also tested the fit of a model (Hood) for structure-function relationship in glaucoma. Results: Structure was significantly related to function for all but the nasal thickness sector. The relationship was strongest for superotemporal RNFL thickness and inferonasal sensitivity (R(2) = 0.314, P < 0.001). The Hood model fitted the data relatively well with 88% of the eyes inside the 95% confidence interval predicted by the model. Conclusions: RNFL thinning measured by the Cirrus SDOCT was associated with correspondent visual field loss in glaucoma.

The local power of the gradient test

The asymptotic expansion of the distribution of the gradient test statistic is derived for a composite hypothesis under a sequence of Pitman alternative hypotheses converging to the null hypothesis at rate n(-1/2), n being the sample size. Comparisons of the local powers of the gradient, likelihood ratio, Wald and score tests reveal no uniform superiority property. The power performance of all four criteria in one-parameter exponential family is examined.

Clinical correlates of social adjustment in patients with obsessive-compulsive disorder

Background: Patients with obsessive-compulsive disorder (OCD) frequently show poor social adjustment, which has been associated with OCD severity. Little is known about the effects that age at symptom onset, specific OCD symptoms, and psychiatric comorbidities have on social adjustment. The objective of this study was to investigate the clinical correlates of social functioning in OCD patients. Methods: Cross-sectional study involving 815 adults with a primary DSM-IV diagnosis of OCD participating in the Brazilian Research Consortium on Obsessive-Compulsive Spectrum Disorders. Patients were assessed with the Social Adjustment Scale, the Medical Outcomes Study 36-item Short-Form Health Survey, the Yale-Brown Obsessive-Compulsive Scale, the Dimensional Yale-Brown Obsessive-Compulsive Scale, and the Structured Clinical Interview for DSM-IV Axis I Disorders. Clinical correlates of social adjustment were assessed with generalized linear models with gamma distribution. Results: Poor overall social functioning was associated with greater OCD severity (p = 0.02); hoarding symptoms (p = 0.004); sexual/religious obsessions (p = 0.005); current major depressive disorder (p = 0.004); current post-traumatic stress disorder (p = 0.002); and current eating disorders (p = 0.02). Poor social adjustment was also associated with impaired quality of life. Conclusions: Patients with OCD have poor social functioning in domains related to personal relationships and professional performance. Hoarding symptoms and sexual/religious obsessions seem to have the strongest negative effects on social functioning. Early age at OCD symptom onset seems to be associated with professional and academic underachievement and impairment within the family unit, whereas current psychiatric comorbidity worsen overall social functioning. In comparison with quality of life, social adjustment measures seem to provide a more comprehensive overview of the OCD-related burden. (C) 2012 Elsevier Ltd. All rights reserved.