Machine learning allows robust classification of lung neoplasm tissue using an electronic biopsy through minimally-invasive electrical impedance spectroscopy
Diagnosing lung diseases often requires tissue characterization, necessitating lung samples for an accurate final diagnosis. Current imaging methods can guide the diagnosis but are limited as they do not provide real-time guidance for collecting lung tissue samples. The studies found in literature , lack in consistency regarding the sensitivity of the biopsy in bronchoscopy procedures for lung neoplasm tissue sampling. The studies found state that the sensitivity for this procedure is between 60% and 88%. In addition, another factor for sensitivity variation would be the biopsy method implemented. Recently, new bronchoscopy techniques like radial probe endobronchial ultrasound and electromagnetic navigation bronchoscopy have been developed for real-time sampling characterization. However, their high cost limits availability in Interventional Pulmonology units, and their use remains suboptimal , . Table summarizes the pros and cons of the existing methods for tissue sampling. To enhance current diagnostic methods, affordable real-time guidance techniques are needed for tissue sampling. Combining electrical impedance spectroscopy (EIS) with artificial intelligence algorithms could enable electronic biopsies, allowing classification and identification of neoplasm tissue. This would ensure accurate tissue sampling and reduce the occurrence of negative and unnecessary biopsies. Bioimpedance (Z) is the tissue’s opposition to electrical current flow, which varies with frequency when alternating current is applied. This measurement across a wide range of frequencies is known as electrical impedance spectroscopy (EIS). Bioimpedance consists of resistance (R), related to extracellular and intracellular mediums, and reactance (Xc), related to cell membrane capacitance. From these, two parameters are derived: the impedance modulus (|Z|) defined as \documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym} 
				\usepackage{amsfonts} 
				\usepackage{amssymb} 
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\:\sqrt{{R}^{2}+X{c}^{2}}$$\end{document} and the phase angle (PA) defined as \documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym} 
				\usepackage{amsfonts} 
				\usepackage{amssymb} 
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\:{\text{tan}}^{-1}(\frac{Xc}{R})$$\end{document} . At low frequencies, current flows through the extracellular medium only, while at high frequencies, it penetrates cell membranes, flowing through both intra- and extracellular mediums – . Thus, changes in bioimpedance values are expected depending on tissue characteristics. Previous studies have explored using bioimpedance for lung tissue characterization – . Meroni et al. . developed an impedance meter for living tissues to test if electrical impedance spectroscopy was helpful in providing information about the structure and the properties of tissues. To validate the instrument, they performed ex-vivo impedance measurements in 3 different rabbits from 6 different tissue types finding statistical significance for the discrimination among the multiple tissues. Toso et al. . evaluated the distribution of the impedance vectors obtained at 50 kHz of frequency from 63 adult male patients with lung cancer and compared the results against 56 healthy subjects obtaining significant differences between cancer patients and control subjects due to a reduced reactance component. Baarends et al. . predicted total body water (TBW) and extracellular water (ECW) in patients with chronic obstructive pulmonary disease (COPD) using bioelectrical impedance spectroscopy (BIS). They concluded that predicted TBW using BIS was comparable to actual TBW, but presented no improvement of the prediction of TBW using bioelectrical impedance analysis (BIA) at 50 kHz. They also found that prediction of ECW had still limitations. These three previous studies focused on the application of impedance measurements for other applications different than our study. Regarding the application of electrical impedance spectroscopy for lung neoplasm differentiation two studies are found. Baghbani et al. . , constructed an electrical bioimpedance sensor with a biopsy forceps shape for measuring electrical conductivity of the tissue inside the body. They obtained and verified the relation between electrical conductivity of the tissue and measured electrical potential with COMSOL software. In addition, they designed and experimentally validated a prototype of the sensor. Furthermore, they measured the impedance of pulmonary tissues in three different samples of tissue founding that the sensor could be potentially beneficial to discriminate tumoral tissues from healthy ones in biopsy process. Baghbani et al. . introduced a method to localize in-depth pulmonary nodules intraoperatively by building a bioimpedance probe with four spherical electrodes. They collected in-vitro bioimpedance data of 286 lung tissue samples and applied principal component analysis (PCA) followed by classification algorithms (support vector machine (SVM), linear discriminant analysis (LDA), and K-nearest neighbors (KNN)) to localize the pulmonary nodules by the bioimpedance spectrum of the lung tissue. Apart from the last two studies introduced, to the extent of the author’s knowledge no studies have applied minimally invasive EIS for lung tissue differentiation, except those by their research group. The first one from Baghbani et al. . , obtained the measures from a biopsy sample while the second applied the classification algorithms from in-vitro samples. Measurements performed by our research group consist on measuring the lung samples directly, during a bronchoscopy process to help in sampling location before performing a biopsy. The first study of our research group was performed by Sanchez et al. . where a bioimpedance device was designed and validated for performing minimally-invasive bioimpedance measurements through bronchoscopy. Later, research focused on validating the best electrode configuration, implementing a calibration method to reduce data variability, and statistically differentiating between lung tissue types – . To the extent of the authors knowledge there are not current reports describing the application of Machine Learning classification algorithms to classify neoplasm lung tissue by using electrical impedance spectroscopy measurements for the implementation of an electronic biopsy measurement method to complement the actual guidance systems for a bronchoscopy procedure. The application of ML algorithms for clinical applications has raise importance in recent years. They present an opportunity to predict outcomes and develop new methods of diagnostic as well as improve prognostics . Current studies regarding the implementation of Artificial Intelligence in pneumology, different studies – performed in lung cancer applying classification algorithms, focused on the early diagnosis of the disease by using genetic data, computed tomography images and dosimetric features. The majority of the studies are focused on the disease prediction based on medical images. With the above mentioned, the aim of this study is to compare different Machine Learning algorithms for neoplasm lung tissue classification using minimally-invasive electrical impedance spectroscopy measurements for the implementation of a low-cost guidance method during bronchoscopy aiding in precise biopsy region detection. 

