News from compositions, the R package

The R-package “compositions”is a tool for advanced compositional analysis. Its basicfunctionality has seen some conceptual improvement, containing now some facilitiesto work with and represent ilr bases built from balances, and an elaborated subsys-tem for dealing with several kinds of irregular data: (rounded or structural) zeroes,incomplete observations and outliers. The general approach to these irregularities isbased on subcompositions: for an irregular datum, one can distinguish a “regular” sub-composition (where all parts are actually observed and the datum behaves typically)and a “problematic” subcomposition (with those unobserved, zero or rounded parts, orelse where the datum shows an erratic or atypical behaviour). Systematic classificationschemes are proposed for both outliers and missing values (including zeros) focusing onthe nature of irregularities in the datum subcomposition(s).To compute statistics with values missing at random and structural zeros, a projectionapproach is implemented: a given datum contributes to the estimation of the desiredparameters only on the subcompositon where it was observed. For data sets withvalues below the detection limit, two different approaches are provided: the well-knownimputation technique, and also the projection approach.To compute statistics in the presence of outliers, robust statistics are adapted to thecharacteristics of compositional data, based on the minimum covariance determinantapproach. The outlier classification is based on four different models of outlier occur-rence and Monte-Carlo-based tests for their characterization. Furthermore the packageprovides special plots helping to understand the nature of outliers in the dataset.Keywords: coda-dendrogram, lost values, MAR, missing data, MCD estimator,robustness, rounded zeros

Total hemoglobin mass--a new parameter to detect blood doping?

PURPOSE: All kinds of blood manipulations aim to increase the total hemoglobin mass (tHb-mass). To establish tHb-mass as an effective screening parameter for detecting blood doping, the knowledge of its normal variation over time is necessary. The aim of the present study, therefore, was to determine the intraindividual variance of tHb-mass in elite athletes during a training year emphasizing off, training, and race seasons at sea level. METHODS: tHb-mass and hemoglobin concentration ([Hb]) were determined in 24 endurance athletes five times during a year and were compared with a control group (n = 6). An analysis of covariance was used to test the effects of training phases, age, gender, competition level, body mass, and training volume. Three error models, based on 1) a total percentage error of measurement, 2) the combination of a typical percentage error (TE) of analytical origin with an absolute SD of biological origin, and 3) between-subject and within-subject variance components as obtained by an analysis of variance, were tested. RESULTS: In addition to the expected influence of performance status, the main results were that the effects of training volume (P = 0.20) and training phases (P = 0.81) on tHb-mass were not significant. We found that within-subject variations mainly have an analytical origin (TE approximately 1.4%) and a very small SD (7.5 g) of biological origin. CONCLUSION: tHb-mass shows very low individual oscillations during a training year (<6%), and these oscillations are below the expected changes in tHb-mass due to Herythropoetin (EPO) application or blood infusion (approximately 10%). The high stability of tHb-mass over a period of 1 year suggests that it should be included in an athlete's biological passport and analyzed by recently developed probabilistic inference techniques that define subject-based reference ranges.

Response to methylphenidate by adult and pediatric patients with attention-deficit/hyperactivity disorder: the Spanish multicenter DIHANA study.

BACKGROUND The purpose of this multicenter Spanish study was to evaluate the response to immediate-release methylphenidate by children and adults diagnosed with attention-deficit/hyperactivity disorder (ADHD), as well as to obtain information on current therapy patterns and safety characteristics. METHODS This multicenter, observational, retrospective, noninterventional study included 730 patients aged 4-65 years with a diagnosis of ADHD. Information was obtained based on a review of medical records for the years 2002-2006 in sequential order. RESULTS The ADHD predominantly inattentive subtype affected 29.7% of patients, ADHD predominantly hyperactive-impulsive was found in 5.2%, and the combined subtype in 65.1%. Overall, a significant lower Clinical Global Impression (CGI) score and mean number of DSM-IV TR (Diagnostic and Statistical Manual of Mental Disorders Fourth Edition, Text Revision) symptoms by subtype were found after one year of treatment with immediate-release methylphenidate; CGI decreased from 4.51 to 1.69, symptoms of inattention from 7.90 to 4.34, symptoms of hyperactivity from 6.73 to 3.39, and combined subtype symptoms from 14.62 to 7.7. Satisfaction with immediate-release methylphenidate after one year was evaluated as "very satisfied" or "satisfied" by 86.90% of the sample; 25.75% of all patients reported at least one adverse effect. At the end of the study, 41.47% of all the patients treated with immediate-release methylphenidate were still receiving it, with a mean time of 3.80 years on therapy. CONCLUSION Good efficacy and safety results were found for immediate-release methylphenidate in patients with ADHD.

