Multiscale distribution entropy analysis of short epileptic EEG signals

Dae Hyeon Kim; Jin-Oh Park; Dae-Young Lee; Young-Seok Choi; Dae Hyeon Kim; Jin-Oh Park; Dae-Young Lee; Young-Seok Choi

doi:10.3934/mbe.2024245

Mathematical Biosciences and Engineering

2024, Volume 21, Issue 4: 5556-5576. doi: 10.3934/mbe.2024245

Previous Article Next Article

Research article Special Issues

Multiscale distribution entropy analysis of short epileptic EEG signals

Department of Electronics and Communications Engineering, Kwangwoon University, Seoul 01897, South Korea
† These authors contributed equally to the work.

Received: 10 February 2024 Revised: 22 March 2024 Accepted: 26 March 2024 Published: 02 April 2024

This paper proposes an information-theoretic measure for discriminating epileptic patterns in short-term electroencephalogram (EEG) recordings. Considering nonlinearity and nonstationarity in EEG signals, quantifying complexity has been preferred. To decipher abnormal epileptic EEGs, i.e., ictal and interictal EEGs, via short-term EEG recordings, a distribution entropy (DE) is used, motivated by its robustness on the signal length. In addition, to reflect the dynamic complexity inherent in EEGs, a multiscale entropy analysis is incorporated. Here, two multiscale distribution entropy (MDE) methods using the coarse-graining and moving-average procedures are presented. Using two popular epileptic EEG datasets, i.e., the Bonn and the Bern-Barcelona datasets, the performance of the proposed MDEs is verified. Experimental results show that the proposed MDEs are robust to the length of EEGs, thus reflecting complexity over multiple time scales. In addition, the proposed MDEs are consistent irrespective of the selection of short-term EEGs from the entire EEG recording. By evaluating the Man-Whitney U test and classification performance, the proposed MDEs can better discriminate epileptic EEGs than the existing methods. Moreover, the proposed MDE with the moving-average procedure performs marginally better than one with the coarse-graining. The experimental results suggest that the proposed MDEs are applicable to practical seizure detection applications.

Keywords:

Citation: Dae Hyeon Kim, Jin-Oh Park, Dae-Young Lee, Young-Seok Choi. Multiscale distribution entropy analysis of short epileptic EEG signals[J]. Mathematical Biosciences and Engineering, 2024, 21(4): 5556-5576. doi: 10.3934/mbe.2024245

Related Papers:

[1]	Abhishek Mallick, Atanu Mondal, Somnath Bhattacharjee, Arijit Roy . Application of nature inspired optimization algorithms in bioimpedance spectroscopy: simulation and experiment. AIMS Biophysics, 2023, 10(2): 132-172. doi: 10.3934/biophy.2023010
[2]	Arijit Roy, Abhishek Mallick, Surajit Das, Abhijit Aich . An experimental method of bioimpedance measurement and analysis for discriminating tissues of fruit or vegetable. AIMS Biophysics, 2020, 7(1): 41-53. doi: 10.3934/biophy.2020004
[3]	Svetlana Kashina, Andrea Monserrat del Rayo Cervantes-Guerrero, Francisco Miguel Vargas-Luna, Gonzalo Paez, Jose Marco Balleza-Ordaz . Tissue-specific bioimpedance changes induced by graphene oxide ex vivo: a step toward contrast media development. AIMS Biophysics, 2025, 12(1): 54-68. doi: 10.3934/biophy.2025005
[4]	Tomasz Rok, Artur Kacprzyk, Eugeniusz Rokita, Dorota Kantor, Grzegorz Tatoń . Quantitative assessment of thermal effects on the auricle region caused by mobile phones operating in different modes. AIMS Biophysics, 2024, 11(4): 427-444. doi: 10.3934/biophy.2024023
[5]	Takayuki Yoshida, Hiroyuki Kojima . Subcutaneous sustained-release drug delivery system for antibodies and proteins. AIMS Biophysics, 2025, 12(1): 69-100. doi: 10.3934/biophy.2025006
[6]	Riyaz A. Mir, Jeff Lovelace, Nicholas P. Schafer, Peter D. Simone, Admir Kellezi, Carol Kolar, Gaelle Spagnol, Paul L. Sorgen, Hamid Band, Vimla Band, Gloria E. O. Borgstahl . Biophysical characterization and modeling of human Ecdysoneless (ECD) protein supports a scaffolding function. AIMS Biophysics, 2016, 3(1): 195-210. doi: 10.3934/biophy.2016.1.195
[7]	Neelabh Datta . A review of molecular biology detection methods for human adenovirus. AIMS Biophysics, 2023, 10(1): 95-120. doi: 10.3934/biophy.2023008
[8]	Min Zhang, Shijun Lin, Wendi Xiao, Danhua Chen, Dongxia Yang, Youming Zhang . Applications of single-cell sequencing for human lung cancer: the progress and the future perspective. AIMS Biophysics, 2017, 4(2): 210-221. doi: 10.3934/biophy.2017.2.210
[9]	Mohamed A. Elblbesy . The refractive index of human blood measured at the visible spectral region by single-fiber reflectance spectroscopy. AIMS Biophysics, 2021, 8(1): 57-65. doi: 10.3934/biophy.2021004
[10]	Shivakumar Keerthikumar . A catalogue of human secreted proteins and its implications. AIMS Biophysics, 2016, 3(4): 563-570. doi: 10.3934/biophy.2016.4.563

Abstract

1. Introduction

The prime interest of measuring bioimpedance is to assess the physiological condition of the biological body under question. In the field of electrical bioimpedance measurement, the Cole impedance model [1] is widely used in the characterization of biological tissue. Due to the tremendous applicability of bioimpedance in medical sciences, bioimpedance measurement remains an attractive area of research in the past couple of decades. Research on biological bodies using Cole model includes: assessment of human body composition [2],[3], estimation of the thickness of human skin [4], monitoring of thoracic fluid content [4], characterization of breast and prostate cancer [5],[6], characterization of liver tissue [7], assessment of muscle injury and porcine spleen tissue [8],[9] etc. and many more which can be found in recently published review articles [10]–[13].

