Compositional analysis of bivariate discrete probabilities

A joint distribution of two discrete random variables with finite support can be displayed as a two way table of probabilities adding to one. Assume that this table hasn rows and m columns and all probabilities are non-null. This kind of table can beseen as an element in the simplex of n · m parts. In this context, the marginals areidentified as compositional amalgams, conditionals (rows or columns) as subcompositions. Also, simplicial perturbation appears as Bayes theorem. However, the Euclideanelements of the Aitchison geometry of the simplex can also be translated into the tableof probabilities: subspaces, orthogonal projections, distances.Two important questions are addressed: a) given a table of probabilities, which isthe nearest independent table to the initial one? b) which is the largest orthogonalprojection of a row onto a column? or, equivalently, which is the information in arow explained by a column, thus explaining the interaction? To answer these questionsthree orthogonal decompositions are presented: (1) by columns and a row-wise geometric marginal, (2) by rows and a columnwise geometric marginal, (3) by independenttwo-way tables and fully dependent tables representing row-column interaction. Animportant result is that the nearest independent table is the product of the two (rowand column)-wise geometric marginal tables. A corollary is that, in an independenttable, the geometric marginals conform with the traditional (arithmetic) marginals.These decompositions can be compared with standard log-linear models.Key words: balance, compositional data, simplex, Aitchison geometry, composition,orthonormal basis, arithmetic and geometric marginals, amalgam, dependence measure,contingency table

Intragenomic variation in the second internal transcribed spacer of the ribosomal DNA of species of the genera Culex and Lutzia (Diptera: Culicidae)

Culex is the largest genus of Culicini and includes vectors of several arboviruses and filarial worms. Many species of Culex are morphologically similar, which makes their identification difficult, particularly when using female specimens. To aid evolutionary studies and species distinction, molecular techniques are often used. Sequences of the second internal transcribed spacer (ITS2) of ribosomal DNA (rDNA) from 16 species of the genus Culex and one of Lutzia were used to assess their genomic variability and to verify their applicability in the phylogenetic analysis of the group. The distance matrix (uncorrected p-distance) that was obtained revealed intragenomic and intraspecific variation. Because of the intragenomic variability, we selected ITS2 copies for use in distance analyses based on their secondary structures. Neighbour-joining topology was obtained with an uncorrected p-distance. Despite the heterogeneity observed, individuals of the same species were grouped together and correlated with the current, morphology-based classification, thereby showing that ITS2 is an appropriate marker to be used in the taxonomy of Culex.

Genetic and morphological characterisation of a new species of the genus Hysterothylacium (Nematoda) from Paralichthys isosceles Jordan, 1890 (Pisces: Teleostei) of the Neotropical Region, state of Rio de Janeiro, Brazil

Taking into account the difficulties of taxonomic identification of larval anisakid nematodes based on morphological characters, genetic analyses were performed, together with those usually applied, in order to identify anisakid larvae found in the flounder Paralichthys isosceles from the littoral of the state of Rio de Janeiro, Brazil. The analysis of 1,820 larvae revealed a new species, similar to Hysterothylacium MD, Hysterothylacium 2, Hysterothylacium KB and Hysterothylacium sp regarding the absence of the larval tooth, an excretory pore situated below the nerve ring level, and slender lateral alae. Moreover, the new species differs from Hysterothylacium fortalezae and Hysterothylacium reliquens with regard to the number and size of spines present on the tail end and from Hysterothylacium patagonicus by the absence of interlabia. The maximum parsimony and neighbour joining tree topologies based on the 18S ribosomal DNA gene, complete internal transcribed spacer region and cytochrome oxidase 2 (COII) gene demonstrated that the Brazilian larvae belong to Raphidascarididae and represent a unique genetic entity, confirmed as a new Hysterothylacium species. Furthermore, the new species presents COII genetic signatures and shares polymorphisms with Raphidascarididae members. This is the first description of a new anisakid species from Brazil through the integration of morphological and molecular taxonomy data.

Aplicació de visió per dispositius mòbils

