Maximum likelihood estimation of marginal pairwise associations with multiple source predictors

Researchers interested in the association of a predictor with an outcome will often collect information about that predictor from more than one source. Standard multiple regression methods allow estimation of the effect of each predictor on the outcome while controlling for the remaining predictors. The resulting regression coefficient for each predictor has an interpretation that is conditional on all other predictors. In settings in which interest is in comparison of the marginal pairwise relationships between each predictor and the outcome separately (e.g., studies in psychiatry with multiple informants or comparison of the predictive values of diagnostic tests), standard regression methods are not appropriate. Instead, the generalized estimating equations (GEE) approach can be used to simultaneously estimate, and make comparisons among, the separate pairwise marginal associations. In this paper, we consider maximum likelihood (ML) estimation of these marginal relationships when the outcome is binary. ML enjoys benefits over GEE methods in that it is asymptotically efficient, can accommodate missing data that are ignorable, and allows likelihood-based inferences about the pairwise marginal relationships. We also explore the asymptotic relative efficiency of ML and GEE methods in this setting.     

Bayesian structural time series for biomedical sensor data: A flexible modeling framework for evaluating interventions

The development of mobile-health technology has the potential to revolutionize personalized medicine. Biomedical sensors (e.g., wearables) can assist with determining treatment plans for individuals, provide quantitative information to healthcare providers, and give objective measurements of health, leading to the goal of precise phenotypic correlates for genotypes. Even though treatments and interventions are becoming more specific and datasets more abundant, measuring the causal impact of health interventions requires careful considerations of complex covariate structures, as well as knowledge of the temporal and spatial properties of the data. Thus, interpreting biomedical sensor data needs to make use of specialized statistical models. Here, we show how the Bayesian structural time series framework, widely used in economics, can be applied to these data. This framework corrects for covariates to provide accurate assessments of the significance of interventions. Furthermore, it allows for a time-dependent confidence interval of impact, which is useful for considering individualized assessments of intervention efficacy. We provide a customized biomedical adaptor tool, MhealthCI, around a specific implementation of the Bayesian structural time series framework that uniformly processes, prepares, and registers diverse biomedical data. We apply the software implementation of MhealthCI to a structured set of examples in biomedicine to showcase the ability of the framework to evaluate interventions with varying levels of data richness and covariate complexity and also compare the performance to other models. Specifically, we show how the framework is able to evaluate an exercise intervention’s effect on stabilizing blood glucose in a diabetes dataset. We also provide a future-anticipating illustration from a behavioral dataset showcasing how the framework integrates complex spatial covariates. Overall, we show the robustness of the Bayesian structural time series framework when applied to biomedical sensor data, highlighting its increasing value for current and future datasets.     

Seagrass Quality Index (SQI), a Water Framework Directive compliant tool for the assessment of transitional and coastal intertidal areas

This work describes a new method for assessing the ecological quality of intertidal seagrass in estuaries and coastal systems, the Seagrass Quality Index (SQI). The design of the SQI aims to fulfil the Water Framework Directive requirements in terms of compliance (e.g., metrics, sampling procedure, pressure relationship, uncertainty of misclassification and comparability to other methodologies in terms of concept). The index includes three common and easy-to-measure structural parameters of seagrass (i.e., the no. of taxa, bed extent and shoot density) combined in a calculation rule that allows the index to report all five of the quality classes (i.e., high, good, moderate, poor and bad). The present study contains analyses of the relationships between the ecosystem-quality results produced by the index and the pressures measured in the system as well as the relationships between the SQI and the seagrass parameters composing it (both the correlation between the SQI and metrics and the SQI sensitivity to the individual variation of each metric). These relationships were tested using a Spearman rank-correlation analysis, producing significant correlations between the biological metrics and the index results as well as between the index results and the environmental quality-pressure category (i.e., the concentration of winter DIN and turbidity). In terms of management, it is possible to apply the methodology on a broad geographical scale in systems where the reference condition for the number of taxa is even higher than one (for the Mondego studied here, the reference value was one species). The tool fulfilled the WFD requirements, had a robust sampling design and proved to be able to track the inertia that usually exists from the moment the pressure is alleviated as well as the biological response that characterises the recovery phase in systems under restoration.     

Identifying Thresholds for Ecosystem-Based Management

Background

One of the greatest obstacles to moving ecosystem-based management (EBM) from concept to practice is the lack of a systematic approach to defining ecosystem-level decision criteria, or reference points that trigger management action.

Methodology/Principal Findings