Participants Minimally invasive EIS measurements were performed between November 2021 and August 2022 in 102 patients (Age: 66 ± 14 year; Weight: 74.5 ± 17.2 kg; BMI: 26.8 ± 4.3 kgm-2) with a bronchoscopy prescribed at the “Hospital de la Santa Creu i Sant Pau” of Barcelona. A total number of 116 samples were obtained divided in 29 samples of lung tissue neoplasm and 87 samples of other lung tissue types (emphysema ( N = 23), healthy lung tissue ( N = 30), pneumonia ( N = 22) and fibrosis ( N = 12)). EIS measurements Minimally-invasive EIS using the 3-electrode method bioimpedance measurements are obtained through the injection of a multisine current signal (from 1 kHz to 1000 kHz) between a distal tetrapolar catheter electrode and a skin electrode during a bronchoscopy procedure. The voltage induced by the injected current is measured between the distal electrode and a second skin electrode. Impedance signal acquisition time was 12 s using a sample frequency of 60 spectra per second. The complete description of the impedance device as well as the calibration of the measurements is at Company-Se et al. . Measurements with abnormally large impedance values due to lose of contact between the catheter electrodes and the tissue samples were discarded for the analysis. Radiological images (CT or PET/CT) are taken in each patient before bronchoscopy following the diagnostic process. The catheter used to obtain the bioimpedance data is inserted through the working channel of the bronchoscope. Patients are placed in a supine position during the bioimpedance acquisition with the upper airways anaesthetized. Moreover, intravenous sedation is also provided. Biopsy was obtained to confirm the neoplasm diagnosis. Prior to the bioimpedance and biopsy acquisition saline solution is injected to clean and homogenize the tissue conditions among the different patients. Ethical statement Ethics approval was obtained from the Hospital de la Santa Creu i Sant Pau (CEIC-73/2020) according to principles of the Declaration of Helsinki for experiments with human beings. The patients/participants provided their written informed consent to participate in this study. Variables included in the study The data features used to apply the classification algorithms is constituted by the 12 s mean averaged spectra of the bioimpedance |Z| and PA obtained from 15 frequencies ranged from 15 kHz to 307 kHz. Given the 29 samples of neoplasm and the 87 samples of other lung tissue types, together with the 15 measures of |Z| and the 15 measures of PA, the total dimensionality of the original dataset is 116 samples * 30 features. Data preprocessing Synthetic Minority Oversampling Technique (SMOTE) is applied twice . First, SMOTE is applied to balance the dataset. From the 29 original neoplasm samples we created 58 synthetic neoplasm cases with a final first step sample size of 87 samples from lung neoplasm and 87 samples from other types of tissue. SMOTE is again applied to increase the sample size in a 50% more, thus obtaining from the 174 samples a final size of 262 samples (adding 44 synthetic samples per each class). Synthetic data creation has been performed using 5 Nearest Neighbors. After data augmentation, the total dimension of the dataset is 262 samples * 30 features. Data analysis The normality of all the features was assessed using the Kolmogorov-Smirnov test. Variables, non-parametric distributed, are described as median (interquartile range, IQR) and (minimum – maximum). Mann-Whitney U test, was used to assess non-normally distributed statistical significance between neoplasm tissue and the other group. The statistical significance was set as P < 0.05. Classification models Decision Tree, Support Vector Machines (SVM), Ensemble Method, K-Nearest Neighbors (KNN), Naïve Bayes and Discriminant Analysis classifiers were evaluated over the impedance modulus and phase angle of 15 frequencies distributed between 15 kHz and 307 kHz. Each dataset was normalized, as the range of values of the impedance modulus is different than the range of values of the impedance phase angle, in order to improve model performance and training efficiency. A grid search analysis was conducted to determine the optimal hyperparameter (HP) configuration for each model, employing a 5-fold cross-validation approach, which has been commonly utilized in prior medical research , . During every cross-validation, each dataset was partitioned into training (≈ 80% of the data) and test (≈ 20% of the data) and the model was trained and evaluated with each set of partitions. The HP optimized and the range of optimization are described in Table . The software MATLAB version: 23.2.0.2485118 (R2023b), Natick, Massachusetts: The MathWorks Inc.; 2023 has been used to implement the algorithms using the Classification Learner App. Classification models interpretability and performance assessment To evaluate the performance of each implemented model, a receiver operating characteristic (ROC) graph was generated. ROC graphs serve as valuable tools for analyzing the effectiveness of classification models by examining their true positive rate in relation to their false positive rate . The diagonal line in a ROC graph represents random guessing, and models positioned below this diagonal are considered less effective than random chance. An ideal classifier is depicted in the top-left corner of the graph, achieving a true positive rate of 1 and a false positive rate of 0. The area under the ROC curve (AUC) is calculated to quantitatively summarize the model’s classification performance. In addition to the ROC graph and AUC, the evaluation also includes precision, which is defined as the ratio of correctly identified positive instances (true positives) to the total predicted positive instances (true positives + false positives). Recall, representing the proportion of correctly identified positive instances (true positives) relative to the total actual positive instances (true positives + false negatives), is also considered. Furthermore, the F1-score, a metric that balances precision and recall to provide a comprehensive measure of classification performance, is computed as part of the assessment . Finally, the implementation time for each of the algorithms is also obtained. The overall process from data acquisition to model evaluation for classification to aid in a precise biopsy region detection is shown in the schematic diagram represented in Fig. . 

