Asymptotic stability of deterministic and stochastic prey-predator models with prey herd immigration

1.   Introduction

Gastric Cancer (GC) is one of the most commonly-occurring cancerous types across the globe. However, its diagnosis rate is comparatively poor and has become a major threat globally ^[1]. Malignant tissues may grow in any part of the Gastrointestinal (GI) tract since it gets directly and constantly exposed to carcinogens in the gut environment via food ingestion. Various authors have identified several hotspots for malevolent conversion of the tissues that get continuously exposed to carcinogens. GI tumors cause 35% of the cancer-related mortality and 26% of tumor incidence globally ^[2]. The current increase in the mortality and incidence rates can be connected to the increasing occurrence of adaptable risk components like unhealthy diet, sedentary lifestyle, other metabolic abnormalities and obesity ^[3]. The outcomes from the recent analysis have exposed a substantial surge in GI cancer occurrence among the individuals aged between 25 and 49 years, thus warning the public about the evolution of medical complications ^[4]. The five main GI tumors lead to substantial global problems, such as hepatocellular carcinoma, esophageal squamous cell carcinoma, colorectal cancer, pancreatic cancer and gastric adenocarcinoma ^[5]. Unfortunately, early diagnosis of the GI-based tumors is frequently missed due to the absence of particular indications. In most of the cases, the medical attention is sought only when non-specific indications turn out to be unbearable ^[6,7].

With the fast expansion of Computer Vision (CV) technologies, particularly, the development of medical image classifiers, any type of electron microscopic images can be inspected in a rapid and an efficient manner ^[8]. Therefore, it presents a chance to realize a solution that is relevant to the diagnosis of GC. Particularly, the image classification process serves a significant role in computer-aided diagnostics ^[9]. The outcomes attained from the image classification approaches are predominantly utilized as significant references by doctors in terms of differentiating the malignant tumors from the benign ones, differentiating the differentiation phases of cancers and distinguishing the subtypes in cancer ^[10,11]. Moreover, with the progression of medical image classification technology, the main aim is to have a high anti-interference ability and attain the maximum precision.

Computer-related analysis of the tissue images is gaining attention among the researchers in digital pathology ^[12]. In recent times, Deep Learning (DL) techniques have been leveraged to classify and distinguish different types of cancers including the GC. DL utilizing the Convolutional Neural Network (CNN) systems includes visual recognition capability since it can discern the features directly from large training datasets, thus outpacing the human beings ^[13]. With the current adoption of digital pathology methods such as the Whole-Slide Imaging (WSI) technique for key prognosis, there is an explosive growth experienced in the digitalized tissue slides. This in turn renders massive volumes of digital pathological images ^[14,15]. But, the application of DL in digital histopathology is still in a nascent stage. So, the accuracy and efficacy of the pathology analysis are enhanced by merging the routine digitalization of DL and WSI.

The current research article introduces an Anas Platyrhynchos Optimizer (APO) with DL-based GC classification (APODL-GCC) technique for endoscopic images. The presented APODL-GCC technique employs a contrast enhancement technique. Next, the feature extraction process is performed with the help of Neural Architectural Search network (NASNet) model to generate a collection of feature vectors. For hyperparameter optimization, the APO algorithm is used which enhances the classification performance. Finally, the GC classification process is performed using the Deep Belief Network (DBN) method. The proposed APODL-GCC technique was simulated using the medical images and the outcomes were achieved.

2.   Related works

In literature ^[16], a new publicly-available gastric histopathology sub-size image database (GasHisSDB) was published to determine the performance of the classifiers. In order to prove that the erstwhile image classification techniques have differences on the GasHisSDB, the authors selected various classifiers for evaluation. A new transformer-based classifier, seven conventional machine learning techniques and three CNN classifiers were chosen to validate the image classification tasks. In the study conducted earlier ^[17], an RF deep feature selection technique was modelled. By compiling the Copy Number Variation (CNV) data and gene expression data, the dimensions of the multi-omics data were minimized and the classification accuracy was improved using a DNN technique and the RF. Li et al. ^[18] projected a DL-related structure called GT-Net for automatic segmentation of the GC. The devised GT-Net implemented diverse structures for deep and shallow layers for a superior extraction of the features.

Zhu et al. ^[19] built a CNN-Computer-Aided Diagnostic (CAD) mechanism related to endoscopic images to screen the patients and determine the depth of invasion using endoscopic resection. An Artificial Intelligence (AI)-related CNN-CAD system was advanced in this study with the help of TL through an existing pre-trained CNN architecture i.e., ResNet50. Lee et al. ^[20] presented a CADx mechanism that differentiates and categorizes the gastric tumors from other pre-cancerous circumstances like gastritis, gastric polyps, bleeding and gastric ulcers. This mechanism leveraged the DL method i.e., Xception which involves depthwise separable convolution for the categorization of non-cancers and cancer. The presented technique had stages such as the fast and robust FCM method for image segmentation during preprocessing stage. These techniques produce a possible technique for differentiation and the classification of tumors from other gastric diseases.

