Next Article in Journal
Vinpocetine Ameliorates Metabolic-Syndrome-Associated Bladder Overactivity in Fructose-Fed Rats by Restoring Succinate-Modulated cAMP Levels and Exerting Anti-Inflammatory Effects in the Bladder Detrusor Muscle
Next Article in Special Issue
Rapid Convolutional Neural Networks for Gram-Stained Image Classification at Inference Time on Mobile Devices: Empirical Study from Transfer Learning to Optimization
Previous Article in Journal
Plasmacytoma Variant Translocation 1 (PVT1) Gene as a Potential Novel Target for the Treatment of Diabetic Nephropathy
Previous Article in Special Issue
How Resilient Are Deep Learning Models in Medical Image Analysis? The Case of the Moment-Based Adversarial Attack (Mb-AdA)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Artificial Intelligence Driven Biomedical Image Classification for Robust Rheumatoid Arthritis Classification

by
Marwa Obayya
1,
Mohammad Alamgeer
2,
Jaber S. Alzahrani
3,
Rana Alabdan
4,
Fahd N. Al-Wesabi
5,*,
Abdullah Mohamed
6 and
Mohamed Ibrahim Alsaid Hassan
7
1
Department of Biomedical Engineering, College of Engineering, Princess Nourah bint Abdulrahman University, P.O. Box 84428, Riyadh 11671, Saudi Arabia
2
Department of Information Systems, College of Science & Art at Mahayil, King Khalid University, Abha 62217, Saudi Arabia
3
Department of Industrial Engineering, College of Engineering at Alqunfudah, Umm Al-Qura University, Mecca 92254, Saudi Arabia
4
Department of Information Systems, College of Computer and Information Science, Majmaah University, Al-Majmaah 11952, Saudi Arabia
5
Department of Computer Science, College of Science & Art at Mahayil, King Khalid University, Abha 62217, Saudi Arabia
6
Research Centre, Future University in Egypt, New Cairo 11845, Egypt
7
Department of Computer and Self Development, Preparatory Year Deanship, Prince Sattam bin Abdulaziz University, AlKharj 16242, Saudi Arabia
*
Author to whom correspondence should be addressed.
Biomedicines 2022, 10(11), 2714; https://doi.org/10.3390/biomedicines10112714
Submission received: 28 September 2022 / Revised: 16 October 2022 / Accepted: 24 October 2022 / Published: 26 October 2022
(This article belongs to the Special Issue Artificial Intelligence in Biological and Biomedical Imaging 2.0)

Abstract

:
Recently, artificial intelligence (AI) including machine learning (ML) and deep learning (DL) models has been commonly employed for the automated disease diagnosis process. AI in biological and biomedical imaging is an emerging area and will be a future trend in the field. At the same time, biomedical images can be used for the classification of Rheumatoid arthritis (RA) diseases. RA is an autoimmune illness that affects the musculoskeletal system causing systemic, inflammatory and chronic effects. The disease frequently becomes progressive and decreases physical function, causing articular damage, suffering, and fatigue. After a time, RA causes harm to the cartilage of the joints and bones, weakens the tendons and joints, and finally causes joint destruction. Sensors (thermal infrared camera sensor, accelerometers and wearable sensors) are more commonly employed to collect data for RA. This study develops an Automated Rheumatoid Arthritis Classification using an Arithmetic Optimization Algorithm with Deep Learning (ARAC-AOADL) model. The goal of the presented ARAC-AOADL technique lies in the classification of health disorders depending upon RA and orthopaedics. Primarily, the presented ARAC-AOADL technique pre-processes the input images by median filtering (MF) technique. Then, the ARAC-AOADL technique uses AOA with an enhanced capsule network (ECN) model to produce feature vectors. For RA classification, the ARAC-AOADL technique uses a multi-kernel extreme learning machine (MKELM) model. The experimental result analysis of the ARAC-AOADL technique on a benchmark dataset reported a maximum accuracy of 98.57%. Therefore, the ARAC-AOADL technique can be employed for accurate and timely RA classification.

1. Introduction

