Next Article in Journal
Physical Therapists’ Opinion of E-Health Treatment of Chronic Low Back Pain
Previous Article in Journal
Dietary Quality in Vegetarian and Omnivorous Female Students in Germany: A Retrospective Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Application of a Hybrid Intelligent System for Diagnosing Primary Headaches

by
Svetlana Simić
1,
José R. Villar
2,
José Luis Calvo-Rolle
3,
Slobodan R. Sekulić
1,
Svetislav D. Simić
4 and
Dragan Simić
4,*
1
Faculty of Medicine, University of Novi Sad, 21000 Novi Sad, Serbia
2
Faculty of Geology, Campus de Llamaquique, University of Oviedo, 33005 Oviedo, Spain
3
Department of Industrial Engineering, University of A Coruña, 15405 Ferrol-A Coruña, Spain
4
Faculty of Technical Sciences, University of Novi Sad, 21000 Novi Sad, Serbia
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2021, 18(4), 1890; https://doi.org/10.3390/ijerph18041890
Submission received: 14 December 2020 / Revised: 10 February 2021 / Accepted: 11 February 2021 / Published: 16 February 2021

Abstract

:
(1) Background: Modern medicine generates a great deal of information that stored in medical databases. Simultaneously, extracting useful knowledge and making scientific decisions for diagnosis and treatment of diseases becomes increasingly necessary. Headache disorders are the most prevalent of all the neurological conditions. Headaches have not only medical but also great socioeconomic significance. The aim of this research is to develop an intelligent system for diagnosing primary headache disorders. (2) Methods: This research applied various mathematical, statistical and artificial intelligence techniques, among which the most important are: Calinski-Harabasz index, Analytical Hierarchy Process, and Weighted Fuzzy C-means Clustering Algorithm. These methods, techniques and methodologies are used to create a hybrid intelligent system for diagnosing primary headache disorders. The proposed intelligent diagnostic system is tested with original real-world data set with different metrics. (3) Results: First at all, nine of 20 attributes – features from International Headache Society (IHS) criteria are selected, and then only five most important attributes from IHS criteria are selected. The calculation result based on the Calinski–Harabasz index value (178) for the optimal number of clusters is three, and they present three classes of headaches: (i) migraine, (ii) tension-type headaches (TTHs), and (iii) other primary headaches (OPHs). The proposed hybrid intelligent system shows the following quality metrics: Accuracy 75%; Precision 67% for migraine, 74% for TTHs, 86% for OPHs, and Average Precision 77%; Recall 86% for migraine, 73% for TTHs, 67% for OPHs, Average Recall 75%; F1 score 75% for migraine, 74% for TTHs, 75% for OPHs, and Average F1 score 75%. (4) Conclusions: The hybrid intelligent system presents qualitative and respectable experimental results. The implementation of existing diagnostics systems and the development of new diagnostics systems in medicine is necessary in order to help physicians make quality diagnosis and decide the best treatments for the patients.

1. Introduction

A great deal of information in the field of modern medicine is stored in medical databases. These databases have also become necessary for the extraction of useful knowledge and scientific decision making for the diagnosis and treatment of diseases. The medical field is concerned with patient care activity primarily and research only secondarily. Thus, the justification for the collection of medical data is established on the benefits for individual patients.
Headache disorders can be considered as the predominant neurological condition and one of the most frequent medical complaints in general practice. More than 90% of the general population report suffering from a headache during any given year, which can be regarded as a lifetime history of head pain [1], while over 10% have had at least one migraine headache [2]. Headache is not a disease which typically shortens one’s life, but it can be a serious social as well as a health problem. A lot of people with headache are unaware of the type of headache they are suffering from.
Headache is one of the most frequent complaints in medicine in general and a frequent disorder in the working population that significantly affects absenteeism and loss of productivity [3,4,5]. The impact of headache disorders is a problem of enormous proportions, both for individuals and the society. Different diseases, including headaches, which may develop into disabilities, were subjects to investigations mostly in the developed countries, regarding the effects of illness on absenteeism and the reduced productivity [6,7]. The cost–benefit and risk–benefit ratios are two of the most relevant items in ongoing health-organization procedures [8]. In the field of neurological diseases, migraine was estimated to cost a total of €27 billion per year for the loss through reduced work productivity in the European Community. The second most costly neurological disease was stroke, totaling €22 billion, followed by epilepsy, Parkinson’s disease and multiple sclerosis with the costs of €16, €11 and €9 billion, respectively [9]. Migraine is highly disabling, associated with substantial economic burden. Approximately 13 billion dollars per year are lost through reduced work productivity in the United States [10].
Therefore, physicians frequently find themselves in situations where they have to consider the economic aspects of a treatment, even though they have not been trained for it. Healthcare requires making decisions about diagnosis and treatment that are compatible with current medical opinions and practices along with taking care of its availability and patient compliance to treatment [11,12]. On the other side, people do not always visit physicians for their complaints, and when they do, physicians do not always make a correct diagnosis. The education of both health professionals and patients is also significant. Better education means better treatment results and less economic loss [13]. The epidemiological research [14] has heavily depended on the diagnostic criteria introduced by the International Headache Society (IHS), and some automatic methods, expert systems, and knowledge-based systems, including some tools that help physicians make diagnoses, have been developed and based on these criteria.
Therefore, according to the previous facts, this research is focused on computerized-based application of a hybrid intelligent system for diagnosing primary headaches. In this research, the clustering method in general and the fuzzy c-means clustering approach in particular are used. The main reason is the fact that in the real-world setting, the physicians do not know the type of primary headache in advance. Therefore, it is necessary to use some artificial intelligent technique which satisfies this objective. Likewise, some other mathematics, statistics, and decision-making techniques are applied to create an intelligent system for diagnosing primary headache disorders. The proposed intelligent system is tested on the real-world data set collected from the Clinical Centre of Vojvodina, Serbia. This paper presents a research that extends from and is applied on the basis of our previous research [15] and on the research in computer-assisted diagnosis methods [16,17,18], as well as on the industrial applications developed for clustering methods [19,20].

2. Materials and Methods

2.1. Primary Headache Classification