En aquest projecte es presenta l’aplicació per a dispositius mòbils Doppelganger. La seva funció és, a partir d’una fotografia, detectar la cara i mostrar la persona famosa de la nostra base de dades que més s’assembla a la persona en la fotografia. Per la implementació s’han utilitzat algoritmes de visió per computador i d’aprenentatge automàtic com per exemple el PCA i el K-Nearest Neighbor, tot utilitzant llibreries gratuïtes com són les OpenCV.

Influence of the ankle arthropathies on the 3D kinematics of the lower limbs during gait: analysis with an ambulatory system.

Introduction: Ankle arthropathy is associated with a decreased motion of the ankle-hindfoot during ambulation. Ankle arthrodesis was shown to result in degeneration of the neighbour joints of the foot. Inversely, total ankle arthroplasty conceptually preserves the adjacent joints because of the residual mobility of the ankle but this has not been demonstrated yet in vivo. It has also been reported that degenerative ankle diseases, and even arthrodesis, do not result in alteration of the knee and hip joints. We present the preliminary results of a new approach of this problem based on ambulatory gait analysis. Patients and Methods: Motion analysis of the lower limbs was performed using a Physilog® (BioAGM, CH) system consisting of three-dimensional (3D) accelerometer and gyroscope, coupled to a magnetic system (Liberty©, Polhemus, USA). Both systems have been validated. Three groups of two patients were included into this pilot study and compared to healthy subjects (controls) during level walking: patients with ankle osteoarthritis (group 1), patients treated by ankle arthrodesis (group 2), patients treated by total ankle prosthesis (group 3). Results: Motion patterns of all analyzed joints over more than 20 gait cycles in each subject were highly repeatable. Motion amplitude of the ankle-hindfoot in control patients was similar to recently reported results. Ankle arthrodesis limited the motion of the ankle-hindfoot in the sagittal and horizontal planes. The prosthetic ankle allowed a more physiologic movement in the sagittal plane only. Ankle arthritis and its treatments did not influence the range of motion of the knee and hip joint during stance phase, excepted for a slight decrease of the hip flexion in groups 1 and 2. Conclusion: The reliability of the system was shown by the repeatability of the consecutive measurements. The results of this preliminary study were similar to those obtained through laboratory gait analysis. However, our system has the advantage to allow ambulatory analysis of 3D kinematics of the lower limbs outside of a gait laboratory and in real life conditions. To our knowledge this is a new concept in the analysis of ankle arthropathy and its treatments. Therefore, there is a potential to address specific questions like the difficult comparison of the benefits of ankle arthroplasty versus arthrodesis. The encouraging results of this pilot study offer the perspective to analyze the consequences of ankle arthropathy and its treatments on the biomechanics of the lower limbs ambulatory, in vivo and in daily life conditions.

Effective diagnosis of Alzheimer's disease by means of large margin-based methodology.

BACKGROUND Functional brain images such as Single-Photon Emission Computed Tomography (SPECT) and Positron Emission Tomography (PET) have been widely used to guide the clinicians in the Alzheimer's Disease (AD) diagnosis. However, the subjectivity involved in their evaluation has favoured the development of Computer Aided Diagnosis (CAD) Systems. METHODS It is proposed a novel combination of feature extraction techniques to improve the diagnosis of AD. Firstly, Regions of Interest (ROIs) are selected by means of a t-test carried out on 3D Normalised Mean Square Error (NMSE) features restricted to be located within a predefined brain activation mask. In order to address the small sample-size problem, the dimension of the feature space was further reduced by: Large Margin Nearest Neighbours using a rectangular matrix (LMNN-RECT), Principal Component Analysis (PCA) or Partial Least Squares (PLS) (the two latter also analysed with a LMNN transformation). Regarding the classifiers, kernel Support Vector Machines (SVMs) and LMNN using Euclidean, Mahalanobis and Energy-based metrics were compared. RESULTS Several experiments were conducted in order to evaluate the proposed LMNN-based feature extraction algorithms and its benefits as: i) linear transformation of the PLS or PCA reduced data, ii) feature reduction technique, and iii) classifier (with Euclidean, Mahalanobis or Energy-based methodology). The system was evaluated by means of k-fold cross-validation yielding accuracy, sensitivity and specificity values of 92.78%, 91.07% and 95.12% (for SPECT) and 90.67%, 88% and 93.33% (for PET), respectively, when a NMSE-PLS-LMNN feature extraction method was used in combination with a SVM classifier, thus outperforming recently reported baseline methods. CONCLUSIONS All the proposed methods turned out to be a valid solution for the presented problem. One of the advances is the robustness of the LMNN algorithm that not only provides higher separation rate between the classes but it also makes (in combination with NMSE and PLS) this rate variation more stable. In addition, their generalization ability is another advance since several experiments were performed on two image modalities (SPECT and PET).

COI barcode versus morphological identification of Culex ( Culex ) (Diptera: Culicidae) species: a case study using samples from Argentina and Brazil

Sequences of the cytochrome c oxidase subunit I (COI) mitochondrial gene from adults of 22 Culex ( Culex ) species from Argentina and Brazil were employed to assess species identification and to test the usefulness of COI for barcoding using the best close match (BCM) algorithm. A pairwise Kimura two-parameter distance matrix including the mean intra and interspecific distances for 71 COI barcode sequences was constructed. Of the 12 COI lineages recovered in the Neighbour-joining topology, five confirmed recognised morphological species ( Cx. acharistus , Cx. chidesteri , Cx. dolosus , Cx. lygrus and Cx. saltanensis ) with intraspecific divergences lower than 1.75%. Cx. bilineatus is formally resurrected from the synonymy of Cx. dolosus . Cx. maxi , Cx. surinamensis and the Coronator group species included were clustered into an unresolved lineage. The intraspecific distance of Cx. pipiens (3%) was almost twice the interspecific between it and Cx. quinquefasciatus (1.6%). Regarding the BCM criteria, the COI barcode successfully identified 69% of all species. The rest of the sequences, approximately 10%, 18% and 3%, remained as ambiguously, mis and unidentified, respectively. The COI barcode does not contain enough information to distinguish Culex ( Cux. ) species.

Classification of sags according to their origin based on the waveform similarity

A statistical method for classification of sags their origin downstream or upstream from the recording point is proposed in this work. The goal is to obtain a statistical model using the sag waveforms useful to characterise one type of sags and to discriminate them from the other type. This model is built on the basis of multi-way principal component analysis an later used to project the available registers in a new space with lower dimension. Thus, a case base of diagnosed sags is built in the projection space. Finally classification is done by comparing new sags against the existing in the case base. Similarity is defined in the projection space using a combination of distances to recover the nearest neighbours to the new sag. Finally the method assigns the origin of the new sag according to the origin of their neighbours

AIDS defining opportunistic infections in patients with high CD4 counts in the combination antiretroviral therapy (cART) era: things ain't what they used to be.

INTRODUCTION: According to reports from observational databases, classic AIDS-defining opportunistic infections (ADOIs) occur in patients with CD4 counts above 500/µL on and off cART. Adjudication of these events is usually not performed. However, ADOIs are often used as endpoints, for example, in analyses on when to start cART. MATERIALS AND METHODS: In the database, Swiss HIV Cohort Study (SHCS) database, we identified 91 cases of ADOIs that occurred from 1996 onwards in patients with the nearest CD4 count >500/µL. Cases of tuberculosis and recurrent bacterial pneumonia were excluded as they also occur in non-immunocompromised patients. Chart review was performed in 82 cases, and in 50 cases we identified CD4 counts within six months before until one month after ADOI and had chart review material to allow an in-depth review. In these 50 cases, we assessed whether (1) the ADOI fulfilled the SHCS diagnostic criteria (www.shcs.ch), and (2) HIV infection with CD4 >500/µL was the main immune-compromising condition to cause the ADOI. Adjudication of cases was done by two experienced clinicians who had to agree on the interpretation. RESULTS: More than 13,000 participants were followed in SHCS in the period of interest. Twenty-four (48%) of the chart-reviewed 50 patients with ADOI and CD4 >500/µL had an HIV RNA <400 copies/mL at the time of ADOI. In the 50 cases, candida oesophagitis was the most frequent ADOI in 30 patients (60%) followed by pneumocystis pneumonia and chronic ulcerative HSV disease (Table 1). Overall chronic HIV infection with a CD4 count >500/µL was the likely explanation for the ADOI in only seven cases (14%). Other reasons (Table 1) were ADOIs occurring during primary HIV infection in 5 (10%) cases, unmasking IRIS in 1 (2%) case, chronic HIV infection with CD4 counts <500/µL near the ADOI in 13 (26%) cases, diagnosis not according to SHCS diagnostic criteria in 7 (14%) cases and most importantly other additional immune-compromising conditions such as immunosuppressive drugs in 14 (34%). CONCLUSIONS: In patients with CD4 counts >500/ µL, chronic HIV infection is the cause of ADOIs in only a minority of cases. Other immuno-compromising conditions are more likely explanations in one-third of the patients, especially in cases of candida oesophagitis. ADOIs in HIV patients with high CD4 counts should be used as endpoints only with much caution in studies based on observational databases.