Biomedical imaging acts as a vital part in the domain of biology and biomedicine, offering data related to the structural and functional mechanism of cells and the human body. Biological and biomedical imaging comprises microscopy, molecular imaging, pathological imaging, optical coherence tomography, nuclear medicine, ultrasound imaging, X-ray radiography, computed tomography, magnetic resonance imaging, and so on. The typical medical phenotype is labelled as non-infectious persistent polyarticular swelling, specifically of small joints, that causes progressive joint destruction and deformity and bone erosion, i.e., rheumatoid arthritis (RA) [1]. However, persistent periarticular synovitis, regardless of the immunopathogenesis, was linked to joint erosion and destruction, so regardless of the immune trigger a similar medical phenotype arises [2]. A logical extension of this view is that clinically different cases are linked with the medical RA phenotype. An unfulfilled necessity occurs in the translational setting for developing a robust pattern for assessing, diagnosing and prognosing patients affected with polyarthritis characteristic of early RA, particularly in recent times where the key role of autoinflammation or innate immunity is effectively detected in other chronic inflammatory illnesses [3]. Here, a new classification is proposed for the full medical disease spectrum of RA by utilizing the paradigmatic shift that occurred, in addition to the description of autoimmunity against citrullinated antigens in numerous RA cases. This immunological disease continuum method of inflammation in RA has consequences for therapeutic techniques [4,5].
Particularly viable reasons for applications utilizing the ensemble ML technique include hospital-based applications, controlling smart homes, information on-request systems, monitoring systems, outpatient care and mobile games communication interface [6]. Additionally, machine learning (ML) can be used in the development and evaluation of electroencephalogram (EEG)-related brain activities to measure using a biosensor [7]. With the developments in on-body wearable sensors and sensor technology, the ML technique will perform effectively in RA disease categorization. The author used a wearable glove-related sensor IoT for identifying RA diseases. Certainly, thermal structure-related camera sensors can be used to monitor temperature disparities in finger joints in several analyses [8]. ML includes different algorithms, procedures and techniques for finding limited associations within particular data and to produce tools constituting prescription, prediction or description combinations. ML has been initiated with many clinical domains and has been exemplified as very accurate in classifying and identifying different illnesses [9]. ML is a commonly used method for enhancing medical services and disease diagnosis with medical data development in various medical fields. Studies of effective applications of artificial intelligence (AI), which includes deep learning (DL) and ML methods. have seen an exponential growth in healthcare and medical fields [10]. Such techniques are critical in offering high-quality care to patients with RA.
This study develops an Automated Rheumatoid Arthritis Classification using an Arithmetic Optimization Algorithm with Deep Learning (ARAC-AOADL) model. The goal of the presented ARAC-AOADL technique lies in the classification of health disorders depending upon RA and orthopaedics. Primarily, the presented ARAC-AOADL technique pre-processes the input images by median filtering (MF) technique. Then, the ARAC-AOADL technique uses AOA with an enhanced capsule network (ECN) model to produce feature vectors. For RA classification, the ARAC-AOADL technique uses a multi-kernel extreme learning machine (MKELM) model. The experimental results analysis of the ARAC-AOADL technique is tested on two medical datasets and the results are examined under several aspects.

2. Related Works

Lim et al. [11] modelled a new feature engineering technique compiling potentially functional coding haplotypes (pfcHap), including ML feature selection to detect biologically meaningful, probably causative genetic factors, that considers effective SNP–SNP interactions in the pfcHap to optimally forecast the methotrexate (MTX) response in RA patients. Ahalya et al. [12], by utilizing modified pre-trained CNN techniques, produced automated patch-related classification of hand Xray images, and then for for automated classification and feature extraction of hand Xray images and, for comparing the efficiency of CNN techniques with linear and non-linear kernels, a customized CNN technique was developed; they finally classified the normal and RA by employing ML methods and framing the hand-crafted feature fusion (SIFT and Customized CNN features).
Yang et al. [13] introduced a grading technique to estimate and detect the texture and geometric features of bone erosion and synovium thickening. This study will use the metrics and texture features of ROI in a dissimilar way to previous studies in this area. The segmented outcomes were examined for the extraction of three quantitative geometric variables, which are integrated with GLCM statistic texture features to describe the ultrasonic image of metacarpophalangeal RA. Tang et al. [14] modelled an automated RA grading technique leveraging DCNN for assistance in medical assessment. Here, the input is the Gray-scale ultrasound images of finger joints whereas the output is the RA grading outcomes. The authors executed data augmentation for increasing the training samples count. The authors pre-trained the GoogLeNet on ImageNet as an extracting feature and then fine-tuned them.
In [15], the classification of clinical disorders related to RA and orthopaedics dataset utilizing Ensemble techniques was conferred. The RA data was collected from the study of WBC classification by utilizing features derived from the lymphocyte image obtained through a digital microscope. The orthopaedic datasets are a benchmark dataset for this work, since they imposed the same classifier issue with some numerical features. In this study, three ensemble techniques, random subspace, bagging and Adaboost, were used. Such ensemble techniques use RF and kNN as the base learners of the ensemble techniques. In [16], for enhancement of disease risk evaluation, ML and matrix factorization methods were combined to find significant and implicit risk factors. A new structure was modelled that can successfully evaluate early disease risks and RA was employed as a case study. This structure has three main phases: 1early disease risk assessment, data preprocessing, and risk factor optimization. This was the first study compiling ML and matrix factorization for disease risk evaluation implemented with nation-wide and longitudinal medical diagnostic databases. Andreu-Perez et al. [17] formulated a technique that could produce optimally grained actigraphies for capturing the effect of the disease on the daily actions of patients. A study of processing techniques related to ML and DL was offered.

