Next Article in Journal
Comparison of Non-Uniform Image Quality Caused by Anode Heel Effect between Two Digital Radiographic Systems Using a Circular Step-Wedge Phantom and Mutual Information
Previous Article in Journal
Numerical Study on Heat-Transfer Characteristics of Convection Melting in Metal Foam under Sinusoidal Temperature Boundary Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Heart Failure Evolution Model Based on Anomalous Diffusion Theory

by
Andrzej Augustyn Walczak
Faculty of Cybernetics, Military University of Technology, Gen. Kaliskiego St. 2, 00-908 Warsaw, Poland
Entropy 2022, 24(12), 1780; https://doi.org/10.3390/e24121780
Submission received: 13 September 2022 / Revised: 24 November 2022 / Accepted: 30 November 2022 / Published: 6 December 2022

Abstract

:
The unexpectable variations of the diagnosed disease symptoms are quite often observed during medical diagnosis. In stochastics, such behavior is called “grey swan” or “black swan” as synonyms of sudden, unpredictable change. Evolution of the disease’s symptoms is usually described by means of Markov processes, where dependency on process history is neglected. The common expectation is that such processes are Gaussian. It is demonstrated here that medical observation can be described as a Markov process and is non-Gaussian. Presented non-Gaussian processes have “fat tail” probability density distribution (pdf). “Fat tail” permits a slight change of probability density distribution and triggers an unexpectable big variation of the diagnosed parameter. Such “fat tail” solution is delivered by the anomalous diffusion model applied here to describe disease evolution and to explain the possible presence of “swans” mentioned above. The proposed model has been obtained as solution of the Fractal Fokker–Planck equation (FFPE). The paper shows a comparison of the results of the theoretical model of anomalous diffusion with experimental results of clinical studies using bioimpedance measurements in cardiology. This allows us to consider the practical usefulness of the proposed solutions.

1. Introduction

In the history of a disease, there is often talk of an unpredictable reaction of the patient’s body. In the language of statistics, we would expect a precise indication of whether to consider sudden changes in the evolution of a disease as the equivalent of a “black swan” or a “grey swan” [1]. In terms of disease evolution, the observed change of disease symptoms is assumed. The basic classification into “grey and black” swans is based on two elements—the predictability of their occurrence and the consequences of their occurrence. The effects of the occurrence of a “grey swan” are known and are assumed to be predictable with an unknown value of variation. The consequences of a “black swan” occurrence are classified as catastrophic and unpredictable [1].
What each of these cases has in common is that the variation of the quantity characterizing the variation of the observed parameter is a power function of time
x 2 ( t ) = K α t α
or, in Hurst exponent notation
x 2 ( t ) = K H t 2 H
with exponent α > 0 [2].
Assuming that H is not equal to 1/2, it is interpreted in the literature as the presence of long-time correlations as in Madelbrot work [3]. When H is not equal to 1/2, but there are no autocorrelations, then it is a Levy process with infinite variation [4].
More generally, this is captured by the density of the probability distribution characterized by so-called “fat tails”, which means that even a slight change in the value of the probability density function corresponds to a significant change in the variation of the observed variable [1]. Thus, the grey and black swans are governed by a stochastic process, in which a small change in the value of the probability density function is associated with an unexpectedly large variance of the parameter because the probability distribution of the observed parameter has a “fat tail,” that is, it is not described by a Gaussian process.
Medical cases, due to their relative predictability of effects, can be classified (if they occur) as “grey swans.” However, it remains an open question whether stochastic processes that characterize the behavior of diagnosed medical parameters can have “fat tails.” If so, the occurrence of “grey swans” would be expected and frequent or even constantly present.
To study the nature of the stochastic process governing changes in medical parameters by observing their variance, one needs to have a large population of observations recorded under the same clinical conditions [5]. This is an extremely long and expensive process even with the current state of IT (meaning computer science) solutions in medicine.
We have proposed in earlier publications a different analysis using heart failure cases as an example [6,7]. Since we do not have a very large population of measurements of the selected parameter to be diagnosed (a few thousand measurements made under the same conditions of clinical observation of the patient, which guarantees their statistical utility but is a relatively small population), we checked whether, in a smaller population, the time series of observations is or is not well approximated by a Markov process. A Markov process, if it approximates well the results of a clinical observation, can indicate that the recorded observable is subject to a Gaussian stochastic process. This is a very rough estimate. Bassaler and co-authors showed that even when process variation would indicate the existence of “fat tails” in the probability density distribution, the processes can be, and most often are, Markov processes. This means that there is no need to consider the effect of history on the current value of the recorded observable [8]. The aforementioned authors proved that the determination of the value of the Hurst exponent, without further demonstration that the increments of the recorded observable are stationary or that the observed dynamics of the process directly depend on the history, is not a basis for claiming the existence of autocorrelation in the process. The admissibility of a Markov process with a Hurst exponent not equal to ½ implies the existence of non-stationary increments in the Markov process. This would mean that a good fit of the experimental results to the Markov process does not rule out non-stationary changes in the process and does not rule out that the observed process is nevertheless a process with “fat tails.”
Our aim is to find out if observed processes in disease evolution have “fat tails” or, in other words, if the Hurst exponent for registered observations is different from H = 1/2. This is because our experiment was conducted in the regime of exacting, multicenter clinical observations.
In the following, we propose to extend the analysis of disease evolution based on the fractal F–P equation and its analytical solutions. Fitting to experimental data allows us to verify whether the non-Gaussian process can be an effective model of disease evolution. The analysis presented here consists of experimental and theoretical parts. In the experimental part of the project, 605 patients with heart failure were recruited. The study was a multi-site, randomized, controlled group trial conducted in multiple centers in parallel (ClinicalTrials.gov Identifier: NCT03476590). The study included hemodynamic parameters measured by bioimpedance techniques and electrocardiography in parallel. The same experimental results were used in earlier work, which allows a common reference to the relationship between theoretical models and clinical observations.