Consumption of vegetables, fruit and other plant foods in the European Prospective Investigation into Cancer and Nutrition (EPIC) cohorts from 10 European countries.

OBJECTIVE To describe and compare the consumption of the main groups and sub-groups of vegetables and fruits (V&F) in men and women from the centres participating in the European Prospective Investigation into Cancer and Nutrition (EPIC). DESIGN Cross-sectional analysis. Dietary intake was assessed by means of a 24-hour dietary recall using computerised interview software and standardised procedures. Crude and adjusted means were computed for the main groups and sub-groups of V&F by centre, separately for men and women. Adjusted means by season, day of the week and age were estimated using weights and covariance analysis. SETTING Twenty-seven centres in 10 European countries participating in the EPIC project. SUBJECTS In total, 35 955 subjects (13 031 men and 22 924 women), aged 35-74 years, randomly selected from each EPIC cohort. RESULTS The centres from southern countries had the highest consumption of V&F, while the lowest intake was seen in The Netherlands and Scandinavia for both genders. These differences were more evident for fruits, particularly citrus. However, slightly different patterns arose for some sub-groups of vegetables, such as root vegetables and cabbage. Adjustment for body mass index, physical activity, smoking habits and education did not substantially modify the mean intakes of vegetables and fruits. CONCLUSIONS Total vegetable and fruit intake follows a south-north gradient in both genders, whereas for several sub-groups of vegetables a different geographic distribution exists. Differences in mean intake of V&F by centre were not explained by lifestyle factors associated with V&F intake.

Markov chain montecarlo method applied to rounding zeros of compositional data: first approach

As stated in Aitchison (1986), a proper study of relative variation in a compositional data set should be based on logratios, and dealing with logratios excludes dealing with zeros. Nevertheless, it is clear that zero observations might be present in real data sets, either because the corresponding part is completelyabsent –essential zeros– or because it is below detection limit –rounded zeros. Because the second kind of zeros is usually understood as “a trace too small to measure”, it seems reasonable to replace them by a suitable small value, and this has been the traditional approach. As stated, e.g. by Tauber (1999) and byMartín-Fernández, Barceló-Vidal, and Pawlowsky-Glahn (2000), the principal problem in compositional data analysis is related to rounded zeros. One should be careful to use a replacement strategy that does not seriously distort the general structure of the data. In particular, the covariance structure of the involvedparts –and thus the metric properties– should be preserved, as otherwise further analysis on subpopulations could be misleading. Following this point of view, a non-parametric imputation method isintroduced in Martín-Fernández, Barceló-Vidal, and Pawlowsky-Glahn (2000). This method is analyzed in depth by Martín-Fernández, Barceló-Vidal, and Pawlowsky-Glahn (2003) where it is shown that thetheoretical drawbacks of the additive zero replacement method proposed in Aitchison (1986) can be overcome using a new multiplicative approach on the non-zero parts of a composition. The new approachhas reasonable properties from a compositional point of view. In particular, it is “natural” in the sense thatit recovers the “true” composition if replacement values are identical to the missing values, and it is coherent with the basic operations on the simplex. This coherence implies that the covariance structure of subcompositions with no zeros is preserved. As a generalization of the multiplicative replacement, in thesame paper a substitution method for missing values on compositional data sets is introduced

Existential behavioural therapy for informal caregivers of palliative patients: a randomised controlled trial.