The uniform terminology, as well as the consistent operational diagnostic criteria applied to a wide range of the headache disorders around the world [21] are established on The International Classification of Headache Disorders–The Third Edition (ICHD-3). The ICHD-3 offers a hierarchy of diagnoses specifying a number of varying degrees. Headache disorders are usually identified using three-digit codes, and sometimes using five-digit codes, which is displayed in detail in Table 1 with the short identification for only two important digit codes.
All headache disorders are classified into two major groups: (A) Primary headaches, from ICHD-3 code 1. to 4., and (B) Secondary headaches, ICHD-3 code from 5. to 12. The first digit specifies the major diagnostic categories (i.e. migraine). The second digit indicates a disorder within the category (i.e. migraine without aura). Each category is then subdivided into groups, types, subtypes and sub-forms. Subsequent digits permit more specific diagnosis for some headache types [21].

2.2. Data Extraction

Our previous study was completed on a sample of 1022 employed subjects, who filled in the questionnaire in their workplace. The statistical analyses of that research are presented in [15]. The experimental results show that, out of the 1022 subjects who completed the survey, 579 (56.6%) reported headaches; 169 (16.5%) had migraines, 224 (21.9%) had tension-type headaches (TTHs), and 186 (18.1%) had other headache types. The category of other headache types included patients with a headache that, according to their answers, could not be classified as a migraine or TTHs.
The previous research [15] is continued with the idea to create a computerized-based intelligent system for diagnosing primary headaches to assist and help physicians make correct diagnosis for the patients who suffer from primary headache disorders and their treatments. This research was approved by the Ethics Committee of the Clinical Center of Vojvodina on 4 June 2020 (no. 00-214). This research includes only subjects with headache disorders, who were further analyzed. This experiment in total used the data set from 579 instances (subject–patients) and 20 features–attributes (answers to questions in the survey). The data set is divided in to three classes (groups of headaches), which included the following types of primary headache: Migraine, Tension Type Headache (TTH), and Other primary headaches (Other); where the missing data does not exist.
The previous study examined a number of studies based on IHS recommendations [22]. Different approaches were presented to attribute to the selection based on automatic methods, knowledge–based systems, expert systems, as well as physicians’ expert knowledge, as displayed in Table 2. The selection of features could be separated into Stochastic and no-Stochastic Feature Selection methodology, which presents a refinement of the initially used stochastic feature selection task with a no-stochastic method in order to further reduce the subset of features that were to be retained [23]. The previous study demonstrates that the most important features are the following: A4, A5, A6, A7, A8), A10, A12), A13 and A15; these are presented in Table 2 in black bold font. The emphasized nine features are used in the remain of the research [22].

2.3. Early Literature Research

During the last several decades, a sufficient number of approaches have been proposed to explain the clustering problem in medical data with the aim of helping physicians make decisions related to patient’s illness and their future treatments. One of the first studies, which contained 726 headache patient analyses, dates back to 1982 and it is presented in [24]. In that study, the cluster analysis was used as a method to find groups (clusters) of patients with similar symptoms. When only two clusters were required, the best two clusters were TTHs and migraine-like. However, in that research, eight clusters could also be distinguished, and the migraine group then became very small.
Based on previously mentioned, it can be concluded that the first step in the clustering problem, defining the number of clusters, is a very complicated question with the uncertain answer. After forming the clusters, the frequency of every symptom in every cluster was tabulated. Twelve physicians (eight internists, a cardiologist, a pediatrician, a neurologist, and a pathologist) were asked to use these symptom frequencies in order to label each cluster and then to prescribe what would a typical therapy be for every cluster. The “IF–THEN” rule algorithm was developed and subsequently tested, proving that 92.3% headache sufferers were placed in the correct cluster.

2.4. Optimal Number of Clusters

The Calinski-Harabasz index, as the concept of dense and well-separated clusters, is used to estimate the optimal number of clusters. Two measures were used, the Variance Ratio Criterion and the Total within Sum of Squares, in order to select the suitable c, i.e., a number of clusters. For creating the Calinski-Harabasz index, it was required to define the inter cluster dispersion [25] first. When N, i.e., the total number of observations (data points), and c number of clusters with their relative centroids and the global centroid, are known, the inter-cluster dispersion B(c) (between cluster variations) can be defined as follows:
B ( c )   =   t = 1 C n t   ( μ     μ t ) T   ( μ     μ t ) ,
In Equation (1), nt is the number of elements that belongs to the cluster c, µ is the global centroid, µt is the centroid of the cluster t.
The intra-cluster dispersion W(c), within cluster variation, is defined as:
W ( c ) =   t = 1 C x C N N ( x μ t ) T ( x μ t ) ,
The Calinski-Harabasz index is defined as the ratio between B(c) and W(c):
C H ( c )   =   N c c 1   ×   B ( c ) W ( c ) ,
The Calinski–Harabasz index is based on the comparison of the weighted ratio of the between cluster sum of squares (the cluster separation measure) and the within cluster sum of squares (the measure presenting how points are tightly packed within a cluster). The number of clusters maximizing this index needs to be recovered in order to determine low intra-cluster dispersion and high inter-cluster dispersion. Ideally, the clusters should be well separated, and the between cluster sum of squares value should be large. On the other hand, points within a cluster should be as close to one another as possible, which would result in smaller values of the within cluster sum of squares measure [25]. The decision that a point will be assigned to a cluster depends only on its features and sometimes on the position of a set of other points. Likewise, there are diverse algorithms based on alternative strategies to solve this problem, which can yield diverse results [26].

2.5. Analytical Hierarchy Process

Analytic Hierarchy Process (AHP) is one of multi criteria decision making methods that was originally developed and first time presented in [27]. It is a method to derive ratio scales from paired comparisons. The AHP starts by decomposing a complex, multi-criteria problem into a hierarchy where each level consists of a few manageable elements (group criteria Ci) which are then decomposed into another set of elements (criteria Cij). Later, these criteria (Cij) are mutually compared in order to get the priority of each criterion in hierarchy. Finally, all alternatives are compared in relation to the set of criteria (Cij) and in this way the comparison of alternatives (Ai) is obtained. The detail description of AHP flowchart is shown in Figure 1 [28].
In order to get the final result and in order to get the gradation of influence for each criterion, first it is necessary to compare criteria mutually. The comparison of any two criteria Ci and Cj with respect to the goal is made using the questions of the type: of the two criteria Ci and Cj which is more important and how much. In reality, however, an individual cannot give estimates such that they would conform to perfect consistency. It means that the final result might be accepted when consistency ratio (CR) is <0.1, less than 10%. Otherwise, the problem has to be revised [29].