Yoon et al. ^[21] developed an optimized method for depth prediction and EGC detection. In this study, the authors inspected the elements that hinder the outcomes of AI-based diagnosis. The authors used the Visual Geometry Group (VGG)-16 method to classify the endoscopic images as either non-EGC or EGC (T1a or T1b). To construct a method for the identification of EGC areas during training, the authors devised a new loss function that concurrently measures the localization and classification errors. In the study conducted earlier ^[22], TL was implemented via a fine-tuned deep CNN technique for mechanical categorization of the M-NBI images into two groups such as the EGC images and normal gastric images. Furthermore, this study also explored the impact of TL upon the classifier's outcome under four aspects such as the size of the network input image, training dataset, basic structures of the deep CNN and the number of optimally-tuned network layers; this study even presented a certain set of guidelines to be followed in this domain for later studies.

3.   Proposed model

In the current research work, a new APODL-GCC approach has been modelled for GC classification on endoscopic images. The presented APODL-GCC technique focuses on the identification of GC using the CV and DL concepts. It comprises a sequence of operations such as the contrast enhancement, feature extraction using the NASNet, APO-based hyperparameter tuning and the DBN classification. Figure 1 illustrates the workflow of the APODL-GCC approach.

3.1.   Image preprocessing

Primarily, the APODL-GCC technique exploits the contrast enhancement technique. CLAHE is different from the Adaptive Histogram Equalization (AHE) method it overcomes the over-amplification of the contrast. CLAHE is applied to smaller regions of the image that are termed as tiles, instead of the whole image. The adjacent tiles are subsequently compiled through bilinear interpolation so as to remove the artificial boundary. This method is executed to enrich the image contrast.

3.2.   Feature extraction using the optimal NASNet model

At this stage, the feature extraction process is performed with the help of NASNet model to produce a collection of feature vectors. Zoph and Le proposed a neural structure search network, named as the NASNet model. This model uses Reinforcement Learning and Recurrent Neural Network (RNN) technique for training so as to attain the most precise parameters ^[23]. Constructing a CNN model requires extensive computational time, in case if the content is larger as in the case of ImageNet data. The CNN model searches for a better structure from the smaller data and transmits it so that it can be used for training using the larger dataset. This phenomenon is referred to as learning transferable structure. The NASNet structure is scaled based on the quantity of information.

Here, N time denotes the time required to identify and segregate the number of cancerous cells from the normal cells. In this stage, the RNN approach is implemented as a search technique. An example operation for detection and classification of the cancerous and normal cells is given below. This is achieved by searching the controller RNN for a proper architecture.

● 1 × 7 then 7 × 1 convolutional layer

● 3 × 3 average pooling layer

● 5 × 5 max pooling layer

● 1 × 1 convolutional layer

● 3 × 3 depth-wise separable convolutional layer

● 7 × 7 depth-wise separable convolutional layer

● 1 × 3 then 3 × 1 convolutional layer

● 3 × 3 dilated convolutional layer

● 3 × 3 max pooling layer

● 7 × 7 max pooling layer

● 3 × 3 convolutional layer

● 5 × 5 depth-wise separable convolutional layer

For hyperparameter optimization, the APO algorithm is used which in turn enhances the classification performance. The population in the APO is established based on the Anas platyrhynchos performance method as given below ^[24].

In Eq (1), $Po{p}_{i}$ represents the ${i}^{th}$ population, $N$ denotes the population size, $low$ and $up$ indicate the lower and upper boundaries of the search space respectively and $rand$ denotes the arbitrary values selected in the range of [0, 1]. With regards to warning performance, the fly-in-danger functions are recognized by the selected probability $Pc$. The basic method of warning classification is established using the following equation:

Step1: Evaluate the probability of distress $Pc$ using the following expression.

In Eq (2), $fit\left(Po{p}_{i}\right)$ indicates the fitness value of $Po{p}_{i}$, and $rank\left(fit\left(Po{p}_{i}\right)\right)$ is considered as a rank of individual $Po{p}_{i}$ values among other individuals in the population.

Step2: Once the probability $Pc$ is attained, a novel individual is constructed based on the equation given below.

Now, $t$ represents the present iteration, $Po{p}_{est}$ shows the lead duck, ${\alpha }_{0} > 0$ denotes the step length scale factor and $sign$ signifies the sign function. Levy Flight (LF) is responsible for arbitrary walks and the distribution equation is formulated as given below.