2. Materials and Methods

2.1. Experimental Data

The experimental data used in this study were sourced from the AMULET research project [9,10,11,12], under which a multicenter, prospective, randomized, open-label and controlled, parallel group trial was conducted (ClinicalTrials.gov Identifier: NCT03476590) at nine locations throughout Poland. In total, 605 patients with heart failure were recruited. In our study, we examined electrocardiogram (ECG—electrocardiogram) and impedance (ICG— impedance cardiogram) curves recorded with the use of an ICG device (Cardioscreen 2000, Medis, Illmenau, Germany). This non-invasive diagnostic method allows one to collect a set of specific hemodynamic parameters, such as heart rate (HR—heart rate), stroke volume (SV—stroke volume), and thoracic fluid content (TFC—thoraric fluid content). In our analysis, TFC has been used as an illustrative example. The value of TFC is the inverse of chest impedance with unit (1/Ohm) and strongly depends on fluid content of the chest. The TFC is basic information used in the diagnosis of heart failure; because of fluid overload in the pathogenesis of HF exacerbation, early detection of fluid retention is of key importance in preventing emergency admissions. For such reason, our analysis took TFC into account as the most important parameter. TFC values were measured during a clinical examination of 605 patients, performed in a relaxed, seated position. The total number of measurements recorded here is 2860, with the number of individual patient observations differing for each patient within the measurement set. The results registered are illustrated in Figure 1. At least half of the patients were investigated several times, with the observation period lasting for up to 12 months per patient. All intervals of the investigated TFC value have been divided arbitrarily into 15 discrete value slots, as shown in Figure 1. The assumed number of slots must be greater than five to allow statistical analysis (especially χ2 tests) but not too high to avoid complexity of calculations. With computer calculation, we estimate that the amount of 15 slots allows for an acceptable accuracy of calculations. The time series was registered as follows: Due to the irregular flow of the measured data per patient, we adopted a registration period that was 31 days long for each slot. Each slot includes hits counted during measurements registered over a five-month period. Longer registration does not carry new results in TFC distribution. Hits registered during each single measurement period are illustrated in Figure 1.
One can observe that, for the increasing ∆t, the total number of TFC hits registered tends to have the same stationary distribution of TFC, as shown in Figure 1. We assumed that the period of 5 months during which the measurements were made was sufficient for proper estimation of stationary distribution.

