Improving generalized regression analysis for the spatial prediction of forest communities

Aim This study used data from temperate forest communities to assess: (1) five different stepwise selection methods with generalized additive models, (2) the effect of weighting absences to ensure a prevalence of 0.5, (3) the effect of limiting absences beyond the environmental envelope defined by presences, (4) four different methods for incorporating spatial autocorrelation, and (5) the effect of integrating an interaction factor defined by a regression tree on the residuals of an initial environmental model. Location State of Vaud, western Switzerland. Methods Generalized additive models (GAMs) were fitted using the grasp package (generalized regression analysis and spatial predictions, http://www.cscf.ch/grasp). Results Model selection based on cross-validation appeared to be the best compromise between model stability and performance (parsimony) among the five methods tested. Weighting absences returned models that perform better than models fitted with the original sample prevalence. This appeared to be mainly due to the impact of very low prevalence values on evaluation statistics. Removing zeroes beyond the range of presences on main environmental gradients changed the set of selected predictors, and potentially their response curve shape. Moreover, removing zeroes slightly improved model performance and stability when compared with the baseline model on the same data set. Incorporating a spatial trend predictor improved model performance and stability significantly. Even better models were obtained when including local spatial autocorrelation. A novel approach to include interactions proved to be an efficient way to account for interactions between all predictors at once. Main conclusions Models and spatial predictions of 18 forest communities were significantly improved by using either: (1) cross-validation as a model selection method, (2) weighted absences, (3) limited absences, (4) predictors accounting for spatial autocorrelation, or (5) a factor variable accounting for interactions between all predictors. The final choice of model strategy should depend on the nature of the available data and the specific study aims. Statistical evaluation is useful in searching for the best modelling practice. However, one should not neglect to consider the shapes and interpretability of response curves, as well as the resulting spatial predictions in the final assessment.

Descriptive Sensory Analysis and Acceptance of Stingless Bee Honey

Stingless bee honey samples were evaluated by sensory descriptive analysis using free choice profile methodology. Appearance, flavor, and aroma were described and the data were treated with Generalized Procrustes Analysis. Individual descriptive terms ranged from 8 to 20. Plotting the samples in a bidimensional plan indicated that appearance attributes (color and viscosity) and sweet, sour and acid flavor were strongly correlated with x axis (Dimension 1) while coconut, wood, acid, sour, and sweet flavor aroma attributes were correlated with y axis (Dimension 2). The affective test was also performed and with the exception of the Melipona scutellaris honey, all the other samples showed good acceptance. Honeys that were described as sweeter and less acid were preferred by nontrained assessors, indicating that the regular consumer recognizes honey produced by Apis mellifera bee as a standard.

Caracterización de poblaciones nativas de maíz mediante Análisis de Procrustes Generalizado y Análisis Factorial Múltiple

El establecimiento de relaciones entre taxones es un paso esencial en el proceso de catalogación y evaluación del material conservado en un Banco de Germoplasma. Existen distintos métodos de evaluación en función del tipo de caracteres estudiados. Cuando el registro de caracteres se repite en el tiempo y en distintos ambientes, se debe separar la variabilidad intrínsecamente genética entre los taxones de aquella que se debe al ambiente, y más aún, de la posible variabilidad debida a la interacción genotipo*ambiente para el posterior establecimiento de relaciones puramente filogenéticas. En el presente trabajo se estudia comparativamente la factibilidad de aplicación de dos estrategias de análisis estadístico para dar solución a este problema. La primera corresponde al análisis tradicional donde se realiza un Análisis de Componentes Principales sobre los caracteres promedios a lo largo de los diferentes ambientes; y la segunda son métodos más complejos en los cuales cada dato es originado por tres modos: individuos, variables y condiciones ambientales, tales como el Análisis Factorial Múltiple y el Análisis de Procrustes Generalizado. Si bien las configuraciones resultantes fueron todas equivalentes, los métodos de tres vías permiten la interpretación de la interacción genotipo*ambiente.

Phenotypic evolution in a venerid bivalve species lineage from the late Middle Miocene Central Paratethys Sea: a multi-approach morphometric analysis

A morphometric analysis was performed for the late Middle Miocene bivalve species lineage of Polititapes tricuspis (Eichwald, 1829) (Veneridae: Tapetini). Specimens from various localities grouped into two stratigraphically successive biozones, i.e. the upper Ervilia Zone and the Sarmatimactra Zone, were investigated using a multi-method approach. A Generalized Procrustes Analysis was computed for fifteen landmarks, covering characteristics of the hinge, muscle scars, and pallial line. The shell outline was separately quantified by applying the Fast Fourier Transform, which redraws the outline by fitting in a combination of trigonometric curves. Shell size was calculated as centroid size from the landmark configuration. Shell thickness, as not covered by either analysis, was additionally measured at the centroid. The analyses showed significant phenotypic differentiation between specimens from the two biozones. The bivalves become distinctly larger and thicker over geological time and develop circular shells with stronger cardinal teeth and a deeper pallial sinus. Data on the paleoenvironmental changes in the late Middle Miocene Central Paratethys Sea suggest the phenotypic shifts to be functional adaptations. The typical habitats for Polititapes changed to extensive, very shallow shores exposed to high wave action and tidal activity. Caused by the growing need for higher mechanical stability, the bivalves produced larger and thicker shells with stronger cardinal teeth. The latter are additionally shifted towards the hinge center to compensate for the lacking lateral teeth and improve stability. The deepening pallial sinus is related to a deeper burrowing habit, which is considered to impede being washed out in the new high-energy settings.

Constant fluctuating asymmetry but not directional asymmetry along the geographic distribution of Drosophila antonietae (Diptera, Drosophilidae)

ABSTRACT The population dynamics of a species tends to change from the core to the periphery of its distribution. Therefore, one could expect peripheral populations to be subject to a higher level of stress than more central populations (the center–periphery hypothesis) and consequently should present a higher level of fluctuating asymmetry. To test these predictions we study asymmetry in wing shape of five populations of Drosophila antonietae collected throughout the distribution of the species using fluctuating asymmetry as a proxy for developmental instability. More specifically, we addressed the following questions: (1) what types of asymmetry occur in populations of D. antonietae? (2) Does the level of fluctuating asymmetry vary among populations? (3) Does peripheral populations have a higher fluctuating asymmetry level than central populations? We used 12 anatomical landmarks to quantify patterns of asymmetry in wing shape in five populations of D. antonietae within the framework of geometric morphometrics. Net asymmetry – a composite measure of directional asymmetry + fluctuating asymmetry – varied significantly among populations. However, once net asymmetry of each population is decomposed into directional asymmetry and fluctuating asymmetry, most of the variation in asymmetry was explained by directional asymmetry alone, suggesting that populations of D. antonietae have the same magnitude of fluctuating asymmetry throughout the geographical distribution of the species. We hypothesize that larval development in rotting cladodes might play an important role in explaining our results. In addition, our study underscores the importance of understanding the interplay between the biology of a species and its geographical patterns of asymmetry.

Analyse de la variation inter-populationnelle du canal obstétrique au sein de trois populations génétiquement différentes

La dystocie représente une des causes principales de la mortalité maternelle et fœtale lors de l’accouchement. Une relation trop étroite entre les dimensions maternelles et fœtales résultent souvent en une césarienne d’urgence. Les manuels obstétricaux actuels proposent qu’une forme gynoïde de la cavité obstétricale est normale car elle favorise une descente fœtale sans difficulté lors de la parturition parce qu’elle est arrondie. Ce mémoire de maîtrise porte sur l’étude de la variation inter-populationnelle de la forme du canal obstétrique chez trois populations génétiquement différentes pour tester si ces populations ont toutes une forme gynoïde ou si elles sont variables morphologiquement. Les cavités pelviennes analysées sont constituées de 59 bassins d’Américaines d’origine européenne, 74 bassins d’Américaines d’origine africaine et 19 bassins de Sadlermiuts. Vingt points anatomiques sur la cavité pelvienne ont été digitalisés à partir du bassin articulé pour permettre une visualisation tridimensionnelle de la forme du canal obstétrique. Une analyse Procrustes suivie d’une analyse en composantes principales ont été faites pour faire ressortir les composantes qui permettent de faire des distinctions entre les trois populations. L’existence d’une relation d’une part entre la taille et/ou la masse et d’autre part la morphologie du canal obstétrique a aussi été explorée. Résultats : La morphologie qui décrit les Sadlermiuts est la forme platypelloïde, soit un diamètre médio-latéral plus grand que le diamètre antéro-postérieur, et ce à tous les détroits. Les Américaines africaines présentent une morphologie qui tend vers une forme gynoïde aux détroits supérieur et moyen, et les Américaines européennes démontrent une morphologie intermédiaire entre les formes gynoïde et androïde au détroit supérieur. Seule la taille du bassin (taille centroïde) est corrélée à des changements de la forme du canal obstétrique pour les détroits supérieur et moyen. Discussion: La forme gynoïde n’est pas représentative de la cavité obstétricale chez les trois populations. Il est possible que la forme platypelloïde des Sadlermiuts soit une conséquence de comportements induisant des changements morphologiques lors du développement squelettique. D’autres études sont requises pour évaluer cette hypothèse.

Determinação do perfil sensorial de amostras comerciais de cachaça

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

A two-centre study on facial morphology in patients with complete bilateral cleft lip and palate at nine years of age

The aim of this study was to compare craniofacial morphology and soft tissue profiles in patients with complete bilateral cleft lip and palate at 9 years of age, treated in two European cleft centres with delayed hard palate closure but different treatment protocols. The cephalometric data of 83 consecutively treated patients were compared (Gothenburg, N=44; Nijmegen, N=39). In total, 18 hard tissue and 10 soft tissue landmarks were digitized by one operator. To determine the intra-observer reliability 20 cephalograms were digitized twice with a monthly interval. Paired t-test, Pearson correlation coefficients and multiple regression models were applied for statistical analysis. Hard and soft tissue data were superimposed using the Generalized Procrustes Analysis. In Nijmegen, the maxilla was protrusive for hard and soft tissue values (P=0.001, P=0.030, respectively) and the maxillary incisors were retroclined (P<0.001), influencing the nasolabial angle, which was increased in comparison with Gothenburg (P=0.004). In conclusion, both centres showed a favourable craniofacial form at 9-10 years of age, although there were significant differences in the maxillary prominence, the incisor inclination and soft tissue cephalometric values. Follow-up of these patients until facial growth has ceased, may elucidate components for outcome improvement.

A longitudinal three-center study of craniofacial morphology at 6 and 12 years of age in patients with complete bilateral cleft lip and palate

In this longitudinal study, the craniofacial morphology and evaluated soft tissue profile changes, at 6 and 12 years of age in patients with complete bilateral cleft lip and palate (CBCLP) were compared. Lateral cephalograms from 148 patients with CBCLP, treated consecutively at three European cleft centers, Gothenburg (n (A) = 37), Nijmegen (n (B) = 26), and Oslo (n (C) = 85), were evaluated. Eighteen hard tissue and ten soft tissue landmarks were digitized. Paired t test, Pearson's correlation coefficients, and multiple regression models were applied for statistical analysis. ANOVA and Tukey-B, as a post hoc test, were used to evaluate the increments and compare centers. Hard and soft tissue data were superimposed using the generalized Procrustes analysis. For Nijmegen, the increments of the variables SNA, ANB, SN-NL, SN-ML, NL-ML, Snss, and Snpg were significantly different than the two other centers (p = 0.041 to <0.001). SNPg increments were significantly different between Nijmegen and Oslo (p = 0.002). The three cleft centers followed different treatment protocols, but the main differences in craniofacial morphology until 12 years of age were the growth pattern and the maxillary and upper incisor variables. Follow-up of these patients until facial growth has ceased, which may elucidate components for improving treatment outcome.

Genetic diversity among alfalfa genotypes (Medicago sativa L.) of non-dormant cultivars using SSR markers and agronomic traits

The aim of this study was to assess genetic diversity among 40 alfalfa (Medicago sativa L.) genotypes of different non-dormant (FD=8) cultivars. Biomass yield, regrowth speed and reaction to spring black stem, lepto leaf spot, and rust were evaluated. Analyses of variances were performed using a mixed model to examine the agronomic variation among individuals. A principal component analysis on standardized agronomic data was performed. Agronomic data were also used to calculate Gower's distance and UPGMA algorithm. For the molecular analysis, six SSR markers were evaluated and 84 alleles were identified. The genetic distance was estimated using standard Nei's distance. Average standard genetic diversity was 0.843, indicating a high degree of variability among genotypes. Finally, a generalized procrustes analysis was performed to calculate the correlation between molecular and agronomic distance, indicating a 65.4% of consensus. This value is likely related to the low number of individuals included in the study, which might have underestimated the real phenotypic variability among genotypes. Despite the low number of individuals and SSR markers analyzed, this study provides a baseline for future diversity studies to identify genetically distant alfalfa individuals or cultivars.

Ordenaciones de material genético a partir de información multidimensional

Nuevas biotecnologías permiten obtener información para caracterizar materiales genéticos a partir de múltiples marcadores, ya sean éstos moleculares y/o morfológicos. La ordenación del material genético a través de la exploración de patrones de variabilidad multidimensionales se aborda mediante diversas técnicas de análisis multivariado. Las técnicas multivariadas de reducción de dimensión (TRD) y la representación gráfica de las mismas cobran sustancial importancia en la visualización de datos multivariados en espacios de baja dimensión ya que facilitan la interpretación de interrelaciones entre las variables (marcadores) y entre los casos u observaciones bajo análisis. Tanto el Análisis de Componentes Principales, como el Análisis de Coordenadas Principales y el Análisis de Procrustes Generalizado son TRD aplicables a datos provenientes de marcadores moleculares y/o morfológicos. Los Árboles de Mínimo Recorrido y los biplots constituyen técnicas para lograr representaciones geométricas de resultados provenientes de TRD. En este trabajo se describen estas técnicas multivariadas y se ilustran sus aplicaciones sobre dos conjuntos de datos, moleculares y morfológicos, usados para caracterizar material genético fúngico.

Generalized canonical correlation analysis of matrices with different row and column orders

A Method is offered that makes it possible to apply generalized canonicalcorrelations analysis (CANCOR) to two or more matrices of different row and column order. The new method optimizes the generalized canonical correlationanalysis objective by considering only the observed values. This is achieved byemploying selection matrices. We present and discuss fit measures to assessthe quality of the solutions. In a simulation study we assess the performance of our new method and compare it to an existing procedure called GENCOM,proposed by Green and Carroll. We find that our new method outperforms the GENCOM algorithm both with respect to model fit and recovery of the truestructure. Moreover, as our new method does not require any type of iteration itis easier to implement and requires less computation. We illustrate the methodby means of an example concerning the relative positions of the political parties inthe Netherlands based on provincial data.

Some results on holomorphic generalized functions

In this note, we present three independent results within generalized complex analysis (in the Colombeau sense). The first of them deals with non-removable singularities; we construct a generalized function u on an open subset Omega of C(n), which is not a holomorphic generalized function on Omega but it is a holomorphic generalized function on Omega\S, where S is a hypersurface contained in Omega. The second result shows the existence of a holomorphic generalized function with prescribed values in the zero-set of a classical holomorphic function. The last result states the existence of a compactly supported solution to the (partial derivative) over bar operator.

Multidimensional Scaling Visualization using Parametric Similarity Indices

In this paper, we apply multidimensional scaling (MDS) and parametric similarity indices (PSI) in the analysis of complex systems (CS). Each CS is viewed as a dynamical system, exhibiting an output time-series to be interpreted as a manifestation of its behavior. We start by adopting a sliding window to sample the original data into several consecutive time periods. Second, we define a given PSI for tracking pieces of data. We then compare the windows for different values of the parameter, and we generate the corresponding MDS maps of ‘points’. Third, we use Procrustes analysis to linearly transform the MDS charts for maximum superposition and to build a global MDS map of “shapes”. This final plot captures the time evolution of the phenomena and is sensitive to the PSI adopted. The generalized correlation, the Minkowski distance and four entropy-based indices are tested. The proposed approach is applied to the Dow Jones Industrial Average stock market index and the Europe Brent Spot Price FOB time-series.