BACKGROUND: Existential behavioural therapy (EBT) was developed to support informal caregivers of palliative patients in the last stage of life and during bereavement as a manualised group psychotherapy comprising six sessions. We tested the effectiveness of EBT on mental stress and quality of life (QOL). METHODS: Informal caregivers were randomly assigned (1:1) to EBT or a treatment-as-usual control group using computer-generated numbers in blocks of 10. Primary outcomes were assessed with the Brief Symptom Inventory (subscales somatisation, anxiety and depression), the Satisfaction with Life Scale (SWLS), the WHOQOL-BREF and a numeric rating scale for QOL (QOL-NRS, range 0-10). Data were collected at baseline, pre-treatment, post-treatment and follow-ups after 3 and 12âeuro0/00months. Treatment effects were assessed with a multivariate analysis of covariance. RESULTS: Out of 160 relatives, 81 were assigned to EBT and 79 to the control group. Participants were 54.5âeuro0/00±âeuro0/0013.2âeuro0/00years old; 69.9% were female. The multivariate model was significant for the pre-/post-comparison (pâeuro0/00=âeuro0/000.005) and the pre-/12-month comparison (pâeuro0/00=âeuro0/000.05) but not for the pre-/3-month comparison. Medium to large effects on anxiety and QOL (SWLS, WHOQOL-BREF, QOL-NRS) were found at post-treatment; medium effects on depression and QOL (QOL-NRS) emerged in the 12-month follow-up. No adverse effects of the intervention were observed. CONCLUSION: Existential behavioural therapy appears to exert beneficial effects on distress and QOL of informal caregivers of palliative patients. Further longitudinal evidence is needed to confirm these findings. Copyright © 2013 John Wiley & Sons, Ltd.

Brain tissue properties differentiate between motor and limbic basal ganglia circuits.

Despite advances in understanding basic organizational principles of the human basal ganglia, accurate in vivo assessment of their anatomical properties is essential to improve early diagnosis in disorders with corticosubcortical pathology and optimize target planning in deep brain stimulation. Main goal of this study was the detailed topological characterization of limbic, associative, and motor subdivisions of the subthalamic nucleus (STN) in relation to corresponding corticosubcortical circuits. To this aim, we used magnetic resonance imaging and investigated independently anatomical connectivity via white matter tracts next to brain tissue properties. On the basis of probabilistic diffusion tractography we identified STN subregions with predominantly motor, associative, and limbic connectivity. We then computed for each of the nonoverlapping STN subregions the covariance between local brain tissue properties and the rest of the brain using high-resolution maps of magnetization transfer (MT) saturation and longitudinal (R1) and transverse relaxation rate (R2*). The demonstrated spatial distribution pattern of covariance between brain tissue properties linked to myelin (R1 and MT) and iron (R2*) content clearly segregates between motor and limbic basal ganglia circuits. We interpret the demonstrated covariance pattern as evidence for shared tissue properties within a functional circuit, which is closely linked to its function. Our findings open new possibilities for investigation of changes in the established covariance pattern aiming at accurate diagnosis of basal ganglia disorders and prediction of treatment outcome.

PACIC Instrument: disentangling dimensions using published validation models.

OBJECTIVE: To better understand the structure of the Patient Assessment of Chronic Illness Care (PACIC) instrument. More specifically to test all published validation models, using one single data set and appropriate statistical tools. DESIGN: Validation study using data from cross-sectional survey. PARTICIPANTS: A population-based sample of non-institutionalized adults with diabetes residing in Switzerland (canton of Vaud). MAIN OUTCOME MEASURE: French version of the 20-items PACIC instrument (5-point response scale). We conducted validation analyses using confirmatory factor analysis (CFA). The original five-dimension model and other published models were tested with three types of CFA: based on (i) a Pearson estimator of variance-covariance matrix, (ii) a polychoric correlation matrix and (iii) a likelihood estimation with a multinomial distribution for the manifest variables. All models were assessed using loadings and goodness-of-fit measures. RESULTS: The analytical sample included 406 patients. Mean age was 64.4 years and 59% were men. Median of item responses varied between 1 and 4 (range 1-5), and range of missing values was between 5.7 and 12.3%. Strong floor and ceiling effects were present. Even though loadings of the tested models were relatively high, the only model showing acceptable fit was the 11-item single-dimension model. PACIC was associated with the expected variables of the field. CONCLUSIONS: Our results showed that the model considering 11 items in a single dimension exhibited the best fit for our data. A single score, in complement to the consideration of single-item results, might be used instead of the five dimensions usually described.

