RESUMEN

Theta oscillations, ranging from 4-8 Hz, play a significant role in spatial learning and memory functions during navigation tasks. Frontal theta oscillations are thought to play an important role in spatial navigation and memory. Electroencephalography (EEG) datasets are very complex, making any changes in the neural signal related to behaviour difficult to interpret. However, multiple analytical methods are available to examine complex data structures, especially machine learning-based techniques. These methods have shown high classification performance, and their combination with feature engineering enhances their capability. This paper proposes using hidden Markov and linear mixed effects models to extract features from EEG data. Based on the engineered features obtained from frontal theta EEG data during a spatial navigation task in two key trials (first, last) and between two conditions (learner and non-learner), we analysed the performance of six machine learning methods on classifying learner and non-learner participants. We also analysed how different standardisation methods used to pre-process the EEG data contribute to classification performance. We compared the classification performance of each trial with data gathered from the same subjects, including solely coordinate-based features, such as idle time and average speed. We found that more machine learning methods perform better classification using coordinate-based data. However, only deep neural networks achieved an area under the ROC curve higher than 80% using the theta EEG data alone. Our findings suggest that standardising the theta EEG data and using deep neural networks enhances the classification of learner and non-learner subjects in a spatial learning task.

RESUMEN

This paper focuses on inferring a general class of hidden Markov models (HMMs) using data acquired from experts. Expert-acquired data contain decisions/actions made by humans/users for various objectives, such as navigation data reflecting drivers' behavior, cybersecurity data carrying defender decisions, and biological data containing the biologist's actions (e.g., interventions and experiments). Conventional inference methods rely on temporal changes in data without accounting for expert knowledge. This paper incorporates expert knowledge into the inference of HMMs by modeling expert behavior as an imperfect reinforcement learning agent. The proposed method optimally quantifies experts' perceptions about the system model, which, alongside the temporal changes in data, contributes to the inference process. The proposed inference method is derived through a combination of dynamic programming and optimal recursive Bayesian estimation. The applicability of this method is demonstrated to a wide range of inference criteria, such as maximum likelihood and maximum a posteriori. The performance of the proposed method is investigated through a comprehensive numerical experiment using a benchmark problem and biological networks.

RESUMEN

Single-particle tracking is a powerful technique to investigate the motion of molecules or particles. Here, we review the methods for analyzing the reconstructed trajectories, a fundamental step for deciphering the underlying mechanisms driving the motion. First, we review the traditional analysis based on the mean squared displacement (MSD), highlighting the sometimes-neglected factors potentially affecting the accuracy of the results. We then report methods that exploit the distribution of parameters other than displacements, e.g., angles, velocities, and times and probabilities of reaching a target, discussing how they are more sensitive in characterizing heterogeneities and transient behaviors masked in the MSD analysis. Hidden Markov Models are also used for this purpose, and these allow for the identification of different states, their populations and the switching kinetics. Finally, we discuss a rapidly expanding field-trajectory analysis based on machine learning. Various approaches, from random forest to deep learning, are used to classify trajectory motions, which can be identified by motion models or by model-free sets of trajectory features, either previously defined or automatically identified by the algorithms. We also review free software available for some of the analysis methods. We emphasize that approaches based on a combination of the different methods, including classical statistics and machine learning, may be the way to obtain the most informative and accurate results.

Asunto(s)

RESUMEN

Determining movement parameters for pest insects such as tephritid fruit flies is critical to developing models which can be used to increase the effectiveness of surveillance and control strategies. In this study, harmonic radar was used to track wild-caught male Queensland fruit flies (Qflies), Bactrocera tryoni, in papaya fields. Experiment 1 continuously tracked single flies which were prodded to induce movement. Qfly movements from this experiment showed greater mean squared displacement than predicted by both a simple random walk (RW) or a correlated random walk (CRW) model, suggesting that movement parameters derived from the entire data set do not adequately describe the movement of individual Qfly at all spatial scales or for all behavioral states. This conclusion is supported by both fractal and hidden Markov model (HMM) analysis. Lower fractal dimensions (straighter movement paths) were observed at larger spatial scales (> 2.5 m) suggesting that Qflies have qualitatively distinct movement at different scales. Further, a two-state HMM fit the observed movement data better than the CRW or RW models. Experiment 2 identified individual landing locations, twice a day, for groups of released Qflies, demonstrating that flies could be tracked over longer periods of time.