Risk factors associated with Trypanosoma cruziexposure in domestic dogs from a rural community in Panama

Chagas disease, caused by Trypanosoma cruzi infection, is a zoonosis of humans, wild and domestic mammals, including dogs. In Panama, the main T. cruzi vector is Rhodnius pallescens, a triatomine bug whose main natural habitat is the royal palm, Attalea butyracea. In this paper, we present results from three T. cruzi serological tests (immunochromatographic dipstick, indirect immunofluorescence and ELISA) performed in 51 dogs from 24 houses in Trinidad de Las Minas, western Panama. We found that nine dogs were seropositive (17.6% prevalence). Dogs were 1.6 times more likely to become T. cruziseropositive with each year of age and 11.6 times if royal palms where present in the peridomiciliary area of the dog’s household or its two nearest neighbours. Mouse-baited-adhesive traps were employed to evaluate 12 peridomestic royal palms. All palms were found infested with R. pallescens with an average of 25.50 triatomines captured per palm. Of 35 adult bugs analysed, 88.6% showed protozoa flagellates in their intestinal contents. In addition, dogs were five times more likely to be infected by the presence of an additional domestic animal species in the dog’s peridomiciliary environment. Our results suggest that interventions focused on royal palms might reduce the exposure to T. cruzi infection.

A Factor analysis of hidrochemical composition of Llobregat river basin

Hydrogeological research usually includes some statistical studies devised to elucidate mean background state, characterise relationships among different hydrochemical parameters, and show the influence of human activities. These goals are achieved either by means of a statistical approach or by mixing modelsbetween end-members. Compositional data analysis has proved to be effective with the first approach, but there is no commonly accepted solution to the end-member problem in a compositional framework.We present here a possible solution based on factor analysis of compositions illustrated with a case study.We find two factors on the compositional bi-plot fitting two non-centered orthogonal axes to the most representative variables. Each one of these axes defines a subcomposition, grouping those variables thatlay nearest to it. With each subcomposition a log-contrast is computed and rewritten as an equilibrium equation. These two factors can be interpreted as the isometric log-ratio coordinates (ilr) of three hiddencomponents, that can be plotted in a ternary diagram. These hidden components might be interpreted as end-members.We have analysed 14 molarities in 31 sampling stations all along the Llobregat River and its tributaries, with a monthly measure during two years. We have obtained a bi-plot with a 57% of explained totalvariance, from which we have extracted two factors: factor G, reflecting geological background enhanced by potash mining; and factor A, essentially controlled by urban and/or farming wastewater. Graphicalrepresentation of these two factors allows us to identify three extreme samples, corresponding to pristine waters, potash mining influence and urban sewage influence. To confirm this, we have available analysisof diffused and widespread point sources identified in the area: springs, potash mining lixiviates, sewage, and fertilisers. Each one of these sources shows a clear link with one of the extreme samples, exceptfertilisers due to the heterogeneity of their composition.This approach is a useful tool to distinguish end-members, and characterise them, an issue generally difficult to solve. It is worth note that the end-member composition cannot be fully estimated but only characterised through log-ratio relationships among components. Moreover, the influence of each endmember in a given sample must be evaluated in relative terms of the other samples. These limitations areintrinsic to the relative nature of compositional data

Model d'adequació de l'hàbitat del llop a Catalunya

