Testing of two-dimensional Gaussian processes by sample cross-covariance function.

In this paper, we address the issue of testing two-dimensional Gaussian processes with a defined cross-dependency structure. Multivariate Gaussian processes are widely used in various applications; therefore, it is essential to identify the theoretical model that accurately describes the data. While it is relatively straightforward to do so in a one-dimensional case, analyzing multi-dimensional vectors requires considering the dependency between the components, which can significantly affect the efficiency of statistical methods. The testing methodology presented in this paper is based on the sample cross-covariance function and can be considered a natural generalization of the approach recently proposed for testing one-dimensional Gaussian processes based on the sample autocovariance function. We verify the efficiency of this procedure on three classes of two-dimensional Gaussian processes: Brownian motion, fractional Brownian motion, and two-dimensional autoregressive discrete-time process. The simulation results clearly demonstrate the effectiveness of the testing methodology, even for small sample sizes. The theoretical and simulation results are supported by analyzing two-dimensional real-time series that describe the main risk factors of a mining company, namely, copper price and exchange rates (USDPLN). We believe that the introduced methodology is intuitive and relatively simple to implement, and thus, it can be applied in many real-world scenarios where multi-dimensional data are examined.

Goodness-of-fit test for stochastic processes using even empirical moments statistic.

In this paper, we introduce a novel framework that allows efficient stochastic process discrimination. The underlying test statistic is based on even empirical moments and generalizes the time-averaged mean-squared displacement framework; the test is designed to allow goodness-of-fit statistical testing of processes with stationary increments and a finite-moment distribution. In particular, while our test statistic is based on a simple and intuitive idea, it enables efficient discrimination between finite- and infinite-moment processes even if the underlying laws are relatively close to each other. This claim is illustrated via an extensive simulation study, e.g., where we confront α-stable processes with stability index close to 2 with their standard Gaussian equivalents. For completeness, we also show how to embed our methodology into the real data analysis by studying the real metal price data.

Morphological and molecular characteristics of Geocenamus longus and the first report of G. brevidens from a karst cave (Nematoda: Merliniidae Siddiqi, 1971).

The taxonomic status of stunt nematodes remains under investigation. Many nematode species belonging to the subfamily Merliniinae still await for more detailed studies, especially for the analyzes on the molecular level. In the presented work, two nematode species belonging to this subfamily were investigated. Characteristics of Geocenamus longus were expanded on the morphological and the rDNA- and mtCOI-derived molecular data. The mitochondrial sequences are for the first time presented for Merliniinae. Additionally, a bisexual population of G. brevidens was reported and characterized from a new, unique underground environment.

Time-averaged mean squared displacement ratio test for Gaussian processes with unknown diffusion coefficient.

The time-averaged mean squared displacement (TAMSD) is one of the most common statistics used for the analysis of anomalous diffusion processes. Anomalous diffusion is manifested by non-linear (mostly power-law) characteristics of the process in contrast to normal diffusion where linear characteristics are expected. One can distinguish between sub- and super-diffusive processes. We consider Gaussian anomalous diffusion models and propose a new approach used for their testing. This approach is based on the TAMSD ratio statistic for different time lags. Similar to the TAMSD, this statistic exhibits a specific behavior in the anomalous diffusion regime. Through its structure, it is independent of the diffusion coefficient, which, in general, does not influence anomalous diffusion behavior. Thus, the TAMSD ratio-based approach does not require preliminary knowledge of the diffusion coefficient's value, in contrast to the TAMSD-approach, where this value is crucial in the testing procedure. Based on the quadratic form representation of the TAMSD ratio, we calculate its main characteristics and propose a step-by-step testing procedure that can be applied for any Gaussian process. For the anomalous diffusion model used here, namely, the fractional Brownian motion, we demonstrate the effectiveness of the proposed methodology. We show that the new approach outperforms the TAMSD-based one, especially for small sample sizes. Finally, the methodology is applied to the real data from the financial market.

Discriminating Gaussian processes via quadratic form statistics.

Gaussian processes are powerful tools for modeling and predicting various numerical data. Hence, checking their quality of fit becomes a vital issue. In this article, we introduce a testing methodology for general Gaussian processes based on a quadratic form statistic. We illustrate the methodology on three statistical tests recently introduced in the literature, which are based on the sample autocovariance function, time average mean-squared displacement, and detrended moving average statistics. We compare the usefulness of the tests by taking into consideration three very important Gaussian processes: the fractional Brownian motion, which is self-similar with stationary increments (SSSIs), scaled Brownian motion, which is self-similar with independent increments (SSIIs), and the Ornstein-Uhlenbeck (OU) process, which is stationary. We show that the considered statistics' ability to distinguish between these Gaussian processes is high, and we identify the best performing tests for different scenarios. We also find that there is no omnibus quadratic form test; however, the detrended moving average test seems to be the first choice in distinguishing between same processes with different parameters. We also show that the detrended moving average method outperforms the Cholesky method. Based on the previous findings, we introduce a novel procedure of discriminating between Gaussian SSSI, SSII, and stationary processes. Finally, we illustrate the proposed procedure by applying it to real-world data, namely, the daily EURUSD currency exchange rates, and show that the data can be modeled by the OU process.

