A novel fMRI group data analysis method based on data-driven reference extracting from group subjects

Exploring brain effective connectivity of early MCI with GRU_GC model on resting-state fMRI

BACKGROUND

Investigating the functional interactions between different brain regions and revealing the transmission of information by computing brain connectivity have great potential and significance in the diagnosis of early Mild Cognitive Impairment (EMCI).

METHODS

The Granger causality with Gate Recurrent Unit (GRU_GC) model is a suitable method that allows the detection of a nonlinear causal relationship and solves the limitation of fixed time lag, which cannot be detected by the classical Granger method. The model can transmit time series signals with any transmission delay length, and the time series can be screened and learned through the gate model.

RESULTS

The classification experiment of 89 EMCI and 73 neurologically healthy controls (HC) shows that the accuracy reached 87.88%. Compared with multivariate variables GC (MVGC) and Long Short-Term Memory-based GC (LSTM_GC), the GRU_GC significantly improved the estimation of brain connectivity communication. Constructing a difference network to explore the brain effective connectivity between EMCI and HC.

CONCLUSIONS

The GRU_GC can discover the abnormal brain regions, including the parahippocampal gyrus, the posterior cingulate gyrus. The method can be used in clinical applications as an effective brain connectivity analysis tool and provides auxiliary means for the medical diagnosis of EMCI.

The effective connectivity analysis of fMRI based on asymmetric detection of transfer brain entropy

It is important to explore causal relationships in functional magnetic resonance imaging study. However, the traditional effective connectivity analysis method is easy to produce false causality, and the detection accuracy needs to be improved. In this paper, we introduce a novel functional magnetic resonance imaging effective connectivity method based on the asymmetry detection of transfer entropy, which quantifies the disparity in predictive information between forward and backward time, subsequently normalizing this disparity to establish a more precise criterion for detecting causal relationships while concurrently reducing computational complexity. Then, we evaluate the effectiveness of this method on the simulated data with different level of nonlinearity, and the results demonstrated that the proposed method outperforms others methods on the detection of both linear and nonlinear causal relationships, including Granger Causality, Partial Granger Causality, Kernel Granger Causality, Copula Granger Causality, and traditional transfer entropy. Furthermore, we applied it to study the effective connectivity of brain functional activities in seafarers. The results showed that there are significantly different causal relationships between different brain regions in seafarers compared with non-seafarers, such as Temporal lobe related to sound and auditory information processing, Hippocampus related to spatial navigation, Precuneus related to emotion processing as well as Supp_Motor_Area associated with motor control and coordination, which reflects the occupational specificity of brain function of seafarers.

Sparse Independence Component Analysis for Competitive Endogenous RNA Co-Module Identification in Liver Hepatocellular Carcinoma

OBJECTIVE

Long non-coding RNAs (lncRNAs) have been shown to be associated with the pathogenesis of different kinds of diseases and play important roles in various biological processes. Although numerous lncRNAs have been found, the functions of most lncRNAs and physiological/pathological significance are still in its infancy. Meanwhile, their expression patterns and regulation mechanisms are also far from being fully understood.

METHODS

In order to reveal functional lncRNAs and identify the key lncRNAs, we develop a new sparse independence component analysis (ICA) method to identify lncRNA-mRNA-miRNA expression co-modules based on the competitive endogenous RNA (ceRNA) theory using the sample-matched lncRNA, mRNA and miRNA expression profiles. The expression data of the three RNA combined together is approximated sparsely to obtain the corresponding sparsity coefficient, and then it is decomposed by using ICA constraint optimization to obtain the common basis and modules. Subsequently, affine propagation clustering is used to perform cluster analysis on the common basis under multiple running conditions to obtain the co-modules for the selection of different RNA elements.

RESULTS

We applied sparse ICA to Liver Hepatocellular Carcinoma (LIHC) dataset and the experiment results demonstrate that the proposed sparse ICA method can effectively discover biologically functional expression common modules.

CONCLUSION

It may provide insights into the function of lncRNAs and molecular mechanism of LIHC. Clinical and Translational Impact Statement-The results on LIHC dataset demonstrate that the proposed sparse ICA method can effectively discover biologically functional expression common modules, which may provide insights into the function of IncRNAs and molecular mechanism of LIHC.

The integrative functional connectivity analysis between seafarer’s brain networks using functional magnetic resonance imaging data of different states