3. The Proposed Model

In this study, we have developed a new ARAC-AOADL technique for accurate RA classification, which helps to identify the health disorders depending upon RA and orthopaedics. The presented ARAC-AOADL technique encompasses MF-based noise removal, ECN feature extraction, AOA hyperparameter tuning and MKELM classification. The working of the 1ARAC-AOADL technique is depicted in Figure 1.

3.1. Noise Filtering Technique

Primarily, the presented ARAC-AOADL technique pre-processes the input images by the MF technique. Initially, the input image was preprocessed by means of the MF algorithm for getting rid of the noise within them. MF based on specificity is one of the applications from medical image noise extraction. The major concept behindhand MF is to present a process for assembling each neighborhood from the increasing order, which selects the median values of arranged numbers and replaces the central pixel as follows:
y i ,   j = m e d i a n x i , j ,   i ,   j C ,  
where C indicates the central neighborhood location of the image. In such cases, the MF was executed by digital noise extraction from the input images, whereby the filter mask with size 3 × 3 was applied.

3.2. Feature Extraction Using Optimal ECN Model

In this study, the ARAC-AOADL technique uses AOA with the ECN models to produce the feature vector. In the presented technique, the splitted pixel set of the images can be labelled as a collection of nerve cells corresponding to the capsule [18]. Consider Y i [healthy, tumor] as i - th output capsules, and w e i j signifies the weighted matrices in the following:
y ^ i j = w e i j y i j
In Equation (5), y ^ i j describes the detection vector that diagnoses output parent j - th capsules with i - th capsule, and pixel range was applied to evaluate the weight quantity. The quantity for the weight was improved if the value was decreased or the pixel involves the positive group. The softmax method is exploited by the preceding layer capsule and the potential parent capsule as a coefficient is encoded c i j wherein major logits b i j show the log preceding probability of i - th routing capsules in the preceding layer to j - th capsules in the succeeding layer. Generally, the “routing-by-agreement” methodology was executed by logit for the capsule in all the layers:
c i j = e b i j i e b i j
The previous layer demonstrates key elements to compute the input of j - th parent capsules as follows:
s j = i c i j y ^ ( i | j )
The compressed pixel vector can be defined in (0, 1) by a non-linear method called squashing and it is computed by the following expression:
v a j = s j | | 2 1 + s j 2 × s j ε + | | s j 2
where ε = 10 7 . The subsequent layer capsule was attained using:
a i j = v a j × y ^ i j
The entire capsule classifier is regarded as margin loss L o s s k in the class capsule k for capsule network based on the loss:
L o s s k = T k   max   ( 0 ,   m + v a k ) 2 + λ 1 T k   max   ( 0 ,   v a k m ) 2
For the hyperparameter selection process, the AOA is exploited. AOA mainly replicates the use of the arithmetical operator during the arithmetical problem-solving method. Arithmetic is the part of mathematics that is exploited to handle the property of operation and numbers. An arithmetical operator is an operation symbol that executes fundamental arithmetic, i.e., a symbol applied to four processes. During optimization, these operators are utilized for selecting the best solution from the candidates. The optimization technique is based on two main processes: exploration and development. Initially, the search space of candidate solution is expansively protected for breaking the deadlock of the methodology within the search stagnation. Then, the performance for the solution is searched more deeply.
The series of possible solutions was arbitrarily generated in the initial stage of the optimization method of AOA, as given below.
X = x 1 , 1 x 1 , j x 1 , n 1 x 1 , n x 2 , 1 x 2 , j x 2 , n x N 1 ,   1 x N 1 , j x N 1 , n   x N , 1 x N , j x N , n 1 x N , n              
Formerly, the AOA implements the optimized methodology, necessitating the resolution of the search process in accordance with the value of the Math Optimizer Accelerated (MOA) process that is calculated by the following equation.
M O A C _ I t e r = Min + C _ I t e r × Max   Min   M I t e r
In Equation (9), M O A C _ I t e r indicates the function value in C _ I t e r iteration; C _ I t e r represent the present iteration; M _ I t e r denotes the maximal iterations count; Min and Max are accelerated function minimal and maximal values.
The exploration phase in the AOA method is realized generally by Division (D) and Multiplication (M) operators [19]. During mathematical computation, these two operators accomplish distributed value for a wide-ranging coverage of candidate solutions. The location of candidate solution is upgraded considerably in the exploration procedure, as follows:
x i , j C _ I t e r + 1 = b e s t x j / M O P + ε × U B j L B j × μ + L B j r 2 < 0.5 b e s t x j × M O P × U B j L B j × μ + L B j o t h e r w i s e
In Equation (16), χ i , j C _ I t e r + 1 denotes the jth location of ith solution in C _ I t e r + 1 th iteration; ε indicates the small values; U B and L B indicates the upper and lower bounds of the location of candidate solution; μ is employed for regulating the exploration stage set to 0.5; the MOP signifies the math optimizer probability of AOA that is described in the following:
M O P   C _ I t e r = 1 C I t e r 1 α M I t e r 1 α
Now, α defines the accuracy of exploitation on the iteration, α = 5 . The execution of exploitation method depends mainly on Subtraction (S) and Addition (A) operators that are easier to cause minimal dispersion, for the candidate solution is executed by a deep searching with larger probability of estimating the optimal solution [20]. In the growth step, the candidate solution was upgraded as follows:
x i , j C I t e r + 1 = b e s t x j M O P × ( U B j L B j × μ + L B j , r 3 < 05 b e s t x j + M 0 P × ( U B j L B j × μ + L B j , o t h e r w i s e
The adaptive conversions among the exploration and exploitation stages are supported by the AOA method which defines an optimal solution and continues with a diversity of possible solutions to conduct a wide-ranging search as illusetrated in Algorithm 1.
Algorithm 1. Pseudocode of AOA
Initialization of the parameter pop-size (N) and maximal iteration (T)
Initialization of the location of every search agent X i i = 1 , 2 ,   ,   N
Set the parameters Min and Max
While t T
    Evaluate the fitness of all the search agent Upgrade bestFitness, X b
    Evaluate the MOP
    Evaluate the M O A
    For every search agent
        If r a n d > M O A
            Upgrade position
        Else
            Upgrade position
        End if
    End for
                    t = t + 1
End While
Return best Fitness, X b

3.3. RA Classification Model

For RA classification, the ARAC-AOADL technique employed the MKELM model. The standard KELM is a single kernel based model. The structure of standard ELM is shown in Figure 2. Since a distinct kernel function provided a similar measure to the sample point, the efficacy of the kernel function could be based considerably to the related dataset. The input signal has features of considerable amount, irregular distribution of instances generated using imbalance, and maximal dimension feature space. Utilizing a single kernel to process the input dataset could not solve the problems effectively. The kernel function was regarded as global or local kernel function depending on rotation or translation invariances [21]. The global kernel function was higher at removing global features, and the local kernel function was better at eliminating the local feature of instance. In multi-kernel learning, an optimal kernel was regarded as linear integration of the group of base kernels, and the better linear integration coefficient and the classification parameter are learned equally using the margin maximization. The polynomial kernel and RBF are global and local kernel functions with optimal efficacy.
In order to balance the integration of the classifier efficiency and generalized capability, an MKELM has been generated using linear integration of the polynomial kernels and RBF and they are defined in the following equation:
K m i x = λ k r b f u , u i + 1 λ k p o l y u , u j )
where λ ( 0 < λ < 1 ) denotes the weighted coefficient of linear integration.
k r b f u , u j = exp u u i 2 2 σ 2 ,
k p o l y u , u i = u · u i + 1 d ,
From the expression, d indicates fixed to two, as the dimension of polynomial space refers to n d ; once the sample size corresponded to thousands and the index corresponded to three, the dimension can be accomplished as 1 billion, and the computation of the inner product generates a dimension disaster [22].
Eventually, the resulting objective of MKELM was defined by the following expression:
f x = K m i x u , u 1 K m i x u , u N ( 1 C + M ) 1   T .

4. Experimental Validation

The proposed model is simulated using Python 3.6.5 tool on PC i5-8600k, GeForce 1050Ti 4GB, 16GB RAM, 250GB SSD, and 1TB HDD. The parameter settings are given as follows: learning rate: 0.01, dropout: 0.5, batch size: 5, epoch count: 50, and activation: ReLU. The experimental validation of the ARAC-AOADL model is performed on a benchmark dataset from the Kaggle repository [23]. The dataset holds 310 samples with three classes, as given in Table 1.
Figure 3 shows the confusion matrices of the ARAC-AOADL model on various training (TR) and testing (TS) data. The figure highlights that the ARAC-AOADL model has reached effective RA classification results.
Table 2 offers a detailed classifier results of the ARAC-AOADL model on 80% of TR data and 20% of TS data. The results reveal the RA classification results of the ARAC-AOADL model on 80% of TR data. The ARAC-AOADL model has recognized Hernia class samples with a c c u y of 96.77%, p r e c n of 93.88%, r e c a l of 90.20%, F s c o r e of 92%, and A U C s c o r e of 94.34%. In addition, the ARAC-AOADL model has categorized Normal class samples with a c c u y of 97.58%, p r e c n of 96.15%, r e c a l of 96.15%, F s c o r e of 96.15%, and A U C s c o r e of 97.19%. Moreover, the ARAC-AOADL model has reached average a c c u y of 97.04%, p r e c n of 95.30%, r e c a l of 94.61%, F s c o r e of 94.94%, and A U C s c o r e of 96.11%. The ARAC-AOADL method has recognized Hernia class samples with a c c u y of 98.39%, p r e c n of 90%, r e c a l of 100.00%, F s c o r e of 94.74%, and A U C s c o r e of 99.06%. Furthermore, the ARAC-AOADL technique has categorized Normal class samples with a c c u y of 98.39%, p r e c n of 100%, r e c a l of 95.45%, F s c o r e of 97.67%, and A U C s c o r e of 97.73%. In addition, the ARAC-AOADL approach has reached average a c c u y of 97.85%, p r e c n of 95.59%, r e c a l of 97.41%, F s c o r e of 96.40%, and A U C s c o r e of 97.85%.
Table 3 shows the overall RA classification results on 70% of TR data and 30% of TS data. Figure 4 demonstrates the RA classification outcomes of the ARAC-AOADL technique on 70% of TR data. The ARAC-AOADL system has recognized Hernia class samples with a c c u y of 98.62%, p r e c n of 97.67%, r e c a l of 95.45%, F s c o r e of 96.55%, and A U C s c o r e of 97.44%. Furthermore, the ARAC-AOADL method has categorized Normal class samples with a c c u y of 94.47%, p r e c n of 89.19%, r e c a l of 94.29%, F s c o r e of 91.67%, and A U C s c o r e of 94.42%. Besides, the ARAC-AOADL technique has obtained average a c c u y of 96.01%, p r e c n of 94.29%, r e c a l of 94.31%, F s c o r e of 94.27%, and A U C s c o r e of 95.57%.
A brief training accuracy (TRAC) and validation accuracy (VAAC) of the ARAC-AOADL model is given in Figure 4. The results inferred that the ARAC-AOADL model has reached maximum TRAC and VAAC values. It is obvious that the VAAC is superior to TRAC.
In Figure 5, a clear training loss (TRAL) and validation loss (VALL) of the ARAC-AOADL model is reported. The figure reported that the ARAC-AOADL model has reached minimal values of TRAL and VALL.
To highlight the enhanced performance of the ARAC-AOADL model, a comparison study is made in Table 4 and Figure 6. An extensive comparison study of the presented ARAC-AOADL model with existing ML models in terms of a c c u y and F s c o r e is provided in Figure 6. The experimental values inferred that the ARAC-AOADL model has shown effective classification performance. For instance, based on a c c u y , the ARAC-AOADL model has offered higher a c c u y of 98.57% whereas the Bagging, Adaboost, DT, and Subspace-k-NN random models have attained lower a c c u y of 94.89%, 89.37%, 92.64%, and 97.50% respectively. Furthermore, based on F s c o r e , the ARAC-AOADL technique has provided maximum a c c u y of 97.67% while the Bagging, Adaboost, DT, and Subspace-k-NN random techniques have accomplished lower a c c u y of 95.23%, 89.21%, 94.68%, and 96.82%, correspondingly.
A comprehensive analysis of the proposed ARAC-AOADL methodology with current ML models with respect to p r e c n and r e c a l is given in Figure 7. The experimental value demonstrates that the ARAC-AOADL approach has demonstrated effective classification performance. For example, based on p r e c n , the ARAC-AOADL technique has provided maximum p r e c n of 98.22% while the Bagging, Adaboost, DT, and Subspace-k-NN random models have accomplished lower p r e c n of 94.89%, 89.37%, 92.64%, and 97.50% correspondingly. Furthermore, based on r e c a l , the ARAC-AOADL technique has given maximum r e c a l of 97.67% while the Bagging, Adaboost, DT, and Subspace-k-NN random approaches have accomplished lower r e c a l of 95.40%, 90.01%, 94.73%, and 97.02%, correspondingly.
After examining the detailed results, the proposed model has gained enhanced performance with maximum a c c u y of 98.57%, p r e c n of 98.22%, r e c a l of 97.21%, and F s c o r e of 97.67%. The enhanced performacne of the proposed model is due to the inclusion of the AOA based hyperparameter tuning process. Since the trial and error hyperparmaeter selection is not an effective process, the optimal hyperparmater tuning AOA helps to accomplish enhanced RA classification performance. Therefore, the proposed model can be employed for precise RA classification, which enables detection of health disorders based on RA and orthopaedics.

5. Conclusions

In this study, we have developed a new ARAC-AOADL technique for accurate RA classification, which helps to identify health disorders depending upon RA and orthopaedics. Primarily, the presented ARAC-AOADL technique pre-processes the input images by the MF technique. Then, the ARAC-AOADL technique uses AOA with the ECN model to produce feature vectors. For RA classification, the ARAC-AOADL technique employed the MKELM model. The experimental result analysis of the ARAC-AOADL technique is tested on two medical datasets and the results are inspected under several aspects. The simulation results ensured the enhancements of the ARAC-AOADL technique in terms of different measures. In future, we can extend the ARAC-AOADL technique by hybrid DL classification models.

Author Contributions

Conceptualization, M.O. and F.N.A.-W.; methodology, M.A.; software, A.M.; validation, M.O., J.S.A. and R.A.; formal analysis, A.M.; investigation, F.N.A.-W.; resources, A.M.; data curation, M.I.A.H.; writing—original draft preparation, M.O., F.N.A.-W., J.S.A. and R.A.; writing—review and editing, M.A. and M.I.A.H.; visualization, A.M.; supervision, M.O.; project administration, F.N.A.-W.; funding acquisition, M.O., M.A and J.S.A. All authors have read and agreed to the published version of the manuscript.

Funding

The authors extend their appreciation to the Deanship of Scientific Research at King Khalid University for funding this work through Large Groups Project under grant number (180/43). Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2022R203), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia. The authors would like to thank the Deanship of Scientific Research at Umm Al-Qura University for supporting this work by Grant Code: 22UQU4340237DSR51.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing not applicable to this article as no datasets were generated during the current study.

Conflicts of Interest

The authors declare that they have no conflict of interest. The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript.

References

  1. Myasoedova, E.; Athreya, A.P.; Crowson, C.S.; Davis, J., III; Warrington, K.J.; Walchak, R.C.; Carlson, E.; Kalari, K.R.; Bongartz, T.; Tak, P.P.; et al. Toward Individualized Prediction of Response to Methotrexate in Early Rheumatoid Arthritis: A Pharmacogenomics-Driven Machine Learning Approach. Arthritis Care Res. 2022, 74, 879–888. [Google Scholar] [CrossRef] [PubMed]
  2. Maarseveen, T.D.; Meinderink, T.; Reinders, M.J.T.; Knitza, J.; Huizinga, T.W.J.; Kleyer, A.; Simon, D.; Akker, E.B.V.D.; Knevel, R. Machine learning electronic health record identification of patients with rheumatoid arthritis: Algorithm pipeline development and validation study. JMIR Med. Inform. 2020, 8, e23930. [Google Scholar] [CrossRef] [PubMed]
  3. Goñi, M.; Basu, N.; Murray, A.D.; Waiter, G.D. Brain predictors of fatigue in rheumatoid arthritis: A machine learning study. PLoS ONE 2022, 17, e0269952. [Google Scholar] [CrossRef] [PubMed]
  4. Prasad, B.; McGeough, C.; Eakin, A.; Ahmed, T.; Small, D.; Gardiner, P.; Pendleton, A.; Wright, G.; Bjourson, A.J.; Gibson, D.S.; et al. ATRPred: A machine learning based tool for clinical decision making of anti-TNF treatment in rheumatoid arthritis patients. PLoS Comput. Biol. 2022, 18, e1010204. [Google Scholar] [CrossRef] [PubMed]
  5. Yoosuf, N.; Maciejewski, M.; Ziemek, D.; Jelinsky, S.A.; Folkersen, L.; Müller, M.; Sahlström, P.; Vivar, N.; Catrina, A.; Berg, L.; et al. Early prediction of clinical response to anti-TNF treatment using multi-omics and machine learning in rheumatoid arthritis. Rheumatology 2022, 61, 1680–1689. [Google Scholar] [CrossRef] [PubMed]
  6. Lötsch, J.; Alfredsson, L.; Lampa, J. Machine-learning-based knowledge discovery in rheumatoid arthritis-related registry data to identify predictors of persistent pain. Pain 2019, 161, 114–126. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Crowson, C.S.; Gunderson, T.M.; Davis, J.M., III; Myasoedova, E.; Kronzer, V.L.; Coffey, C.M.; Atkinson, E.J. Using Unsupervised Machine Learning Methods to Cluster Comorbidities in a Population-based Cohort of Patients with Rheumatoid Arthritis. Arthritis Care Res. 2022. [Google Scholar] [CrossRef] [PubMed]
  8. Singh, U.V.; Gupta, E.; Choudhury, T. Detection of rheumatoid arthritis using machine learning. In Proceedings of the 2019 International Conference on Computational Intelligence and Knowledge Economy (ICCIKE), Dubai, United Arab Emirates, 11–12 December 2019. [Google Scholar]
  9. Sharon, H.; Elamvazuthi, I.; Lu, C.; Parasuraman, S.; Natarajan, E. Classification of Rheumatoid Arthritis using Machine Learning Algorithms. In Proceedings of the 2019 IEEE Student Conference on Research and Development (SCOReD), Bandar Seri Iskandar, Malaysia, 15–17 October 2019; pp. 245–250. [Google Scholar] [CrossRef]
  10. Kedra, J.; Davergne, T.; Braithwaite, B.; Servy, H.; Gossec, L. Machine learning approaches to improve disease management of patients with rheumatoid arthritis: Review and future directions. Expert Rev. Clin. Immunol. 2021, 17, 1311–1321. [Google Scholar] [CrossRef] [PubMed]
  11. Lim, A.J.; Lim, L.J.; Ooi, B.N.; Koh, E.T.; Tan, J.W.L.; Chong, S.S.; Khor, C.C.; Tucker-Kellogg, L.; Leong, K.P.; Lee, C.G.; et al. Functional coding haplotypes and machine-learning feature elimination identifies predictors of Methotrexate Response in Rheumatoid Arthritis patients. eBioMedicine 2022, 75, 103800. [Google Scholar] [CrossRef] [PubMed]
  12. Ahalya, R.K.; Umapathy, S.; Krishnan, P.T.; Raj, A.N.J. Automated evaluation of rheumatoid arthritis from hand radiographs using Machine Learning and deep learning techniques. Proc. Inst. Mech. Eng. Part H J. Eng. Med. 2022, 236, 1238–1249. [Google Scholar] [CrossRef] [PubMed]
  13. Yang, T.; Zhu, H.; Gao, X.; Zhang, Y.; Hui, Y.; Wang, F. Grading of metacarpophalangeal rheumatoid arthritis on ultrasound images using machine learning algorithms. IEEE Access 2020, 8, 67137–67146. [Google Scholar] [CrossRef]
  14. Tang, J.; Jin, Z.; Zhou, X.; Chu, H.; Yuan, J.; Wu, M.; Cheng, Q.; Wang, X. Grading of rheumatoid arthritis on ultrasound images with deep convolutional neural network. In Proceedings of the 2018 IEEE International Ultrasonics Symposium (IUS), Kobe, Japan, 22–25 October 2018; pp. 1–4. [Google Scholar]
  15. Sharon, H.; Elamvazuthi, I.; Lu, C.-K.; Parasuraman, S.; Natarajan, E. Development of rheumatoid arthritis classification from electronic image sensor using ensemble method. Sensors 2019, 20, 167. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Chin, C.-Y.; Hsieh, S.-Y.; Tseng, V.S. eDRAM: Effective early disease risk assessment with matrix factorization on a large-scale medical database: A case study on rheumatoid arthritis. PLoS ONE 2018, 13, e0207579. [Google Scholar] [CrossRef] [PubMed]
  17. Andreu-Perez, J.; Garcia-Gancedo, L.; McKinnell, J.; Van Der Drift, A.; Powell, A.; Hamy, V.; Keller, T.; Yang, G.-Z. Developing fine-grained actigraphies for rheumatoid arthritis patients from a single accelerometer using machine learning. Sensors 2017, 17, 2113. [Google Scholar] [CrossRef] [PubMed]
  18. Zhang, Z.; Ye, S.; Liao, P.; Liu, Y.; Su, G.; Sun, Y. Enhanced capsule network for medical image classification. In Proceedings of the 2020 42nd Annual International Conference of the IEEE Engineering in Medicine & Biology Society (EMBC), Montreal, QC, Canada, 20–24 July 2020; pp. 1544–1547. [Google Scholar]
  19. Abualigah, L.; Diabat, A.; Mirjalili, S.; Abd Elaziz, M.; Gandomi, A.H. The arithmetic optimization algorithm. Comput. Methods Appl. Mech. Eng. 2021, 376, 113609. [Google Scholar] [CrossRef]
  20. Kaveh, A.; Hamedani, K.B. January. Improved arithmetic optimization algorithm and its application to discrete structural optimization. In Structures; Elsevier: Amsterdam, The Netherlands, 2021; Volume 35, pp. 748–764. [Google Scholar]
  21. Bisoi, R.; Dash, P.K.; Das, P.P. Short-term electricity price forecasting and classification in smart grids using optimized multikernel extreme learning machine. Neural Comput. Appl. 2018, 32, 1457–1480. [Google Scholar] [CrossRef]
  22. Xing, Y.; Ban, X.; Guo, C. Probabilistic forecasting of traffic flow using multikernel based extreme learning machine. Sci. Program. 2017, 2017, 2073680. [Google Scholar] [CrossRef]
  23. Biomechanical Features of Orthopedic Patients. Available online: https://www.kaggle.com/datasets/uciml/biomechanical-features-of-orthopedic-patients?resource=download (accessed on 12 August 2022).