Minimally invasive EIS measurements were performed between November 2021 and August 2022 in 102 patients (Age: 66 ± 14 year; Weight: 74.5 ± 17.2 kg; BMI: 26.8 ± 4.3 kgm-2) with a bronchoscopy prescribed at the “Hospital de la Santa Creu i Sant Pau” of Barcelona. A total number of 116 samples were obtained divided in 29 samples of lung tissue neoplasm and 87 samples of other lung tissue types (emphysema ( N = 23), healthy lung tissue ( N = 30), pneumonia ( N = 22) and fibrosis ( N = 12)). 

Minimally-invasive EIS using the 3-electrode method bioimpedance measurements are obtained through the injection of a multisine current signal (from 1 kHz to 1000 kHz) between a distal tetrapolar catheter electrode and a skin electrode during a bronchoscopy procedure. The voltage induced by the injected current is measured between the distal electrode and a second skin electrode. Impedance signal acquisition time was 12 s using a sample frequency of 60 spectra per second. The complete description of the impedance device as well as the calibration of the measurements is at Company-Se et al. . Measurements with abnormally large impedance values due to lose of contact between the catheter electrodes and the tissue samples were discarded for the analysis. Radiological images (CT or PET/CT) are taken in each patient before bronchoscopy following the diagnostic process. The catheter used to obtain the bioimpedance data is inserted through the working channel of the bronchoscope. Patients are placed in a supine position during the bioimpedance acquisition with the upper airways anaesthetized. Moreover, intravenous sedation is also provided. Biopsy was obtained to confirm the neoplasm diagnosis. Prior to the bioimpedance and biopsy acquisition saline solution is injected to clean and homogenize the tissue conditions among the different patients. Ethical statement Ethics approval was obtained from the Hospital de la Santa Creu i Sant Pau (CEIC-73/2020) according to principles of the Declaration of Helsinki for experiments with human beings. The patients/participants provided their written informed consent to participate in this study. 