To assist resource managers and policymakers in developing EBM decision criteria, we introduce a quantitative, transferable method for identifying utility thresholds. A utility threshold is the level of human-induced pressure (e.g., pollution) at which small changes produce substantial improvements toward the EBM goal of protecting an ecosystem''s structural (e.g., diversity) and functional (e.g., resilience) attributes. The analytical approach is based on the detection of nonlinearities in relationships between ecosystem attributes and pressures. We illustrate the method with a hypothetical case study of (1) fishing and (2) nearshore habitat pressure using an empirically-validated marine ecosystem model for British Columbia, Canada, and derive numerical threshold values in terms of the density of two empirically-tractable indicator groups, sablefish and jellyfish. We also describe how to incorporate uncertainty into the estimation of utility thresholds and highlight their value in the context of understanding EBM trade-offs.

Conclusions/Significance

For any policy scenario, an understanding of utility thresholds provides insight into the amount and type of management intervention required to make significant progress toward improved ecosystem structure and function. The approach outlined in this paper can be applied in the context of single or multiple human-induced pressures, to any marine, freshwater, or terrestrial ecosystem, and should facilitate more effective management.     

Ichnofossil tiering analysis of a rhythmically bedded chalk-marl sequence in the Upper Cretaceous of Alabama

Ichnofossils are typically well expressed at bed transitions within rhythmically bedded marine sequences. This is because most discrete biogenic structures therein are actively or passively filled with sediment derived wholly or partly from overlying layers; hence they contrast sharply with ambient sediments. These bed transitions, or 'piped zones', provide the opportunity to assess the nature of infaunal tiering using two different approaches: (1) cross-cutting relationships among recurring ichnotaxa; and (2) apparent maximum penetration depths of piped-zone ichnofossils below their intervals of origination (primary strata). Both approaches were applied to piped zones at chalk-marl and marl-chalk transitions within the rhythmically bedded Upper Cretaceous (Campanian) Demopolis Chalk of western Alabama. Quantitative analysis of cross-cutting relationships indicates the existence of a deep Chondrites tier, a shallow Thalassinoides tier, and a diverse intermediate tier characterized by Anconichnus, Planolites, Taenidium, Teichichnus , and Zoophycos. Analyses of penetration depths similarly document the presence of a deep Chondrites tier but indicate that the Thalassinoides producer occupied the same tier as the tracemakers of the remaining five ichnotaxa. The latter approach, which should be broadly applicable to heterolithic sequences of various sorts, appears to yield a more accurate picture of tiering because it provides a more direct reading of the depths of activity of tracemakers within a substrate, and because results are less likely to be influenced by behavioral factors (e.g., phobotaxis and preferential sediment exploitation) or other parameters (e.g., burrow size differences) that can impact cross-cutting relationships. Results of both approaches indicate that the tiering structure in pelagic carbonate substrates may not be as organized or complex as indicated by previous qualitative analyses of cross-cutting relationships in analogous chalk sequences.     

Bio-SCoRes: A Smorgasbord Architecture for Coreference Resolution in Biomedical Text

Coreference resolution is one of the fundamental and challenging tasks in natural language processing. Resolving coreference successfully can have a significant positive effect on downstream natural language processing tasks, such as information extraction and question answering. The importance of coreference resolution for biomedical text analysis applications has increasingly been acknowledged. One of the difficulties in coreference resolution stems from the fact that distinct types of coreference (e.g., anaphora, appositive) are expressed with a variety of lexical and syntactic means (e.g., personal pronouns, definite noun phrases), and that resolution of each combination often requires a different approach. In the biomedical domain, it is common for coreference annotation and resolution efforts to focus on specific subcategories of coreference deemed important for the downstream task. In the current work, we aim to address some of these concerns regarding coreference resolution in biomedical text. We propose a general, modular framework underpinned by a smorgasbord architecture (Bio-SCoRes), which incorporates a variety of coreference types, their mentions and allows fine-grained specification of resolution strategies to resolve coreference of distinct coreference type-mention pairs. For development and evaluation, we used a corpus of structured drug labels annotated with fine-grained coreference information. In addition, we evaluated our approach on two other corpora (i2b2/VA discharge summaries and protein coreference dataset) to investigate its generality and ease of adaptation to other biomedical text types. Our results demonstrate the usefulness of our novel smorgasbord architecture. The specific pipelines based on the architecture perform successfully in linking coreferential mention pairs, while we find that recognition of full mention clusters is more challenging. The corpus of structured drug labels (SPL) as well as the components of Bio-SCoRes and some of the pipelines based on it are publicly available at https://github.com/kilicogluh/Bio-SCoRes. We believe that Bio-SCoRes can serve as a strong and extensible baseline system for coreference resolution of biomedical text.     