The particularity of seafarers’ occupation makes their brain functional activities vulnerable to the influence of working environments, which leads to abnormal functional connectivities (FCs) between brain networks. To further investigate the influences of maritime environments on the seafarers’ functional brain networks, the functional magnetic resonance imaging (fMRI) datasets of 33 seafarers before and after sailing were used to study FCs among the functional brain networks in this paper. On the basis of making full use of the intrinsic prior information from fMRI data, six resting-state brain functional networks of seafarers before and after sailing were obtained by using group independent component analysis with intrinsic reference, and then the differences between the static and dynamic FCs among these six brain networks of seafarers before and after sailing were, respectively, analyzed from both group and individual levels. Subsequently, the potential dynamic functional connectivity states of seafarers before and after sailing were extracted by using the affine propagation clustering algorithm and the probabilities of state transition between them were obtained simultaneously. The results show that the dynamic FCs among large-scale brain networks have significant difference seafarers before and after sailing both at the group level and individual level, while the static FCs between them varies only at the individual level. This suggests that the maritime environments can indeed affect the brain functional activity of seafarers in real time, and the degree of influence is different for different subjects, which is of a great significance to explore the neural changes of seafarer’s brain functional network.

Brain Response of Major Depressive Disorder Patients to Emotionally Positive and Negative Music

Depression is characterized by poor emotion regulation that makes it difficult to escape the effects of emotional pain, but the neuromodulation behind these symptoms is still unclear. This study investigated the neural mechanism of emotional state-related responses during music stimuli in participants with major depressive disorder (MDD) compared to never-depressed (ND) controls. A novel two-level feature selection method, integrating recursive feature elimination based on support vector machine (SVM-RFE) and random forest algorithm (RF), was proposed to screen emotional recognition brain regions (ERBRs). On this basis, the differences of functional connectivity (FC) were systematically analyzed by two-sample t-test. The results demonstrate that ND participants show eight pairs of FCs with a significant difference between positive emotional music stimuli (pEMS) versus negative emotional music stimuli (nEMS) in 15 ERBRs of MDD, but the participants with MDD show one pair of significant difference in FC. The decreased number reflects the fuzzy response to positive and negative emotions in MDD, which appears to arise from obstacle to emotional cognition and regulation. Furthermore, there was no significant difference in FC between MDDs and NDs under pEMS, but a significant difference was detected between the two groups under nEMS (p < 0.01), revealing a 'bias' against the negative state in MDD. The current study may help to better comprehend the abnormal evolution from normal to depression and inform the utilization of pEMS in formal treatment for depression.

Functional Covariance Connectivity of Gray and White Matter in Olfactory-Related Brain Regions in Parkinson’s Disease

Before the onset of motor symptoms, Parkinson’s disease (PD) involves dysfunction of the anterior olfactory nucleus and olfactory bulb, causing olfactory disturbance, commonly resulting in hyposmia in the early stages of PD. Accumulating evidence has shown that blood oxygen level dependent (BOLD) signals in white matter are altered by olfactory disorders and related stimuli, and the signal changes in brain white matter pathways show a certain degree of specificity, which can reflect changes of early olfactory dysfunction in Parkinson’s disease. In this study, we apply the functional covariance connectivity (FCC) method to decode FCC of gray and white matter in olfactory-related brain regions in Parkinson’s disease. Our results show that the dorsolateral prefrontal, anterior entorhinal cortex and fronto-orbital cortices in the gray matter have abnormal connectivity with the posterior corona radiata and superior corona radiata in white matter in patients with Parkinson’s hyposmia. The functional covariance connection strength (FCS) of the right dorsolateral prefrontal cortex and white matter, and the covariance connection strength of the left superior corona radiata and gray matter function have potential diagnostic value. These results demonstrate that alterations in FCC of gray and white matter in olfactory-related brain regions can reflect the change of olfactory function in the early stages of Parkinson’s disease, indicating that it could be a potential neuroimaging marker for early diagnosis.

A Novel Constrained Non-negative Matrix Factorization Method for Group Functional Magnetic Resonance Imaging Data Analysis of Adult Attention-Deficit/Hyperactivity Disorder