Meta-analysis of gene-level associations for rare variants based on single-variant statistics.

Meta-analysis of genome-wide association studies (GWASs) has led to the discoveries of many common variants associated with complex human diseases. There is a growing recognition that identifying "causal" rare variants also requires large-scale meta-analysis. The fact that association tests with rare variants are performed at the gene level rather than at the variant level poses unprecedented challenges in the meta-analysis. First, different studies may adopt different gene-level tests, so the results are not compatible. Second, gene-level tests require multivariate statistics (i.e., components of the test statistic and their covariance matrix), which are difficult to obtain. To overcome these challenges, we propose to perform gene-level tests for rare variants by combining the results of single-variant analysis (i.e., p values of association tests and effect estimates) from participating studies. This simple strategy is possible because of an insight that multivariate statistics can be recovered from single-variant statistics, together with the correlation matrix of the single-variant test statistics, which can be estimated from one of the participating studies or from a publicly available database. We show both theoretically and numerically that the proposed meta-analysis approach provides accurate control of the type I error and is as powerful as joint analysis of individual participant data. This approach accommodates any disease phenotype and any study design and produces all commonly used gene-level tests. An application to the GWAS summary results of the Genetic Investigation of ANthropometric Traits (GIANT) consortium reveals rare and low-frequency variants associated with human height. The relevant software is freely available.

Whole-body protein turnover and resting energy expenditure in obese, prepubertal children.

BACKGROUND: Obesity is becoming more frequent in children; understanding the extent to which this condition affects not only carbohydrate and lipid metabolism but also protein metabolism is of paramount importance. OBJECTIVE: We evaluated the kinetics of protein metabolism in obese, prepubertal children in the static phase of obesity. DESIGN: In this cross-sectional study, 9 obese children (x +/- SE: 44+/-4 kg, 30.9+/-1.5% body fat) were compared with 8 lean (28+/-2 kg ,16.8+/-1.2% body fat), age-matched (8.5+/-0.2 y) control children. Whole-body nitrogen flux, protein synthesis, and protein breakdown were calculated postprandially over 9 h from 15N abundance in urinary ammonia by using a single oral dose of [15N]glycine; resting energy expenditure (REE) was assessed by indirect calorimetry (canopy) and body composition by multiple skinfold-thickness measurements. RESULTS: Absolute rates of protein synthesis and breakdown were significantly greater in obese children than in control children (x +/- SE: 208+/-24 compared with 137+/-14 g/d, P < 0.05, and 149+/-20 compared with 89+/-13 g/d, P < 0.05, respectively). When these variables were adjusted for fat-free mass by analysis of covariance, however, the differences between groups disappeared. There was a significant relation between protein synthesis and fat-free mass (r = 0.83, P < 0.001) as well as between protein synthesis and REE (r = 0.79, P < 0.005). CONCLUSIONS: Obesity in prepubertal children is associated with an absolute increase in whole-body protein turnover that is consistent with an absolute increase in fat-free mass, both of which contribute to explaining the greater absolute REE in obese children than in control children.

Modelling the effect of minor orthopaedic day surgery on patient mood at the early post-operative period: a prospective population-based cohort study.

OBJECTIVE: The effect of minor orthopaedic day surgery (MiODS) on patient's mood. METHODS: A prospective population-based cohort study of 148 consecutive patients with age above 18 and less than 65, an American Society of Anaesthesiology (ASA) score of 1, and the requirement of general anaesthesia (GA) were included. The Medical Outcomes Study - Short Form 36 (SF-36), Beck Anxiety Inventory (BAI) and Beck Depression Inventory (BDI) were used pre- and post-operatively. RESULTS: The mean physical component score of SF-36 before surgery was 45.3 (SD=+/-10.1) and 8 weeks following surgery was 44.9 (SD=+/-11.04) [n=148, p=0.51, 95% CI=(-1.03 to 1.52)]. For the measurement of the changes in mood using BDI, BAI and SF-36, latent construct modelling was employed to increase validity. The covariance between mood pre- and post-operatively (cov=69.44) corresponded to a correlation coefficient, r=0.88 indicating that patients suffering a greater number of mood symptoms before surgery continue to have a greater number of symptoms following surgery. When the latent mood constructs were permitted to have different means the model fitted well with chi(2) (df=1)=0.86 for which p=0.77, thus the null hypothesis that MiODS has no effect on patient mood was rejected. CONCLUSIONS: MiODS affects patient mood which deteriorates at 8 weeks post-operatively regardless of the pre-operative patient mood state. More importantly patients suffering a greater number of mood symptoms before MiODS continue to have a greater number of symptoms following surgery.