2.6. Weighted Fuzzy C-means Clustering Algorithm

Fuzzy c-means (FCM) clustering algorithm is one of the most popular fuzzy clustering techniques, which was proposed in [30] and partially modified in [31]. It is an approach where the data points have their membership values with the cluster centers (centroids), which will be updated iteratively. In this research, the classical FCM algorithm is extended with the weight matrix of the features. The basic steps of the proposed Weight FCM clustering algorithm, employing a Euclidean distance norm, are summarized by the pseudo code shown in Algorithm 1.
Algorithm 1: Weighted fuzzy c-means clustering algorithm
Begin
Step 1:
---Initialization.
X, c, ε > 0, W
Step 2:
---Randomly select V cluster centers.
2 ≤ c ≤ N
Step 3:
---Choose an appropriate level of cluster fuzzinessf.
f [1, ∞], f > 1
Step 4:
---Choose an appropriate membership matrix U with size N × c × M
Uijm ∈ [0, 1] and j = 1 c U i j m   =   1 for a fixed value of m
Step 5:
---Calculate the cluster centers.
Repeat for jth cluster and its mth dimension
C C j m   =   i = 1 N U i j m f   x i m   w m i + 1 N U i j m f
Step 6:
---Calculate the Euclidean distance
D i j m   =   ( x i m   w m   C C j m   )
Step 7:
---Update fuzzy membership matrix U according to Dijm
U i j m   =   1 c = 1 C ( D i j m D i c m ) 2 f     1
Step 8:
Until U ≤ ε
End.
Let us suppose X = {x1, x2, x3,..., xN} be the set of data points, where N is the number of data points, in M-dimensional (M = 1, 2,..., m) data space. Let V = {v1, v2, v3,..., vc} be the set of cluster centers, and c represents the number of cluster centers. Let W be M-dimensional weighted matrix of the features. Let ε be the termination tolerance > 0 where termination criterion is between [0, 1] (Step 1). Assume the number of clusters c to be made, where 2 ≤ cN (Step 2). Choose an appropriate membership matrix U with size N × c × M, and conditions which are presented in Step 4. Determine the cluster centers CCjm, including weighted matrix of the features, for jth cluster and its mth dimension by using the expression in Step 5. Calculate the Euclidean distance Dijm between ith data point and jth cluster center with respect to mth dimension (Step 6). Update fuzzy membership matrix U according to Dijm, if Dijm > 0 (Step 7). Finally (Step 8), Repeat from Step 5 to Step 7 Until the changes in U ≤ ε, ε is a pre-specified termination criterion.

2.7. Implemented System for Diagnosing Primary Headaches

Figure 2 depicts the proposed hybrid intelligent system used for the diagnosis of primary headache disorders implemented in the research. The system comprises of two phases. The first phase is realized in four steps: (1) define and select nine important attributes of input data set based on our summarized research presented in Table 2, called new data set; (2) estimate the optimal and correct number of clusters calculated by Calinski–Harabasz index value; (3) according to previous estimation, select five most important attributes from new data set based on means and standard deviation, called Selected data; (4) apply Analytical Hierarchy Process to determine weighted coefficients for attributes in Selected data. It could be concluded that all activities in the first phase are focused on preparing the original data set with selection attributes and their weighted coefficients for implementation, which is called Define attributes and weights.
In the second phase, FCM clustering method is applied, and called Fuzzy Method and Decision. After fuzzy c-means clustering method, Solution presents three clusters with three cluster centers in five-dimensional data space. Three clusters display Diagnosed Patients with Primary Headache with adjoining types which suffer of: migraines, TTHs, and other primary headaches.

2.8. Metrics

In this research, the focus is on the multi-class problem. The confusion matrix for the multi-class problem is presented in Table 3.
TP (true positive) stands for the number of correctly classified positive instances. The existing elements in the form EXY present true negative, which stands for correctly classified negative examples, or false positive, as the numbers of incorrectly classified objects from classes. On the basis of the confusion matrix, the following metric can be calculated, showing the quality of the correct predictions. Accuracy is one of the most popular thresholds metrics. This is a measurement of the correct predictions between all predicted and actual values. Accuracy for the multi-class problem is presented by Equation (4):
A c c u r a c y   =   i T P i i T P i   +   i j E i j ,
The other popular thresholds metrics, which will be defined for the multi-class problem, are: precision (Equation (5)), recall (another term for Sensitivity) (Equation (6)) and F1 score (Equation (7)):
p r e c i s i o n   i   =     T P i j E j i ,
r e c a l l   i   =     T P i i E i j ,
F 1   s c o r e   =   2   ×   p r e c i s i o n   ×   r e c a l l p r e c i s i o n   +   r e c a l l ,
These appropriate metrics will be used for measuring quality for the proposed hybrid intelligent system for diagnosing primary headaches.

3. Results

