RESUMO
The Burridge-Knopoff model implements an earthquake fault as a mechanical block-spring chain. While numerical studies of the model are abundant, experimental investigations are limited to a two-blocks, analog electronic implementation that was proposed by drawing an analogy between mechanical and electrical quantities. Although elegant, this approach is not versatile, mostly because of its heavy reliance on inductors. Here, we propose an alternative, inductorless implementation of the same system. The experimental characterization of the proposed circuit shows very good agreement with theoretical predictions. Besides periodic oscillations, the circuit exhibits a chaotic regime: the corresponding markers of chaoticity, namely, the correlation dimension and the maximum Lyapunov exponent, were experimentally assessed to be consistent with those provided by numerical simulations. The improved versatility and scalability of the circuit is expected to allow for experimental implementations of the Burridge-Knopoff model with a large number of blocks. In addition, the circuit can be used as the basic element of scalable platforms to investigate the dynamics of networks of oscillators and related phenomena.
RESUMO
In the last decade permutation entropy (PE) has become a popular tool to analyze the degree of randomness within a time series. In typical applications, changes in the dynamics of a source are inferred by observing changes of PE computed on different time series generated by that source. However, most works neglect the crucial question related to the statistical significance of these changes. The main reason probably lies in the difficulty of assessing, out of a single time series, not only the PE value, but also its uncertainty. In this paper we propose a method to overcome this issue by using generation of surrogate time series. The analysis conducted on both synthetic and experimental time series shows the reliability of the approach, which can be promptly implemented by means of widely available numerical tools. The method is computationally affordable for a broad range of users.
RESUMO
The task of identifying and characterizing network structures out of experimentally observed time series is tackled by implementing different solutions, ranging from entropy-based techniques to the evaluation of the significance of observed correlation estimators. Among the metrics that belong to the first class, mutual information is of major importance due to the relative simplicity of implementation and its relying on the crucial concept of entropy. With regard to the second class, a method that allows us to assess the connectivity strength of a link in terms of a time scale of its observability via the significance estimate of measured cross correlation was recently shown to provide a reliable tool to study network structures. In this paper, we investigate the relationship between this last metric and mutual information by simultaneously assessing both metrics on large sets of data extracted from three experimental contexts, human brain magnetoencephalography, human brain electroencephalography, and surface wind measurements carried out on a small regional scale, as well as on simulated coupled, auto-regressive processes. We show that the relationship is well described by a power law and provide a theoretical explanation based on a simple noise and signal model. Besides further upholding the reliability of cross-correlation time scale of observability, the results show that the combined use of this metric and mutual information can be used as a valuable tool to identify and characterize connectivity links in a wide range of experimental contexts.
Assuntos
Eletroencefalografia , Magnetoencefalografia , Encéfalo , Entropia , Humanos , Reprodutibilidade dos TestesRESUMO
We analyze the permutation entropy of deterministic chaotic signals affected by a weak observational noise. We investigate the scaling dependence of the entropy increase on both the noise amplitude and the window length used to encode the time series. In order to shed light on the scenario, we perform a multifractal analysis, which allows highlighting the emergence of many poorly populated symbolic sequences generated by the stochastic fluctuations. We finally make use of this information to reconstruct the noiseless permutation entropy. While this approach works quite well for Hénon and tent maps, it is much less effective in the case of hyperchaos. We argue about the underlying motivations.
RESUMO
The detection of an underlying chaotic behavior in experimental recordings is a longstanding issue in the field of nonlinear time series analysis. Conventional approaches require the assessment of a suitable dimension and lag pair to embed a given input sequence and, thereupon, the estimation of dynamical invariants to characterize the underlying source. In this work, we propose an alternative approach to the problem of identifying chaos, which is built upon an improved method for optimal embedding. The core of the new approach is the analysis of an input sequence on a lattice of embedding pairs whose results provide, if any, evidence of a finite-dimensional, chaotic source generating the sequence and, if such evidence is present, yield a set of equivalently suitable embedding pairs to embed the sequence. The application of this approach to two experimental case studies, namely, an electronic circuit and magnetoencephalographic recordings of the human brain, highlights how it can make up a powerful tool to detect chaos in complex systems.
RESUMO
Many studies in nonlinear science heavily rely on surrogate-based hypothesis testing to provide significance estimations of analysis results. Among the complex data produced by nonlinear systems, spike trains are a class of sequences requiring algorithms for surrogate generation that are typically more sophisticated and computationally demanding than methods developed for continuous signals. Although algorithms to specifically generate surrogate spike trains exist, the availability of open-source, portable implementations is still incomplete. In this paper, we introduce the SpiSeMe (Spike Sequence Mime) software package that implements four algorithms for the generation of surrogate data out of spike trains and more generally out of any sequence of discrete events. The purpose of the package is to provide a unified and portable toolbox to carry out surrogate generation on point-process data. Code is provided in three languages, namely, C++, Matlab, and Python, thus allowing straightforward integration of package functions into most analysis pipelines.
Assuntos
Idioma , Software , Potenciais de Ação , Algoritmos , Simulação por Computador , Humanos , Modelos NeurológicosRESUMO
The study of many dynamical systems relies on the analysis of experimentally-recorded sequences of events for which information is encoded in the sequence of interevent intervals. A correct interpretation of the results of the application of analytical techniques to these sequences requires the assessment of statistical significance. In most cases, the corresponding null-hypothesis distribution is unknown, thus forbidding an evaluation of the significance. An alternative solution, which is efficient in the case of continuous signals, is provided by the generation of surrogate data that share statistical and spectral properties with the original dataset. However, in the case of event sequences, the available algorithms for the generation of surrogate data can become cumbersome and computationally demanding. In this work, we present a new method for the generation of surrogate event sequences that relies on the joint distribution of successive interevent intervals. Our method, which was tested on both synthetic and experimental sequences, performs equally well or even better than conventional methods in terms of interevent interval distribution and autocorrelation while abating the computational time by at least one order of magnitude.
RESUMO
Fractal structures pervade nature and are receiving increasing engineering attention towards the realization of broadband resonators and antennas. We show that fractal resonators can support the emergence of high-dimensional chaotic dynamics even in the context of an elementary, single-transistor oscillator circuit. Sierpinski gaskets of variable depth are constructed using discrete capacitors and inductors, whose values are scaled according to a simple sequence. It is found that in regular fractals of this kind, each iteration effectively adds a conjugate pole/zero pair, yielding gradually more complex and broader frequency responses, which can also be implemented as much smaller Foster equivalent networks. The resonators are instanced in the circuit as one-port devices, replacing the inductors found in the initial version of the oscillator. By means of a highly simplified numerical model, it is shown that increasing the fractal depth elevates the dimension of the chaotic dynamics, leading to high-order hyperchaos. This result is overall confirmed by SPICE simulations and experiments, which however also reveal that the non-ideal behavior of physical components hinders obtaining high-dimensional dynamics. The issue could be practically mitigated by building the Foster equivalent networks rather than the verbatim fractals. Furthermore, it is shown that considerably more complex resonances, and consequently richer dynamics, can be obtained by rendering the fractal resonators irregular through reshuffling the inductors, or even by inserting a limited number of focal imperfections. The present results draw attention to the potential usefulness of fractal resonators for generating high-dimensional chaotic dynamics, and underline the importance of irregularities and component non-idealities.
RESUMO
BACKGROUND: During the last decade, many scientific works have concerned the possible use of miRNA levels as diagnostic and prognostic tools for different kinds of cancer. The development of reliable classifiers requires tackling several crucial aspects, some of which have been widely overlooked in the scientific literature: the distribution of the measured miRNA expressions and the statistical uncertainty that affects the parameters that characterize a classifier. In this paper, these topics are analysed in detail by discussing a model problem, i.e. the development of a Bayesian classifier that, on the basis of the expression of miR-205, miR-21 and snRNA U6, discriminates samples into two classes of pulmonary tumors: adenocarcinomas and squamous cell carcinomas. RESULTS: We proved that the variance of miRNA expression triplicates is well described by a normal distribution and that triplicate averages also follow normal distributions. We provide a method to enhance a classifiers' performance by exploiting the correlations between the class-discriminating miRNA and the expression of an additional normalized miRNA. CONCLUSIONS: By exploiting the normal behavior of triplicate variances and averages, invalid samples (outliers) can be identified by checking their variability via chi-square test or their displacement by the respective population mean via Student's t-test. Finally, the normal behavior allows to optimally set the Bayesian classifier and to determine its performance and the related uncertainty.
Assuntos
Adenocarcinoma/classificação , Teorema de Bayes , Carcinoma de Células Escamosas/classificação , Marcadores Genéticos , Neoplasias Pulmonares/classificação , MicroRNAs/genética , Modelos Estatísticos , Adenocarcinoma/genética , Carcinoma de Células Escamosas/genética , Regulação Neoplásica da Expressão Gênica , Humanos , Neoplasias Pulmonares/genética , Prognóstico , SoftwareRESUMO
Alpha-synuclein is a natively unfolded protein widely expressed in neurons at the presynaptic level. It is linked to Parkinson's disease by two lines of evidence: amyloid fibrils of the protein accumulate in patients' brains and three genetic mutants cause autosomal dominant forms of the disease. The biological role of the protein and the mechanisms involved in the etiopathogenesis of Parkinson's disease are still unknown. Membrane binding causes the formation of an amphipathic alpha-helix, which lies on the surface without crossing the bilayer. Recent observations however reported that the application of a voltage induces a pore-like activity of alpha-synuclein. This study aims to characterize the pore forming activity of the protein starting from its monomeric form. In particular, experiments with planar lipid membranes allowed recording of conductance activity bursts with a defined and reproducible fingerprint. Additional experiments with deletion mutants and covalently bound alpha-synuclein dimers were performed to understand both pore assembly and stoichiometry. The information acquired allowed formulation of a model for pore formation at different conductance levels.
Assuntos
Membrana Celular , alfa-Sinucleína/metabolismo , Membrana Celular/metabolismo , Dimerização , Bicamadas LipídicasRESUMO
BACKGROUND: Several studies investigated the role of inflammation in the etiopathogenesis of mood disorders. The aim of our cross-sectional study is evaluating baseline high-sensitivity C-reactive-protein (hsCRP) levels in a cohort of unipolar and bipolar depressive inpatients, in relation with psychopathological, temperamental and chronotype features. METHODS: Among 313 screened inpatients, we retrospectively recruited 133 moderate-to-severe depressive patients who were assessed for hsCRP levels, chronotype with Morningness-Eveningness Questionnaire (MEQ) and affective temperament with Temperament Evaluation of Memphis, Pisa, Paris and San Diego (TEMPS). LIMITATIONS: The cross-sectional and retrospective design of the study, the small sample size, the exclusion of hypomanic, maniac and euthymic bipolar patients. RESULTS: hsCRP levels were significantly higher among those with previous suicide attempt (p = 0.05), death (p = 0.018) and self-harm/self-injury thoughts (p = 0.011). Linear regression analyses, adjusted for all covariates, demonstrated that higher scores at the TEMPS-M depressive, while lower scores at the hyperthymic and irritable affective temperaments [F = 88.955, R2 = 0.710, p < 0.001] and lower MEQ scores [F = 75.456, R2 = 0.405, p < 0.001] statistically significantly predicted higher hsCRP. CONCLUSION: Eveningness chronotype and a depressive affective temperament appeared to be associated with higher hsCRP levels during moderate-to-severe unipolar and bipolar depression. Further longitudinal and larger studies should better characterise patients with mood disorders by investigating the influence of chronotype and temperament.
Assuntos
Transtorno Bipolar , Humanos , Transtorno Bipolar/psicologia , Temperamento , Proteína C-Reativa , Estudos Transversais , Estudos Retrospectivos , Cronotipo , Inquéritos e Questionários , Inventário de PersonalidadeRESUMO
Background: Dementia is one of the most common diseases in elderly people and hundreds of thousand new cases per year of Alzheimer's disease (AD) are estimated. While the recent decade has seen significant advances in the development of novel biomarkers to identify dementias at their early stage, a great effort has been recently made to identify biomarkers able to improve differential diagnosis. However, only few potential candidates, mainly detectable in cerebrospinal fluid (CSF), have been described so far. Methods: We searched for miRNAs regulating MAPT translation. We employed a capture technology able to find the miRNAs directly bound to the MAPT transcript in cell lines. Afterwards, we evaluated the levels of these miRNAs in plasma samples from FTD (n = 42) and AD patients (n = 33) and relative healthy controls (HCs) (n = 42) by using qRT-PCR. Results: Firstly, we found all miRNAs that interact with the MAPT transcript. Ten miRNAs have been selected to verify their effect on Tau levels increasing or reducing miRNA levels by using cell transfections with plasmids expressing the miRNAs genes or LNA antagomiRs. Following the results obtained, miR-92a-3p, miR-320a and miR-320b were selected to analyse their levels in plasma samples of patients with FTD and AD respect to HCs. The analysis showed that the miR-92a-1-3p was under-expressed in both AD and FTD compared to HCs. Moreover, miR-320a was upregulated in FTD vs. AD patients, particularly in men when we stratified by sex. Respect to HC, the only difference is showed in men with AD who have reduced levels of this miRNA. Instead, miR-320b is up-regulated in both dementias, but only patients with FTD maintain this trend in both genders. Conclusions: Our results seem to identify miR-92a-3p and miR-320a as possible good biomarkers to discriminate AD from HC, while miR-320b to discriminate FTD from HC, particularly in males. Combining three miRNAs improves the accuracy only in females, particularly for differential diagnosis (FTD vs. AD) and to distinguish FTD from HC.
RESUMO
In 1632, Galileo described his experience of motion aboard a large ship and exposed in detail the invariance principle, which was then rightly named after him. I suggest that more than three centuries earlier, in The Divine Comedy, his fellow countryman Dante Alighieri intuitively grasped what Galileo was later to establish as one of the pillars of modern science.
RESUMO
The inference of Shannon entropy out of sample histograms is known to be affected by systematic and random errors that depend on the finite size of the available data set. This dependence was mostly investigated in the multinomial case, in which states are visited in an independent fashion. In this paper the asymptotic behavior of the distribution of the sample Shannon entropy, also referred to as plug-in estimator, is investigated in the case of an underlying finite Markov process characterized by a regular stochastic matrix. As the size of the data set tends to infinity, the plug-in estimator is shown to become asymptotically normal, though in a way that substantially deviates from the known multinomial case. The asymptotic behavior of bias and variance of the plug-in estimator are expressed in terms of the spectrum of the stochastic matrix and of the related covariance matrix. Effects of initial conditions are also considered. By virtue of the formal similarity with Shannon entropy, the results are directly applicable to the evaluation of permutation entropy.
RESUMO
The statistical analysis of data stemming from dynamical systems, including, but not limited to, time series, routinely relies on the estimation of information theoretical quantities, most notably Shannon entropy. To this purpose, possibly the most widespread tool is provided by the so-called plug-in estimator, whose statistical properties in terms of bias and variance were investigated since the first decade after the publication of Shannon's seminal works. In the case of an underlying multinomial distribution, while the bias can be evaluated by knowing support and data set size, variance is far more elusive. The aim of the present work is to investigate, in the multinomial case, the statistical properties of an estimator of a parameter that describes the variance of the plug-in estimator of Shannon entropy. We then exactly determine the probability distributions that maximize that parameter. The results presented here allow one to set upper limits to the uncertainty of entropy assessments under the hypothesis of memoryless underlying stochastic processes.
RESUMO
The amount of information exchanged per unit of time between two dynamic processes is an important concept for the analysis of complex systems. Theoretical formulations and data-efficient estimators have been recently introduced for this quantity, known as the mutual information rate (MIR), allowing its continuous-time computation for event-based data sets measured as realizations of coupled point processes. This work presents the implementation of MIR for point process applications in Network Physiology and cardiovascular variability, which typically feature short and noisy experimental time series. We assess the bias of MIR estimated for uncoupled point processes in the frame of surrogate data, and we compensate it by introducing a corrected MIR (cMIR) measure designed to return zero values when the two processes do not exchange information. The method is first tested extensively in synthetic point processes including a physiologically-based model of the heartbeat dynamics and the blood pressure propagation times, where we show the ability of cMIR to compensate the negative bias of MIR and return statistically significant values even for weakly coupled processes. The method is then assessed in real point-process data measured from healthy subjects during different physiological conditions, showing that cMIR between heartbeat and pressure propagation times increases significantly during postural stress, though not during mental stress. These results document that cMIR reflects physiological mechanisms of cardiovascular variability related to the joint neural autonomic modulation of heart rate and arterial compliance.
RESUMO
The divergence rate method, which is used to determine the maximum Lyapunov exponent out of time series, is based on the evaluation of the time-dependent divergence exponent. For chaotic systems and in the small time regime, this exponent grows linearly in time. The asymptotic regime is instead characterized by a time-independent behavior due to the system eventually losing its memory of the starting conditions. The amplitude of this "plateau"-like divergence exponent depends both on the choice of the embedding dimension and lag and on the maximum distance of nearby starting trajectories in a way that is characteristic of the underlying dynamical system. In this paper, upon introducing the basic mathematical tools, we address the plateau evaluation for two classes of time series, those generated by a white noise source and those generated by a finite-dimensional chaotic system. The different behavior provides a novel tool to distinguish purely stochastic sources from deterministic ones, as well as to provide a precise estimate of the correlation dimension in the latter case. The method is also sensitive to correlated noise sources.
RESUMO
Physical connections between nodes in a complex network are constrained by limiting factors, such as the cost of establishing links and maintaining them, which can hinder network capability in terms of signal propagation speed and processing power. Trade-off mechanisms between cost constraints and performance requirements are reflected in the topology of a network and, ultimately, on the dependence of connectivity on geometric distance. This issue, though rarely addressed, is crucial in neuroscience, where physical links between brain regions are associated with a metabolic cost. In this work we investigate brain connectivity-estimated by means of a recently developed method that evaluates time scales of cross-correlation observability-and its dependence on geometric distance by analyzing resting state magnetoencephalographic recordings collected from a large set of healthy subjects. We identify three regimes of distance each showing a specific behavior of connectivity. This identification makes up a new tool to study the mechanisms underlying network formation and sustainment, with possible applications to the investigation of neuroscientific issues, such as aging and neurodegenerative diseases.
RESUMO
Lung cancer is still one of the leading cause of death worldwide. The clinical variability of lung cancer is high and drives treatment decision. In this context, correct discrimination of pulmonary neuroendocrine tumors is still of critical relevance. The spectrum of neuroendocrine tumors is various, and each type has molecular and phenotypical differences. In order to advance in the discrimination of neuroendocrine from non-neuroendocrine lung tumors, we tested a series of 95 surgically resected and formalin-fixed paraffin embedded lung cancer tissues, and we analyzed the expression of miR205-5p and miR375-3p via TaqMan RT-qPCR. Via a robust mathematical approach, we excluded technical outliers increasing the data reproducibility. We found that miR375-3p levels are higher in low-grade neuroendocrine lung tumor samples compared to non-neuroendocrine lung tumors. However, miR375-3p is not able to distinguish among different types of neuroendocrine lung tumors. In this work, we provide a new molecular marker for distinguishing non-neuroendocrine from low-grade neuroendocrine lung tumors samples establishing an easy miRNA score to be used in clinical settings, enabling the pathologist to classify more accurately lung tumors biopsies, which may be ambiguously cataloged in routine examination.
RESUMO
In any network, the dependence of connectivity on physical distance between nodes is a direct consequence of trade-off mechanisms between costs of establishing and sustaining links, processing rates, propagation speed of signals between nodes. Despite its universality, there are still few studies addressing this issue. Here we apply a recently-developed method to infer links between nodes, and possibly subnetwork structures, to determine connectivity strength as a function of physical distance between nodes. The model system we investigate is brain activity reconstructed on the cortex out of magnetoencephalography recordings sampled on a set of healthy subjects in resting state. We found that the dependence of the time scale of observability of a link on its geometric length follows a power-law characterized by an exponent whose extent is inversely proportional to connectivity. Our method provides a new tool to highlight and investigate networks in neuroscience.