Attention-deficit/hyperactivity disorder (ADHD) is a common childhood psychiatric disorder that often persists into adulthood. Extracting brain networks from functional magnetic resonance imaging (fMRI) data can help explore neurocognitive disorders in adult ADHD. However, there is still a lack of effective methods to extract large-scale brain networks to identify disease-related brain network changes. Hence, this study proposed a spatial constrained non-negative matrix factorization (SCNMF) method based on the fMRI real reference signal. First, non-negative matrix factorization analysis was carried out on each subject to select the brain network components of interest. Subsequently, the available spatial prior information was mined by integrating the interested components of all subjects. This prior constraint was then incorporated into the NMF objective function to improve its efficiency. For the sake of verifying the effectiveness and feasibility of the proposed method, we quantitatively compared the SCNMF method with other classical algorithms and applied it to the dynamic functional connectivity analysis framework. The algorithm successfully extracted ten resting-state brain functional networks from fMRI data of adult ADHD and healthy controls and found large-scale brain network changes in adult ADHD patients, such as enhanced connectivity between executive control network and right frontoparietal network. In addition, we found that older ADHD spent more time in the pattern of relatively weak connectivity. These findings indicate that the method can effectively extract large-scale functional networks and provide new insights into understanding the neurobiological mechanisms of adult ADHD from the perspective of brain networks.

Internet-like brain hierarchical network model: Alzheimer's disease study as an example

BACKGROUND AND OBJECTIVE

The modular structure and hierarchy are important topological characteristics in real complex networks such as brain networks on temporal scale. However, there are few studies investigating the hierarchical structure at the spatial scale of brain networks, the application of which still remains to be further studied.

METHODS

In this study, a novel model of brain hierarchical network based on the hierarchical characteristic of Internet topology is proposed for the first time, which is called Internet-like brain hierarchical network (IBHN). In this model, the whole brain network is partitioned into multiple levels: brain wide area network (Brain-WAN), brain metropolitan network (Brain-MAN), and brain local area network (Brain-LAN). A Brain-MAN is formed by the interconnection of multiple Brain-LANs, and the interconnection of multiple Brain-MANs forms a Brain-WAN. A multivariate analysis method is employed to measure overall functional connectivity between two brain networks at the same network level rather than detecting the change of each node pair's functional connection. Furthermore, we demonstrate the utility of IBHN model with application to a practical case-control study involving 64 patients with Alzheimer's disease and 75 healthy controls.

RESULTS

The proposed model identified enhanced functional connectivity (P-value<0.05) at Brain-WAN level and reduced functional connectivity (P-value=0.004) at Brain-LAN level of Alzheimer's disease patients, which can be used as a multi-dimension functional reference for AD's diagnosis.

CONCLUSIONS

This study not only provides insight into AD pathophysiology, but also further proves the effectiveness of the proposed IBHN model. In addition, the IBHN model makes it possible to explore the brain's functional organization from multiple dimensions and offers a multi-level perspective for the research of complex brain network.

The Brain Alteration of Seafarer Revealed by Activated Functional Connectivity Mode in fMRI Data Analysis

As a special occupational group, the working and living environments faced by seafarers are greatly different from those of land. It is easy to affect the psychological and physiological activities of seafarers, which inevitably lead to changes in the brain functional activities of seafarers. Therefore, it is of great significance to study the neural activity rules of seafarers' brain. In view of this, this paper studied the seafarers' brain alteration at the activated voxel level based on functional magnetic resonance imaging technology by comparing the differences in functional connectivities (FCs) between seafarers and non-seafarers. Firstly, the activated voxels of each group were obtained by independence component analysis, and then the distribution of these voxels in the brain and the common activated voxels between the two groups were statistically analyzed. Next, the FCs between the common activated voxels of the two groups were calculated and obtained the FCs that had significant differences between them through two-sample T-test. Finally, all FCs and FCs with significant differences (DFCs) between the common activated voxels were used as the features for the support vector machine to classify seafarers and non-seafarers. The results showed that DFCs between the activated voxels had better recognition ability for seafarers, especially for Precuneus_L and Precuneus_R, which may play an important role in the classification prediction of seafarers and non-seafarers, so that provided a new perspective for studying the specificity of neurological activities of seafarers.

Extraction and Analysis of Dynamic Functional Connectome Patterns in Migraine Sufferers: A Resting-State fMRI Study