The Input Data Set consists of 579 instances (patients); 169 of them have migraines, 224 have TTHs, and 186 have another primary headache type. Following the flowchart for the hybrid intelligent system for diagnosing primary headache presented in Figure 2, the first step in the first phase is defining nine important attributes. The important features are: A4, A5, A6, A7, A8, A10, A12, A13 and A15 presented in Table 2 and discussed in details in [22].
The next step is calculating Calinski-Harabasz index value presented in Figure 3. The Calinski-Harabasz index value is as follows: for two clusters it is 189; for three clusters it is 178; for four clusters it is 165; for five clusters it is 148; and for six clusters it is 123. As mentioned before, whole data set will be clustered in three classes; therefore, it is decided for Calinski-Harabasz index value to be 178 and the optimal number of clusters to be three (Figure 3).
In the next step in phase one, five most important attributes based on the diversity means and standard deviation values in three clusters are selected, and presented in Table 4. It could be concluded, after reducing the attributes, that the most important attributes are: A4, A6, A7, A8, and A15. These attributes describe the following characteristics of headache types: A4-length of headache attacks; A6-intensity of the pain; A7-quality of the pain; A8-headache is worsen after physical activities; andA15-tension or heightened tenderness of head during a headache attack.
In the last step of the first phase, AHP method is included. The pairwise comparison matrix, which compares the priority of the attributes, is presented in Table 5.
The consistency ratio (CR) is CR = 0.0174, which is less than 0.1, and pairwise comparison matrix is satisfied. Then, the priority weights for the attributes are: Attribute (4) is 0.2260; Attribute (5) is 0.1530; Attribute (6) is 0.2260; Attribute (8) is 0.1975; and Attribute (15) is 0.1975:
Priority Weights = [0.2260; 0.1530; 0.2260; 0.1975; 0.1975]
In the second phase, weighted phase fuzzy c-means clustering algorithm is applied. The confusion matrix for actual classes and predicted classes (Migraine, TTHs, Other) headaches is presented in Table 6.
Accuracy = (146 + 163 + 125)/579 = 0.7496,
Now, the Accuracy for the whole Hybrid Intelligent System for Diagnosing Primary Headaches can be calculated, and it is presented with Equation 4. Accuracy for the Hybrid Intelligent System for Diagnosing Primary Headaches is 0.7496 Equation (8), which is 75.0%, presenting a qualitative and respectable result.
The quality measurement for the proposed hybrid intelligent system for diagnosing primary headache is presented in Table 7. The appropriate metrics measuring qualities are: precision, recall and F1 score. In the multi-class problem, precision, recall, and F1 score are presented for every attribute–feature separately, and average values for every metrics are calculated. It is a different situation compared to metrics in two-class problem, when only one value for every metrics, precision, recall, and F1 score, is defined.
The implementation of the hybrid intelligent system for diagnosing primary headaches presents qualitative and respectable experimental results. The average of all metrics values is about 75%, which presents a qualitative result.

4. Discussion

The discussion in this section will be mostly focused on the usage of different clustering algorithms and decision support systems and their application in the diagnosis of primary headaches. The total number of 150 headache-prone subjects responded as follows: there were 49 migraines and 101 TTHs examined [32]. A cluster analysis was introduced, and the adjectives were thus grouped into seven clusters, which included five sensory and two affective sub-groups. Headache was most commonly described in relation to clusters for reflecting discomfort and aching pain sensations. There was not any marked difference in pain quality referred by migraine and TTH sufferers; however, the intensity of pain differentiated between these groups. In order to examine the relationship between descriptors, a complete linkage cluster analysis as one of several methods of agglomerative hierarchical clustering, was performed. Experimental results demonstrated the correspondence between migraine and TTH in two new cluster groups which were derived after using the cluster analysis: migraine, 93 of 101 (92.1%); TTH, 49 of 49 (100%); and misclassified, eight (5.3%) patients; thus, the diagnosis average accuracy is 94.7%.
A new proposed migraine analysis method was based on using electroencephalography signals under flash stimulation in time domain. These types of signals are normally pre-processed before the analysis procedure, where the pre-processing techniques affect the analysis results. Histogram differences in the case of flash stimulation are calculated and the results are used as features for the healthy subjects and migraine patients. These features are also applied to a k-means clustering algorithm in order to be able to observe clustering results of the proposed technique. Silhouette clustering results demonstrate that a good clustering performance is evaluated as 86.6% of the correct clustering rate in migraine patients [33].
The ant colony optimization classification algorithm is used for the diagnosis of primary headaches by using a website questionnaire expert system; it has shown the overall diagnosis accuracy as 96.9% [34]. In order to evaluate diagnosis accuracy of artificial immune system, algorithms for the classification of migraine, TTH and cluster headache have the maximum accuracy of 71% [35].
A hybrid fuzzy clustering approach for diagnosing primary headache disorder which combines fuzzy partitioning method and maximum likelihood estimation clustering algorithm is presented in [36]. Calinski–Harabasz index value is used to estimate the optimal and correct number of clusters. The proposed hybrid system is tested on the data set which has 579 patients who suffer from headaches, and it is facilitated by the application of the IHS criteria for diagnosing primary headache disorder. The headaches are clarified in four classes: migraine without aura (MWoA), migraine with aura (MWA), TTHs, and Other primary headaches. There are pairwise comparisons used between classes, where the minimum accuracy value is 78.1% and the maximum accuracy value is 86.4%. The Average Accuracy of that applied fuzzy clustering system is 82.9%.
Decision Support System (DSS) application determines the type of Headache using Fuzzy Multi Attribute Decision Making (FMADM) Method with AHP, as one of the methods that assists decision makers in making decisions; it is considered in [37]. A triangular membership function in fuzzy process is selected and used. The fuzzy compatibility index (FCI) is calculated for three types of headaches: FCI = 62.5% for migraines, FCI = 77.5% for TTHs, FCI = 57.5% for cluster headache.
The diagnostics system that was developed in Italy in 1998 is known as “IHS Diagnostic Criteria for Primary Headache”. The system was developed strictly based on the IHS operational diagnostic criteria. The system was tested on the computerized structured record in the cases of 500 patients in nine headache centers in Italy and the experimental results are depicted in [38]. The rate of concordance was calculated between the diagnosis provided by the computerized structured record and the one reported by physicians. The concordance between two diagnoses was recognized in 345 of 500 cases examined (69%). Further software implementation from 2004 was focused on 200 patients who reported primary chronic headaches. A certain diagnosis of a chronic form and a probable form of headaches was obtained in 50.8% cases [39].
In [40], a computerized program designed to diagnose primary headache was developed in Iran and tested in Tehran (Iran). It was based on International Classification of Headache Disorders, 2nd edition (ICHD-II) criteria which was to be used by physicians. A total of 80 patients were tested with the developed system. The software tool was able to come up with correct results in 78 out of 80 cases. Migraine headache accounted for 71 cases, five patients had TTHs, and two had cluster headaches; all were correctly diagnosed by the software. Two cases were not concordant with the neurologist’s diagnosis which represents an Average Accuracy of 97.5%.
The experimental results from this research are comparable with researches discussed in this Section, as well as with other researches. The Average Accuracy in this research is 75.0%, which is in the middle compared to other studies where Minimum Accuracy is 57.5% and Maximum Accuracy is 94.6%. On the other side, it is important to mention that in many papers methodologies are presented, but in fewer papers the detailed steps are provided, and only in a very small number of papers, are detailed steps of the methodology and intermediate experimental results presented, as is the case in our study. In other research papers, final experimental results are mainly presented without intermediate experimental results; therefore, it may be a bit difficult to follow the complete applied methodology.
Conversely, a number of applied statistical methods in the diagnosis of primary headache disorders is used for analyzing different features and variables. For example, descriptive statistics with categorical data presented using frequencies and percentages can be used for presenting continuous descriptive data that uses means and standard deviations, which can in turn be compared using the Chi-square test or Fisher’s exact test in order to compare two groups of categorical variables. Then, the Student ‘t’ test may be used, as well as Mann–Whitney ‘U’ test for continuous variables, or logistic regression analysis to determine binary dependent variables [41,42,43,44].
Regarding the socioeconomic impact on society, the studies have evaluated both the indirect costs of migraine as well as the direct costs [45,46,47]. The indirect cost of pain may be estimated with the number of missed work days. Indirect costs include the aggregate effects of migraine on productivity at work (paid employment), on household work, and in other roles [48]. Working people with migraine dropped out of work 14 days per year, which adds up to 6.8 million working days per year in Austria [49]. This remains a substantial economic factor. Migraine is estimated to account for about 1,154,336 lost working days each year in Flanders and Brussels [50]. The largest component of indirect costs is the productivity loss caused by absenteeism and reduced productivity; it is estimated that productivity losses due to migraine cost American employers 13 billion dollars per year [45]. A European study estimated that 5.7 working days were lost per year for every working or student migraineur (although the most disabled 10% of migraineurs accounted for 85% of the total), which projects to a loss of 25 million days from work or school each year in England [51]. A realistic, rational, and economical approach with adherence to the current guidelines for diagnosis and treatment is almost an impossible task for a practicing physician. In less developed countries, such as Serbia, this problem is even greater [13]. Therefore, it is necessary to research new automatic methods, intelligent and expert systems, or knowledge–based systems to help and improve physicians work to make better diagnoses and treatment for the patients.
Weaknesses in diagnosis occur in everyday clinical work, since atypical disease forms exist, and such cases are certainly present in our research as well. Likewise, the characteristics of the disease - which are presented as attributes in this study, do not mean that each observed case meets all the characteristics in 100%, as established in The International Classification of Headache Disorders–Third Edition (ICHD-3) [21]. The difference in clinical presentations in some cases leads to deviations from the recommended criteria in ICHD-3, which is a weakness of the recommended hybrid intelligent system for diagnosing primary headaches. On the other hand, the patient may suffer from several types of headaches that occur at different times. Therefore, the obtained answers, i.e. attributes, can be a combination of two types of headaches, leading to discrepancies in the proposed diagnosis provided by the intelligent system and the diagnosis by a neurologist, representing their expert assessment. The occurrence of two types of headaches with the same patient can be interpreted as the mean value of the observed attributes, having similar values for different types of headaches, as displayed in Table 4 for attributes A5, A10, A12, and A13 for TTHs and Other headache types.
The advantage of this system is the identification of the attributes that are dominant for different types of headaches. In this hybrid intelligent system for diagnosing primary headaches, their significance, in relation to other attributes, is increased by weighting coefficients in weighted fuzzy c-means clustering algorithm. The goal of our research is to help physicians have a reliable tool that would help them make the better diagnosis of primary headaches.
It can be concluded that the main contributions of our paper are:
  • The application of the Hybrid Intelligent System for Diagnosing Primary Headaches and computational complexity calculation.
  • Experimental evaluation of the application of the Hybrid Intelligent System for Diagnosing Primary Headaches and comparison with the selected state-of-the-art methods.