In Cole model, a lump tissue is represented by its equivalent circuit and is studied as a single dispersion material in the frequency domain. Cole model is basically a modified Debye model by the replacement of an ideal capacitor with a more general Constant Phase Element (CPE) and Cole model provides better fitting to the experimental data, especially for biological tissue. Therefore, for effective characterization of skin [14], we have considered Cole bioimpedance model in this work.

On the other hand, the effect of moisturizing cream on human skin is studied from various perspectives, such as effectiveness, improvement of skin barrier and its repair, maintaining the integrity and appearance of skin by acting as occlusive, emollients and humectants, each with its own mechanism of actions [15]. Basically, the function of a moisturizer is to improve hydration and water content of the stratum corneum. By moisturizing the skin, a cover of small skin fissures is formed which provides a soothing protection film and protects the skin from friction. Skin care products such as a moisturizing cream do not stay inactive on the upper layer of the skin, instead it penetrates the upper layer and reaches the underlying layer thereby influencing the structure of the skin and its function [16]. In the market, different kinds (difference in chemical composition) of skin moisturizers are available. Similarly, as there are differences among the variety of skin types, a certain structure of a given activity may be ideal for one skin but unsuitable for another [17]. Electrical susceptance-based measurement on human skin is also studied to obtain the effect of moisturizer cream [18]. These studies are clinical in nature and are neither cost effective nor time efficient. Therefore, bioimpedance appeals as a promising tool for characterizing the effect of moisturizing cream on human skin. To apply bioimpedance as a characterization tool for human skin, it is essential that it should be capable to distinguish treated and untreated types of skin. The method of bioimpedance with the discrimination capability of living tissue using Cole bioimpedance model is an established method [1] and recently this discrimination method for biological tissues is described by Roy et al. [19]. Similar procedure of discrimination capability is indeed necessary to study the bioimpedance of human skin and such a procedure is presented in this paper.

The effect of moisturizing cream on the human skin by means of non-invasive bioimpedance measurement is studied here. A novel electrode-pair is fabricated and an LCR meter is used to measure the bioimpedance as a function of frequency. The discrimination capability of our measurement process to distinguish various skins is verified through experiments. The bioimpedance is measured for before treatment, during treatment and post-treatment period of skin. From these measurements, the Cole parameters are extracted and the relaxation-time of the tissues under investigation is estimated. Moisturizing cream is applied on the skin of a particular area with the dosage of twice a day for 22 days and the measurement is conducted at regular intervals. After 22 days, the application of cream is stopped but the measurement is continued to see the changes in relaxation-time. In addition to this, measurement on skin at untreated area is also performed for comparison. Two different types of moisturizing cream with study duration of 47 days are considered. Skin of both male and female are considered in this study.

2. Theory and model

In this work, human skin is modeled according to Cole bioimpedance model. The model and its equivalent circuit are schematically shown in Figure 1. Basically, Cole model considers the tissue under investigation as a dielectric material and provides an equivalent electrical circuit of that dielectric material. To characterize a particular tissue using the Cole model, it is essential to determine the four parameters, namely: R₀, R_∞, C, and α (see Figure 1 for details). These parameters are not dimension independent. Therefore, each of these parameters is expected to vary depending on various parameters of the measurement method. For example, variation in R₀ is obvious if the length of the tissue (gap between the electrodes) changes. This simply shows that uniqueness in the measurement process is desirable. A high degree of repeatability in measuring these parameters is required. The relaxation-time (τ) is a fundamental property of dielectric material and the Cole model also enables us to estimate the relaxation-time in terms of these four parameters. When an identical process of measurement is imposed to estimate the four Cole parameters, a slight variation in fixing the electrode would result in variability in these parameters and these variabilities cancel out in such a way that τ remains invariant. This is indeed found true experimentally [19]. To extract R₀ and R_∞, the tissue impedance is measured by the LCR meter. The other two parameters (namely C and α) are determined by fitting the experimental data with the theoretical Cole impedance model [19].

The Cole bioimpedance is given by [19],

$Z = R_{\infty} + \frac{(R_{0} - R_{\infty})}{1 + {(jωτ)}^{α}}$ (1)

Let, $j^{α} = (j \sin \frac{απ}{2} + \cos \frac{απ}{2})$

Therefore, $Z = R_{\infty} + \frac{(R_{0} - R_{\infty})}{1 + ω^{α} τ^{α} (j \sin \frac{απ}{2} + \cos \frac{απ}{2})}$

Or, $Z = R_{\infty} + \frac{(R_{0} - R_{\infty})}{(1 + ω^{α} τ^{α} \cos \frac{απ}{2}) + j (ω^{α} τ^{α} \sin \frac{απ}{2})}$

Finally,

$| Z | = \sqrt{{(R_{\infty} + \frac{(R_{0} - R_{\infty}) (1 + ω^{α} τ^{α} \cos \frac{απ}{2})}{{(1 + ω^{α} τ^{α} \cos \frac{απ}{2})}^{2} + {(ω^{α} τ^{α} \sin \frac{απ}{2})}^{2}})}^{2} + {(\frac{(R_{0} - R_{\infty}) (ω^{α} τ^{α} \sin \frac{απ}{2})}{{(1 + ω^{α} τ^{α} \cos \frac{απ}{2})}^{2} + {(ω^{α} τ^{α} \sin \frac{απ}{2})}^{2}})}^{2}}$ (2)

The expression for the relaxation-time is given by,

$τ = ([(R_{0} - R_{\infty}) C]^{1 / α})$ (3)

Where, Z = Complex impedance

R₀ = Resistance at very low frequency