Migraine seriously affects the physical and mental health of patients because of its recurrence and the hypersensitivity to the environment that it causes. However, the pathogenesis and pathophysiology of migraine are not fully understood. We addressed this issue in the present study using an autodynamic functional connectome model (A-DFCM) with twice-clustering to compare dynamic functional connectome patterns (DFCPs) from resting-state functional magnetic resonance imaging data from migraine patients and normal control subjects. We used automatic localization of segment points to improve the efficiency of the model, and intergroup differences and network metrics were analyzed to identify the neural mechanisms of migraine. Using the A-DFCM model, we identified 17 DFCPs-including 1 that was specific and 16 that were general-based on intergroup differences. The specific DFCP was closely associated with neuronal dysfunction in migraine, whereas the general DFCPs showed that the 2 groups had similar functional topology as well as differences in the brain resting state. An analysis of network metrics revealed the critical brain regions in the specific DFCP; these were not only distributed in brain areas related to pain such as Brodmann area 1/2/3, basal ganglia, and thalamus but also located in regions that have been implicated in migraine symptoms such as the occipital lobe. An analysis of the dissimilarities in general DFCPs between the 2 groups identified 6 brain areas belonging to the so-called pain matrix. Our findings provide insight into the neural mechanisms of migraine while also identifying neuroimaging biomarkers that can aid in the diagnosis or monitoring of migraine patients.

Dynamic visual cortical connectivity analysis based on functional magnetic resonance imaging

BACKGROUND

Studies of brain functional connectivity (FC) and effective connectivity (EC) using the functional magnetic resonance imaging (fMRI) have advanced our understanding of functional organization on visual cortex of human brain. The current studies mainly focus on static or dynamic connectivity, while the relationships between them have not been well characterized especially for static EC (sEC) and dynamic EC (dEC), as well as the consistency characteristics of changing trend of dFCs and dECs, which is of great importance to reveal the neural information processing mechanism in visual cortex region.

METHOD

In this study, we explore these relationships among several subareas of human visual cortex (V1-V5) by calculating the connection intensity and information flow among them over time by sliding window method, which are defined by Pearson correlation coefficient and Granger causality analysis, respectively, in each window.

RESULTS

The results demonstrate that there are extensive connections existing in human visual network, which are time-varying both in resting and task-related states. sFC intensity is negatively correlated with the variance of dFC, while sEC intensity is positively correlated with the variance of dEC. Furthermore, we also find that dFC within visual cortex at rest shows more consistency, while dEC shows less compared with task state in changing trend.

CONCLUSION

Therefore, this study provides novel findings about dynamics of connectivity in human visual cortex from the perspective of functional and effective connectivity.

The Identification of Alzheimer's Disease Using Functional Connectivity Between Activity Voxels in Resting-State fMRI Data

Background: Alzheimer’s disease (AD) is a common neurodegenerative disease occurring in the elderly population. The effective and accurate classification of AD symptoms by using functional magnetic resonance imaging (fMRI) has a great significance for the clinical diagnosis and prediction of AD patients. Methods: Therefore, this paper proposes a new method for identifying AD patients from healthy subjects by using functional connectivities (FCs) between the activity voxels in the brain based on fMRI data analysis. Firstly, independent component analysis is used to detect the activity voxels in the fMRI signals of AD patients and healthy subjects; Secondly, the FCs between the common activity voxels of the two groups are calculated, and then the FCs with significant differences are further identified by statistical analysis between them; Finally, the classification of AD patients from healthy subjects is realized by using FCs with significant differences as the feature samples in support vector machine. Results: The results show that the proposed identification method can obtain higher classification accuracy, and the FCs between activity voxels within prefrontal lobe as well as those between prefrontal and parietal lobes play an important role in the prediction of AD patients. Furthermore, we also find that more brain regions and much more voxels in some regions are activity in AD group compared with health control group. Conclusion: It has a great potential value for the AD pathogenesis mechanism study.

SCTICA: Sub-packet constrained temporal ICA method for fMRI data analysis

