MOTIVATION: The quantification of long-range correlation of electroencephalogram (EEG) signals is an important research direction for its relevance in helping understanding the brain activity. Epileptic seizures have been studied in the past years where different non-linear statistical approaches have been employed to understand the relationship between the EEG signal and the epileptic discharge. One of the most widely used method for to analyse long memory processes is the detrended fluctuation analysis (DFA). However, no objective and pragmatic methods have been developed to detect crossover points and reference channels in DFA. RESULTS: In this article, we propose: (i) two automatic approaches that successfully detect crossover points in DFA related methods on the log-log plot and (ii) a criteria to choose the reference channel for the log-amplitude function. Moreover, the DFA is applied to EEG signals of 10 epileptic patients collected from the CHB-MIT database, being the log-amplitude function used to compare the different brain hemispheres by making use of the methodology proposed in the article. The existence of long-range power-law correlations is demonstrated and indicates that the EEG signals of epileptic patients present three well-defined regions with the first region showing a 1/f noise (pink noise) for seven subjects and a random walk behaviour for three subjects. The second and third regions show anti-persistence behaviour. Moreover, the results of the log-amplitude function were divided in two groups: the first, including seven subjects, where the increase in the scales results in an increase in the fluctuation in the frontal channels and the second, included three subjects, where the fluctuation for large scales are greater for the parietal channels. AVAILABILITY AND IMPLEMENTATION: The functions used in this article are available in the R package DFA (Mesquita et al., 2020). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

MOTIVATION: In genetic association studies, linear mixed models (LMMs) are used to test for associations between phenotypes and candidate single nucleotide polymorphisms (SNPs). These same models are also used to estimate heritability, which is central not only to evolutionary biology but also to the prediction of the response to selection in plant and animal breeding, as well as the prediction of disease risk in humans. However, when one or more of the underlying assumptions are violated, the estimation of variance components may be compromised and therefore so may the estimates of heritability and any other functions of these. Considering that datasets obtained from real life experiments are prone to several sources of contamination, which usually induce the violation of the assumption of the normality of the errors, a robust derivative-free restricted-maximum likelihood framework (DF-REML) together with a robust coefficient of determination are proposed for the LMM in the context of genetic studies of continuous traits. RESULTS: The proposed approach, in addition to the robust estimation of variance components and robust computation of the coefficient of determination, allows in particular for the robust estimation of SNP-based heritability by reducing the bias and increasing the precision of its estimates. The performance of both classical and robust DF-REML approaches is compared via a Monte Carlo simulation study. Additionally, three examples of application of the methodologies to real datasets are given in order to validate the usefulness of the proposed robust approach. Although the main focus of this article is on plant breeding applications, the proposed methodology is applicable to both human and animal genetic studies. AVAILABILITY AND IMPLEMENTATION: Source code implemented in R is available in the Supplementary Material. CONTACT: [email protected]. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

MOTIVATION: One of the most widely used models to analyse genotype-by-environment data is the additive main effects and multiplicative interaction (AMMI) model. Genotype-by-environment data resulting from multi-location trials are usually organized in two-way tables with genotypes in the rows and environments (location-year combinations) in the columns. The AMMI model applies singular value decomposition (SVD) to the residuals of a specific linear model, to decompose the genotype-by-environment interaction (GEI) into a sum of multiplicative terms. However, SVD, being a least squares method, is highly sensitive to contamination and the presence of even a single outlier, if extreme, may draw the leading principal component towards itself resulting in possible misinterpretations and in turn lead to bad practical decisions. Since, as in many other real-life studies the distribution of these data is usually not normal due to the presence of outlying observations, either resulting from measurement errors or sometimes from individual intrinsic characteristics, robust SVD methods have been suggested to help overcome this handicap. RESULTS: We propose a robust generalization of the AMMI model (the R-AMMI model) that overcomes the fragility of its classical version when the data are contaminated. Here, robust statistical methods replace the classic ones to model, structure and analyse GEI. The performance of the robust extensions of the AMMI model is assessed through a Monte Carlo simulation study where several contamination schemes are considered. Applications to two real plant datasets are also presented to illustrate the benefits of the proposed methodology, which can be broadened to both animal and human genetics studies. AVAILABILITY AND IMPLEMENTATION: Source code implemented in R is available in the supplementary material under the function r-AMMI. CONTACT: [email protected] SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.