Fractional Dynamics Identification via Intelligent Unpacking of the Sample Autocovariance Function by Neural Networks.

Many single-particle tracking data related to the motion in crowded environments exhibit anomalous diffusion behavior. This phenomenon can be described by different theoretical models. In this paper, fractional Brownian motion (FBM) was examined as the exemplary Gaussian process with fractional dynamics. The autocovariance function (ACVF) is a function that determines completely the Gaussian process. In the case of experimental data with anomalous dynamics, the main problem is first to recognize the type of anomaly and then to reconstruct properly the physical rules governing such a phenomenon. The challenge is to identify the process from short trajectory inputs. Various approaches to address this problem can be found in the literature, e.g., theoretical properties of the sample ACVF for a given process. This method is effective; however, it does not utilize all of the information contained in the sample ACVF for a given trajectory, i.e., only values of statistics for selected lags are used for identification. An evolution of this approach is proposed in this paper, where the process is determined based on the knowledge extracted from the ACVF. The designed method is intuitive and it uses information directly available in a new fashion. Moreover, the knowledge retrieval from the sample ACVF vector is enhanced with a learning-based scheme operating on the most informative subset of available lags, which is proven to be an effective encoder of the properties inherited in complex data. Finally, the robustness of the proposed algorithm for FBM is demonstrated with the use of Monte Carlo simulations.

Probabilistic properties of detrended fluctuation analysis for Gaussian processes.

Detrended fluctuation analysis (DFA) is one of the most widely used tools for the detection of long-range dependence in time series. Although DFA has found many interesting applications and has been shown to be one of the best performing detrending methods, its probabilistic foundations are still unclear. In this paper, we study probabilistic properties of DFA for Gaussian processes. Our main attention is paid to the distribution of the squared error sum of the detrended process. We use a probabilistic approach to derive general formulas for the expected value and the variance of the squared fluctuation function of DFA for Gaussian processes. We also get analytical results for the expected value of the squared fluctuation function for particular examples of Gaussian processes, such as Gaussian white noise, fractional Gaussian noise, ordinary Brownian motion, and fractional Brownian motion. Our analytical formulas are supported by numerical simulations. The results obtained can serve as a starting point for analyzing the statistical properties of DFA-based estimators for the fluctuation function and long-memory parameter.

Identifying diffusive motions in single-particle trajectories on the plasma membrane via fractional time-series models.

In this paper we show that an autoregressive fractionally integrated moving average time-series model can identify two types of motion of membrane proteins on the surface of mammalian cells. Specifically we analyze the motion of the voltage-gated sodium channel Nav1.6 and beta-2 adrenergic receptors. We find that the autoregressive (AR) part models well the confined dynamics whereas the fractionally integrated moving average (FIMA) model describes the nonconfined periods of the trajectories. Since the Ornstein-Uhlenbeck process is a continuous counterpart of the AR model, we are also able to calculate its physical parameters and show their biological relevance. The fitted FIMA and AR parameters show marked differences in the dynamics of the two studied molecules.

Statistical properties of the anomalous scaling exponent estimator based on time-averaged mean-square displacement.

The most common way of estimating the anomalous scaling exponent from single-particle trajectories consists of a linear fit of the dependence of the time-averaged mean-square displacement on the lag time at the log-log scale. We investigate the statistical properties of this estimator in the case of fractional Brownian motion (FBM). We determine the mean value, the variance, and the distribution of the estimator. Our theoretical results are confirmed by Monte Carlo simulations. In the limit of long trajectories, the estimator is shown to be asymptotically unbiased, consistent, and with vanishing variance. These properties ensure an accurate estimation of the scaling exponent even from a single (long enough) trajectory. As a consequence, we prove that the usual way to estimate the diffusion exponent of FBM is correct from the statistical point of view. Moreover, the knowledge of the estimator distribution is the first step toward new statistical tests of FBM and toward a more reliable interpretation of the experimental histograms of scaling exponents in microbiology.

Mean-squared-displacement statistical test for fractional Brownian motion.

Anomalous diffusion in crowded fluids, e.g., in cytoplasm of living cells, is a frequent phenomenon. A common tool by which the anomalous diffusion of a single particle can be classified is the time-averaged mean square displacement (TAMSD). A classical mechanism leading to the anomalous diffusion is the fractional Brownian motion (FBM). A validation of such process for single-particle tracking data is of great interest for experimentalists. In this paper we propose a rigorous statistical test for FBM based on TAMSD. To this end we analyze the distribution of the TAMSD statistic, which is given by the generalized chi-squared distribution. Next, we study the power of the test by means of Monte Carlo simulations. We show that the test is very sensitive for changes of the Hurst parameter. Moreover, it can easily distinguish between two models of subdiffusion: FBM and continuous-time random walk.

Elucidating distinct ion channel populations on the surface of hippocampal neurons via single-particle tracking recurrence analysis.