Optical trait indicators for remote sensing of plant species composition: Predictive power and seasonal variability

Most plant species feature similar biochemical compositions and thus similar spectral signals. Still, empirical evidence suggests that the spectral discrimination of species and plant assemblages is possible. Success depends on the presence or absence of faint but detectable differences in biochemical (e.g., pigments, leaf water and dry matter content) and structural properties (e.g., leaf area, angle, and leaf structure), i.e., optical traits. A systematic analysis of the contributions and spatio-temporal variability of optical traits for the remote sensing of organismic vegetation patterns has not yet been conducted. We thus use time series of optical trait values retrieved from the reflectance signal using physical models (optical trait indicators, OTIs) to answer the following questions: How are optical traits related among patterns of floristic composition and reflectance? How variable are these relations in space and time? Are OTIs suitable predictors of plant species composition?We conducted a case study of three temperate open study sites with semi-natural vegetation. The canopy reflectance of permanent vegetation plots was measured on multiple dates over the vegetation period using a field spectrometer. We recorded the cover fractions of all plant species found in the vegetation plots and extracted gradients of species composition from these data. The physical PROSAIL leaf and canopy optical properties model was inverted with random forest regression models to retrieve time series of OTIs for each plot from the reflectance spectra. We analyzed these data sets using correlation analyses. This approach allowed us to assess the distribution of optical traits across gradients of species composition. The predictive performance of OTIs was tested in relation to canopy reflectance using random forest models.OTIs showed pronounced relationships with floristic patterns in all three study sites. These relationships were subject to considerable temporal variability. Such variability was driven by short-term vegetation dynamics introduced by local resource stress. In 72% of all cases OTIs out-performed the original canopy reflectance spectra as indicators of plant species composition. OTIs are also easier to interpret in an ecological sense than spectral bands or features. We thus conclude that optical traits retrieved from reflectance data have a high indicative value for ecological research and applications.     

Analyzing the sequence-structure relationship of a library of local structural prototypes

We present a thorough analysis of the relation between amino acid sequence and local three-dimensional structure in proteins. A library of overlapping local structural prototypes was built using an unsupervised clustering approach called “hybrid protein model” (HPM). The HPM carries out a multiple structural alignment of local folds from a non-redundant protein structure databank encoded into a structural alphabet composed of 16 protein blocks (PBs). Following previous research focusing on the HPM protocol, we have considered gaps in the local structure prototype. This methodology allows to have variable length fragments. Hence, 120 local structure prototypes were obtained. Twenty-five percent of the protein fragments learnt by HPM had gaps.An investigation of tight turns suggested that they are mainly derived from three PB series with precise locations in the HPM. The amino acid information content of the whole conformational classes was tackled by multivariate methods, e.g., canonical correlation analysis. It points out the presence of seven amino acid equivalence classes showing high propensities for preferential local structures. In the same way, definition of “contrast factors” based on sequence-structure properties underline the specificity of certain structural prototypes, e.g., the dependence of Gly or Asn-rich turns to a limited number of PBs, or, the opposition between Pro-rich coils to those enriched in Ser, Thr, Asn and Glu. These results are so useful to analyze the sequence-structure relationships, but could also be used to improve fragment-based method for protein structure prediction from sequence.     

LB3D: a protein three-dimensional substructure search program based on the lower bound of a root mean square deviation value

Searching for protein structure-function relationships using three-dimensional (3D) structural coordinates represents a fundamental approach for determining the function of proteins with unknown functions. Since protein structure databases are rapidly growing in size, the development of a fast search method to find similar protein substructures by comparison of protein 3D structures is essential. In this article, we present a novel protein 3D structure search method to find all substructures with root mean square deviations (RMSDs) to the query structure that are lower than a given threshold value. Our new algorithm runs in O(m + N/m(0.5)) time, after O(N log N) preprocessing, where N is the database size and m is the query length. The new method is 1.8-41.6 times faster than the practically best known O(N) algorithm, according to computational experiments using a huge database (i.e., >20,000,000 C-alpha coordinates).     

Nomograms Aid Interpretation of Complex Regression Models

ABSTRACT Ecologists often develop complex regression models that include multiple categorical and continuous variables, interactions among predictors, and nonlinear relationships between the response and predictor variables. Nomograms, which are graphical devices for presenting mathematical functions and calculating output values, can aid biologists in interpreting and presenting these complex models. To illustrate benefits of nomograms, we developed a logistic regression model of elk (Cervus elaphus) resource selection. With this model, we demonstrated how a nomogram helps scientists and managers interpret interactions among variables, compare the relative biological importance of variables, and examine predicted shapes of relationships (e.g., linear vs. nonlinear) between response and predictor variables. Although our example focused on logistic regression, nomograms are equally useful for other linear and nonlinear models. Regardless of the approach used for model development, nomograms and other graphical summaries can help scientists and managers develop, interpret, and apply statistical models.     

Exploring and linking biomedical resources through multidimensional semantic spaces

BACKGROUND: The semantic integration of biomedical resources is still a challenging issue which is required for effective information processing and data analysis. The availability of comprehensive knowledge resources such as biomedical ontologies and integrated thesauri greatly facilitates this integration effort by means of semantic annotation, which allows disparate data formats and contents to be expressed under a common semantic space. In this paper, we propose a multidimensional representation for such a semantic space, where dimensions regard the different perspectives in biomedical research (e.g., population, disease, anatomy and protein/genes). RESULTS: This paper presents a novel method for building multidimensional semantic spaces from semantically annotated biomedical data collections. This method consists of two main processes: knowledge and data normalization. The former one arranges the concepts provided by a reference knowledge resource (e.g., biomedical ontologies and thesauri) into a set of hierarchical dimensions for analysis purposes. The latter one reduces the annotation set associated to each collection item into a set of points of the multidimensional space. Additionally, we have developed a visual tool, called 3D-Browser, which implements OLAP-like operators over the generated multidimensional space. The method and the tool have been tested and evaluated in the context of the Health-e-Child (HeC) project. Automatic semantic annotation was applied to tag three collections of abstracts taken from PubMed, one for each target disease of the project, the Uniprot database, and the HeC patient record database. We adopted the UMLS Meta-thesaurus 2010AA as the reference knowledge resource. CONCLUSIONS: Current knowledge resources and semantic-aware technology make possible the integration of biomedical resources. Such an integration is performed through semantic annotation of the intended biomedical data resources. This paper shows how these annotations can be exploited for integration, exploration, and analysis tasks. Results over a real scenario demonstrate the viability and usefulness of the approach, as well as the quality of the generated multidimensional semantic spaces.     

A machine‐learning approach for extending classical wildlife resource selection analyses

Resource selection functions (RSFs) are tremendously valuable for ecologists and resource managers because they quantify spatial patterns in resource utilization by wildlife, thereby facilitating identification of critical habitat areas and characterizing specific habitat features that are selected or avoided. RSFs discriminate between known‐use resource units (e.g., telemetry locations) and available (or randomly selected) resource units based on an array of environmental features, and in their standard form are performed using logistic regression. As generalized linear models, standard RSFs have some notable limitations, such as difficulties in accommodating nonlinear (e.g., humped or threshold) relationships and complex interactions. Increasingly, ecologists are using flexible machine‐learning methods (e.g., random forests, neural networks) to overcome these limitations. Herein, we investigate the seasonal resource selection patterns of mule deer (Odocoileus hemionus) by comparing a logistic regression framework with random forest (RF), a popular machine‐learning algorithm. Random forest (RF) models detected nonlinear relationships (e.g., optimal ranges for slope and elevation) and complex interactions which would have been very challenging to discover and characterize using standard model‐based approaches. Compared with standard RSF models, RF models exhibited improved predictive skill, provided novel insights about resource selection patterns of mule deer, and, when projected across a relevant geographic space, manifested notable differences in predicted habitat suitability. We recommend that wildlife researchers harness the strengths of machine‐learning tools like RF in addition to “classical” tools (e.g., mixed‐effects logistic regression) for evaluating resource selection, especially in cases where extensive telemetry data sets are available.     

Guidelines for the design and statistical analysis of experiments using laboratory animals