Independent component analysis (ICA) has become a widely used method for functional magnetic resonance imaging (fMRI) data analysis. However, spatial ICA usually performs better than temporal ICA with regard to the stability and accuracy of functional connectivity detection, and temporal ICA is often not feasible when it is applied to the analysis of real fMRI data of the whole brain because of the excessive spatial dimensions. In this paper, to overcome these problems, we propose a sub-packet constrained temporal ICA (SCTICA) method to take advantage of the a priori information using a multi-objective optimization framework with the Newton iterative algorithm. Moreover, a splitting strategy is presented to improve the feasibility of the temporal ICA for whole brain fMRI data analysis. The experimental results of real data show that the splitting strategy improved the ability of the temporal ICA to analyze whole brain fMRI data. Furthermore, the experimental results also demonstrated that the proposed SCTICA method can not only improve the stability of the temporal ICA, but can also improve the functional connectivity detection ability compared with the classical ICA and ICA with a priori information methods. In brief, the proposed SCTICA method overcomes the problem that prevents temporal ICA from being applied to fMRI data of the whole brain, and the functional connectivity detection performance is greatly improved compared with that of traditional methods.

A New Constrained Spatiotemporal ICA Method Based on Multi-Objective Optimization for fMRI Data Analysis

Compared with independent component analysis (ICA), constrained ICA (CICA) has unique advantages in functional magnetic resonance image (fMRI) data analysis by introducing some priori information into the estimation process. However, there are still some controversies in the current CICA methods, such as how to choose the threshold parameter to restrain the similarity, and how to reduce the accuracy requirements for a priori information. In this paper, we propose a new constrained spatiotemporal ICA (CSTICA) method based on the framework of multi-objective optimization, where the inequality constraint of the traditional CICA method is transformed into the objective optimization function of the CSTICA, and both temporal and spatial a priori information are included simultaneously. The simulated, hybrid, and real fMRI data experiments are designed to evaluate the performance of the proposed CSTICA method in comparison with the classical ICA and CICA methods. Compared with the traditional CICA methods, the CSTICA has circumvented the problem of threshold parameter selection. Furthermore, the experimental results demonstrate that the source recovery ability of the CSTICA has been improved to a certain extent especially in the cases of a priori information with low accuracies. Meanwhile, the results also indicate that the CSTICA reduces dependency on the accuracy of a priori information.

A multi-variate blind source separation algorithm

BACKGROUND AND OBJECTIVE

The study follows the proposal of decomposing a given data matrix into a product of independent spatial and temporal component matrices. A multi-variate decomposition approach is presented, based on an approximate diagonalization of a set of matrices computed using a latent space representation.

METHODS

The proposed methodology follows an algebraic approach, which is common to space, temporal or spatiotemporal blind source separation algorithms. More specifically, the algebraic approach relies on singular value decomposition techniques, which avoids computationally costly and numerically instable matrix inversion. The method is equally applicable to correlation matrices determined from second order correlations or by considering fourth order correlations.

RESULTS

The resulting algorithms are applied to fMRI data sets either to extract the underlying fMRI components or to extract connectivity maps from resting state fMRI data collected for a dynamic functional connectivity analysis. Intriguingly, our algorithm shows increased spatial specificity compared to common approaches, while temporal precision stays similar.

CONCLUSION

The study presents a novel spatiotemporal blind source separation algorithm, which is both robust and avoids parameters that are difficult to fine tune. Applied on experimental data sets, the new method yields highly confined and focused areas with least spatial extent in the retinotopy case, and similar results in the dynamic functional connectivity analyses compared to other blind source separation algorithms. Therefore, we conclude that our novel algorithm is highly competitive and yields results, which are superior or at least similar to existing approaches.

Brain Functional Plasticity Driven by Career Experience: A Resting-State fMRI Study of the Seafarer

The functional connectome derived from BOLD resting-state functional magnetic resonance imaging data represents meaningful functional organizations and a shift between distinct cognitive states. However, the body of knowledge on how the long-term career experience affects the brain's functional plasticity is still very limited. In this study, we used a dynamic functional connectome characterization (DBFCC) model with the automatic target generation process K-Means clustering to explore the functional reorganization property of resting brain states, driven by long-term career experience. Taking sailors as an example, DBFCC generated seventeen reproducibly common atomic connectome patterns (ACP) and one reproducibly distinct ACP, i.e., ACP14. The common ACPs indicating the same functional topology of the resting brain state transitions were shared by two control groups, while the distinct ACP, which mainly represented functional plasticity and only existed in the sailors, showed close relationships with the long-term career experience of sailors. More specifically, the distinct ACP14 of the sailors was made up of four specific sub-networks, such as the auditory network, visual network, executive control network, and vestibular function-related network, which were most likely linked to sailing experience, i.e., continuously suffering auditory noise, maintaining balance, locating one's position in three-dimensional space at sea, obeying orders, etc. Our results demonstrated DBFCC's effectiveness in revealing the specifically functional alterations modulated by sailing experience and particularly provided the evidence that functional plasticity was beneficial in reorganizing brain's functional topology, which could be driven by career experience.