Asunto(s)

RESUMEN

BACKGROUND: Movement plays a key role in allowing animal species to adapt to sudden environmental shifts. Anthropogenic climate and land use change have accelerated the frequency of some of these extreme disturbances, including megafire. These megafires dramatically alter ecosystems and challenge the capacity of several species to adjust to a rapidly changing landscape. Ungulates and their movement behaviors play a central role in the ecosystem functions of fire-prone ecosystems around the world. Previous work has shown behavioral plasticity is an important mechanism underlying whether large ungulates are able to adjust to recent changes in their environments effectively. Ungulates may respond to the immediate effects of megafire by adjusting their movement and behavior, but how these responses persist or change over time following disturbance is poorly understood. METHODS: We examined how an ecologically dominant ungulate with strong site fidelity, Columbian black-tailed deer (Odocoileus hemionus columbianus), adjusted its movement and behavior in response to an altered landscape following a megafire. To do so, we collected GPS data from 21 individual female deer over the course of a year to compare changes in home range size over time and used resource selection functions (RSFs) and hidden Markov movement models (HMMs) to assess changes in behavior and habitat selection. RESULTS: We found compelling evidence of adaptive capacity across individual deer in response to megafire. Deer avoided exposed and severely burned areas that lack forage and could be riskier for predation immediately following megafire, but they later altered these behaviors to select areas that burned at higher severities, potentially to take advantage of enhanced forage. CONCLUSIONS: These results suggest that despite their high site fidelity, deer can navigate altered landscapes to track rapid shifts in encounter risk with predators and resource availability. This successful adjustment of movement and behavior following extreme disturbance could help facilitate resilience at broader ecological scales.

RESUMEN

This study addresses the environmental problem of PET plastic through in silico bioprospecting for the identification and experimental validation of novel PET degrading eukaryotes through the in silico bioprospectingI of PETases, employing a methodology that combines Hidden Markov Models (HMMs), clustering techniques, molecular docking, and dynamic simulations. A total of 424 putative PETase sequences were identified from 219 eukaryotic organisms, highlighting six sequences with low affinity energies. The Aspergillus luchuensis sequence showed the lowest Gibbs free energy and exhibited stability at different temperatures in molecular dynamics assays. Experimental validation, through a plate clearance assay and HPLC, confirmed PETase activity in three wild-type fungal strains, with A. luchuensis showing the highest efficiency. The results obtained demonstrate the effectiveness of combining computational and experimental approaches as proof of concept to discover and validate eukaryotes with PET-degrading capabilities opening new perspectives for the sustainable management of this type of waste and contributing to its environmental mitigation.

Asunto(s)

RESUMEN

Player movement in rugby league is complex, being spatiotemporal and multifaceted. Modeling this complexity to provide robust measures of player activity and load has proved difficult, with important aspects of player movement yet to be considered. These include the influence of time-varying covariates on player activity and the combination of different dimensions of player movement. Few studies have simultaneously categorized player activity into different activity states and investigated factors influencing the transition between states, or compared player activity and load profiles between matches and training. This study applied hidden Markov models (HMMs)-a data-driven, multivariate approach-to rugby league training and match GPS data to i) demonstrate how HMMs can combine multiple variables in a data-driven way to effectively categorize player movement states, ii) investigate the influence of two time-varying covariates, score difference and elapsed match time on player activity states, and iii) compare player activity and load profiles within and between training and match modalities. HMMs were fitted to player GPS, accelerometer and heart rate data of one English Super League team across 60 training sessions and 35 matches. Distinct activity states were detected for both matches and training, with transitions between states in matches influenced by score difference and elapsed time and clear differences in activity and load profiles between training and matches. HMMs can model the complexity of player movement to effectively profile player activity and load in rugby league and have the potential to facilitate new research across several sports.