LF is a distinct type of arbitrary walk, and the distribution probability of the step length follows a heavy‐tailed distribution as given below.

In Eq (5), $s$ indicates the LF step length. Moreover, $\lambda = 1+\beta$; also, ${\alpha }_{0} = 0.01$ and $\beta = 3/2$, as CS is fixed. Additionally, $\mu$ and $\nu$ are selected for the standard distribution $\mu = N\left(0, {\delta }_{\mu }^{2}\right)$ and $\nu = N\left(0, {\delta }_{v}^{2}\right)$.

An intelligent technique is selected for this study to simulate the biological migration process. Similar to the SOMA, the APO essentially concentrates on the Anas platyrhynchos, thus simulating its migration movement. SOMA and APO differ from each other in terms of place upgrade method and the flow of the technique. In this performance, a basic technique is established as follows.

Step1: Once an optimal particle is defined, then other search particles try to shift towards the best particles which is mathematically expressed as follows:

In Eq (7), $A$ and $C$ indicate the co-efficient vectors that are obtained as given below.

Here, $a$ indicates the co-efficient vector that gets linearly decreased with multiple iterations. The value of $a$ is expressed as follows.

In this expression, $T$ indicates the maximal number of iterations.

Step2: Once the novel individual solution is found to be worse than the older solution, then other particle $Po{p}_{and}$ is randomly selected

Step 3: When $Po{p}_{and}$ is better than Pop, then the ${i}^{th}$ individual is transferred to an arbitrary particle $Po{p}_{and}$ as given below.

Now $l$ shows the distance of the arbitrary particles and the ${i}^{th}$ individual.

Step4: When $Po{p}_{and}$ corresponds to $Po{p}_{i}$, it remains the same.

Step5: Once the $Po{p}_{rad}$ is less than $Po{p}_{i}$, then the particle moves towards the ${i}^{th}$ individual in a random manner using the following expression.

Fitness selection is an essential feature in APO system. The solution encoded is exploited to calculate the aptitude (goodness) of the candidate solutions. At this point, the accuracy value is the main criterion, employed in this study, to decide the fitness function.

Here, $TP$ denotes the True Positive values whereas $FP$ stands for the False Positive value.

3.3.   GC classification using DBN

Lastly, the GC classification process is carried out with the help of the DBN model. During the classification process, the normalized feature vector is provided as the input to the DBN mechanism for both detection and the classification of the GC ^[25]. Artificial Neural Network (ANN) is modeled by dissimilar output and hidden unit layers and is named as DL. It encompasses two levels, namely,

● Pre-training phase

● Fine-tuning phase

Pre-training phase

In DBN technique, an essential feed forward network and a deep structure are involved in which the sample is derived from the input to output layers through the highest quantity of hidden nodes that possess additional nodes. Based on the application of the DBN approach, the algorithm generates the activation function according to the hidden unit that differentiates the DBN method. Moreover, the Restricted Boltzmann Machine (RBM) is defined to resolve the problems encountered in potential activation function. Figure 2 showcases the framework of DBN. RBM is a type of Markov subjective field that encompasses the individual layers with stochastic hidden unit.

Step 1: Initializing the rich unit $\nu$ for RBM training.

In Eq (4), ${W}_{kl}$ indicates the symmetric transmission between the ${v}_{k}$ visible and ${h}_{l}$ hidden layers, the term $\alpha \wedge \beta$ indicates the biases and $K$ and $L$ show the quantity of visible and hidden layers respectively. The subordinate of the $log$ possibility, of the preparation vector, with regards to its weight is irregular. From the hidden unit of RBM, there is no unpredicted replies that tend to accomplish the impartial sample from ${\left({V}_{k}, {h}_{l}\right)}_{data}$.

In Eq (16), $\delta \left(x\right)$ indicates the logistic sigmoid function, $\frac{1}{\left(1+exp\left(x\right)\right)}$ and ${v}_{k}and{h}_{l}$ shows the unbiased samples.

Updating phase

The hidden layer is upgraded whereas the visible unit is regarded as a simultaneous element from both visible and the hidden layers. It leads to a complex technique as given below.

In Eq (17), the RBM undergoes training whereas the divergent RBM is constructed with a multiple layer approach. In general, different RBMs are stacked together and the visible input unit is identified as a quality measure. The vector for the units are effectively placed in the RBM layer that is shared over the application of shared model in the present weight and bias. Consequently, the final layer is officially trained to secure the RBM. Thereby, the obtained DNN weights are occupied in the fine-tuning phase.

4.   Experimental validation

The proposed model was simulated using Python 3.6.5 tool on a PC configured with 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. In this section, the experimental outcomes of the proposed APODL-GCC method were validated using the GC dataset ^[26]. This dataset has a total of 1500 images as briefed in Table 1. Figure 3 showcases some of the sample images from the dataset. The dataset holds images under three classes with 500 samples under each class.

The GC classification output of the proposed APODL-GCC model is shown in Figure 4 in the form of confusion matrix. The results infer that the APODL-GCC model accomplished an effectual categorization of the dataset under healthy and GC classes.

Table 2 and Figure 5 portray the GC classification results accomplished by the proposed APODL-GCC model with 80: 20 of TR/TS databases. The outcomes denote that the APODL-GCC model recognized both healthy and the GC cases accurately. On 80% of TR database, the APODL-GCC model gained an average $acc{u}_{y}$ of 99.67%, $pre{c}_{n}$ of 99.50%, $rec{a}_{l}$ of 99.50%, ${F}_{score}$ of 99.50%, and an MCC of 99.25%. Meanwhile, on 20% of TS database, the APODL-GCC method obtained an average $acc{u}_{y}$ of 99.11%, $pre{c}_{n}$ of 98.69%, $rec{a}_{l}$ of 98.67%, ${F}_{score}$ of 98.67%, and an MCC of 98.01%.

Table 3 and Figure 6 showcase the GC classification results achieved by the proposed APODL-GCC approach with 70: 30 of TR/TS databases. The experimental outcomes infer that the APODL-GCC approach recognized both healthy and GC cases accurately. On 70% of TR database, the APODL-GCC method reached an average $acc{u}_{y}$ of 97.84%, $pre{c}_{n}$ of 96.81%, $rec{a}_{l}$ of 96.76%, ${F}_{score}$ of 96.77%, and an MCC of 95.17%. In the meantime, on 30% of TS database, the APODL-GCC method attained an average $acc{u}_{y}$ of 98.67%, $pre{c}_{n}$ of 98.03%, $rec{a}_{l}$ of 98%, ${F}_{score}$ of 97.97%, and an MCC of 97.03%.

The TACC and VACC values, obtained by the proposed APODL-GCC technique in terms of GC classification, are shown in Figure 7. The figure shows that the APODL-GCC approach exhibited an improved performance with increased TACC and VACC values. Notably, the APODL-GCC technique reached the maximum TACC outcomes.

The TLS and VLS values, accomplished by the proposed APODL-GCC approach in terms of GC classification, are portrayed in Figure 8. The figure infers that the APODL-GCC technique achieved better performance with minimal TLS and VLS values. Seemingly, the APODL-GCC approach produced low VLS outcomes.

A clear precision-recall inspection was conducted upon the APODL-GCC approach using the test database and the results are shown in Figure 9. The figure exhibits that the APODL-GCC technique achieved enhanced precision-recall values under all the classes.

The GC classification outcomes of the APODL-GCC technique were compared with that of the results achieved by the existing DL methods and the results are shown in Table 4 and Figure 10 ^[26]. The outcomes found the close performance of the existing models such as SSD, CNN, and Mask R-CNN models with $acc{u}_{y}$ values being 96.22, 96.77, and 96.97% respectively. Concurrently, the U-Net-CNN method and the cascade CNN method obtained moderately improved $acc{u}_{y}$ values such as 98.21 and 97.07% respectively.

On the contrary, the MRFOTL-GCDC technique managed to achieve a considerable GC classifier outcome with an $acc{u}_{y}$ of 99.17. However, the proposed APODL-GCC model obtained the maximum performance, in terms of GC classification, with an $acc{u}_{y}$ of 99.67%. Therefore, the APODL-GCC model can be employed in the near future for accurate GC classification.

5.   Conclusions

In the current research work, a new APODL-GCC approach has been developed for GC classification on endoscopy images. The presented APODL-GCC method focuses on the identification of the GC using both CV and DL concepts. Primarily, the APODL-GCC technique employs a contrast enhancement technique. Next, the feature extraction process is performed with the help of the NASNet model to generate a collection of feature vectors. For hyperparameter optimization, the APO algorithm is used which in turn enhances the classification performance. Finally, the GC classification process is executed by the DBN model. The proposed APODL-GCC technique was simulated using the medical images. The experimental results confirmed that the proposed APODL-GCC technique accomplished an enhanced performance over other models. Therefore, the proposed model can be employed in real-time GC detection and classification tasks, which will help the physicians in decision making process. In the future, the performance of the APODL-GCC method can be boosted by the ensemble learning models.

Acknowledgments

The authors extend their appreciation to the Deputyship for Research & Innovation, Ministry of Education in Saudi Arabia for funding this research work through the project no. (IFKSURG-1742).

Conflict 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.