R_∞ = Resistance at very high frequency

α = A dimensionless exponent (0 < α < 1 for biological tissue)

τ = Relaxation-time

ω = Angular frequency

C = Capacitance of skin tissue

Experimentally |Z| is measured as a function of frequency in the range of 1 Hz to 5 kHz. The parameters R₀ and R_∞ are the measured impedance of the tissue under investigation at the lowest frequency (1 Hz here) and at the highest frequency (5 kHz) respectively [19]. The other two parameters (namely C and α) are determined by fitting the experimental data with the theoretical Cole impedance expressed by Equation 2 [19].

The dielectric property of human skin depends on various factors which include its physical and chemical characteristics, the composition of intra- and extra-cellular fluids [20], psychological condition and physiological changes occurring in the subject and the condition of health or disease of the skin. When suitable parameters are chosen, a great amount of information can be obtained for Dermatological and Cosmetological purposes from dielectrics or bioimpedance measurements.

Figure 1. Schematic of human skin and its equivalent electrical circuit known as Cole bioimpedance model..

DownLoad: Full-Size Img PowerPoint

The skin impedance variation is quite complex and depends on a variety of factors. First of all, Galvanic Skin Reflex (GSR) [21] plays a significant role in the determination of skin impedance. GSR is the sudden variations of skin impedance due to internal or psychological factors. Secondly, we have noticed that for a calmly seated person who suffers no pain or stress, the value of impedance remains fairly the same and can be assessed for physiological purposes. Nevertheless, statistical analysis of bioimpedance data is essential to understand the actual variation if any with certainty. In order to improve the effectiveness of measurement, both the separation between the electrodes and the electrode width (contact area) are to be kept small compared to the thickness of the stratum corneum [22].

3. Experimental protocol

3.1. Electrode setup

Choice of material for the electrodes and the shape of the electrodes are important factors in bioimpedance measurement, because the success of the investigation is highly dependent on these factors. Various electrode configurations are being used in the measurement of bioimpedance. Among the two-, three-, four- and five-electrode configurations, the two-electrode configuration is more attractive because of its simplicity and effectiveness [19],[23].

Human skin near the wrist is investigated in this study. A pair of electrodes is fabricated, then attached to a holder and then the holder is stitched to a glove. The electrode-pair system with the glove is shown in Figure 2. The electrode holder is made of a piece of PCB (2.3 cm × 1.0 cm × 0.1 cm in dimension) stitched to the gloves by using micro-screws. The screws (2.5 mm × 1.4 mm × 8.0 mm in dimension and are made of stainless steel) which are used to tighten the PCB to the gloves also serve as electrodes. The distance between the electrodes is 0.7 mm, which is much greater than 100 µm and because of this difference, more subtle changes are easily accessible in the deeper sections of the stratum corneum [24],[25]. When the hand is placed into the glove and its belt is strapped, the head of the screws are pressed on the skin near the wrist.

Figure 2. Schematic representation of the two-electrode circuit (left) and electrode setup used in bioimpedance measurement (right)..

DownLoad: Full-Size Img PowerPoint

The gloves have a pre-curved form fit design with palm grip for a stable fit on the hand. The glove comes with a cuff wrist strap with hook-and-loop wrist closure for a secure fit. The tail portion of the screws are connected to LCR meter through ordinary multi strand wire.

3.2. Instruments and materials

Electrical impedance is measured using a 1 MHz precision LCR meter (Model: 8101G, GWINSTEK). Before any measurements, the LCR meter is calibrated by the usual open- and short-circuit method. Calibration eliminates stray capacitance and series impedance from the testing fixture. Each time, before measurement, the electrode-pair is cleaned by medicated spirit and dried properly. The meter is set to 2 V ac and frequency sweep from 1 Hz to 5 kHz with a total of 200 steps. The meter is controlled by a PC using dedicated software through an RS-232C serial port. The measured data is transferred to the PC for further analysis. The entire bioimpedance measurement system is shown in Figure 3. After each measurement, the strap of the glove containing the electrode-pair is released and after one minute, the strap is tightened again to repeat the measurement. In this way, the measurement is repeated five times. There is a possibility of slight change in position of electrode-pair in repeating the measurement. This factor may cause a variation in the measured parameters. To reduce this effect, each time five measurement is conducted and using ANOVA, it is found that the variation of estimated relaxation-time due to this factor is negligible. Therefore, average relaxation-time from five replicas is considered for subsequent ANOVA procedures.

Two varieties of cream are used in this investigation. The Cream-1 used is a refreshing soft cream with jojoba oil and vitamin E for invigorated and smooth skin. According to the manufacturer, it provides substantial moisturization. It also contains denatured alcohol which has a drying effect on the skin. The Cream-2 used is a herbal cream. The manufacturer of this cream claims that the cream is made from multiple fruit extracts that act as quality elements to achieve smoother, brighter and well-protected skin. It is enriched with rare fruit and potent plant extracts that penetrate deep into the skin to provide enduring action. It results in fairer, flawless and radiant skin. It is free from harsh chemicals that generally does not have side effects. The ingredients of these two creams are listed in Table 1. Both the moisturizing creams are commercial products and are purchased from the local market.

Table 1. Data of the ingredients of the creams used in the experiment.