We successfully derived player activity and load profiles in both training and match contexts in a data-driven and multivariate way using hidden Markov models.HMMs can be used to investigate the probability of changing between activity states as a function of time-varying covariates, augmenting current activity profiling practice.We discovered key differences between the activity and load profiles between training and matches in rugby league. In particular, a very directed high-speed running state in training that is seldom accessed by players in matches.We demonstrated how visualizing the output of HMMs can provide decision support by facilitating comparisons of activity and load profiles within and between players in matches and training.We posit that the methodology detailed in this paper can become a standardized approach to player activity and load profiling based on player movement data across multiple sports because it is flexible, data-driven, multivariate and statistically robust.

RESUMEN

Profile hidden Markov models (pHMMs) are able to achieve high sensitivity in remote homology search, making them popular choices for detecting novel or highly diverged viruses in metagenomic data. However, many existing pHMM databases have different design focuses, making it difficult for users to decide the proper one to use. In this review, we provide a thorough evaluation and comparison for multiple commonly used profile HMM databases for viral sequence discovery in metagenomic data. We characterized the databases by comparing their sizes, their taxonomic coverage, and the properties of their models using quantitative metrics. Subsequently, we assessed their performance in virus identification across multiple application scenarios, utilizing both simulated and real metagenomic data. We aim to offer researchers a thorough and critical assessment of the strengths and limitations of different databases. Furthermore, based on the experimental results obtained from the simulated and real metagenomic data, we provided practical suggestions for users to optimize their use of pHMM databases, thus enhancing the quality and reliability of their findings in the field of viral metagenomics.

Asunto(s)

RESUMEN

Improving the management of diabetic patients is receiving increasing attention in the health policy agenda due to increasing prevalence in the population and raising pressure on healthcare resources. This paper examines the determinants of healthcare services utilisation in patients with type-2 diabetes, investigating the potential substitution effect of general practice visits on the utilisation of emergency department visits. By using rich longitudinal data from Denmark and a bivariate econometric model, our analysis highlights primary care services that are more effective in preventing emergency department visits and socioeconomic groups of patients with a weak substitution response. Our results suggest that empowering primary care services, such as preventive assessment visits, may contribute to reducing emergency department visits significantly. Moreover, special attention should be devoted to vulnerable groups, such as patients from low socioeconomic background and older patients, who may find more difficult achieving a large substitution response.

Asunto(s)

RESUMEN

A pervasive challenge in neuroscience is testing whether neuronal connectivity changes over time due to specific causes, such as stimuli, events, or clinical interventions. Recent hardware innovations and falling data storage costs enable longer, more naturalistic neuronal recordings. The implicit opportunity for understanding the self-organised brain calls for new analysis methods that link temporal scales: from the order of milliseconds over which neuronal dynamics evolve, to the order of minutes, days, or even years over which experimental observations unfold. This review article demonstrates how hierarchical generative models and Bayesian inference help to characterise neuronal activity across different time scales. Crucially, these methods go beyond describing statistical associations among observations and enable inference about underlying mechanisms. We offer an overview of fundamental concepts in state-space modeling and suggest a taxonomy for these methods. Additionally, we introduce key mathematical principles that underscore a separation of temporal scales, such as the slaving principle, and review Bayesian methods that are being used to test hypotheses about the brain with multiscale data. We hope that this review will serve as a useful primer for experimental and computational neuroscientists on the state of the art and current directions of travel in the complex systems modelling literature.

RESUMEN