5. Conclusions

Headache disorders are the most prevalent of all neurological conditions. The goal of this research is focused on the design and implementation of a computer-based hybrid intelligent system for diagnosing primary headaches. This research applied various mathematical, statistical and artificial intelligence techniques, including decision-making methodology and clustering methods. The proposed hybrid intelligent system for diagnosing primary headaches is tested on real-world data set collected from the Clinical Centre of Vojvodina (Serbia). The Accuracy of the proposed intelligent system is 75.0%, which is a respectable value and a value comparable with similar research. Further research could be focused on creating new and more efficient tools and systems to help and improve physicians’ work and make diagnoses better.

Author Contributions

Conceptualization, S.S. and D.S.; Methodology, D.S. and S.D.S.; Formal Analysis, J.R.V.; S.D.S.; Investigation, S.S and S.R.S.; Resources, S.S. and D.S.; Data Curation, S.D.S.; Writing—Original Draft Preparation, S.S. and D.S.; Writing—Review & Editing, S.S. and D.S.; Visualization, S.R.S. and J.L.C.-R.; Supervision, J.R.V. and J.L.C.-R. All authors have read and agreed to the published version of the manuscript.

Funding

This research has been funded by the Spanish Ministry of Science and Innovation under project MINECO-TIN2017-84804-R and by the Grant FCGRUPIN-IDI/2018/000226 project from the Asturias Regional Government. This research has been supported by the Ministry of Education, Science and Technological Development, Republic of Serbia, through the project no. 451-03-68/2020-14/200156: “Innovative scientific and artistic research from the FTS (activity) domain”.

Institutional Review Board Statement