For ethical and economic reasons, it is important to design animal experiments well, to analyze the data correctly, and to use the minimum number of animals necessary to achieve the scientific objectives---but not so few as to miss biologically important effects or require unnecessary repetition of experiments. Investigators are urged to consult a statistician at the design stage and are reminded that no experiment should ever be started without a clear idea of how the resulting data are to be analyzed. These guidelines are provided to help biomedical research workers perform their experiments efficiently and analyze their results so that they can extract all useful information from the resulting data. Among the topics discussed are the varying purposes of experiments (e.g., exploratory vs. confirmatory); the experimental unit; the necessity of recording full experimental details (e.g., species, sex, age, microbiological status, strain and source of animals, and husbandry conditions); assigning experimental units to treatments using randomization; other aspects of the experiment (e.g., timing of measurements); using formal experimental designs (e.g., completely randomized and randomized block); estimating the size of the experiment using power and sample size calculations; screening raw data for obvious errors; using the t-test or analysis of variance for parametric analysis; and effective design of graphical data.     

Optimization of economic policies and investment projects using a fuzzy logic based cost-effectiveness model of coral reef quality: empirical results for Montego Bay, Jamaica

 For effective mitigation of human impacts, quantitative models are required that facilitate a comprehensive analysis of the effects of human activity on reefs. Fuzzy logic procedures generate a complex dose-response surface that models the relationships among coral abundance and various inputs (e.g., physical damage, sedimentation, nutrient influx), within the context of the abiotic marine environment. This is linked to a nonlinear economic structure incorporating technical interventions (e.g., pollution treatment) and policy interventions (e.g., taxation) in eight economic sectors. Optimization provides insights into the most cost-effective means for protecting coral reefs under different reef quality targets. The research demonstrates that: (1) it is feasible to use fuzzy logic to model complex interactions in coral reef ecosystems; and, (2) conventional economic procedures for modeling cost-effectiveness can result in sub-optimal policy choices when applied to complex systems such as coral reefs. In Montego Bay, Jamaica, up to a 20% increase in coral abundance may be achievable through using appropriate policy measures having a present value cost of US$153 million over 25 years; a 10% increase is achievable at a cost of US$12 million. Accepted: 20 June 1999     

Perspectives on the ecomorphology of bony fishes

Synopsis The field of ecomorphology has a long history with early roots in Europe. In this half of the century the application of ecomorphology to the biology of fishes has developed in the former Soviet Union, Poland and Czechoslovakia, The Netherlands, and in North America. While the specific approaches vary among countries, many North American studies begin by comparing morphological variation with variation in ecological characteristics at the intra or interspecific levels. These initial correlative studies form the ground work for hypotheses that explore the mechanistic underpinnings of the observed ecomorphological associations. Supporting these mechanistic hypotheses are insights from functional studies which demonstrate the limits to potential resource use resulting from a particular morphology; however, the actual resource use is likely to be more limited due to additional constraints provided by internal (e.g., behavior, physiology) and external (e.g., resource abundance, predator distribution) factors. The results from performance studies in the laboratory or field can be used to test specific ecomorphological hypotheses developed from the initial correlational and functional studies. Such studies may, but rarely do, incorporate an ontogenetic analysis of the ecomorphological association to determine their effect on performance. Finally, input from phylogenetic analyses allow an investigator to examine the evolution of specific features and to assess the rates and directionality of character evolution. The structural and ecological diversity of fishes provides a fertile ground to investigate these interactions. The contributions in this volume highlight some of the specific directions for ecomorphological research covering a variety of biological processes in fishes. These include foraging, locomotion, reproduction, respiration, and sensory systems. Running through these papers are new insights into universal ecomorphological issues, i.e., the relationships between form and ecological role and the factors that modify these relationships.     

Development and Validation of the Behavioral Tendencies Questionnaire

At a fundamental level, taxonomy of behavior and behavioral tendencies can be described in terms of approach, avoid, or equivocate (i.e., neither approach nor avoid). While there are numerous theories of personality, temperament, and character, few seem to take advantage of parsimonious taxonomy. The present study sought to implement this taxonomy by creating a questionnaire based on a categorization of behavioral temperaments/tendencies first identified in Buddhist accounts over fifteen hundred years ago. Items were developed using historical and contemporary texts of the behavioral temperaments, described as “Greedy/Faithful”, “Aversive/Discerning”, and “Deluded/Speculative”. To both maintain this categorical typology and benefit from the advantageous properties of forced-choice response format (e.g., reduction of response biases), binary pairwise preferences for items were modeled using Latent Class Analysis (LCA). One sample (n₁ = 394) was used to estimate the item parameters, and the second sample (n₂ = 504) was used to classify the participants using the established parameters and cross-validate the classification against multiple other measures. The cross-validated measure exhibited good nomothetic span (construct-consistent relationships with related measures) that seemed to corroborate the ideas present in the original Buddhist source documents. The final 13-block questionnaire created from the best performing items (the Behavioral Tendencies Questionnaire or BTQ) is a psychometrically valid questionnaire that is historically consistent, based in behavioral tendencies, and promises practical and clinical utility particularly in settings that teach and study meditation practices such as Mindfulness Based Stress Reduction (MBSR).     

miniTUBA: medical inference by network integration of temporal data using Bayesian analysis

MOTIVATION: Many biomedical and clinical research problems involve discovering causal relationships between observations gathered from temporal events. Dynamic Bayesian networks are a powerful modeling approach to describe causal or apparently causal relationships, and support complex medical inference, such as future response prediction, automated learning, and rational decision making. Although many engines exist for creating Bayesian networks, most require a local installation and significant data manipulation to be practical for a general biologist or clinician. No software pipeline currently exists for interpretation and inference of dynamic Bayesian networks learned from biomedical and clinical data. RESULTS: miniTUBA is a web-based modeling system that allows clinical and biomedical researchers to perform complex medical/clinical inference and prediction using dynamic Bayesian network analysis with temporal datasets. The software allows users to choose different analysis parameters (e.g. Markov lags and prior topology), and continuously update their data and refine their results. miniTUBA can make temporal predictions to suggest interventions based on an automated learning process pipeline using all data provided. Preliminary tests using synthetic data and laboratory research data indicate that miniTUBA accurately identifies regulatory network structures from temporal data. AVAILABILITY: miniTUBA is available at http://www.minituba.org.     

Ratios as a size adjustment in morphometrics

Simple ratios in which a measurement variable is divided by a size variable are commonly used but known to be inadequate for eliminating size correlations from morphometric data. Deficiencies in the simple ratio can be alleviated by incorporating regression coefficients describing the bivariate relationship between the measurement and size variables. Recommendations have included: 1) subtracting the regression intercept to force the bivariate relationship through the origin (intercept-adjusted ratios); 2) exponentiating either the measurement or the size variable using an allometry coefficient to achieve linearity (allometrically adjusted ratios); or 3) both subtracting the intercept and exponentiating (fully adjusted ratios). These three strategies for deriving size-adjusted ratios imply different data models for describing the bivariate relationship between the measurement and size variables (i.e., the linear, simple allometric, and full allometric models, respectively). Algebraic rearrangement of the equation associated with each data model leads to a correctly formulated adjusted ratio whose expected value is constant (i.e., size correlation is eliminated). Alternatively, simple algebra can be used to derive an expected value function for assessing whether any proposed ratio formula is effective in eliminating size correlations. Some published ratio adjustments were incorrectly formulated as indicated by expected values that remain a function of size after ratio transformation. Regression coefficients incorporated into adjusted ratios must be estimated using least-squares regression of the measurement variable on the size variable. Use of parameters estimated by any other regression technique (e.g., major axis or reduced major axis) results in residual correlations between size and the adjusted measurement variable. Correctly formulated adjusted ratios, whose parameters are estimated by least-squares methods, do control for size correlations. The size-adjusted results are similar to those based on analysis of least-squares residuals from the regression of the measurement on the size variable. However, adjusted ratios introduce size-related changes in distributional characteristics (variances) that differentially alter relationships among animals in different size classes. © 1993 Wiley-Liss, Inc.     

Networks of Causal Linkage Between Eigenmodes Characterize Behavioral Dynamics of Caenorhabditis elegans

Behavioral phenotyping of model organisms has played an important role in unravelling the complexities of animal behavior. Techniques for classifying behavior often rely on easily identified changes in posture and motion. However, such approaches are likely to miss complex behaviors that cannot be readily distinguished by eye (e.g., behaviors produced by high dimensional dynamics). To explore this issue, we focus on the model organism Caenorhabditis elegans, where behaviors have been extensively recorded and classified. Using a dynamical systems lens, we identify high dimensional, nonlinear causal relationships between four basic shapes that describe worm motion (eigenmodes, also called “eigenworms”). We find relationships between all pairs of eigenmodes, but the timescales of the interactions vary between pairs and across individuals. Using these varying timescales, we create “interaction profiles” to represent an individual’s behavioral dynamics. As desired, these profiles are able to distinguish well-known behavioral states: i.e., the profiles for foraging individuals are distinct from those of individuals exhibiting an escape response. More importantly, we find that interaction profiles can distinguish high dimensional behaviors among divergent mutant strains that were previously classified as phenotypically similar. Specifically, we find it is able to detect phenotypic behavioral differences not previously identified in strains related to dysfunction of hermaphrodite-specific neurons.