Advances in technology have enabled the use of sensors with varied modalities to monitor different parts of systems, each providing diverse levels of information about the underlying system. However, resource limitations and computational power restrict the number of sensors/data that can be processed in real-time in most complex systems. These challenges necessitate the need for selecting/scheduling a subset of sensors to obtain measurements that guarantee the best monitoring objectives. This paper focuses on sensor scheduling for systems modeled by hidden Markov models. Despite the development of several sensor selection and scheduling methods, existing methods tend to be greedy and do not take into account the long-term impact of selected sensors on monitoring objectives. This paper formulates optimal sensor scheduling as a reinforcement learning problem defined over the posterior distribution of system states. Further, the paper derives a deep reinforcement learning policy for offline learning of the sensor scheduling policy, which can then be executed in real-time as new information unfolds. The proposed method applies to any monitoring objective that can be expressed in terms of the posterior distribution of the states (e.g., state estimation, information gain, etc.). The performance of the proposed method in terms of accuracy and robustness is investigated for monitoring the security of networked systems and the health monitoring of gene regulatory networks.

RESUMEN

Step-selection models are widely used to study animals' fine-scale habitat selection based on movement data. Resource preferences and movement patterns, however, often depend on the animal's unobserved behavioral states, such as resting or foraging. As this is ignored in standard (integrated) step-selection analyses (SSA, iSSA), different approaches have emerged to account for such states in the analysis. The performance of these approaches and the consequences of ignoring the states in step-selection analysis, however, have rarely been quantified. We evaluate the recent idea of combining iSSAs with hidden Markov models (HMMs), which allows for a joint estimation of the unobserved behavioral states and the associated state-dependent habitat selection. Besides theoretical considerations, we use an extensive simulation study and a case study on fine-scale interactions of simultaneously tracked bank voles (Myodes glareolus) to compare this HMM-iSSA empirically to both the standard and a widely used classification-based iSSA (i.e., a two-step approach based on a separate prior state classification). Moreover, to facilitate its use, we implemented the basic HMM-iSSA approach in the R package HMMiSSA available on GitHub.

Asunto(s)

RESUMEN

We present BATH, a tool for highly sensitive annotation of protein-coding DNA based on direct alignment of that DNA to a database of protein sequences or profile hidden Markov models (pHMMs). BATH is built on top of the HMMER3 code base, and simplifies the annotation workflow for pHMM-based annotation by providing a straightforward input interface and easy-to-interpret output. BATH also introduces novel frameshift-aware algorithms to detect frameshift-inducing nucleotide insertions and deletions (indels). BATH matches the accuracy of HMMER3 for annotation of sequences containing no errors, and produces superior accuracy to all tested tools for annotation of sequences containing nucleotide indels. These results suggest that BATH should be used when high annotation sensitivity is required, particularly when frameshift errors are expected to interrupt protein-coding regions, as is true with long read sequencing data and in the context of pseudogenes.

RESUMEN

The multilevel hidden Markov model (MHMM) is a promising method to investigate intense longitudinal data obtained within the social and behavioral sciences. The MHMM quantifies information on the latent dynamics of behavior over time. In addition, heterogeneity between individuals is accommodated with the inclusion of individual-specific random effects, facilitating the study of individual differences in dynamics. However, the performance of the MHMM has not been sufficiently explored. We performed an extensive simulation to assess the effect of the number of dependent variables (1-8), number of individuals (5-90), and number of observations per individual (100-1600) on the estimation performance of a Bayesian MHMM with categorical data including various levels of state distinctiveness and separation. We found that using multivariate data generally alleviates the sample size needed and improves the stability of the results. Moreover, including variables only consisting of random noise was generally not detrimental to model performance. Regarding the estimation of group-level parameters, the number of individuals and observations largely compensate for each other. However, only the former drives the estimation of between-individual variability. We conclude with guidelines on the sample size necessary based on the level of state distinctiveness and separation and study objectives of the researcher.

Asunto(s)

RESUMEN