Protein and lipid nanodomains are prevalent on the surface of mammalian cells. In particular, it has been recently recognized that ion channels assemble into surface nanoclusters in the soma of cultured neurons. However, the interactions of these molecules with surface nanodomains display a considerable degree of heterogeneity. Here, we investigate this heterogeneity and develop statistical tools based on the recurrence of individual trajectories to identify subpopulations within ion channels in the neuronal surface. We specifically study the dynamics of the K^{+} channel Kv1.4 and the Na^{+} channel Nav1.6 on the surface of cultured hippocampal neurons at the single-molecule level. We find that both these molecules are expressed in two different forms with distinct kinetics with regards to surface interactions, emphasizing the complex proteomic landscape of the neuronal surface. Further, the tools presented in this work provide new methods for the analysis of membrane nanodomains, transient confinement, and identification of populations within single-particle trajectories.

Identifying ergodicity breaking for fractional anomalous diffusion: Criteria for minimal trajectory length.

In this paper, we study ergodic properties of α-stable autoregressive fractionally integrated moving average (ARFIMA) processes which form a large class of anomalous diffusions. A crucial practical question is how long trajectories one needs to observe in an experiment in order to claim that the analyzed data are ergodic or not. This will be solved by checking the asymptotic convergence to 0 of the empirical estimator F(n) for the dynamical functional D(n) defined as a Fourier transform of the n-lag increments of the ARFIMA process. Moreover, we introduce more flexible concept of the ε-ergodicity.

Estimating the anomalous diffusion exponent for single particle tracking data with measurement errors - An alternative approach.

Accurately characterizing the anomalous diffusion of a tracer particle has become a central issue in biophysics. However, measurement errors raise difficulty in the characterization of single trajectories, which is usually performed through the time-averaged mean square displacement (TAMSD). In this paper, we study a fractionally integrated moving average (FIMA) process as an appropriate model for anomalous diffusion data with measurement errors. We compare FIMA and traditional TAMSD estimators for the anomalous diffusion exponent. The ability of the FIMA framework to characterize dynamics in a wide range of anomalous exponents and noise levels through the simulation of a toy model (fractional Brownian motion disturbed by Gaussian white noise) is discussed. Comparison to the TAMSD technique, shows that FIMA estimation is superior in many scenarios. This is expected to enable new measurement regimes for single particle tracking (SPT) experiments even in the presence of high measurement errors.

Guidelines for the fitting of anomalous diffusion mean square displacement graphs from single particle tracking experiments.

Single particle tracking is an essential tool in the study of complex systems and biophysics and it is commonly analyzed by the time-averaged mean square displacement (MSD) of the diffusive trajectories. However, past work has shown that MSDs are susceptible to significant errors and biases, preventing the comparison and assessment of experimental studies. Here, we attempt to extract practical guidelines for the estimation of anomalous time averaged MSDs through the simulation of multiple scenarios with fractional Brownian motion as a representative of a large class of fractional ergodic processes. We extract the precision and accuracy of the fitted MSD for various anomalous exponents and measurement errors with respect to measurement length and maximum time lags. Based on the calculated precision maps, we present guidelines to improve accuracy in single particle studies. Importantly, we find that in some experimental conditions, the time averaged MSD should not be used as an estimator.

Fractional process as a unified model for subdiffusive dynamics in experimental data.

We show how to use a fractional autoregressive integrated moving average (FARIMA) model to a statistical analysis of the subdiffusive dynamics. The discrete time FARIMA(1,d,1) model is applied in this paper to the random motion of an individual fluorescently labeled mRNA molecule inside live E. coli cells in the experiment described in detail by Golding and Cox [Phys. Rev. Lett. 96, 098102 (2006)] as well as to the motion of fluorescently labeled telomeres in the nucleus of live human cells (U2OS cancer) in the experiment performed by Bronstein et al. [Phys. Rev. Lett. 103, 018102 (2009)]. It is found that only the memory parameter d of the FARIMA model completely detects an anomalous dynamics of the experimental data in both cases independently of the observed distribution of random noises.

Diffusive and subdiffusive dynamics of indoor microclimate: a time series modeling.

The indoor microclimate is an issue in modern society, where people spend about 90% of their time indoors. Temperature and relative humidity are commonly used for its evaluation. In this context, the two parameters are usually considered as behaving in the same manner, just inversely correlated. This opinion comes from observation of the deterministic components of temperature and humidity time series. We focus on the dynamics and the dependency structure of the time series of these parameters, without deterministic components. Here we apply the mean square displacement, the autoregressive integrated moving average (ARIMA), and the methodology for studying anomalous diffusion. The analyzed data originated from five monitoring locations inside a modern office building, covering a period of nearly one week. It was found that the temperature data exhibited a transition between diffusive and subdiffusive behavior, when the building occupancy pattern changed from the weekday to the weekend pattern. At the same time the relative humidity consistently showed diffusive character. Also the structures of the dependencies of the temperature and humidity data sets were different, as shown by the different structures of the ARIMA models which were found appropriate. In the space domain, the dynamics and dependency structure of the particular parameter were preserved. This work proposes an approach to describe the very complex conditions of indoor air and it contributes to the improvement of the representative character of microclimate monitoring.