Figure 1. Working Process of ARAC-AOADL technique.
Figure 1. Working Process of ARAC-AOADL technique.
Biomedicines 10 02714 g001
Figure 2. Structure of ELM model.
Figure 2. Structure of ELM model.
Biomedicines 10 02714 g002
Figure 3. Confusion matrices of ARAC-AOADL model.
Figure 3. Confusion matrices of ARAC-AOADL model.
Biomedicines 10 02714 g003
Figure 4. TRAC and VAAC of ARAC-AOADL model.
Figure 4. TRAC and VAAC of ARAC-AOADL model.
Biomedicines 10 02714 g004
Figure 5. TRAL and VALL of ARAC-AOADL model.
Figure 5. TRAL and VALL of ARAC-AOADL model.
Biomedicines 10 02714 g005
Figure 6. Comparative RA Classification results of ARAC-AOADL model interms of a c c u y and F s c o r e .
Figure 6. Comparative RA Classification results of ARAC-AOADL model interms of a c c u y and F s c o r e .
Biomedicines 10 02714 g006
Figure 7. Comparative RA Classification results of ARAC-AOADL model interms of p r e c n and r e c a l .
Figure 7. Comparative RA Classification results of ARAC-AOADL model interms of p r e c n and r e c a l .
Biomedicines 10 02714 g007
Table 1. Dataset used.
Table 1. Dataset used.
ClassNo. of Instances
Hernia60
Spondylolisthesis150
Normal100
Total Number of Samples310
Table 2. RA Classification results of ARAC-AOADL model on 80:20 of TR/TS data.
Table 2. RA Classification results of ARAC-AOADL model on 80:20 of TR/TS data.
Training/Testing (80:20)
LabelsAccuracyPrecisionRecallF-ScoreAUC Score
Training Phase
Hernia96.7793.8890.2092.0094.34
Spondylolisthesis96.7795.8797.4896.6796.80
Normal97.5896.1596.1596.1597.19
Average97.0495.3094.6194.9496.11
Testing Phase
Hernia98.3990.00100.0094.7499.06
Spondylolisthesis96.7796.7796.7796.7796.77
Normal98.39100.0095.4597.6797.73
Average97.8595.5997.4196.4097.85
Table 3. RA Classification results of ARAC-AOADL model on 70:30 of TR/TS data.
Table 3. RA Classification results of ARAC-AOADL model on 70:30 of TR/TS data.
Training/Testing (70:30)
LabelsAccuracyPrecisionRecallF-ScoreAUC Score
Training Phase
Hernia98.6297.6795.4596.5597.44
Spondylolisthesis94.9396.0093.2094.5894.85
Normal94.4789.1994.2991.6794.42
Average96.0194.2994.3194.2795.57
Testing Phase
Hernia98.92100.0093.7596.7796.88
Spondylolisthesis97.8597.8797.8797.8797.85
Normal98.9296.77100.0098.3699.21
Average98.5798.2297.2197.6797.98
Table 4. Comparative RA Classification results of ARAC-AOADL model.
Table 4. Comparative RA Classification results of ARAC-AOADL model.
MethodsAccuracyPrecisionRecallF-Score
ARAC-AOADL98.5798.2297.2197.67
Bagging94.8994.9795.495.23
Adaboost89.3789.9690.0189.21
DT92.6493.3294.7394.68
Subspace-k-NN Random97.5097.8397.0296.82
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Obayya, M.; Alamgeer, M.; S. Alzahrani, J.; Alabdan, R.; N. Al-Wesabi, F.; Mohamed, A.; Alsaid Hassan, M.I. Artificial Intelligence Driven Biomedical Image Classification for Robust Rheumatoid Arthritis Classification. Biomedicines 2022, 10, 2714. https://doi.org/10.3390/biomedicines10112714

AMA Style

Obayya M, Alamgeer M, S. Alzahrani J, Alabdan R, N. Al-Wesabi F, Mohamed A, Alsaid Hassan MI. Artificial Intelligence Driven Biomedical Image Classification for Robust Rheumatoid Arthritis Classification. Biomedicines. 2022; 10(11):2714. https://doi.org/10.3390/biomedicines10112714

Chicago/Turabian Style

Obayya, Marwa, Mohammad Alamgeer, Jaber S. Alzahrani, Rana Alabdan, Fahd N. Al-Wesabi, Abdullah Mohamed, and Mohamed Ibrahim Alsaid Hassan. 2022. "Artificial Intelligence Driven Biomedical Image Classification for Robust Rheumatoid Arthritis Classification" Biomedicines 10, no. 11: 2714. https://doi.org/10.3390/biomedicines10112714

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