Des de l’any 2000 es té constància de la presencia del llop a Catalunya. Des de llavors, com a mínim 14 llops diferents han entrat i sortit del territori català, encara que cap d’ells s’ha assentat de manera permanent. L’estudi analitza l’entorn català utilitzant GIS, creant un model d’adequació de l’hàbitat tenint en compte les següents variables: la distància a la carretera més propera, la biomassa disponible a la zona, l’altitud i el tipus i tant per cent de recobriment. El model es basa en la informació obtinguda mitjançant la consulta a experts tant del llop com del territori català, així com en una recerca bibliogràfica sobre l’adequació de l’hàbitat del llop. L’enquesta que es dirigí als experts té en compte els valors que cada variable pot prendre dins l’àrea d’estudi, estableix rangs dels valors de cada variable i pregunta als experts com cada rang pot afectar a l’adequació de l’hàbitat pel llop. Els resultats mostren com bona part de la zona Nord de Catalunya té unes condicions adequades perquè el llop pugui arribar a reproduir-s’hi. Es desenvolupa també una anàlisi dels possibles punts de conflicte humà-llop i una superposició dels espais protegits amb les zones més adequades per l’establiment del llop.

The Geographical Scope of Industrial Location Determinants: An Alternative Approach

This paper considers the estimation of the geographical scope of industrial location determinants. While previous studies impose strong assumptions on the weighting scheme of the spatial neighbour matrix, we propose a exible parametrisation that allows for di fferent (distance-based) de finitions of neighbourhood and di fferent weights to the neighbours. In particular, we estimate how far can reach indirect marginal e ffects and discuss how to report them. We also show that the use of smooth transition functions provides tools for policy analysis that are not available in the traditional threshold modelling. Keywords: count data models, industrial location, smooth transition functions, threshold models. JEL-Codes: C25, C52, R11, R30.

Exploring the conditions for a mortality crisis: bringing context back into the debate

This study engages with the debate over the mortality crises in the former Soviet Union and Central and Eastern Europe by 1) considering at length and as complementary to each other the two most prominent explanations for the post-communist mortality crisis, stress and alcohol consumption; 2) emphasizing the importance of context by exploiting systematic similarities and differences across the region. Differential mortality trajectories reveal three country groups that cluster both spatially and in terms of economic transition experiences. The first group are the countries furthest west in which mortality rates increased minimally after the transition began. The second group experienced a severe increase in mortality rates in the early 1990s, but recovered previous levels within a few years. These countries are located peripherally to Russia and its nearest neighbours. The final group consists of countries that experienced two mortality increases or in which mortality levels had not recovered to pre-transition levels well into the 21st century. Cross-sectional time-series data analyses of men’s and women’s age and cause-specific death rates reveal that the clustering of these countries and their mortality trajectories can be partially explained by the economic context, which is argued to be linked to stress and alcohol consumption. Above and beyond many basic differences in the country groups that are held constant—including geographically and historically shared cultural, lifestyle and social characteristics—poor economic conditions account for a remarkably consistent share of excess age-specific and cause-specific deaths.

Local-scale spatial variation in diversity of social wasps in an Amazonian rain forest in Caxiuanã, Pará, Brazil (Hymenoptera, Vespidae, Polistinae)

Polistine wasps are important in Neotropical ecosystems due to their ubiquity and diversity. Inventories have not adequately considered spatial attributes of collected specimens. Spatial data on biodiversity are important for study and mitigation of anthropogenic impacts over natural ecosystems and for protecting species. We described and analyzed local-scale spatial patterns of collecting records of wasp species, as well as spatial variation of diversity descriptors in a 2500-hectare area of an Amazon forest in Brazil. Rare species comprised the largest fraction of the fauna. Close range spatial effects were detected for most of the more common species, with clustering of presence-data at short distances. Larger spatial lag effects could also be identified in some species, constituting probably cases of exogenous autocorrelation and candidates for explanations based on environmental factors. In a few cases, significant or near significant correlations were found between five species (of Agelaia, Angiopolybia, and Mischocyttarus) and three studied environmental variables: distance to nearest stream, terrain altitude, and the type of forest canopy. However, association between these factors and biodiversity variables were generally low. When used as predictors of polistine richness in a linear multiple regression, only the coefficient for the forest canopy variable resulted significant. Some level of prediction of wasp diversity variables can be attained based on environmental variables, especially vegetation structure. Large-scale landscape and regional studies should be scheduled to address this issue.