Ensemble Machine Learning Approach for Parkinson’s Disease Detection Using Speech Signals
Abstract
:1. Introduction
1.1. Contributions
- In this study, we employ a comprehensive dataset comprising voice recordings from PD patients and healthy individuals, allowing for a detailed analysis of speech characteristics and the extraction of relevant features. This dataset encompasses a wide range of clinically significant attributes, including time-frequency fading (TFF) features, Mel frequency cepstral coefficients (MFCCs), vocal fold features (VFF), and tremor waveform quality time (TWQT), among others. By using these diverse features, the proposed model can capture subtle variations in speech patterns associated with PD pathology, thereby enhancing diagnostic accuracy.
- Our study utilizes advanced pre-processing techniques, including synthetic minority oversampling technique (SMOTE) and principal component analysis (PCA), to address imbalanced datasets and reduce dimensionality while preserving essential information. These techniques enhance the robustness and efficiency of our model, allowing for more accurate PD detection.
- The study adopts an ensemble learning approach, specifically AdaBoost, to construct a predictive model for PD detection. By combining multiple weak learners, AdaBoost improves the model’s performance and generalization ability, making it well suited for handling complex relationships in speech datasets. This novel application of ensemble learning techniques to PD detection represents a significant advancement in the field, offering a more effective and reliable diagnostic tool.
- In contrast to conventional diagnostic approaches that often overlook non-motor symptoms, this study explores the diagnostic potential of such symptoms. By incorporating non-motor manifestations, the model aims to enrich the diagnostic process for PD.
- By facilitating the detection of PD, this proposed study strives to mitigate the substantial economic burden associated with delayed diagnosis. Timely identification of the disease can potentially reduce both direct healthcare costs and indirect economic burdens on individuals and healthcare systems.
1.2. Motivations
2. Related Work
3. Materials and Methods
3.1. Retrieval of Dataset
- Time frequency fading (TFF) features: These features combine time and frequency components of voice signals.
- Mel frequency cepstral coefficients (MFCCs): MFCCs capture vital information about the spectral data of speech signals and are obtained using a series of steps like framing, windowing, and Fourier transform.
- Wavelet transform-based features (WTF): These features include pitch, formants, energy, and spectral density of speech signals obtained from wavelet transform functions.
- Vocal fold features (VFF): VFFs characterize the human vocal cord using certain parameters like jitter, shimmer, and fundamental frequency.
- Tremor waveform quality time (TWQT): It measures the quality of the oscillatory movement of voice signals.
- Entropy: Entropy provides information about the disorder in the signal and can help the model identify and adapt to noisy conditions, as in the presence of noise, the signal’s entropy may change [33].
- Detrended fluctuation analysis (DFA) features: To enhance the model’s resilience to short-term noise, DFA features, which are effective for evaluating long-term correlations in the signal, have been leveraged. This approach enables the model to prioritize the underlying structure of the speech signal.
- Quantile, skewness, and kurtosis features: Quantile, skewness, and kurtosis features capture deviations in the voice signals caused by noise, aiding the model in identifying and handling noisy segments.
- Signal-to-noise and noise-to-signal ratios: Integrating these features during model training allows the model to understand the proportion of noise in the signal and distinguish between the actual signal and surrounding noise, promoting robustness in the presence of varying noise levels [34].
3.2. Data Pre-Processing
3.2.1. Feature Weighting
- For the given minority class set X, the k-nearest neighbors of each sample a in X are determined by computing the Euclidean distance between a and all other samples in set X.
- The sampling rate N is determined based on the degree of class imbalance. For every sample in the minority class set X, N instances (a1, a2,…, an) are chosen randomly from their k-nearest neighbors to form the set X1.
- A new example is generated for each instance ak in X1 (where k = 1, 2, 3,…, N) using Equation (1).
3.2.2. Feature Scaling
3.3. Feature Extraction
Principal Component Analysis
- Standardizing the data: To ensure that each variable has the same influence in the analysis, it is important to standardize the data. This involves adjusting the data so that each variable has a mean of zero (0) and a variance of one (1). Without standardization, variables with larger values can have a disproportionate impact on the analysis, leading to biased results.
- Constructing the covariance matrix: The covariance matrix is an N×N matrix that displays the covariance between each pair of variables in the dataset. Two features can have positive covariance, i.e., they tend to vary in the same direction, or negative covariance, i.e., with an increase in one feature, the other feature tends to decrease.
- Decompose the covariance matrix: The covariance matrix is decomposed into Eigen vectors and Eigen values that represent the principal components and their magnitudes, respectively.
- Projecting the data onto a new sub-space: In the last step, the most important k eigenvectors are selected based on their respective eigenvalues (ones with the largest Eigen values are chosen), and those are used to transform the data and project it onto the new subspace.
3.4. Model Selection and Training
- Choosing an appropriate ML model: The process of selecting an ML model involves choosing one model from several potential candidates to learn the complexities in a dataset. In this study, an ensemble modeling approach has been adopted to process the speech dataset for its effectiveness in handling non-linear relationships often exhibited by the speech datasets. Ensemble methods are also recognized for their ability to mitigate the impact of irrelevant information (noise) in the training data, consequently enhancing the proposed model’s robustness against speech data affected by noise.
- Fitting the ML model: Fitting an ML model entails evaluating how effectively the model can extrapolate to comparable data to that for which it was trained. In this study, the dataset was partitioned into training and testing subsets with an 80:20 ratio. Post the train-test split, the ensemble model was fitted to the training subset, allowing it to capture the underlying intricacies of the data. Subsequently, the model’s predictive capability was assessed on the testing subset, which remained unseen during the training phase. Generalizing the model to unseen data is essential for ensuring its ability to make predictions.
Model Building
- Initialization: Initially, each sample in the training dataset is assigned an equal weight.
- Base learning algorithm: A base learning algorithm, typically a decision tree with a single level (also known as a decision stump), is trained on the weighted dataset. The decision stump aims to find the simplest rule that can classify the data.
- Weight update: The misclassified samples are assigned higher weights, while the correctly classified samples are assigned lower weights, emphasizing the importance of the former in subsequent iterations.
- Iterative process: AdaBoost repeats the process of training a new base learner on the weighted dataset and updating the sample weights iteratively. Each subsequent learner focuses more on the instances that were misclassified by the previous learners.
- Weighted voting: After all iterations, the final model combines the predictions of all base learners through a weighted voting mechanism, where the weight of each learner’s prediction is determined by its accuracy on the training data.
- Output: The final prediction is made by aggregating the weighted predictions of all base learners, resulting in a strong classifier capable of accurately classifying instances, even in the presence of noise or complex patterns.
Algorithm 1: AdaBoost for the detection of Parkinson’s Disease | |
; Base algorithm = L; Number of iterations = m. | |
Process: | |
[1] | do |
[2] | |
[3] | end for |
[4] | |
[5] | do |
[6] | |
[7] | set |
[8] | |
[9] | |
[10] | |
[11] | do |
[12] | |
[13] | end for |
[14] | end for |
[15] | Return L; |
[16] |
- Initializing weights: For each sample in the dataset, initially the weight is set to N, where N is the total number of instances in the dataset. Sample the dataset using weights.
- Calculating Gini impurity (GI) for each feature variable: GI for each node can be calculated using Equation (3).
- The GI for a feature can be computed by taking the weighted average of the impurities at each node. The feature variable with the lowest GI is used to create the first DS.
- Determine the amount of say for the newly created DS: This is performed by calculating the total error, which is equal to the sum of the weights of the misclassified samples, as shown in Equation (4).
- Calculate sample weights for the next DS: In this step, the sample weight for misclassified data points is increased and that of correctly classified instances is decreased by the use of Equations (5) and (6).
4. Results and Discussion
4.1. Performance Analysis
4.2. Comparative Analysis
5. Conclusions
Author Contributions
Funding
Data Availability Statement
Acknowledgments
Conflicts of Interest
References
- Parkinson, J. An essay on the shaking palsy. J. Neuropsychiatry Clin. Neurosci. 2002, 14, 223–236. [Google Scholar] [CrossRef] [PubMed]
- Lombardo, J.M.; Lopez, M.A.; Miron, F.; López, M.; León, M.; Arambarri, J.; Álvarez, D. MOBEEZE. Natural interaction technologies, virtual reality and artificial intelligence for gait disorders analysis. Int. J. Interact. Multimed. Artif. Intell. 2019, 5, 54. [Google Scholar]
- Adam, H.; Gopinath, S.C.; Md Arshad, M.K.; Adam, T.; Parmin, N.A.; Husein, I.; Hashim, U. An update on pathogenesis and clinical scenario for Parkinson’s disease: Diagnosis and treatment. 3 Biotech 2023, 13, 142. [Google Scholar] [CrossRef]
- Church, F.C. Treatment options for motor and non-motor symptoms of parkinson’s disease. Biomolecules 2021, 11, 612. [Google Scholar] [CrossRef] [PubMed]
- Park, Y.H.; Suh, J.H.; Kim, Y.W.; Kang, D.R.; Shin, J.; Yang, S.N.; Yoon, S.Y. Machine learning based risk prediction for Parkinson’s disease with nationwide health screening data. Sci. Rep. 2022, 12, 19499. [Google Scholar] [CrossRef] [PubMed]
- Saeed, F.; Al-Sarem, M.; Al-Mohaimeed, M.; Emara, A.; Boulila, W.; Alasli, M.; Ghabban, F. Enhancing Parkinson’s disease prediction using machine learning and feature selection methods. Comput. Mater. Contin. 2022, 71, 5639–5658. [Google Scholar] [CrossRef]
- Pramanik, A.; Sarker, A. Parkinson’s disease detection from voice and speech data using machine learning. In International Joint Conference on Advances in Computational Intelligence; Springer: Singapore, 2021. [Google Scholar]
- Cherubini, A.; Nistico, R.; Novellino, F.; Salsone, M.; Nigro, S.; Donzuso, G.; Quattrone, A. Magnetic resonance support vector machine discriminates essential tremor with rest tremor from tremor-dominant Parkinson disease. Mov. Disord. 2014, 29, 1216–1219. [Google Scholar] [CrossRef] [PubMed]
- Moro-Velazquez, L.; Garcia, J.A.G.; Arias-Londono, J.D.; Dehak, N.; Godino-Llorente, J.I. Advances in Parkinson’s disease detection and assessment using voice and speech: A review of the articulatory and phonatory aspects. Biomed. Signal Process. Control. 2021, 66, 102418. [Google Scholar] [CrossRef]
- Narendra, N.P.; Schuller, B.; Alku, P. The detection of Parkinson’s disease from speech using voice source information. IEEE/ACM Trans. Audio Speech Lang. Process. 2021, 29, 1925–1936. [Google Scholar] [CrossRef]
- Almahadin, G.; Lotfi, A.; Carthy, M.M.; Breedon, P. Task-oriented intelligent solution to measure Parkinson’s disease tremor severity. J. Healthc. Eng. 2021, 4, 9624386. [Google Scholar] [CrossRef]
- Iyer, A.; Kemp, A.; Rahmatallah, Y.; Pillai, L.; Glover, A.; Prior, F.; Larson-Prior, L.; Virmani, T. A machine learning method to process voice samples for identification of Parkinson’s disease. Sci. Rep. 2023, 13, 20615. [Google Scholar] [CrossRef]
- Rahman, S.; Hasan, M.; Sarkar, A.K.; Khan, F. Classification of Parkinson’s disease using speech signal with machine learning and deep learning approaches. Eur. J. Electr. Eng. Comput. Sci. 2023, 7, 20–27. [Google Scholar] [CrossRef]
- Alshammri, R.; Alharbi, G.; Alharbi, E.; Almubark, I. Machine learning approaches to identify Parkinson’s disease using voice signal features. Front. Artif. Intell. 2023, 6, 1084001. [Google Scholar] [CrossRef] [PubMed]
- Govindu, A.; Palwe, S.A. Early detection of Parkinson’s disease using machine learning. Preced. Comput. Sci. 2023, 218, 249–261. [Google Scholar] [CrossRef]
- Alalayah, K.M.; Senan, E.M.; Altam, H.F.; Ahmed, I.A.; Shatnawi, H.S.A. Automatic and early detection of Parkinson’s disease by analyzing acoustic signals using classification algorithms based on recursive feature elimination method. Diagnostics 2023, 13, 1924. [Google Scholar] [CrossRef] [PubMed]
- Mondol, S.R.; Kim, R.; Lee, S. Hybrid machine learning framework for multistage Parkinson’s disease classification using acoustic features of sustained korean vowels. Bioengineering 2023, 10, 984. [Google Scholar] [CrossRef]
- Lamba, R.; Gulati, T.; Alharbi, H.; Jain, A. A hybrid system for Parkinson’s disease diagnosis using machine learning techniques. Int. J. Speech Technol. 2022, 25, 583–593. [Google Scholar] [CrossRef]
- Das, R. A comparison of multiple classification methods for diagnosis of Parkinson disease. Expert Syst. Appl. 2010, 37, 1568–1572. [Google Scholar] [CrossRef]
- Rehman, A.; Saba, T.; Mujahid, M.; Alamri, F.; ElHakim, N. Parkinson’s disease detection using hybrid LSTM-GRU deep learning model. Electronics 2023, 12, 2856. [Google Scholar] [CrossRef]
- Lihore, U.K.; Dalal, S.; Faujdar, N.; Margala, M. Hybrid CNN-LSTM model with efficient hyperparameter tuning for prediction of Parkinson’s disease. Sci. Rep. 2023, 13, 14605. [Google Scholar] [CrossRef]
- Goyal, J.; Khandnor, P.; Aseri, T.C. A Hybrid Approach for Parkinson’s Disease diagnosis with resonance and time-frequency based features from speech signals. Expert Syst. Appl. 2021, 182, 115283. [Google Scholar] [CrossRef]
- Chowdhary, C.L.; Srivatsan, R. Non-invasive detection of Parkinson’s disease using deep learning. Int. J. Image Graph. Signal Process. 2022, 14, 38–46. [Google Scholar] [CrossRef]
- Asuroglu, T.; Ogul, H. A deep learning approach for parkinson’s disease severity assessment. Health Technol. 2022, 12, 943–953. [Google Scholar] [CrossRef]
- Lamba, R.; Gulati, T.; Jain, A. An intelligent system for Parkinson’s diagnosis using hybrid feature selection approach. Int. J. Softw. Innov. 2022, 10, 1–13. [Google Scholar] [CrossRef]
- Liu, X.; Li, W.; Liu, Z.; Du, F.; Zou, Q. A dual-branch model for diagnosis of Parkinson’s disease based on the independent and joint features of the left and right gait. Appl. Intell. 2021, 51, 7221–7232. [Google Scholar] [CrossRef]
- Shabu, S.J.; Sivapriya, V.; Refonaa, J.; Dhamodaran, S.; Poornima. Parkinson’s Disease Detection using Machine Learning Algorithm. In Proceedings of the 2023 8th International Conference on Communication and Electronics Systems (ICCES), Coimbatore, India, 1–3 June 2023; pp. 990–997. [Google Scholar] [CrossRef]
- Latifoğlu, F.; Penekli, S.; Orhanbulucu, F.; Chowdhury, M.E.H. A novel approach for Parkinson’s disease detection using Vold-Kalman order filtering and machine learning algorithms. Neural Comput. Appl. 2024, 36, 9297–9311. [Google Scholar] [CrossRef]
- Tran, C.; Shen, K.; Liu, K.; Ashok, A.; Ramirez-Zamora, A.; Chen, J.; Li, Y.; Fang, R. Deep learning predicts prevalent and incident Parkinson’s disease from UK Biobank fundus imaging. Sci. Rep. 2024, 14, 3637. [Google Scholar] [CrossRef]
- Fenza, G.; Gallo, M.; Loia, V.; Orciuoli, F.; Herrera-Viedma, E. Data set quality in machine learning: Consistency measure based on Group Decision Making. Appl. Soft Comput. 2021, 106, 107366. [Google Scholar] [CrossRef]
- Şakar, C.O.; Serbes, G.; Gunduz, A.; Tunc, H.C. A comparative analysis of speech signal processing algorithms for Parkinson’s disease classification and the use of the tunable Q-factor wavelet transform. Appl. Soft Comput. 2019, 74, 255–263. [Google Scholar] [CrossRef]
- Yin, S.; Liu, C.; Zhang, Z.; Lin, Y.; Wang, D.; Tejedor, J.; Zheng, T.F.; Li, Y. Noisy training for deep neural networks in speech recognition. EURASIP J. Audio Speech Music. Process. 2015, 2015, 2. [Google Scholar] [CrossRef]
- Toh, M.; Togneri, R.; Nordholm, S. Spectral entropy as speech features for speech recognition. Proc. PEECS 2005, 1, 92. [Google Scholar]
- Kumalija, E.; Nakamoto, Y. Performance evaluation of automatic speech recognition systems on integrated noise-network distorted speech. Front. Signal Process. 2022, 2, 999457. [Google Scholar] [CrossRef]
- Chawla, N.V.; Bowyer, K.; Lawrence, H.O.; Kegelmeyer, W.P. SMOTE: Synthetic minority over-sampling technique. J. Artif. Intell. Res. 2002, 16, 321–357. [Google Scholar] [CrossRef]
- Guyon, I.; Elisseeff, A. An introduction to variable and feature selection. J. Mach. Learn. Res. 2003, 3, 1157–1182. [Google Scholar]
- Hubel, D.H.; Wiesel, T.N. Receptive fields, binocular interaction and functional architecture in the cat’s visual cortex. J. Physiol. 1962, 160, 106–154. [Google Scholar] [CrossRef] [PubMed]
- Hotelling, H. Analysis of a complex of statistical variables into principal components. J. Educ. Psychol. 1933, 24, 498–520. [Google Scholar] [CrossRef]
- Dietterich, T.G. Ensemble methods in machine learning. In International Workshop on Multiple Classifier Systems; Springer: Berlin/Heidelberg, Germany, 2000. [Google Scholar]
- Freund, Y.; Schapire, R.E. Experiments with a new boosting algorithm. In Proceedings of the International Conference on Machine Learning, Bari, Italy, 3–6 July 1996. [Google Scholar]
- Bukhari, S.N.H.; Webber, J.; Mehbodniya, A. Decision tree based ensemble machine learning model for the prediction of Zika virus T-cell epitopes as potential vaccine candidates. Sci. Rep. 2022, 12, 7810. [Google Scholar] [CrossRef] [PubMed]
- Şentürk, Z.K. Early diagnosis of Parkinson’s disease using machine learning algorithms. Med. Hypotheses 2020, 138, 109603. [Google Scholar] [CrossRef] [PubMed]
- Bukhari, S.N.H.; Jain, A.; Haq, E.; Mehbodniya, A.; Webber, J. Ensemble machine learning model to predict SARS-CoV-2 t-cell epitopes as potential vaccine targets. Diagnostics 2021, 11, 1990. [Google Scholar] [CrossRef]
- Nour, M.; Şentürk, Ü.K.; Polat, K. Diagnosis and classification of Parkinson’s disease using ensemble learning and 1D-PDCovNN. Comput. Biol. Med. 2023, 161, 107031. [Google Scholar] [CrossRef]
- Anoruo, C.M.; Bukhari, S.N.H.; Nwofor, O.K. Modeling and spatial characterization of aerosols at Middle East AERONET stations. Theor. Appl. Climatol. 2023, 152, 617–625. [Google Scholar] [CrossRef]
- Masoodi, F.; Quasim, M.; Bukhari, S.N.H.; Dixit, S.; Alam, S. Applications of Machine Learning and Deep Learning on Biological Data; CRC Press: Boca Raton, FL, USA, 2023. [Google Scholar]
- Bukhari, S.N.H.; Masoodi, F.; Dar, M.A.; Wani, N.I.; Sajad, A.; Hussain, G. Prediction of erythemato-squamous diseases using machine learning. In Applications of Machine Learning and Deep Learning on Biological Data; CRC Press and Taylor & Francis: Boca Raton, FL, USA, 2023; pp. 87–96. [Google Scholar]
- Bukhari, S.N.H.; Jain, A.; Haq, E. A Novel Ensemble Machine Learning Model for Prediction of Zika Virus T-Cell Epitopes. In Proceedings of Data Analytics and Management; Gupta, D., Polkowski, Z., Khanna, A., Bhattacharyya, S., Castillo, O., Eds.; Springer: Singapore, 2022; pp. 275–292. [Google Scholar]
- Ola, V.; Puri, I.; Goswami, D.; Vibha, D.; Shukla, G.; Goyal, V.; Srivastava, A.; Behari, M. Annual cost of care of Parkinson’s Disease and its determinants in North India—A cost of illness study with patient perspective. Ann. Indian Acad. Neurol. 2022, 25, 660. [Google Scholar] [PubMed]
No. of DT | Acc | Precision | Recall | F1 Score | FNR |
---|---|---|---|---|---|
10 | 0.81 | 0.83 | 0.78 | 0.80 | 0.22 |
50 | 0.88 | 0.90 | 0.85 | 0.97 | 0.15 |
100 | 0.92 | 0.94 | 0.89 | 0.91 | 0.11 |
500 | 0.96 | 0.98 | 0.93 | 0.95 | 0.07 |
Tree Depth | Acc | Precision | Recall | F1 Score | FNR |
---|---|---|---|---|---|
3 | 0.93 | 0.95 | 0.90 | 0.92 | 0.10 |
5 | 0.94 | 0.96 | 0.91 | 0.93 | 0.09 |
7 | 0.96 | 0.98 | 0.93 | 0.95 | 0.07 |
9 | 0.96 | 0.98 | 0.93 | 0.95 | 0.07 |
LR | Acc | Precision | Recall | F1 Score | FNR |
---|---|---|---|---|---|
0.1 | 0.80 | 0.82 | 0.77 | 0.79 | 0.23 |
0.5 | 0.81 | 0.83 | 0.78 | 0.80 | 0.22 |
1.0 | 0.81 | 0.83 | 0.78 | 0.80 | 0.22 |
Performance Metrics | Values |
---|---|
Acc | 0.96 |
Precision | 0.98 |
Recall | 0.93 |
F1 score | 0.95 |
FNR | 0.07 |
AUC score | 0.99 |
Study | Feature Set | Pre-Processing Method | Classification Algorithm | Accuracy | AUC Score |
---|---|---|---|---|---|
Iyer et al., 2023 [12] | Acoustic signals, Spectrogram | Resultant wav files were filtered with floor and ceiling values of 75 decibels (dB) and 300 dB, respectively, for males and 100 dB and 600 dB for females, along with a scaling range of [−1, 1]. | Inception v3 | NA (Only AUC has been used) | 0.97 |
Rahman et al., 2023 [13] | MDVP, NHR, HNR, RPDE, DFA, PPE, etc. | Random over sampler, PCA | DNN2 | 95.41% | 0.96 |
Lihore et al., 2023 [21] | MDVP, NHR, HNR, RPDE, DFA, PPE, etc. | dynamic feature breakdown using CNN and LSTM | CNN, LSTM | 93.51% | 1.0 |
Mondol et al., 2023 [17] | Frequency, jitter shimmer, HNR, MFCC, etc. | SMOTE | MLP | 95.48% | 0.98 |
Lamba et al., 2023 [18] | MDVP, NHR, HNR, RPDE, DFA, PPE, etc. | Genetic algorithm | RF | 95.58% | 0.98 |
Proposed model | MFCC, TFF, NHR, HNR, DFA, WTF, VFF, etc. | SMOTE, PCA | AdaBoost | 96% | 0.98 |
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Bukhari, S.N.H.; Ogudo, K.A. Ensemble Machine Learning Approach for Parkinson’s Disease Detection Using Speech Signals. Mathematics 2024, 12, 1575. https://doi.org/10.3390/math12101575
Bukhari SNH, Ogudo KA. Ensemble Machine Learning Approach for Parkinson’s Disease Detection Using Speech Signals. Mathematics. 2024; 12(10):1575. https://doi.org/10.3390/math12101575
Chicago/Turabian StyleBukhari, Syed Nisar Hussain, and Kingsley A. Ogudo. 2024. "Ensemble Machine Learning Approach for Parkinson’s Disease Detection Using Speech Signals" Mathematics 12, no. 10: 1575. https://doi.org/10.3390/math12101575
APA StyleBukhari, S. N. H., & Ogudo, K. A. (2024). Ensemble Machine Learning Approach for Parkinson’s Disease Detection Using Speech Signals. Mathematics, 12(10), 1575. https://doi.org/10.3390/math12101575