How geographic distance and depth drive ecological variability and isolation of demersal fish communities in an archipelago system (Cape Verde, Eastern Atlantic Ocean)

Cape Verde is a tropical oceanic ecosystem, highly fragmented and dispersed, with islands physically isolated by distance and depth. To understand how isolation affects the ecological variability in this archipelago, we conducted a research project on the community structure of the 18 commercially most important demersal fishes. An index of ecological distance based on species relative dominance (Di) is developed from Catch Per Unit Effort, derived from an extensive database of artisanal fisheries. Two ecological measures of distance between islands are calculated: at the species level, DDi, and at the community level, DD (sum of DDi). A physical isolation factor (Idb) combining distance (d) and bathymetry (b) is proposed. Covariance analysis shows that isolation factor is positively correlated with both DDi and DD, suggesting that Idb can be considered as an ecological isolation factor. The effect of Idb varies with season and species. This effect is stronger in summer (May to November), than in winter (December to April), which appears to be more unstable. Species react differently to Idb, independently of season. A principal component analysis on the monthly (DDi) for the 12 islands and the 18 species, complemented by an agglomerative hierarchical clustering, shows a geographic pattern of island organization, according to Idb. Results indicate that the ecological structure of demersal fish communities of Cape Verde archipelago, both in time and space, can be explained by a geographic isolation factor. The analytical approach used here is promising and could be tested in other archipelago systems.

Validity of the Katz Index to assess activities of daily living by informants in neuropathological studies

Abstract OBJECTIVE To analyze the evidences of construct validity of the Katz Index for the retrospective assessment of activities of daily living (ADL) by informants, to assist neuropathological studies in the elderly. METHOD A cross-sectional study analyzed the functional ability of ADL measure by the Katz Index, of 650 cases randomly selected from the Brazilian Brain Bank of the Ageing Brain Study Group (BBBABSG) database. Sample was divided in two subsamples for the analysis (N=325, each) and then stratified according to cognitive decline assessed by the Clinical Dementia Rating Scale (CDR). Factor analyses with calculations of internal consistency and invariance were performed. RESULTS Factor analysis evidenced a unidimensional instrument with optimal internal consistency, in all subgroups. Goodness of fit indices were obtained after two treatments of covariance, indicating adequacy of the scale for assessing ADL by informants. The scale is invariant to cognitive decline meaning that it can be used for subjects with or without cognitive impairment. CONCLUSION Katz Index is valid for the retrospective assessment of basic ADL by informants, with optimal reliability.

A new method for constructing exact tests without making any assumptions

We present a new method for constructing exact distribution-free tests (and confidence intervals) for variables that can generate more than two possible outcomes.This method separates the search for an exact test from the goal to create a non-randomized test. Randomization is used to extend any exact test relating to meansof variables with finitely many outcomes to variables with outcomes belonging to agiven bounded set. Tests in terms of variance and covariance are reduced to testsrelating to means. Randomness is then eliminated in a separate step.This method is used to create confidence intervals for the difference between twomeans (or variances) and tests of stochastic inequality and correlation.

Weak approximations. A Malliavin calculus approach

We introduce a variation of the proof for weak approximations that issuitable for studying the densities of stochastic processes which areevaluations of the flow generated by a stochastic differential equation on a random variable that maybe anticipating. Our main assumption is that the process and the initial random variable have to be smooth in the Malliavin sense. Furthermore if the inverse of the Malliavin covariance matrix associated with the process under consideration is sufficiently integrable then approximations fordensities and distributions can also be achieved. We apply theseideas to the case of stochastic differential equations with boundaryconditions and the composition of two diffusions.