Ethics approval was obtained from the Hospital de la Santa Creu i Sant Pau (CEIC-73/2020) according to principles of the Declaration of Helsinki for experiments with human beings. The patients/participants provided their written informed consent to participate in this study. 

The data features used to apply the classification algorithms is constituted by the 12 s mean averaged spectra of the bioimpedance |Z| and PA obtained from 15 frequencies ranged from 15 kHz to 307 kHz. Given the 29 samples of neoplasm and the 87 samples of other lung tissue types, together with the 15 measures of |Z| and the 15 measures of PA, the total dimensionality of the original dataset is 116 samples * 30 features. 

Synthetic Minority Oversampling Technique (SMOTE) is applied twice . First, SMOTE is applied to balance the dataset. From the 29 original neoplasm samples we created 58 synthetic neoplasm cases with a final first step sample size of 87 samples from lung neoplasm and 87 samples from other types of tissue. SMOTE is again applied to increase the sample size in a 50% more, thus obtaining from the 174 samples a final size of 262 samples (adding 44 synthetic samples per each class). Synthetic data creation has been performed using 5 Nearest Neighbors. After data augmentation, the total dimension of the dataset is 262 samples * 30 features. 

The normality of all the features was assessed using the Kolmogorov-Smirnov test. Variables, non-parametric distributed, are described as median (interquartile range, IQR) and (minimum – maximum). Mann-Whitney U test, was used to assess non-normally distributed statistical significance between neoplasm tissue and the other group. The statistical significance was set as P < 0.05. 

Decision Tree, Support Vector Machines (SVM), Ensemble Method, K-Nearest Neighbors (KNN), Naïve Bayes and Discriminant Analysis classifiers were evaluated over the impedance modulus and phase angle of 15 frequencies distributed between 15 kHz and 307 kHz. Each dataset was normalized, as the range of values of the impedance modulus is different than the range of values of the impedance phase angle, in order to improve model performance and training efficiency. A grid search analysis was conducted to determine the optimal hyperparameter (HP) configuration for each model, employing a 5-fold cross-validation approach, which has been commonly utilized in prior medical research , . During every cross-validation, each dataset was partitioned into training (≈ 80% of the data) and test (≈ 20% of the data) and the model was trained and evaluated with each set of partitions. The HP optimized and the range of optimization are described in Table . The software MATLAB version: 23.2.0.2485118 (R2023b), Natick, Massachusetts: The MathWorks Inc.; 2023 has been used to implement the algorithms using the Classification Learner App. 

To evaluate the performance of each implemented model, a receiver operating characteristic (ROC) graph was generated. ROC graphs serve as valuable tools for analyzing the effectiveness of classification models by examining their true positive rate in relation to their false positive rate . The diagonal line in a ROC graph represents random guessing, and models positioned below this diagonal are considered less effective than random chance. An ideal classifier is depicted in the top-left corner of the graph, achieving a true positive rate of 1 and a false positive rate of 0. The area under the ROC curve (AUC) is calculated to quantitatively summarize the model’s classification performance. In addition to the ROC graph and AUC, the evaluation also includes precision, which is defined as the ratio of correctly identified positive instances (true positives) to the total predicted positive instances (true positives + false positives). Recall, representing the proportion of correctly identified positive instances (true positives) relative to the total actual positive instances (true positives + false negatives), is also considered. Furthermore, the F1-score, a metric that balances precision and recall to provide a comprehensive measure of classification performance, is computed as part of the assessment . Finally, the implementation time for each of the algorithms is also obtained. The overall process from data acquisition to model evaluation for classification to aid in a precise biopsy region detection is shown in the schematic diagram represented in Fig. . 

Impedance |Z| and PA mean impedance spectrum Figures and show the mean impedance spectra of the modulus and phase angle respectively for neoplasm (black), and other tissue (healthy and pathologic, blue) along all the frequencies (15 kHz to 307 kHz) for the original data (left) and for the dataset augmented after SMOTE application (right). The continuous line represents de mean while the dashed and pointed lines represent the ± SD. It represents the evolution of the impedance values (|Z| and PA) with respect to the frequency for both tissue groups. Analyzed variables Table show the descriptive information expressed as median (IQR) (minimum – maximum) of the variables included in the classification models for neoplasm group and for the group that included different types of lung tissue (Other tissue group). In addition, Table shows the statistic U of Mann Whitney and the statistical significance P. Classification models Table shows the optimal hyperparameter configuration for each of the classification models obtained from the 5-cross-validation grid search analysis, along with the corresponding test accuracy values and the standard deviation of the accuracy obtained from the cross-validation. Figure shows the confusion matrices, representing the relationship between actual and predicted classes, for each of the classification models implemented. Classification models performance assessment Figure show the ROC curves, that represent a visual representation of the models’ performance, obtained for each of the classification models implemented. In addition, the precision, recall and F1-score metrics are also specified. Table shows the implementation time for each of the algorithms implemented. 