This research was approved by the Ethics Committee of the Clinical Center of Vojvodina on 4 June 2020 (no. 00-214).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Data is contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hagen, K.; Åsberg, N.A.; Uhlig, B.L.; Tronvik, E.; Brenner, E.; Stjern, M.; Helde, G.; Gravadahl, G.B.; Sand, T. The epidemiology of headache disorders: A face-to-face interview of participants in HUNT4. J. Headache Pain 2018, 19, 25. [Google Scholar] [CrossRef] [Green Version]
  2. Schwartz, B.S.; Stewart, W.F.; Simon, D.; Lipton, R.B. Epidemiology of tension-type headache. JAMA 1998, 279, 381–383. [Google Scholar] [CrossRef] [Green Version]
  3. Saylor, D.; Steiner, T.J. The global burden of headache. Semin. Neurol. 2018, 38, 182–190. [Google Scholar] [CrossRef]
  4. Burch, R.C.; Loder, S.; Loder, E.; Smitherman, T.A. The prevalence and burden of migraine and severe headache in the United States: Updated statistics from government health surveillance studies. Headache 2015, 55, 21–34. [Google Scholar] [CrossRef] [PubMed]
  5. Takeshima, T.; Wan, Q.; Zhang, Y.; Komon, M.; Stretton, S.; Rajan, N.; Treuer, T.; Ueda, K. Prevalence, burden, and clinical management of migraine in China, Japan, and South Korea: A comprehensive review of the literature. J. Headache Pain 2019, 20, 111. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Mennini, F.S.; Gitto, L.; Martelletti, P. Improving care through health economics analyses: Cost of illness and headache. J. Headache Pain 2008, 9, 199–206. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Agosti, R. Migraine burden of disease: From the patient’s experience to a socio-economic view. Headache 2018, 58, 17–32. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Pini, L.A.; Cainazzo, M.M.; Brovia, D. Risk–Benefit and cost–benefit ratio in headache treatment. J. Headache Pain 2005, 6, 315–318. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Andlin-Sobocki, P.; Jonsson, B.; Wittchen, H.U.; Olsen, J. Cost of disorder of the brain in Europe. J. Eur. Neurol. 2005, 12 (Suppl. 1), 1–27. [Google Scholar] [CrossRef] [PubMed]
  10. Hu, X.H.; Makson, L.E.; Lipton, R.B.; Stewart, W.F.; Berger, M.L. Burden of migraine in the United States: Disability and economic cost. Arch. Intern. Med. 1999, 159, 813–818. [Google Scholar] [CrossRef]
  11. Kernick, D.P. An introduction to the basic principles of health economics for those involved in the development and delivery of headache care. Cephalalgia 2005, 25, 709–714. [Google Scholar] [CrossRef]
  12. Porter, J.; Di Tanna, G.L.; Lipton, R.B.; Sapra, S.; Villa, G. Costs of acute headache medication use and productivity losses among patients with migraine: Insights from three randomized controlled trials. Pharm. Open 2019, 3, 411–417. [Google Scholar] [CrossRef] [Green Version]
  13. Becker, W.J.; Findlay, T.; Moga, C.; Scott, N.A.; Harstall, C.; Taenzer, P. Guideline for primary care management of headache in adults. Can. Fam. Physician 2015, 8, 670–679. [Google Scholar]
  14. Headache Classification Committee of the International Headache Society. The International Classification of Headache Disorders, 3rd edition (beta version). Cephalalgia 2013, 33, 629–808. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Simić, S.; Rabi-Žikić, T.; Villar, J.R.; Calvo-Rolle, J.L.; Simić, D.; Simić, S.D. Impact of individual headache types on the work and work efficiency of headache sufferers. Int. J. Environ. Res. Public Health 2020, 17, 6918. [Google Scholar] [CrossRef]
  16. Simić, S.; Simić, D.; Slankamenac, P.; Simić-Ivkov, M. Computer-Assisted diagnosis of primary headaches. In Lecture Notes in Computer Science; Springer: Berlin/Heidelberg, Germany, 2008; Volume 5271, pp. 314–321. [Google Scholar] [CrossRef]
  17. Simić, S.; Simić, D.; Slankamenac, P.; Simić-Ivkov, M. Rule-Based Fuzzy Logic System for Diagnosing Migraine. In Lecture Notes in Computer Science; Springer: Berlin/Heidelberg, Germany, 2008; Volume 5138, pp. 383–388. [Google Scholar] [CrossRef]
  18. Simić, S.; Banković, Z.; Simić, D.; Simić, S.D. A hybrid clustering approach for diagnosing medical diseases. In Lecture Notes in Computer Science; Springer: Cham, Switzerland, 2018; Volume 10870, pp. 741–752. [Google Scholar] [CrossRef]
  19. Simić, D.; Ilin, V.; Tanackov, I.; Svirčević, V.; Simić, S. A Hybrid Analytic Hierarchy process for clustering and ranking best location for logistics distribution center. In Lecture Notes in Computer Science; Springer: Cham, Switzerland, 2015; Volume 9121, pp. 477–488. [Google Scholar] [CrossRef]
  20. Simić, D.; Ilin, V.; Svirčević, V.; Simić, S. A hybrid clustering and ranking method for best positioned logistics distribution centre in Balkan Peninsula. Log. J. IGPL 2017, 25, 991–1005. [Google Scholar] [CrossRef]
  21. Headache Classification Committee of the International Headache Society. The international classification of headache disorders, 3rd edition. Cephalalgia 2018, 38, 1–211. [Google Scholar] [CrossRef]
  22. Krawczyk, B.; Simić, D.; Simić, S.; Wozniak, M. Automatic diagnosis of primary headaches by machine learning methods. Open Med. 2013, 8, 157–165. [Google Scholar] [CrossRef]
  23. Tallon-Ballesteros, A.J.; Correia, L.; Cho, S.-B. Stochastic and non-stochastic feature selection. In Lecture Notes in Computer Science; Springer: Cham, Switzerland, 2017; Volume 10585, pp. 592–598. [Google Scholar] [CrossRef]
  24. Diehr, P.; Diehr, G.; Koepsell, T.; Wood, R.; Beach, K.; Wolcott, B.; Tompkins, R.K. Cluster analysis to determine headache types. J. Chronic Dis. 1982, 35, 623–633. [Google Scholar] [CrossRef]
  25. Bonaccorso, G. Machine Learning Algorithms, 2nd ed.; Packt Publishing: Birmingham, UK, 2018. [Google Scholar]
  26. Babuška, R.; van der Veen, P.J.; Kaymak, U. Improved covariance estimation for Gustafson-Kessel clustering. In Proceedings of the 2002 IEEE International Conference on Fuzzy Systems, Honolulu, HI, USA, 12–17 May 2002; IEEE: New York, NY, USA, 2002; pp. 1081–1085. [Google Scholar] [CrossRef] [Green Version]
  27. Saaty, T.L. Scaling method for priorities in hierarchical structures. J. Math. Psychol. 1977, 15, 123–281. [Google Scholar] [CrossRef]
  28. Saaty, T.L. Analytic Hierarchy Process, Planning, Priority Setting, Resource Allocation; McGraw Hill Inc.: New York, NY, USA, 1980. [Google Scholar]
  29. Everitt, B.S.; Landau, S.; Leese, M.; Stahl, D. ClusterAnalysis; Wiley: Hoboken, NJ, USA, 2011. [Google Scholar]
  30. Dunn, J.C. A Fuzzy relative of the ISODATA Process and its use in detecting compact well-separated clusters. J. Cybern. 1973, 3, 32–57. [Google Scholar] [CrossRef]
  31. Bezdek, J.C. Pattern Recognition with Fuzzy Objective Function Algorithms; Kluwer Academic Publishers: Norwell, MA, USA, 1981. [Google Scholar]
  32. Hunter, M. The headache scale: A new approach to the assessment of headache pain based on pain descriptions. Pain 1983, 16, 361–373. [Google Scholar] [CrossRef]
  33. Alkan, A.; Akben, S.B. Use of k-means clustering in migraine detection by using EEG records under flash stimulation. Int. J. Phys. Sci. 2011, 6, 641–650. [Google Scholar] [CrossRef]
  34. Celik, U.; Yurtay, N. An ant colony optimization algorithm-based classification for the diagnosis of primary headaches using a website questionnaire expert system. Turk. J. Electr. Eng. Comput. Sci. 2017, 25, 4200–4210. [Google Scholar] [CrossRef]
  35. Celik, U.; Yurtay, N.; Koc, E.R.; Tepe, N.; Gulluoglu, H.; Ertas, M. Diagnostic accuracy comparison of artificial immune algorithms for primary headaches. Comput. Math. Methods Med. 2015, 465192. [Google Scholar] [CrossRef]
  36. Simić, S.; Banković, Z.; Villar, J.R.; Simić, D.; Simić, S.D. A hybrid fuzzy clustering approach for diagnosing primary headache disorder. Log. J. IGPL 2020. [Google Scholar] [CrossRef]
  37. Maseleno, A.; Tang, A.Y.C.; Mahmoud, M.A.; Othman, M.; Saputra, S.; Muslihudin, M. Fuzzy AHP method to determine headache types based on symptoms. Investig. Clin. 2017, 58, 71–82. [Google Scholar]
  38. Gallai, V.; Sarchielli, P.; Alberti, A.; Pedini, M.; Gallai, B.; Rossi, C.; Cittadini, E. Application of the 1988 International Headache Society diagnostic criteria in nine Italian headache centers using a computerized structured record. Headache 2002, 42, 1016–1024. [Google Scholar] [CrossRef] [PubMed]
  39. Sarchielli, P.; Pedini, M.; Alberti, A.; Rossi, C.; Baldi, A.; Corbelli, I.; Calabresi, P. Application of ICHD 2nd edition criteria for primary headaches with the aid of a computerised, structured medical record for the specialist. J. Headache Pain 2005, 6, 205–210. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Eslami, V.; Rouhani-Esfahani, S.; Hafezi-Nejad, N.; Refaeian, F.; Abdi, S.; Togha, M. A computerized expert system for diagnosing primary headache based on International Classification of Headache Disorder (ICHD-II). Springer Plus 2013, 2, 199. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Seo, J.-G.; Kim, K.T.; Moon, H.-J.; Do, J.K.; Kim, S.-Y.; Park, S.-P. Suicidality and its risk factors in tension-type headache patients: A multicenter case-control study. J. Clin. Neurosci. 2019, 69, 21–25. [Google Scholar] [CrossRef] [PubMed]
  42. Moore, C.; Leaver, A.; Sibbritt, D.; Adams, J. The features and burden of headaches within a chiropractic clinical population: A cross-sectional analysis. Complement. Ther. Med. 2020, 48, 102276. [Google Scholar] [CrossRef]
  43. Huang, Y.-S.; Tseng, W.-T.; Chiang, C.-Y.; Zhang, E.-W.; Wu, C.-H.; Cheng, F.-J. The effect of peer influence on the use of CT by emergency physicians for patients with headaches. Am. J. Emerg. Med. 2019, 37, 710–714. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Laskar, S.; Kalita, J.; Misra, U.K. Comparison of chronic daily headache with and without medication overuse headache using ICHD II R and ICHD 3 beta criteria. Clin. Neurol. Neurosurg. 2019, 183, 105382. [Google Scholar] [CrossRef]
  45. Stewart, W.F.; Ricci, J.A.; Chee, E.; Morganstein, D.; Lipton, R. Lost productive time and cost due to common pain conditions in the US workforce. JAMA 2003, 290, 2443–2454. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Duru, G.; Auray, J.P.; Gaudin, A.F.; Dartigues, J.F.; Henry, P.; Lanteri-Minet, M.; Lucas, C.; Pradalier, A.; Chazot, G.; El Hasnaoui, A. Impact of headache on quality of life in general population survey in France (GRIM 2000 Study). Headache 2004, 44, 571–580. [Google Scholar] [CrossRef]
  47. El Hasnaoui, A.; Vray, M.; Blin, P.; Nachit-Ouinekh, F.; Boureau, F.; HEMISFERE study group. Assessment of migraine severity using the MIGSEV scale: Relationship to migraine features and quality of life. Cephalalgia 2004, 24, 262–270. [Google Scholar] [CrossRef]
  48. Darrell, J.; Gaskin, D.J.; Richard, P. The economic costs of pain in the United States. J. Pain 2012, 13, 715–724. [Google Scholar] [CrossRef]
  49. Lampl, C.; Buzath, A.; Baumhackl, U.; Klingler, D. One-Year prevalence of migraine in Austria: A nation-wide survey. Cephalalgia 2003, 23, 280–286. [Google Scholar] [CrossRef] [PubMed]
  50. Moens, G.; Johannik, K.; Verbeek, C.; Bulterys, S. The prevalence and characteristics of migraine among the Belgian working population. ActaNeurol. Belg. 2007, 107, 84–90. [Google Scholar]
  51. Steiner, T.J.; Scher, A.I.; Stewart, W.F.; Kolodner, K.; Liberman, J.; Lipton, R.B. The prevalence and disability burden of adult migraine in England and their relationships to age, gender and ethnicity. Cephalalgia 2003, 23, 519–527. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Analytical Hierarchy Process methodology [28].