There is some evidence for heterogeneity in attentional processes among individuals with social anxiety. However, there is limited work considering how attentional processes may differ as a mechanism in a naturalistic task-based context (e.g., public speaking). In this secondary analysis we tested attentional heterogeneity among individuals diagnosed with social anxiety disorder (N = 21) in the context of a virtual reality exposure treatment study. Participants completed a public speaking challenge in an immersive 360°-video virtual reality environment with eye tracking at pre-treatment, post-treatment, and at 1-week follow-up. Using a Hidden Markov Model (HMM) approach with clustering we tested whether there were distinct profiles of attention pre-treatment and whether there were changes following the intervention. As a secondary aim we tested whether the distinct attentional profiles at pre-treatment predicted differential treatment outcomes. We found two distinct attentional profiles pre-treatment that we characterized as audience-focused and audience-avoidant. However, by the 1-week follow-up the two profiles were no longer meaningfully different. We found a meaningful difference between HMM groups for fear of public speaking at post-treatment b = -8.54, 95% Highest Density Interval (HDI) [-16.00, -0.90], Bayes Factor (BF) = 8.31 but not at one-week follow-up b = -5.83, 95% HDI [-13.25, 1.81], BF = 2.28. These findings provide support for heterogeneity in attentional processes among socially anxious individuals, but our findings indicate that this may change following treatment. Moreover, our results offer preliminary mechanistic evidence that patterns of avoidance may be specifically related to poorer treatment outcomes for virtual reality exposure therapy.

Asunto(s)

RESUMEN

Presence sensing systems are gaining importance and are utilized in various contexts such as smart homes, Ambient Assisted Living (AAL) and surveillance technology. Typically, these systems utilize motion sensors or cameras that have a limited field of view, leading to potential monitoring gaps within a room. However, humans release carbon dioxide (CO2) through respiration which spreads within an enclosed space. Consequently, an observable rise in CO2 concentration is noted when one or more individuals are present in a room. This study examines an approach to detect the presence or absence of individuals indoors by analyzing the ambient air's CO2 concentration using simple Markov Chain Models. The proposed scheme achieved an accuracy of up to 97% in both experimental and real data demonstrating its efficacy in practical scenarios.

RESUMEN

Response process data from computer-based problem-solving items describe respondents' problem-solving processes as sequences of actions. Such data provide a valuable source for understanding respondents' problem-solving behaviors. Recently, data-driven feature extraction methods have been developed to compress the information in unstructured process data into relatively low-dimensional features. Although the extracted features can be used as covariates in regression or other models to understand respondents' response behaviors, the results are often not easy to interpret since the relationship between the extracted features, and the original response process is often not explicitly defined. In this paper, we propose a statistical model for describing response processes and how they vary across respondents. The proposed model assumes a response process follows a hidden Markov model given the respondent's latent traits. The structure of hidden Markov models resembles problem-solving processes, with the hidden states interpreted as problem-solving subtasks or stages. Incorporating the latent traits in hidden Markov models enables us to characterize the heterogeneity of response processes across respondents in a parsimonious and interpretable way. We demonstrate the performance of the proposed model through simulation experiments and case studies of PISA process data.

RESUMEN

Various computational approaches have been developed to annotate epigenomes on a per-position basis by modeling combinatorial and spatial patterns within epigenomic data. However, such annotations are less suitable for gene-based analyses. We present ChromGene, a method based on a mixture of learned hidden Markov models, to annotate genes based on multiple epigenomic maps across the gene body and flanks. We provide ChromGene assignments for over 100 cell and tissue types. We characterize the mixture components in terms of gene expression, constraint, and other gene annotations. The ChromGene method and annotations will provide a useful resource for gene-based epigenomic analyses.

Asunto(s)

RESUMEN