Figures and show the mean impedance spectra of the modulus and phase angle respectively for neoplasm (black), and other tissue (healthy and pathologic, blue) along all the frequencies (15 kHz to 307 kHz) for the original data (left) and for the dataset augmented after SMOTE application (right). The continuous line represents de mean while the dashed and pointed lines represent the ± SD. It represents the evolution of the impedance values (|Z| and PA) with respect to the frequency for both tissue groups. 

Table show the descriptive information expressed as median (IQR) (minimum – maximum) of the variables included in the classification models for neoplasm group and for the group that included different types of lung tissue (Other tissue group). In addition, Table shows the statistic U of Mann Whitney and the statistical significance P. 

Table shows the optimal hyperparameter configuration for each of the classification models obtained from the 5-cross-validation grid search analysis, along with the corresponding test accuracy values and the standard deviation of the accuracy obtained from the cross-validation. Figure shows the confusion matrices, representing the relationship between actual and predicted classes, for each of the classification models implemented. 

Figure show the ROC curves, that represent a visual representation of the models’ performance, obtained for each of the classification models implemented. In addition, the precision, recall and F1-score metrics are also specified. Table shows the implementation time for each of the algorithms implemented. 

The application of Machine Learning classification algorithms are promising tools for the future of the medicine. They present an opportunity to help in the diagnosis of diseases for the tissue characterization. Different pathologies can affect to the respiratory system such as emphysema, neoplasm, fibrosis or pneumonia. Each of these disorders have their own anatomical and histological changes, thus differences in bioimpedance values are expected . Neoplasm is characterized by an increase in cell concentration as well as tissue vascularization which lowers the module impedance and increases the phase angle with respect to other types of tissue (Figs. and ). Regarding data augmentation, large datasets generally enhance classification accuracy, whereas small datasets are prone to overfitting. Data augmentation techniques can mitigate these challenges by generating extra samples for training models, thereby extracting more meaningful insights from limited data. Researchers leverage these methods to increase dataset diversity . Regarding the application of SMOTE, the new data is created from the average of 5 nearest original samples. According to the results of Figs. and , it can be concluded that the algorithm created synthetic data accordingly as the mean impedance spectrum both, for impedance modulus and phase angle does not change from the original data to the dataset with synthetic data. For the implementation of the classification algorithms, the impedance modulus and phase angle from 15 different frequencies from 15 kHz to 307 kHz have been used. Mann-Whitney U test can only test significance for each of the parameters alone. However, classification algorithms can use all the variable to find complex interconnections among variables to find a robust classification between neoplasm lung tissue or other type of lung tissue. As both, modulus and phase angle have obtained high statistical significance ( P < 0.001) for all the frequencies, and the U statistic does not differ much among variables (Table ), all the variables have been used for the classification models implementation. With respect to the classification models, the learning problem consist on finding a complex equation able to classify the samples as better as possible using training data by minimizing the classification error and increasing the accuracy for classification in never seen data. To do that, the different characteristic parameters of each of the algorithms have been optimized with the grid search method. It remains highly effective and ensures the identification of the most optimal combination of model parameters to obtain the higher metrics as possible for the classification task. Thus, the different models implemented are adequately optimized (Table ). In order to minimize as much as possible the overfitting of the algorithm, validation have been performed by using cross-validation with 5 folds . Regarding the accuracy obtained for each of the classification models implemented, all of them obtained an accuracy higher than 0.95 except for the Naïve Bayes, which has obtained an accuracy of 0.92 in addition of the highest SD of the accuracy (0.05). Furthermore, the HP of Table show that K-Nearest Neighbors used 1 neighbor to classify. This leads to the fact that although the bias of the 1-nearest neighbor estimate is often low, the variance is high . The different studies published regarding the use of classification algorithms for medical applications highlight the high performance of Support Vector Machines algorithm, although there has not been found consistence regarding which algorithm could fit best for medical applications – . Therefore, it is important to compare the performance of different algorithms to decide which suits best for the application under study. In our case, all of them, except Naïve Bayes and KNN algorithms could serve to detect neoplasm lung tissue through an electronic biopsy, with special focus on Support Vector Machines. Comparing the performance of the implemented algorithms with the study in the field that also applied classification algorithms , they found that SVM performed slightly better than LDA and KNN. However, the differences between the performance of the three algorithms are minor. In our case, the differences in performance between SVM and LDA are also minor. However, as the literature suggests, we focused especially on SVM since it seems to be the algorithm that performs the best on medical data. The findings in the literature lead to the conclusion that the application of already created classification algorithms are useful and safe for new medical applications or to complement already existing medical approaches. To fully evaluate the performance of the algorithms implemented, the classification report for all the algorithms have been obtained (Fig. ). We have plotted the ROC curve, used to evaluate the performance of a binary classification method for diagnostic, together with the AUC (Area Under Curve), which is a measure used to evaluate the accuracy of the test . All the classification methods reached an AUC higher than 95% which means the test accuracy is excellent. In order to fully assess the performance of the classification tests, the precision, recall and F1-score have also been calculated. While precision takes more importance when false positives are wanted to be avoided, the recall takes importance when false negative are not desired . F1-score offers a trade-off between precision and recall. For our clinical problem, false positives are not desired, as the biopsies of other tissue rather than neoplasm are not of clinical interest. On the other hand, false negative is also not desired, as we want to accurately localize neoplasm tissue to help in sampling location during bronchoscopy. Therefore, we focus on the F1-score to evaluate the performance of the classification methods implemented. Given that the classification method is intended for detection of tissue types during bronchoscopy, the duration of bioimpedance data acquisition and processing is critical, with shorter execution times being a significant advantage. The implementation time of the various classification methods employed are described in Table with shorter times in Decision Tree, SVM and Discriminant Analysis algorithms. According to Fig. , except for Naïve Bayes algorithm, all the algorithms obtained an F1-score higher than 0.95, which, together with the accuracy displayed in Table and the execution times shown in Table , leads to the conclusion that Decision Tree, SVM and Discriminant Analysis are suitable for our clinical application. In future studies the performance of Decision Tree, Discriminant Analysis and, specially, SVM, will be tested in real-time. In addition, by using an electromagnetic navigation bronchoscopy the usefulness of the classification algorithms to detect peripheral nodules will be evaluated. Clinical significance This study demonstrated the usefulness of machine learning algorithms for detecting neoplasm lung tissue during a bronchoscopy by performing an electronic biopsy and measuring the bioimpedance of the tissue. Indeed, thanks to machine learning, healthcare professionals can take advantage of tissue electrical properties variations which may remain unmeasured by actual navigation systems. In clinical practice, this minimally-invasive sampling localizer could be employed in the interventional pulmonology unit of hospitals for accurate biopsy sampling location during bronchoscopy enabling the decrease the negative biopsies due to sampling errors. Limitations Lung nodules, according to their location, can be classified as central nodules and peripheric nodules. During bronchoscopy, patients receive sedation, not anesthesia which implies that movement of the patients are frequent. In order to ensure and minimize the risk of not piercing the pleura, this study only included central nodules, which are analyzed with a conventional bronchoscopy. 

This study demonstrated the usefulness of machine learning algorithms for detecting neoplasm lung tissue during a bronchoscopy by performing an electronic biopsy and measuring the bioimpedance of the tissue. Indeed, thanks to machine learning, healthcare professionals can take advantage of tissue electrical properties variations which may remain unmeasured by actual navigation systems. In clinical practice, this minimally-invasive sampling localizer could be employed in the interventional pulmonology unit of hospitals for accurate biopsy sampling location during bronchoscopy enabling the decrease the negative biopsies due to sampling errors. 

Lung nodules, according to their location, can be classified as central nodules and peripheric nodules. During bronchoscopy, patients receive sedation, not anesthesia which implies that movement of the patients are frequent. In order to ensure and minimize the risk of not piercing the pleura, this study only included central nodules, which are analyzed with a conventional bronchoscopy. 

The results obtained after the comparison among all the algorithms implemented for neoplasm lung tissue classification using minimally-invasive electrical impedance spectroscopy measurements show that Decision Tree, Discriminant Analysis and Support Vector Machines algorithms, with special emphasis in the last one, are suitable for the implementation of a low-cost guidance method during bronchoscopy and that a new tool could be designed as new guidance tool. 