Cream Number	Ingredients
Cream-1	Water, Myristyl Alcohol, Mineral Oil, Glycerin, Butylene Glycol, Alcohol Denatured, Stearic Acid, Petrolatum, Myristyl Myristate, Glyceryl Stearate, Hydrogenated Coco Glycerides, Di methicone, Simmondsia Chinensis Seed Oil, Tocopheryl Acetate, Polyglyceryl 2 Caprate, Phenoxyethanol, Lanolin Alcohol, Fragrance, Carbomer, Sodium Hydroxide, Cetyl Phosphate, Methylparaben, Ethylparaben, Propylparaben, Butylparaben and Isobutylparaben.
Cream-2	Water, Almond Oil, Octyl Methoxy Cinnamate, Glyceryl Mono Hydroxy Stearate, Triticum Vulgare (Wheat Germ) oil, Light Liquid Paraffin, Butyl Methoxydibenzoylmethane, Arctostaphylos Uva-Ursi (Bearberry) Extract, Camellia Sinensis (Green Tea) Extract, Lactic Acid, Glycerin, Stearic Acid, Milk Enzymes, Sodium Lactate, Propylene Glycol, Titanium Dioxide, Ceto Stearyl Alcohol, Cetyl Alcohol, Triethanolamine, Triter Penoids, Phenyl Flavones (Mulberry) Extract, Epigallocathechin (Saxifraga) Extract, Vitis Vinifera (Grape) Extract, 2-Phenoxyethanol, Tocopheryl Acetate, Methyl Paraben, Alpha-Arbutin, Fragrance.

Remarks: All the above data are claimed by the manufacturer.

| Show Table

DownLoad: CSV

Figure 3. Schematic representation of the measurement setup..

DownLoad: Full-Size Img PowerPoint

A total of three volunteers are considered for this study, which includes two males and one female, ageing 22 (mean) years. In this report, we label the two male volunteers as M1 (weight = 81 kg; height = 180 cm) and M2 (weight = 70 kg; height = 177.8 cm), whereas the female volunteer is labelled as F (weight = 60.5 kg; height = 157.48 cm). It is to be noted here that this kind of experiment tends to divert from the desired outcome if the experiment is not conducted in a controlled environment [21]. As a preventive measure to this issue, we have conducted all the measurements in an air-conditioned room with fixed temperature (at 25 °C) and humidity (around 63%–64%).

The cream is applied twice a day (once at 9 am and again at 9 pm), every day for 22 days near the wrist of the right arm (treated wrist). Personal care products which include shampoo, soaps or any other beauty products are not used by the volunteers till the last day of this study. In case of Cream-1, bioimpedance is measured six times during the treatment period and seven times during the post-treatment period. The measurement matrix for Cream-1 is schematically shown in Figure 4.

Figure 4. Period of application and measurement matrix for Cream-1..

DownLoad: Full-Size Img PowerPoint

Similar to Cream-1, the effect of Cream-2 is also investigated. The study matrix for Cream-2 is shown in Figure 5. In order to understand the effect of the application of the moisturizing cream on the skin, whenever bioimpedance is measured for the treated wrist, measurement is also conducted for the untreated wrist.

Figure 5. Period of application and measurement matrix for Cream-2..

DownLoad: Full-Size Img PowerPoint

The impedance as a function of frequency is measured for the skin near the wrist in the range of 1 Hz-5 kHz. For a particular measurement, R₀ and R_∞ are taken along with the measured impedance at 1 Hz and 5 kHz respectively [19]. The other two parameters i.e. C and α are extracted from the fitting of experimental data.

The Cole impedance parameters are usually extracted by using the direct method which has its own limitations [26]. In fact, the lower and higher frequency gain measurements used for fitting affect the accuracy of the direct method. Nonlinear least-square fit (NLSF) is used to overcome the limitations of direct method [27]. In NLSF, every datapoint is used for optimization in order to extract the parameters for the best fit. In order to fit the experimental data with Equation 2, NLSF is performed using commercial software package MS-Excel 2016 version. From the curve fitting, the other two parameters i.e. C and α are determined. In this way, all the Cole parameters are experimentally measured. The relaxation-time, τ is then estimated using Equation 3. The entire measurement is then repeated four times and thus in this study, each measurement contains five repetitions. Analysis of variance (ANOVA) is performed on the average value of τ using the same software package. More specifically, two-factor ANOVA without replication is a data analysis tool which is extensively used to check if there is any significant difference in the parameters of a given set of data.

Typical measured impedance data as a function of frequency is fitted with Equation 2 and is shown in Figure 6. Here the curve shows the data for Day-0 (before treatment of Cream-1) to Day-22 (treatment period of Cream-1). This experimental data justifies the suitability of adopting the Cole impedance model to describe the frequency-dependent behavior of skin tissues.

In this study, the frequency dependent impedance is measured 250 times for each volunteer. Therefore, a total of 750 measurement is conducted. Each measurement contains 200 data points. Here, a large number of data is analyzed by using ANOVA to obtain experimental results.

Figure 6. Measured impedance is fitted as a function of frequency (for M2)..

DownLoad: Full-Size Img PowerPoint

4. Results and discussion

The moisturizing cream is applied on the upper layer of the skin known as stratum corneum. Stratum corneum is composed of dead cells called keratinocytes. The three main functions of the stratum corneum are to defend against dehydration, toxins and bacteria. The main function of a moisturizer is to prevent the skin from drying out and to hold moisture inside the skin. The electrical parameters of biological tissues [20],[22] are partially related to the water content in them. Thus, change in the water content in the tissue results in a little change in its electrical parameters. However, with the advancement of measurement equipment, little change in the electrical parameters that yields due to the change in moisture content in tissues can be captured. It is observed that as the moisture level increases, the cell size [11] increases correspondingly which implies that a longer path is to be traversed by the charged particles (ions), so there is a gross increase in impedance of the skin [28]. Conversely, the water content increases the relaxation-time τ.

The variation of relaxation-time as a function of time (in the unit of days) for all the three volunteers and for both the creams is shown in Figure 7. From this figure, it is observed that at Day-0 (before treatment), the relaxation-time is minimum, then it increases sharply during the application period of the cream and decreases during the period when the application of cream is stopped. In order to understand the effect of a cream, say Cream-1 on an individual volunteer, the finer details are shown in Figures 8, 9 and 10. In these figures, the mean of each dataset is used to construct the respective curves. Note that the variation of relaxation-time is very similar for all the three volunteers. The nature of variation of relaxation-time is found to be very similar to that of skin susceptance [18].

Figure 7. Variation of τ with the days of the experiment for all volunteers..

DownLoad: Full-Size Img PowerPoint

Figure 8. Boxplot of relaxation-time throughout the experiment (for Cream-1)..

DownLoad: Full-Size Img PowerPoint

Figure 9. Boxplot of relaxation-time throughout the experiment (for Cream-1)..

DownLoad: Full-Size Img PowerPoint

Figure 10. Boxplot of relaxation-time throughout the experiment (for Cream-1)..

DownLoad: Full-Size Img PowerPoint

In order to find the effectiveness of the creams, the variation of relaxation-time for treated and untreated skin for each moisturizing cream are shown in Figures 11, 12 and 13 for the three volunteers. It can be seen that Cream-2 is more effective than the Cream-1 as the relaxation-time is higher for Cream-2 for all the volunteers. Further, Cream-1 takes time more time than Cream-2 in the increment of relaxation-time, or in other words, Cream-2 is more effective in moisturization than Cream-1. After the treatment period, the relaxation-time drops and Cream-2 takes longer time than Cream-1 in the decrement of relaxation-time. Thus, Cream-2 is more effective as a moisturizer. Note that the relaxation-time for untreated skin remains almost constant. On treating the skin by a moisturizing cream, the relaxation-time of the skin can be increased to as high as 7 times.

Figure 11. Variation of τ with respect to time for Volunteer-M1..

DownLoad: Full-Size Img PowerPoint

Figure 12. Variation of τ with respect to time for Volunteer-F..

DownLoad: Full-Size Img PowerPoint

Figure 13. Variation of τ with respect to time for Volunteer-M2..

DownLoad: Full-Size Img PowerPoint

The ANOVA analysis (based on τ) before starting the treatment is shown in Table 2. From the Table 2, it is observed that the value of F_critical is greater than the value of F. Thus, we can conclude that the volunteers are indistinguishable and hence the skin of every volunteer cannot be identified statistically.

Table 2. ANOVA results before the treatment period of cream.

Case	F value	F_critical value	P value	Conclusion
Case-1 (M1 vs F)	0.13303787	7.70864742	0.73376864	Indistinguishable
Case-2 (F vs M2)	1.00507646	7.70864742	0.372814	Indistinguishable
Case-3 (M1 vs M2)	0.08096312	7.70864742	0.79011956	Indistinguishable

| Show Table

DownLoad: CSV

In order to understand the effect of moisturizing cream on the skin, ANOVA is performed extensively for the data obtained throughout the experiment. From ANOVA, we also can determine the gender distinguishability for a given cream with a certain level of confidence. The changes in human skin with the treatment of the above moisturizer is clearly visible as shown in Figures 11, 12 and 13. Table 3 shows the ANOVA results of Day-7 vs Day-22 for each volunteer. From Table 3, it can be seen that the skin of individual volunteer is distinguishable at a confidence level of 95% or more. Table 4 shows ANOVA considering gender as a factor. From this table, it can be seen that gender is distinguishable during the treatment period while the same was indistinguishable before the treatment of the skin (see Table 2). This fact clearly indicates that the degree of changes in the skin due to the moisturizing cream is different for male and female. Table 3 depicts the gender distinguishability during the treatment period of Cream-1. It shows that the Day-7 vs Day-22 skin for female volunteer is distinguishable at confidence level more than 95% while the same the for male volunteers is more than 99%. This fact indicates that the Cream-1 is more effective for male during Day-7 to Day-22. When the effectiveness of Cream-2 is tested in a similar manner, it is found that Cream-2 is more effective for female.

Table 3. ANOVA results for the treatment period of the Cream-1.

Period	Volunteer	F value	F_critical value	P value	Conclusion	Confidence Level
Application Period	Volunteer-M1 (Day-7 vs Day-22)	27.7221976	25.9735903	0.00623248	Distinguishable	99.3%
	Volunteer-F (Day-7 vs Day-22)	8.25093556	8.17090943	0.0453568	Distinguishable	95.4%
	Volunteer-M2 (Day-7 vs Day-22)	49.2440205	41.4072734	0.00217182	Distinguishable	99.7%

| Show Table

DownLoad: CSV

The ANOVA result for gender sensitivity is shown in Table 4. From the Table 4, it is observed that F value is greater than the value of F_critical. Therefore, the volunteers of different gender are distinguishable statistically and hence the effect of moisturizing cream is clearly visible. On the other hand, volunteers of same gender are not distinguishable in treatment period. These facts clearly indicate that the effect of treatment is gender sensitive.

Table 4. ANOVA results of the application period of the Cream-1 for different genders.

Case	F value	F_critical value	P value	Conclusion	Confidence Level
Case-1 (M1 vs F)	1.56176403	1.55824737	0.27953649	Distinguishable	72.0%
Case-2 (F vs M2)	4.20154543	4.1923271	0.1097193	Distinguishable	89.0%
Case-3 (M1 vs M2)	2.01251907	2.01834673	0.52844496	Indistinguishable	-

| Show Table

DownLoad: CSV

Post-treatment ANOVA is also performed and that for Cream-1 is shown in Table 5. From Table 5, it is observed that F value is greater than the value of F_critical and p value is less than the value of α (this α should not be confused with the Cole parameter designated by α) in statistical analysis. Thus, we can state that the volunteers are distinguishable. So, the effect of moisturizing cream on different gender in the post-treatment period too.

Table 5. ANOVA results of the Cream-1 after stopping its use.

Period	Volunteer	F value	F_critical value	P value	Conclusion	Level of Confidence
Decay Period	Volunteer-M1 (Day-24 vs Day-45)	8.99191252	8.99116941	0.03999511	Distinguishable	96%
	Volunteer-F (Day-24 vs Day-45)	10.5203056	10.427692	0.03156914	Distinguishable	96.8%
	Volunteer-M2 (Day-24 vs Day-45)	8.17493833	8.1790943	0.04596729	Distinguishable	95.4%

| Show Table

DownLoad: CSV

Among the four Cole parameters, interestingly, significant variation of skin impedance at low frequencies as an effect of moisturizing cream is noticed. The skin impedance at very low frequency (1 Hz in our study) increases rapidly [29] with the application of the cream on each successive day. The variation of the parameter R₀ with the period of the experiment is shown in Figure 14.

Figure 14. Variation of resistance at low frequency (R₀) during experiment period..

DownLoad: Full-Size Img PowerPoint

From Figure 14, it can be seen that the low frequency impedance increases monotonically, shows a peak and then decreases monotonically. The variation in R₀ is found to be as high as 500% in contrast to the variation of α and C (less than 1% and 40% respectively). This fact suggests that the DC resistance of skin can be a good indicator for characterizing the moisturizing study. A systemic investigation is needed to establish a DC resistance-based characterization.

5. Conclusion

The bioimpedance of human skin is measured and the effect of moisturizing cream on skin is studied in this paper. An efficient electrode-pair is fabricated and an experimental setup is demonstrated for measuring human skin impedance of a small area. Experimental study for a long period of time (47 days) is presented here. A large number of measurements is conducted and ANOVA is performed for statistical relevance. It is observed that relaxation-time increases due to the treatment by moisturizing cream and decreases during the post-treatment period. This is confirmed by measuring the relaxation-time of the untreated area of skin. The method demonstrated here is very efficient in comparing the effectiveness of various moisturizing creams. The gender distinguishability during the treatment of the skin with a moisturizer also can be assessed. In short, an entire experimental setup with measurement procedure and data extraction to study the effect of moisturizing cream on human skin (of a small area) is demonstrated.

References

[1]	A. K. Ngugi, C. Bottomley, I. Kleinschmidt, J. W. Sander, C. R. Newton, Estimation of the burden of active and life-time epilepsy: A meta-analytic approach, Epilepsia, 51 (2010), 883–890. https://doi.org/10.1111/j.1528-1167.2009.02481.x doi: 10.1111/j.1528-1167.2009.02481.x
[2]	W. A. Hauser, E. Beghi, First seizure definitions and worldwide incidence and mortality, Epilepsia, 49 (2008), 8–12.
[3]	R. S. Fisher, C. Acevedo, A. Arzimanoglou, A. Bogacz, J. H. Cross, C. E. Elger, et al., ILAE official report: A practical clinical definition of epilepsy, Epilepsia, 55 (2014), 475–482. https://doi.org/10.1111/epi.12550 doi: 10.1111/epi.12550
[4]	O. Devinsky, T. Spruill, D. Thurman, D. Friedman, Recognizing and preventing epilepsy-related mortality: A call for action, Neurology, 86 (2016), 779–786.
[5]	H. O. Lüders, I. Najm, D. Nair, P. Widdess-Walsh, W. Bingman, The epileptogenic zone: general principles, Epileptic Disord., 8 (2006), 1–2. https://doi.org/10.1684/j.1950-6945.2006.tb00152.x doi: 10.1684/j.1950-6945.2006.tb00152.x
[6]	Y. Paul, Various epileptic seizure detection techniques using biomedical signals: A review, Brain Inf., 5 (2018). https://doi.org/10.1186/s40708-018-0084-z doi: 10.1186/s40708-018-0084-z
[7]	S. Noachtar, J. Rémi, The role of EEG in epilepsy: A critical review, Epilepsy Behav., 15 (2009), 22–33. https://doi.org/10.1016/j.yebeh.2009.02.035 doi: 10.1016/j.yebeh.2009.02.035
[8]	U. Seneviratne, M. Cook, W. D'Souza, Brainwaves beyond diagnosis: Wider applications of electroencephalography in idiopathic generalized epilepsy, Epilepsia, 63 (2022), 22–41. https://doi.org/10.1111/epi.17119 doi: 10.1111/epi.17119
[9]	R. Cherian, E. G. Kanaga, Theoretical and methodological analysis of EEG based seizure detection and prediction: An exhaustive review, J. Neurosci. Methods, 369 (2022), 109483. https://doi.org/10.1016/j.jneumeth.2022.109483 doi: 10.1016/j.jneumeth.2022.109483
[10]	S. J. J. Jui, R. C. Deo, P. D. Barua, A. Devi, J. Soar, U. R. Acharya, Application of entropy for automated detection of neurological disorders with electroencephalogram signals: A review of the last decade (2012–2022), IEEE Access, 11 (2023), 71905–71924. https://doi.org/10.1109/ACCESS.2023.3294473 doi: 10.1109/ACCESS.2023.3294473
[11]	A. Ulate-Campos, F. Coughlin, M. Gaínza-Lein, I. S. Fernández, P. L. Pearl, T. Loddenkemper, Automated seizure detection systems and their effectiveness for each type of seizure, Seizure, 40 (2016), 88–101. https://doi.org/10.1016/j.seizure.2016.06.008 doi: 10.1016/j.seizure.2016.06.008
[12]	A. Shoeibi, N. Ghassemi, R. Alizadehsani, M. Rouhani, H. Hosseini-Nejad, A. Khosravi, et al., A comprehensive comparison of handcrafted features and convolutional autoencoders for epileptic seizures detection in EEG signals, Expert Syst. Appl., 163 (2021), 113788. https://doi.org/10.1109/ACCESS.2023.3294473 doi: 10.1109/ACCESS.2023.3294473
[13]	D. Y. Lee, Y. S. Choi, Multiscale distribution entropy analysis of short-term heart rate variability, Entropy, 20 (2018), 952. https://doi.org/10.3390/e20120952 doi: 10.3390/e20120952
[14]	V. Srinivasan, C. Eswaran, N. Sriraam, Approximate entropy-based epileptic EEG detection using artificial neural networks, IEEE Trans. Inf. Technol. Biomed., 11 (2007), 288–295. https://doi.org/10.1109/TITB.2006.884369 doi: 10.1109/TITB.2006.884369
[15]	H. Ocak, Automatic detection of epileptic seizures in EEG using discrete wavelet transform and approximate entropy, Expert Syst. Appl., 36 (2009), 2027–2036. https://doi.org/10.1016/j.eswa.2007.12.065 doi: 10.1016/j.eswa.2007.12.065
[16]	J. S. Richman, J. R. Moorman, Physiological time-series analysis using approximate entropy and sample entropy, Am. J. Physiol. Heart Circ. Physiol., 278 (2000), H2039–2049. https://doi.org/10.1152/ajpheart.2000.278.6.H2039 doi: 10.1152/ajpheart.2000.278.6.H2039
[17]	Y. Kumar, M. L. Dewal, R. S. Anand, Epileptic seizure detection using DWT based fuzzy approximate entropy and support vector machine, Neurocomputing, 133 (2014), 271–279. https://doi.org/10.1016/j.neucom.2013.11.009 doi: 10.1016/j.neucom.2013.11.009
[18]	Y. S. Choi, K. Hyun, J. Y. Choi, Assessing multiscale permutation entropy for short electroencephalogram recordings, Cluster Comput., 19 (2016), 2305–2314. https://doi.org/10.1007/s10586-016-0648-8 doi: 10.1007/s10586-016-0648-8
[19]	M. Costa, A. L. Goldberger, C. K. Peng, Multiscale entropy analysis of complex physiologic time series, Phys. Rev. Lett., 89 (2002), 068102. https://doi.org/10.1103/PhysRevLett.89.068102 doi: 10.1103/PhysRevLett.89.068102
[20]	S. D. Wu, C. W. Wu, K. Y. Lee, S. G. Lin, Modified multiscale entropy for short-term time series analysis, Phys. A, 392 (2013), 5865–5873. https://doi.org/10.1016/j.physa.2013.07.075 doi: 10.1016/j.physa.2013.07.075
[21]	Y. Zhang, P. Shang, Refined composite multiscale weighted-permutation entropy of financial time series, Phys. A, 496 (2018), 189–199. https://doi.org/10.1016/j.physa.2013.07.075 doi: 10.1016/j.physa.2013.07.075
[22]	Y. Li, B. Tang, S. Jiao, Q. Su, Snake optimization-based variable-step multiscale single threshold slope entropy for complexity analysis of signals, IEEE Trans. Instrum. Meas., 72 (2023), 6505313. https://doi.org/10.1109/TIM.2023.3317908 doi: 10.1109/TIM.2023.3317908
[23]	Y. Li, B. Tang, S. Jiao, Y. Zhou, Optimized multivariate multiscale slope entropy for nonlinear dynamic analysis of mechanical signals, Chaos Solitons Fractals, 179 (2024), 114436. https://doi.org/10.1016/j.chaos.2023.114436 doi: 10.1016/j.chaos.2023.114436
[24]	P. Li, C. Liu, K. Li, D. Zheng, C. Liu, Y. Hou, Assessing the complexity of short-term heartbeat interval series by distribution entropy, Med. Biol. Eng. Comput., 53 (2015), 77–87. https://doi.org/10.1016/j.chaos.2023.114436 doi: 10.1016/j.chaos.2023.114436
[25]	R. K. Udhayakumar, C. Karmakar, P. Li, M. Palaniswami, Effect of data length and bin numbers on distribution entropy (DistEn) measurement in analyzing healthy aging, in 2015 37th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), (2015), 7877–7880. https://doi.org/10.1109/EMBC.2015.7320218
[26]	P. Li, C. Karmakar, C. Yan, M. Palaniswami, C. Liu, Classification of 5-S Epileptic EEG recordings using distribution entropy and sample entropy, Front. Physiol., 7 (2016), 136. https://doi.org/10.3389/fphys.2016.00136 doi: 10.3389/fphys.2016.00136
[27]	D. Y. Lee, Y. S. Choi, Multiscale distribution entropy analysis of heart rate variability using differential inter-beat intervals, IEEE Access, 8 (2020), 48761–48773. https://doi.org/10.1109/ACCESS.2020.2978930 doi: 10.1109/ACCESS.2020.2978930
[28]	S. Supriya, S. Siuly, H. Wang, Y. Zhang, Epilepsy detection from eeg using complex network techniques: A Review, IEEE Rev. Biomed. Eng., 16 (2023), 292–306. https://doi.org/10.1109/RBME.2021.3055956 doi: 10.1109/RBME.2021.3055956
[29]	É. Lemoine, J. N. Briard, B. Rioux, R. Podbielski, B. Nauche, D. Toffa, et al., Computer-assisted analysis of routine electroencephalogram to identify hidden biomarkers of epilepsy: Protocol for a systematic review, BMJ Open, 13 (2023), e066932. https://doi.org/10.1136/bmjopen-2022-066932 doi: 10.1136/bmjopen-2022-066932
[30]	R. G. Andrzejak, K. Lehnertz, F. Mormann, C. Rieke, P. David, C. E. Elger, Indications of nonlinear deterministic and finite-dimensional structures in time series of brain electrical activity: Dependence on recording region and brain state, Phys. Rev. E, 64 (2001), 061907. https://doi.org/10.1103/PhysRevE.64.061907 doi: 10.1103/PhysRevE.64.061907
[31]	R. G. Andrzejak, K. Schindler, C. Rummel, Nonrandomness, nonlinear dependence, and nonstationarity of electroencephalographic recordings from epilepsy patients, Phys. Rev. E, 86 (2012), 046206. https://doi.org/10.1103/PhysRevE.86.046206 doi: 10.1103/PhysRevE.86.046206
[32]	T. Gautama, D. P. Mandic, M. M. Van Hulle, A differential entropy based method for determining the optimal embedding parameters of a signal, in 2003 IEEE International Conference on Acoustics, Speech, and Signal Processing, (2023).
[33]	M. R. Islam, X. Zhao, Y. Miao, H. Sugano, T. Tanaka, Epileptic seizure focus detection from interictal electroencephalogram: A survey, Cogn. Neurodyn., 17 (2023), 1–23. https://doi.org/10.1007/s11571-022-09816-z doi: 10.1007/s11571-022-09816-z
[34]	L. Ilias, D. Askounis, J. Psarras, Multimodal detection of epilepsy with deep neural networks, Expert Syst. Appl., 213 (2023), 119010. https://doi.org/10.1016/j.eswa.2022.119010 doi: 10.1016/j.eswa.2022.119010
[35]	M. Sunkara, S. R. Reeja, Tri-SeizureDualNet: A novel multimodal brain seizure detection using triple stream skipped feature extraction module entrenched dual parallel attention transformer, Biomed. Signal Process Control, 88 (2024), 105593. https://doi.org/10.1016/j.bspc.2023.105593 doi: 10.1016/j.bspc.2023.105593
[36]	F. Dong, Z. Yuan, D. Wu, L. Jiang, J. Liu, W. Hu, Novel seizure detection algorithm based on multi-dimension feature selection, Biomed. Signal Process Control, 84 (2023), 104747. https://doi.org/10.1016/j.bspc.2023.104747 doi: 10.1016/j.bspc.2023.104747
[37]	L. Jiang, J. He, H. Pan, D. Wu, T. Jiang, J. Liu, Seizure detection algorithm based on improved functional brain network structure feature extraction, Biomed. Signal Process Control, 79 (2023), 104053. https://doi.org/10.1016/j.bspc.2022.104053 doi: 10.1016/j.bspc.2022.104053
[38]	A. Einizade, S. Nasiri, M. Mozafari, S. H. Sardouie, G.D. Clifford, Explainable automated seizure detection using attentive deep multi-view networks, Biomed. Signal Process Control, 79 (2023), 104076. https://doi.org/10.1016/j.bspc.2022.104076 doi: 10.1016/j.bspc.2022.104076
[39]	X. Qiu, F. Yan, H. Liu, A difference attention ResNet-LSTM network for epileptic seizure detection using EEG signal, Biomed. Signal Process Control, 83 (2023), 104652. https://doi.org/10.1016/j.bspc.2023.104652 doi: 10.1016/j.bspc.2023.104652
[40]	D. Lee, B. Kim, T. Kim, I. Joe, J. Chong, K. Min, K. Jung, A ResNet-LSTM hybrid model for predicting epileptic seizures using a pretrained model with supervised contrastive learning, Sci. Rep., 14 (2024), 1319. https://doi.org/10.1038/s41598-023-43328-y doi: 10.1038/s41598-023-43328-y
[41]	T. K. K. Ho, N. Armanfard, Self-supervised learning for anomalous channel detection in EEG graphs: application to seizure analysis, in Proceedings of the AAAI Conference on Artificial Intelligence, (2023), 7866–7874. https://doi.org/10.1609/aaai.v37i7.25952

This article has been cited by:

1.	Lucian Pîslaru-Dănescu, George-Claudiu Zărnescu, Gabriela Telipan, Victor Stoica, Design and Manufacturing of Equipment for Investigation of Low Frequency Bioimpedance, 2022, 13, 2072-666X, 1858, 10.3390/mi13111858
2.	Jiawen Chen, Seyed Mohammad Mir, Maria R. Hudock, Meghan R. Pinezich, Panpan Chen, Matthew Bacchetta, Gordana Vunjak-Novakovic, Jinho Kim, Opto-electromechanical quantification of epithelial barrier function in injured and healthy airway tissues, 2023, 7, 2473-2877, 016104, 10.1063/5.0123127
3.	Nour Ammar, Cherif Ouni, Ahmed Yahia Kallel, Ahmed Fakhfakh, Nabil Derbel, Olfa Kanoun, 2024, Parameters Estimation of Cole-Cole Bio-Impedance Model with a Minimum Number of Frequencies, 979-8-3503-8090-3, 1, 10.1109/I2MTC60896.2024.10561098
4.	Ana Arché-Núñez, Peter Krebsbach, Bara Levit, Daniel Possti, Aaron Gerston, Thorsten Knoll, Thomas Velten, Chen Bar-Haim, Shani Oz, Shira Klorfeld-Auslender, Gerardo Hernandez-Sosa, Anat Mirelman, Yael Hanein, Bio-potential noise of dry printed electrodes: physiology versus the skin-electrode impedance, 2023, 44, 0967-3334, 095006, 10.1088/1361-6579/acf2e7
5.	Merwa Al‐Rasheed, Emily Lam, Mohammad Jambar, Jean Paul Ilogon, Sandra Gardner, Ladan Eskandarian, Amirali Toossi, Industry‐Scalable Reusable Textile Electrodes for Neurostimulation Applications, 2024, 2192-2640, 10.1002/adhm.202401642
6.	Pablo Dutra da Silva, Pedro Bertemes Filho, Switched CMOS current source compared to enhanced Howland circuit for bio-impedance applications, 2024, 15, 1891-5469, 145, 10.2478/joeb-2024-0017
7.	Nour Ammar, Cherif Ouni, Ahmed Yahia Kallel, Hanen Nouri, Bilel Ben Atitallah, Nabil Derbel, Olfa Kanoun, 2023, Cole-Cole Bio-Impedance Parameters Estimation From Sinewave Excitation Signal with a Minimum Number of Frequencies, 979-8-3503-5895-7, 22, 10.1109/IWIS61214.2023.10302771

Reader Comments

Your name:*

Email:*
© 2024 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)