BACKGROUND AND OBJECTIVES: Computerized Cardiotocography (cCTG) allows to analyze the Fetal Heart Rate (FHR) objectively and thoroughly, providing valuable insights on fetal condition. A challenging but crucial task in this context is the automatic identification of fetal activity and quiet periods within the tracings. Different neural mechanisms are involved in the regulation of the fetal heart, depending on the behavioral states. Thereby, their correct identification has the potential to increase the interpretability and diagnostic capabilities of FHR quantitative analysis. Moreover, the most common pathologies in pregnancy have been associated with variations in the alternation between quiet and activity states. METHODS: We address the problem of fetal states clustering by means of an unsupervised approach, resorting to the use of a multivariate Hidden Markov Models (HMM) with discrete emissions. A fixed length sliding window is shifted on the CTG traces and a small set of features is extracted at each slide. After an encoding procedure, these features become the emissions of a multivariate HMM in which quiet and activity are the hidden states. After an unsupervised training procedure, the model is used to automatically segment signals. RESULTS: The achieved results indicate that our developed model exhibits a high degree of reliability in identifying quiet and activity states within FHR signals. A set of 35 CTG signals belonging to different pregnancies were independently annotated by an expert gynecologist and segmented using the proposed HMM. To avoid any bias, the physician was blinded to the results provided by the algorithm. The overall agreement between the HMM's predictions and the clinician's interpretations was 90%. CONCLUSIONS: The proposed method reliably identified fetal behavioral states, the alternance of which is an important factor in the fetal development. One key strength of our approach lies in the ease of interpreting the obtained results. By utilizing a small set of parameters that are already used in cCTG and possess clear intrinsic meanings, our method provides a high level of explainability. Another significant advantage of our approach is its fully unsupervised learning process. The states identified by our model using the Baum-Welch algorithm are associated with the "Active" and "Quiet" states only after the clustering process, removing the reliance on expert annotations. By autonomously identifying the clusters based solely on the intrinsic characteristics of the signal, our method achieves a more objective evaluation that overcomes the limitations of subjective interpretations. Indeed, we believe it could be integrated in cCTG systems to obtain a more complete signal analysis.

Asunto(s)

RESUMEN

BACKGROUND: Freshwater ecosystems are some of the most affected by biological invasions due, in part, to the introduction of invasive carp worldwide. Where carp have become established, management programs often seek to limit further range expansion into new areas by reducing their movement through interconnected rivers and waterways. Lock and dams are important locations for non-physical deterrents, such as carbon dioxide (CO2), to reduce unwanted fish passage without disrupting human use. The purpose of this study was to evaluate the behavioral responses of common carp (Cyprinus carpio) to non-physical deterrents within a navigation structure on the Fox River, Wisconsin. Acoustic telemetry combined with hidden Markov models (HMMs) was used to analyze variation in carp responses to treatments. Outcomes may inform CO2 effectiveness at preventing invasive carp movement through movement pinch-points. METHODS: Carbon dioxide (CO2) was recently registered as a pesticide in the United States for use as a deterrent to invasive carp movement. As a part of a multi-component study to test a large-scale CO2 delivery system within a navigation lock, we characterized the influence of elevated CO2 and forced water circulation in the lock chamber on carp movements and behavior. Through time-to-event analyses, we described the responses of acoustic-tagged carp to experimental treatments including (1) CO2 injection in water with forced water circulation, (2) forced water circulation without CO2 and (3) no forced water circulation or CO2. We then used hidden Markov models (HMMs) to define fine-scale carp movement and evaluate the relationships between carp behavioral states and CO2 concentration, forced water circulation, and temperature. RESULTS: Forced water circulation with and without CO2 injection were effective at expelling carp from the lock chamber relative to null treatments where no stimulus was applied. A portion of carp exposed to forced water circulation with CO2 transitioned from an exploratory to an encamped behavioral state with shorter step-lengths and a unimodal distribution in turning angles, resulting in some carp remaining in the lock chamber. Whereas carp exposed to forced water circulation only remained primarily in an exploratory behavioral state, resulting in all carp exiting the lock chamber. CONCLUSION: Our findings illustrate the potential of forced water circulation, alone, as a non-physical deterrent and the efficacy of CO2 injection with forced water circulation in expelling carp from a navigation lock. Results demonstrate how acoustic telemetry and HMMs in an experimental context can describe fish behavior and inform management strategies.