A Novel Feature-Map Based ICA Model for Identifying the Individual, Intra/Inter-Group Brain Networks across Multiple fMRI Datasets

Independent component analysis (ICA) has been widely used in functional magnetic resonance imaging (fMRI) data analysis to evaluate functional connectivity of the brain; however, there are still some limitations on ICA simultaneously handling neuroimaging datasets with diverse acquisition parameters, e.g., different repetition time, different scanner, etc. Therefore, it is difficult for the traditional ICA framework to effectively handle ever-increasingly big neuroimaging datasets. In this research, a novel feature-map based ICA framework (FMICA) was proposed to address the aforementioned deficiencies, which aimed at exploring brain functional networks (BFNs) at different scales, e.g., the first level (individual subject level), second level (intragroup level of subjects within a certain dataset) and third level (intergroup level of subjects across different datasets), based only on the feature maps extracted from the fMRI datasets. The FMICA was presented as a hierarchical framework, which effectively made ICA and constrained ICA as a whole to identify the BFNs from the feature maps. The simulated and real experimental results demonstrated that FMICA had the excellent ability to identify the intergroup BFNs and to characterize subject-specific and group-specific difference of BFNs from the independent component feature maps, which sharply reduced the size of fMRI datasets. Compared with traditional ICAs, FMICA as a more generalized framework could efficiently and simultaneously identify the variant BFNs at the subject-specific, intragroup, intragroup-specific and intergroup levels, implying that FMICA was able to handle big neuroimaging datasets in neuroscience research.

An improved multi-objective optimization-based CICA method with data-driver temporal reference for group fMRI data analysis

Group independent component analysis (GICA) has been successfully applied to study multi-subject functional magnetic resonance imaging (fMRI) data, and the group independent component (GIC) represents the commonality of all subjects in the group. However, some studies show that the performance of GICA can be improved by incorporating a priori information, which is not always considered when looking for GICs in existing GICA methods. In this paper, we propose an improved multi-objective optimization-based constrained independent component analysis (CICA) method to take advantage of the temporal a priori information extracted from all subjects in the group by incorporating it into the computational process of GICA for group fMRI data analysis. The experimental results of simulated and real data show that the activated regions and the time course detected by the improved CICA method are more accurate in some sense. Moreover, the GIC computed by the improved CICA method has a higher correlation with the corresponding independent component of each subject in the group, which means that the improved CICA method with the temporal a priori information extracted from the group can better reflect the commonality of the subjects. These results demonstrate that the improved CICA method has its own advantages in fMRI data analysis.

SCGICAR: Spatial concatenation based group ICA with reference for fMRI data analysis

BACKGROUND AND OBJECTIVE

With the rapid development of big data, the functional magnetic resonance imaging (fMRI) data analysis of multi-subject is becoming more and more important. As a kind of blind source separation technique, group independent component analysis (GICA) has been widely applied for the multi-subject fMRI data analysis. However, spatial concatenated GICA is rarely used compared with temporal concatenated GICA due to its disadvantages.

METHODS

In this paper, in order to overcome these issues and to consider that the ability of GICA for fMRI data analysis can be improved by adding a priori information, we propose a novel spatial concatenation based GICA with reference (SCGICAR) method to take advantage of the priori information extracted from the group subjects, and then the multi-objective optimization strategy is used to implement this method. Finally, the post-processing means of principal component analysis and anti-reconstruction are used to obtain group spatial component and individual temporal component in the group, respectively.

RESULTS

The experimental results show that the proposed SCGICAR method has a better performance on both single-subject and multi-subject fMRI data analysis compared with classical methods. It not only can detect more accurate spatial and temporal component for each subject of the group, but also can obtain a better group component on both temporal and spatial domains.

CONCLUSIONS

These results demonstrate that the proposed SCGICAR method has its own advantages in comparison with classical methods, and it can better reflect the commonness of subjects in the group.