Figure 1. Analytical Hierarchy Process methodology [28].
Ijerph 18 01890 g001
Figure 2. Flowchart of the Hybrid Intelligent System for Diagnosing Primary Headaches.
Figure 2. Flowchart of the Hybrid Intelligent System for Diagnosing Primary Headaches.
Ijerph 18 01890 g002
Figure 3. Calinski-Harabasz Index values for the estimation of the optimal number of clusters.
Figure 3. Calinski-Harabasz Index values for the estimation of the optimal number of clusters.
Ijerph 18 01890 g003
Table 1. The International Classification of Headache Disorders–the Third Edition [21].
Table 1. The International Classification of Headache Disorders–the Third Edition [21].
ICHD-3 CodeDiagnosis
A Primary Headache Disorders
1. Migraine
1.1Migraine without aura (MWoA)
1.2Migraine with aura (MWA)
¦¦
1.6Episodic syndromes that may be associated with migraine
2. Tension-type headache (TTH)
2.1Infrequent episodic tension-type headache
¦¦
2.4Probable tension-type headache
3. Trigeminal autonomic cephalalgias (TACs)
4. Other primary headache disorders
B5. Secondary headache disorders
¦
12.
Table 2. Comparison of selected attributes for primary headache on the basis of IHS diagnostic criteria: 1. consistency measure filter, 2. relief greedy, 3. relief top10, 4. genetic algorithm wrapper, 5. ACO-based classification algorithm, 6. rule-based fuzzy logic system, 7. physician’s expert choice, 8. column RES—final decision for the important attribute selection [22].
Table 2. Comparison of selected attributes for primary headache on the basis of IHS diagnostic criteria: 1. consistency measure filter, 2. relief greedy, 3. relief top10, 4. genetic algorithm wrapper, 5. ACO-based classification algorithm, 6. rule-based fuzzy logic system, 7. physician’s expert choice, 8. column RES—final decision for the important attribute selection [22].
Attributes1234567RES
1Sex
2How old were you when the headache occurred for the first time?
3How often do you have headache attacks?
4How long do the headache attacks last?
5Where is the headache located?
6How intense is the pain?
7What is the quality of the pain you experience?
8Do your headaches worsen after physical activities such as walking?
9Do you avoid routine physical activities because you fear they might trigger your headache?
10Are the headaches accompanied by? a) Nausea
11Are the headaches accompanied by? b) Vomiting
12Are the headaches accompanied by? c) Photophobia
13Are the headaches accompanied by? d) Phonophobia
14Do you have temporary visual, sensory or speech disturbance?
15Do you, during a headache attack, have tension and/or heightened tenderness of head or neck muscles? TTH
16Do you have any body numbness or weakness?
17Do you have any indications of oncoming headache?
18Headache is usually triggered by: Menstrual periods
19In the half or my visual field, lasting 5 minutes to an hour, along with the headache attack or an hour before.
20Along with the headache attack or an hour before one I have sensory symptoms.
Table 3. Confusion matrix for the multi-class problem.
Table 3. Confusion matrix for the multi-class problem.
Actual Class
Predicted class MigraineTTHOther
MigraineTPMETMEOM
TTHEMTTPTEOT
OtherEMOETOTPO
Table 4. Matrix of Means and standard deviation values per classes.
Table 4. Matrix of Means and standard deviation values per classes.
ClassesAttributes
(4)(5)(6)(7)(8)(10)(12)(13)(15)
Migraine2.0 ± 0.22.0 ± 1.0 2.8 ± 0.61.3 ± 0.51.4 ± 0.51.2 ± 0.41.3 ± 0.41.2 ± 0.41.3 ± 0.4
TTHs2.3 ± 0.62.7 ± 1.0 2.0 ± 0.61.8 ± 0.51.7 ± 0.51.7 ± 0.51.7 ± 0.51.6 ± 0.51.1 ± 0.4
Other1.5 ± 0.82.7 ± 0.72.0 ± 0.71.7 ± 0.61.7 ± 0.41.7 ± 0.51.6 ± 0.51.6 ± 0.51.4 ± 0.5
Bold numbers: the most important attributes are: A4, A6, A7, A8, and A15.
Table 5. Pairwise comparison matrix.
Table 5. Pairwise comparison matrix.
Attribute (4)Attribute (6)Attribute (7)Attribute (8)Attribute (15)
Attribute (4)12111
Attribute (6)½1½11
Attribute (7)12111
Attribute (8)11111
Attribute (15)11111
Table 6. Confusion matrix for actual classes and predicted classes.
Table 6. Confusion matrix for actual classes and predicted classes.
Actual Classes
MigraineTTHsOtherΣ
Predicted classesMigraine1464427217
TTH2116334218
Other217125144
Total 169224186579
Table 7. Quality measurement for the Hybrid Intelligent System for Diagnosing Primary Headache.
Table 7. Quality measurement for the Hybrid Intelligent System for Diagnosing Primary Headache.
MigraineTTHsOtherAverage
Precision %67.374.886.876.3
Recall %86.472.767.275.4
F1score %75.673.775.775.0
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Simić, S.; Villar, J.R.; Calvo-Rolle, J.L.; Sekulić, S.R.; Simić, S.D.; Simić, D. An Application of a Hybrid Intelligent System for Diagnosing Primary Headaches. Int. J. Environ. Res. Public Health 2021, 18, 1890. https://doi.org/10.3390/ijerph18041890

AMA Style

Simić S, Villar JR, Calvo-Rolle JL, Sekulić SR, Simić SD, Simić D. An Application of a Hybrid Intelligent System for Diagnosing Primary Headaches. International Journal of Environmental Research and Public Health. 2021; 18(4):1890. https://doi.org/10.3390/ijerph18041890

Chicago/Turabian Style

Simić, Svetlana, José R. Villar, José Luis Calvo-Rolle, Slobodan R. Sekulić, Svetislav D. Simić, and Dragan Simić. 2021. "An Application of a Hybrid Intelligent System for Diagnosing Primary Headaches" International Journal of Environmental Research and Public Health 18, no. 4: 1890. https://doi.org/10.3390/ijerph18041890

APA Style

Simić, S., Villar, J. R., Calvo-Rolle, J. L., Sekulić, S. R., Simić, S. D., & Simić, D. (2021). An Application of a Hybrid Intelligent System for Diagnosing Primary Headaches. International Journal of Environmental Research and Public Health, 18(4), 1890. https://doi.org/10.3390/ijerph18041890

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop