Latent residual analysis in binary regression with skewed link

Model diagnostics is an integral part of model determination and an important part of the model diagnostics is residual analysis. We adapt and implement residuals considered in the literature for the probit, logistic and skew-probit links under binary regression. New latent residuals for the skew-probit link are proposed here. We have detected the presence of outliers using the residuals proposed here for different models in a simulated dataset and a real medical dataset.

Trends of the genetic connectedness measures among Nelore beef cattle herds

Validity of comparisons between expected breeding values obtained from best linear unbiased prediction procedures in genetic evaluations is dependent on genetic connectedness among herds. Different cattle breeding programmes have their own particular features that distinguish their database structure and can affect connectedness. Thus, the evolution of these programmes can also alter the connectedness measures. This study analysed the evolution of the genetic connectedness measures among Brazilian Nelore cattle herds from 1999 to 2008, using the French Criterion of Admission to the group of Connected Herds (CACO) method, based on coefficients of determination (CD) of contrasts. Genetic connectedness levels were analysed by using simple and multiple regression analyses on herd descriptors to understand their relationship and their temporal trends from the 19992003 to the 20042008 period. The results showed a high level of genetic connectedness, with CACO estimates higher than 0.4 for the majority of them. Evaluation of the last 5-year period showed only a small increase in average CACO measures compared with the first 5 years, from 0.77 to 0.80. The percentage of herds with CACO estimates lower than 0.7 decreased from 27.5% in the first period to 16.2% in the last one. The connectedness measures were correlated with percentage of progeny from connecting sires, and the artificial insemination spread among Brazilian herds in recent years. But changes in connectedness levels were shown to be more complex, and their complete explanation cannot consider only herd descriptors. They involve more comprehensive changes in the relationship matrix, which can be only fully expressed by the CD of contrasts.

Effects of tannery sludge application on physiological and fatty acid profiles of the soil microbial community

The impact of tannery sludge application on soil microbial community and diversity is poorly understood. We studied the microbial community in an agricultural soil following two applications (2006 and 2007) of tannery sludge with annual application rates of 0.0,2.3 and 22.6 Mg ha(-1). The soil was sampled 12 and 271 days after the second (2007) application. Community structure was assessed via a phospholipid fatty acid analysis, and the physiological profile of the soil microbial community via the Biolog method. Tannery sludge application changed soil chemical properties, increasing the soil pH and electrical conductivity as well as available P and mineral N concentrations. The higher sludge application rate changed the community structure and the physiological profile of the microbial community at both sampling dates. However, there is no clear link between community structure and carbon substrate utilization. According to the Distance Based Linear Models Analysis, the fatty acids 16:0 and 117:0 together contributed 84% to the observed PLFA patterns, whereas the chemical properties available P, mineral N, and Ca, and pH together contributed 54%. At 12 days, tannery sludge application increased the average well color development from 0.46 to 0.87 after 48 h, and reduced the time elapsed before reaching the midpoint carbon substrate utilization (s) from 71 to 44 h, an effect still apparent nine months after application of the higher sludge application rate. The dominant signature fatty acids and kinetic parameters (r and s) were correlated to the concentrations of available P. Ca, mineral N, pH and EC. (c) 2012 Elsevier B.V. All rights reserved.

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.

Effectiveness of a long-term use of a minimalist footwear versus habitual shoe on pain, function and mechanical loads in knee osteoarthritis: a randomized controlled trial

Background: Recent studies have shown an important reduction of joint overload during locomotion in elderly women with knee osteoarthritis (OA) after short- term use of minimalist shoes. Our aim is to investigate the chronic effect of inexpensive and minimalist footwear on the clinical and functional aspects of OA and gait biomechanics of elderly women with knee OA. Methods/Design: Fifty-six elderly women with knee OA grade 2 or 3 (Kellgren and Lawrence) are randomized into blocks and allocated to either the intervention group, which will use flexible, non-heeled shoes-Moleca (R)-for six months for at least six hours daily, or the control group, which could not use these shoes. Neither group is undergoing physical therapy treatment throughout the intervention period. Moleca (R) is a women's double canvas, flexible, flat walking shoe without heels, with a 5-mm anti-slip rubber sole and a 3-mm internal wedge of ethylene vinyl acetate. Both groups will be followed for six months and will be assessed at baseline condition, after three months, and after six months (end of intervention). All the assessments will be performed by a physiotherapist that is blind to the group allocation. The primary outcome is the pain Western Ontario and McMaster Universities Osteoarthritis (WOMAC) score. The secondary outcomes are global WOMAC score; joint stiffness and disability WOMAC scores; knee pain with a visual analogue scale; walking distance in the six-minute walk test; Lequesne score; amount and frequency (number of days) of paracetamol (500 mg) intake over six months; knee adduction moment during gait; global medical assessment score; and global patient auto-assessment score. At baseline, all patients receive a diary to record the hours of daily use of the footwear intervention; every two weeks, the same physiotherapist makes phone calls to all patients in order to verify adherence to treatment. The statistical analysis will be based on intention to treat analysis, as well as general linear models of analysis of variance for repeated measure to detect treatment-time interactions (alpha = 5%). Discussion: This is the first randomized, clinical trial protocol to assess the chronic effect of minimalist footwear on the clinical and functional aspects and gait biomechanics of elderly women with knee osteoarthritis. We expect that the use of Moleca (R) shoes for six months will provide pain relief, reduction of the knee adduction moment when walking, and improve joint function in elderly women with knee OA, and that the treatment, thus, can be considered another inexpensive and easy-to-use option for conservative OA treatment.

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.

Longitudinal familial analysis of blood pressure involving parametric (co)variance functions

Abstract Background For analyzing longitudinal familial data we adopted a log-linear form to incorporate heterogeneity in genetic variance components over the time, and additionally a serial correlation term in the genetic effects at different levels of ages. Due to the availability of multiple measures on the same individual, we permitted environmental correlations that may change across time. Results Systolic blood pressure from family members from the first and second cohort was used in the current analysis. Measures of subjects receiving hypertension treatment were set as censored values and they were corrected. An initial check of the variance and covariance functions proposed for analyzing longitudinal familial data, using empirical semi-variogram plots, indicated that the observed trait dispersion pattern follows the assumptions adopted. Conclusion The corrections for censored phenotypes based on ordinary linear models may be an appropriate simple model to correct the data, ensuring that the original variability in the data was retained. In addition, empirical semi-variogram plots are useful for diagnosis of the (co)variance model adopted.

Inequalities in mortality of men by oral and pharyngeal cancer in Barcelona, Spain and São Paulo, Brazil, 1995–2003

Abstract Background Large inequalities of mortality by most cancers in general, by mouth and pharynx cancer in particular, have been associated to behaviour and geopolitical factors. The assessment of socioeconomic covariates of cancer mortality may be relevant to a full comprehension of distal determinants of the disease, and to appraise opportune interventions. The objective of this study was to compare socioeconomic inequalities in male mortality by oral and pharyngeal cancer in two major cities of Europe and South America. Methods The official system of information on mortality provided data on deaths in each city; general censuses informed population data. Age-adjusted death rates by oral and pharyngeal cancer for men were independently assessed for neighbourhoods of Barcelona, Spain, and São Paulo, Brazil, from 1995 to 2003. Uniform methodological criteria instructed the comparative assessment of magnitude, trends and spatial distribution of mortality. General linear models assessed ecologic correlations between death rates and socioeconomic indices (unemployment, schooling levels and the human development index) at the inner-city area level. Results obtained for each city were subsequently compared. Results Mortality of men by oral and pharyngeal cancer ranked higher in Barcelona (9.45 yearly deaths per 100,000 male inhabitants) than in Spain and Europe as a whole; rates were on decrease. São Paulo presented a poorer profile, with higher magnitude (11.86) and stationary trend. The appraisal of ecologic correlations indicated an unequal and inequitably distributed burden of disease in both cities, with poorer areas tending to present higher mortality. Barcelona had a larger gradient of mortality than São Paulo, indicating a higher inequality of cancer deaths across its neighbourhoods. Conclusion The quantitative monitoring of inequalities in health may contribute to the formulation of redistributive policies aimed at the concurrent promotion of wellbeing and social justice. The assessment of groups experiencing a higher burden of disease can instruct health services to provide additional resources for expanding preventive actions and facilities aimed at early diagnosis, standardized treatments and rehabilitation.

The impact of asymptomatic vertebral fractures on quality of life in older community-dwelling women: the São 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.

Local power properties of some asymptotic tests in symmetric linear regression models

In this paper we obtain asymptotic expansions, up to order n(-1/2) and under a sequence of Pitman alternatives, for the nonnull distribution functions of the likelihood ratio, Wald, score and gradient test statistics in the class of symmetric linear regression models. This is a wide class of models which encompasses the t model and several other symmetric distributions with longer-than normal tails. The asymptotic distributions of all four statistics are obtained for testing a subset of regression parameters. Furthermore, in order to compare the finite-sample performance of these tests in this class of models, Monte Carlo simulations are presented. An empirical application to a real data set is considered for illustrative purposes. (C) 2011 Elsevier B.V. All rights reserved.

Design-based random permutation models with auxiliary information

We extend the random permutation model to obtain the best linear unbiased estimator of a finite population mean accounting for auxiliary variables under simple random sampling without replacement (SRS) or stratified SRS. The proposed method provides a systematic design-based justification for well-known results involving common estimators derived under minimal assumptions that do not require specification of a functional relationship between the response and the auxiliary variables.

Assessment of variance components in nonlinear mixed-effects elliptical models

The issue of assessing variance components is essential in deciding on the inclusion of random effects in the context of mixed models. In this work we discuss this problem by supposing nonlinear elliptical models for correlated data by using the score-type test proposed in Silvapulle and Silvapulle (1995). Being asymptotically equivalent to the likelihood ratio test and only requiring the estimation under the null hypothesis, this test provides a fairly easy computable alternative for assessing one-sided hypotheses in the context of the marginal model. Taking into account the possible non-normal distribution, we assume that the joint distribution of the response variable and the random effects lies in the elliptical class, which includes light-tailed and heavy-tailed distributions such as Student-t, power exponential, logistic, generalized Student-t, generalized logistic, contaminated normal, and the normal itself, among others. We compare the sensitivity of the score-type test under normal, Student-t and power exponential models for the kinetics data set discussed in Vonesh and Carter (1992) and fitted using the model presented in Russo et al. (2009). Also, a simulation study is performed to analyze the consequences of the kurtosis misspecification.

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.

The concept of quasi-integrability for modified non-linear Schrodinger models

We consider modifications of the nonlinear Schrodinger model (NLS) to look at the recently introduced concept of quasi-integrability. We show that such models possess an in finite number of quasi-conserved charges which present intriguing properties in relation to very specific space-time parity transformations. For the case of two-soliton solutions where the fields are eigenstates of this parity, those charges are asymptotically conserved in the scattering process of the solitons. Even though the charges vary in time their values in the far past and the far future are the same. Such results are obtained through analytical and numerical methods, and employ adaptations of algebraic techniques used in integrable field theories. Our findings may have important consequences on the applications of these models in several areas of non-linear science. We make a detailed numerical study of the modified NLS potential of the form V similar to (vertical bar psi vertical bar(2))(2+epsilon), with epsilon being a perturbation parameter. We perform numerical simulations of the scattering of solitons for this model and find a good agreement with the results predicted by the analytical considerations. Our paper shows that the quasi-integrability concepts recently proposed in the context of modifications of the sine-Gordon model remain valid for perturbations of the NLS model.

Degree of multicollinearity and variables involved in linear dependence in additive-dominant models

The objective of this work was to assess the degree of multicollinearity and to identify the variables involved in linear dependence relations in additive-dominant models. Data of birth weight (n=141,567), yearling weight (n=58,124), and scrotal circumference (n=20,371) of Montana Tropical composite cattle were used. Diagnosis of multicollinearity was based on the variance inflation factor (VIF) and on the evaluation of the condition indexes and eigenvalues from the correlation matrix among explanatory variables. The first model studied (RM) included the fixed effect of dam age class at calving and the covariates associated to the direct and maternal additive and non-additive effects. The second model (R) included all the effects of the RM model except the maternal additive effects. Multicollinearity was detected in both models for all traits considered, with VIF values of 1.03 - 70.20 for RM and 1.03 - 60.70 for R. Collinearity increased with the increase of variables in the model and the decrease in the number of observations, and it was classified as weak, with condition index values between 10.00 and 26.77. In general, the variables associated with additive and non-additive effects were involved in multicollinearity, partially due to the natural connection between these covariables as fractions of the biological types in breed composition.