2.2. Theoretical Models

As mentioned, two theoretical models—the Markov process [6] and the analytical solution of the Fokker–Planck equation (FPE) of diffusion containing a drift component for modeling the presence of therapy [7]—were previously studied in parallel.
The results described in papers [6,7] indicated good agreement between experimental observations and theoretical models, confirming to some extent that the observed processes can be considered Gaussian processes. In the following, we examine anomalous diffusion to describe the distribution of observed hemodynamic parameters to show that the process may nevertheless have “fat tails.”
Anomalous diffusion is found in a wide variety of such systems, with non-linear growth of the mean squared displacement, observed over the course of time, being its hallmark. The above means that process x(t) complies with the rule, which is common in a diverse number of complex systems [13,14,15,16,17,18,19,20,21,22,23].
In order to obtain a more general and adequate class of FPE solutions, it is necessary to use a fractional FP equation (FFPE) to describe a complex system [24,25,26,27,28,29,30,31,32]. This is particularly true in the context of continuous time random walk (CTRW) and anomalous diffusion. The procedure requires the inclusion of a non-Debye relaxation in the process of changing from one slot of the measured value to another. Debye relaxation fulfils the following rule [33,34,35].
d f ( t ) d t = 1 τ   f ( t )   , t > 0
The following is an integral form of (3):
f ( t ) f 0 = 1 τ d 1 f ( t ) d t 1 = 1 τ 0 t d t   f ( t   )
Generalization from an integer to a rational derivative form requires the following transformation:
f ( t ) f 0 = 1 τ α d α f ( t ) d t α = { 1 Γ E u l e r   ( α ) 0 t d t f ( t ) ( t t ) ,   α > 0 f ( t ) ,       α = 0
In the case of non-Debye processes, the Fractal Fokker–Planck Equation (FFPE) is of the form:
P ( x , t ) t = D t 1 α K α 2 x 2   P ( x , t )
In (6), the Riemann Liouville fractional differintegration is of the same form as in the non-Debye diffusion model above:
D t 1 α = t D t α 1 Γ ( α ) t 0 t d t 1 ( t t ) 1 α
Finally, we should write:
α f ( t ) t α = 1 Γ ( α n ) n t n 0 t d t f ( t ) ( t t ) ( 1 + α n )
With n −1 ≤ α < n, and the propagator (Green function) of (6), we obtain the following form:
P α , K α ( x , t , x 0 ) = 1 2 Γ ( α + 1 ) K α 2 f β , γ ( t ,   | x x 0 | Γ ( α + 1 ) K α 2 , )
f β , γ ( t ,   g ) = 1 t 1 + β k = 0 1 k !   Γ ( k γ β ) ( g t γ ) k
β = α 2 1   , γ = α 2 , g = | x x 0 | Γ ( α + 1 ) K α 2
The solution of FFPE (9) allows value classes of functions P(x, t) that exceed 1 and that may also be lower than 0. So, for P(x, t) to be interpreted as a probability density function (pdf), one must adopt specific constraints for α and Kα to ensure proper pdf interpretation. The constraint applies just to solution classes in Equation (9). Now, one may simulate the obtained solution (9–11), which must obey the predefined constraint 0 ≤ Pα, (x,t) ≤ 1 for a set of pairs {α, Kα}. The general constraint of 0 < α < 2 applies, as a rule, for fractional FFPE. The condition that P(x, t) is a decreased function of both arguments also imposes a constraint on the permissible pairs {α, Kα}. Figure 2 presents an example of the simulated Pα, (x,t) with arbitrary values of α and Kα.
We use the Green’s function of the FPE Equation (9) to calculate the time evolution of the probability density distribution for TFC starting from the selected stage of the stationary distribution obtained during the recording of the results and shown in Figure 1.
P α ,   K α ( x , t ) = 0 P α ,   K α ( x , t ,   x )   T F C   ( x ) d x

3. Results

The results of the fit between the theoretical model and clinical observations are presented in Figure 3.
Values of α and Kα placed inside Figure 3 have been calculated by means of root mean square error (RMSE) between measured experimental data and the theoretical “fat tail” distribution presented in Equation (12), closed for a time of observation as long as 5 months. One can see that the Hurst exponent is equal to α/2 = H = 0.235.

4. Discussion

The observed experimental results may depend on the ongoing therapy that regulates TFC. In other words, subjecting a patient to therapy does not rule out the existence of a low or very low dependence of the observed parameters on their history, that is, it does not rule out a possible temporal correlation. However, the effect of history on the observation performed during therapy, if any, appears to be negligibly small. This is indicated by the low value of the Hurst exponent (H = 0.235).
Importantly, it is both the indicated previous work based on Markov process models and the present model based on anomalous diffusion that effectively approximate the course of disease evolution. Thus, we do not settle on which process—Gaussian or non-Gaussian—is a better approximation. In view of the work of Besseler and co-authors as well as results presented here, it should be assumed that the process of disease evolution is a process with “fat tails”, but the use of Markov models is perfectly acceptable and effective, and the existing temporal correlation in the observed experiment is weak. This may be caused by the therapy conducted during the observation of patients, which significantly affects the temporal correlations in the evolution of the disease.
One of the peculiarities of observing disease symptoms is that we do not analyze the evolution of the disease as its effect on the patient’s body. There are always three elements: the therapy used, the patient’s body responding to the therapy, and the evolution of the therapy-treated disease in the patient’s body. The evolution of any disease described in the medical literature usually has such an entangled relationship with both the applied therapy and the patient’s body’s reaction to the disease and therapy simultaneously. When we talk about the evolution of a disease, we are describing and measuring the effects of the interaction of the aforementioned three elements, and the stochastic process we are analyzing concerns the temporal evolution of the combination of these three elements. We do suppose that such complicated relations cause unexpected jumps of registered disease symptoms. In fact, disease observation is a registration of system dynamics composed with three interacting parts. So, non-Gaussian processes seem to be expected in such systems.
The results obtained, which indicate that the stochastic process describing the evolution of a disease can have “fat tails,” show at the same time that “grey swans” in the observation of a patient’s condition are a natural phenomenon and can occur often. If the process were Gaussian, then their occurrence would be rather unexpected, and it is known from medical practice that they are a fairly common part of clinical observations. Finally, we achieve the aim of the investigation and prove that “fat tail” probability distribution is present in the disease symptoms’ observation.

Funding

The experimental part of this work was supported by the National Centre for Research and Development and prepared within the framework of the scientific project “A new Model of medical care with Use of modern methods of non-invasive clinical assessment and Telemedicine in patients with heart failure” (AMULET) (STRATEGMED3/305274/8/NCBR/2017).

Conflicts of Interest

The author declare no conflict of interest.

References

  1. Taleb, N.N. The Black Swan—The Impact of the Highly Improbable; Random House: New York, NY, USA, 2007. [Google Scholar]
  2. Hurst, H.E. Long-Term Storage Capacity of Reservoirs. Trans. Am. Soc. Civ. Eng. 1951, 116, 770–808. [Google Scholar] [CrossRef]
  3. Mandelbrot, B.; Wallis, J.R. Noah, Joseph, and operational hydrology. Water Resour. Res. 1968, 4, 909–918. [Google Scholar] [CrossRef] [Green Version]
  4. Stanley, H.E.; Mantegna, R.N. An Introduction to Econophysics; Cambridge University Press: Cambridge, UK, 2000. [Google Scholar]
  5. do Valle, I.F.; Roweth, H.G.; Malloy, M.W.; Moco, S.; Barron, D.; Battinelli, E.; Loscalzo, J.; Barabási, A.L. Network medicine framework shows that proximity of polyphenol targets and disease proteins predicts therapeutic effects of polyphenols. Nat. Food 2021, 2, 143–155. [Google Scholar] [CrossRef]
  6. Moszczynski, P.; Walczak, A.; Krzesinski, P. Prediction of cardiological parameters evolution with the Markov model. In Proceedings of the 37th IBIMA Conference, Cordoba, Spain, 30–31 May 2021. [Google Scholar]
  7. Walczak, A.; Moszczyński, P.; Krzesiński, P. Evolution of Hemodynamic Parameters Simulated by Means of Diffusion Models. Appl. Sci. 2021, 11, 11412. [Google Scholar] [CrossRef]
  8. Bassler, K.; Gemunu, E.; Gunaratne, H.; McCauley, J.L. Markov Processes, Hurst Exponents, and Nonlinear Diffusion Equations with application to finance. Phys. A Stat. Mech. Its Appl. 2006, 369, 343–353. [Google Scholar] [CrossRef] [Green Version]
  9. Krzesinski, P.; Sobotnicki, A.; Gacek, A.; Siebert, J.; Walczak, A.; Murawski, P.; Gielerak, G. Noninvasive Bioimpedance Methods From the Viewpoint of Remote Monitoring in Heart Failure. JMIR mHealth uHealth 2021, 9, e25937. [Google Scholar] [CrossRef] [PubMed]
  10. Krzesiński, P.; Jankowska, E.A.; Siebert, J.; Galas, A.; Piotrowicz, K.; Stańczyk, A.; Siwołowski, P.; Gutknecht, P.; Chrom, P.; Murawski, P.; et al. Effects of an outpatient intervention comprising nurse-led non-invasive assessments, telemedicine support and remote cardiologists’ decisions in patients with heart failure (AMULET study): A randomised controlled trial. Eur. J. Hear. Fail. 2021, 24, 565–577. [Google Scholar] [CrossRef] [PubMed]
  11. Krzesiński, P.; Siebert, J.; Jankowska, E.A.; Banasiak, W.; Piotrowicz, K.; Stańczyk, A.; Galas, A.; Walczak, A.; Murawski, P.; Chrom, P.; et al. Rationale and design of the AMULET study: A new Model of telemedical care in patients with heart failure. ESC Hear. Fail. 2021, 8, 2569–2579. [Google Scholar] [CrossRef]
  12. Krzesiński, P.; Siebert, J.; Jankowska, E.A.; Galas, A.; Piotrowicz, K.; Stańczyk, A.; Siwołowski, P.; Gutknecht, P.; Chrom, P.; Murawski, P.; et al. Nurse-led ambulatory care supported by non-invasive haemodynamic assessment after acute heart failure decompensation. ESC Hear. Fail. 2021, 8, 1018–1026. [Google Scholar] [CrossRef]
  13. de Groot, S.R.; Mazur, P. Non-Equilibrium Thermodynamics; North-Holland: Amsterdam, The Netherlands, 1969. [Google Scholar]
  14. Van Kampen, N.G. Stochastic Processes in Physics and Chemistry, 3rd ed.; North Holland: Amsterdam, The Netherlands, 2007. [Google Scholar]
  15. Risken, H. The Fokker Planck Equation; Springer: Berlin, Germany, 1989. [Google Scholar]
  16. Coley, W.T.; Kalmykov, Y.P.; Waldron, J.T. The Langevin Equation; World Scientific: Singapore, 1996. [Google Scholar]
  17. Hughes, B.D. Random Walks and Random Environments, Volume 1: Random Walks; Oxford University Press: Oxford, UK, 1995. [Google Scholar]
  18. Weiss, G.H. Aspects and Applications of the Random Walk; North-Holland: Amsterdam, The Netherlands, 1994. [Google Scholar]
  19. Zaslavsky, G.M.; Benkadda, S. Kinetics and Nonlinear Dynamics in Fluids and Plasmas, Chaos; Springer: Berlin, Germany, 1998. [Google Scholar]
  20. Barabási, A.-L. The origin of bursts and heavy tails in human dynamics. Nature 2005, 435, 207–211. [Google Scholar] [CrossRef]
  21. Oliveira, J.G.; Barabási, A.-L. Darwin and Einstein correspondence patterns. Nature 2005, 437, 1251. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Dezsö, Z.; Almaas, E.; Lukács, A.; Rácz, B.; Szakadát, I.; Barabási, A.-L. Dynamics of information access on the web. Phys. Rev. E 2006, 73, 066132. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Vázquez, A. Exact Results for the Barabási Model of Human Dynamics. Phys. Rev. Lett. 2005, 95, 248701. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Metzler, R.; Klafter, J. The random walk’s guide to anomalous diffusion: A fractional dynamics approach. Phys. Rep. 2000, 339, 1–77. [Google Scholar] [CrossRef]
  25. Metzler, R.; Klafter, J. The restaurant at the end of the random walk: Recent developments in the description of anomalous transport by fractional dynamics. J. Phys. A Math. Gen. 2004, 37, R161–R208. [Google Scholar] [CrossRef]
  26. Barkai, E.; Fleurov, V. Simulating Brownian Type of Motion—The Rescaling Velocity Approach vs Langevin Approach. Phys. Rev. E 1995, 52, 137. [Google Scholar] [CrossRef]
  27. Barkai, E.; Fleurov, V. Brownian type of motion of a randomly kicked particle far from and close to the diffusion limit. Phys. Rev. E 1995, 52, 1558–1570. [Google Scholar] [CrossRef]
  28. Barkai, E.; Eisenberg, R.S.; Schuss, Z. Bi-Directional Shot Noise in a Singly Occupied Channel. Phys. Rev. E 1996, 54, 1161. [Google Scholar] [CrossRef]
  29. Barkai, E.; Fleurov, V. Dissipation and fluctuation for a randomly kicked particle: Normal and anomalous diffusion. Chem. Phys. 1996, 212, 69–88. [Google Scholar] [CrossRef]
  30. Barkai, E.; Klafter, J. Crossover from Dispersive to Regular Transport in Biased Maps. Phys. Rev. Lett. 1997, 79, 2245–2248. [Google Scholar] [CrossRef]
  31. Flomenbom, O. Mathematical treatments that solve single molecule. Biophys. Rev. Lett. 2013, 8, 109–136. [Google Scholar] [CrossRef]
  32. Schlesinger, M.F.; Zaslavsky, G.M.; Klafter, J. Strange Kinetics. Nature 1993, 363, 31–37. [Google Scholar] [CrossRef]
  33. Hill, R.M.; Dissado, L.A. Debye, and non-Debye relaxation. J. Phys. C Solid State Phys. 1985, 18, 3829–3836. [Google Scholar] [CrossRef]
  34. Lukichev, A. Relaxation function for the non-Debye relaxation spectra description. Chem. Phys. 2014, 428, 29–33. [Google Scholar] [CrossRef]
  35. Fronczak, A.; Fronczak, P.; Hołyst, J.A. Microscopic explanation of non-Debye relaxation for heat transfer. Phys. A Stat. Mech. Its Appl. 2007, 375, 571–576. [Google Scholar] [CrossRef]
Figure 1. Measured number of hits of TFC values for each period of registration. Items from t1 to t5 mean one month of hits registration indicated by bar color in the figure.
Figure 1. Measured number of hits of TFC values for each period of registration. Items from t1 to t5 mean one month of hits registration indicated by bar color in the figure.
Entropy 24 01780 g001
Figure 2. Theoretical P(x, t) shape for selected α and Kα.
Figure 2. Theoretical P(x, t) shape for selected α and Kα.
Entropy 24 01780 g002
Figure 3. Comparison of TFC clinical observations (blue) and theoretical results (grey) for parameters α and Kα minimizing RMSE error value (RMSE = 0.003).
Figure 3. Comparison of TFC clinical observations (blue) and theoretical results (grey) for parameters α and Kα minimizing RMSE error value (RMSE = 0.003).
Entropy 24 01780 g003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Walczak, A.A. Heart Failure Evolution Model Based on Anomalous Diffusion Theory. Entropy 2022, 24, 1780. https://doi.org/10.3390/e24121780

AMA Style

Walczak AA. Heart Failure Evolution Model Based on Anomalous Diffusion Theory. Entropy. 2022; 24(12):1780. https://doi.org/10.3390/e24121780

Chicago/Turabian Style

Walczak, Andrzej Augustyn. 2022. "Heart Failure Evolution Model Based on Anomalous Diffusion Theory" Entropy 24, no. 12: 1780. https://doi.org/10.3390/e24121780

APA Style

Walczak, A. A. (2022). Heart Failure Evolution Model Based on Anomalous Diffusion Theory. Entropy, 24(12), 1780. https://doi.org/10.3390/e24121780

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop