Intelligent Big Data system based on scientific machine learning of cyber‐physical systems of medical and biological processes Vasyl Martsenyuka, Aleksandra Klos-Witkowskaa, Andriy Sverstiukb, Oksana Bahrii-Zaiatsb, Marcin Bernasa, Krzysztof Witosa a University of Bielsko-Biala, Willowa St., 2, Bielsko-Biala, 43-300, Poland b I. Horbachevsky Ternopil National Medical University, Maidan Voli St., 1, Ternopil, 46002, Ukraine Abstract The work focuses on developing a Big Data system allowing us to process medical data intelligently, which are ingested from different sources and of various types. It aims to design biosensor devices with desired qualitative characteristics, namely their operational stability. The article suggests the mathematical representation of the discrete population dynamics and the dynamic logic of the studied models of the most advanced biosensors as biopixels arrays which are used at the analytics stage of scientific machine learning. Application software for the intelligent Big Data system for investigating the stability of cyber-physical systems for medical purposes using the R package has been developed. The software package for the intelligent Big Data system for investigating the stability of cyber-physical systems for medical applications consists of the following main software modules: unit of ingesting and processing Big Data used to identify input parameters of the model of cyber-physical systems; the software module of research of dynamic behavior of cyber-physical systems; the software module of research of dynamic logic of cyber-physical systems; the block of decision-making on the stability of cyber-physical systems; the block of visualization, lattice images of macrophage/monoclonal antibody data source, a database that receives lattice images of macrophage binding to monoclonal antibodies and a database that receives images of fluorescent biopixels. Using the software package, the results of computer simulation of mathematical models of cyber-physical systems of medical and biological processes in the form of images of macrophage, monoclonal antibodies, connections of macrophages with monoclonal antibodies, fluorescent pixels, and an electrical signal from the converter are obtained. Keywords 1 cyber-physical system, Big Data, stability of the model, differential equations, difference equations, rectangular lattice, hexagonal lattice 1. Introduction The most popular areas of use of Big Data in medicine are diagnostics using medical images in ophthalmology [1-3], oncology [4, 5], dermatology [6], cardiology [7]. These are the branches of modern medicine where most images are used for diagnosis. Over the last decade, we have seen a rapid growth of digital developments, including a variety of sensors, microcontrollers, communication systems to meet societal needs. Thanks to the Internet of Things (IoT) and the use of biosensors, digital electronic health services are increasingly used every year. Biosensors play a critical role in the Internet of Things when it comes to eHealth. There are a CMIS-2021: The Fourth International Workshop on Computer Modeling and Intelligent Systems, April 27, 2021, Zaporizhzhia, Ukraine EMAIL: vmartsenyuk@ath.bielsko.pl (V. Martsenyuk); awitkowska@ath.bielsko.pl (A. Klos-Witkowska); sverstyuk@tdmu.edu.ua (A. Sverstiuk); bagrijzayats@tdmu.edu.ua (O. Bahrii-Zaiats); mbernas@ath.bielsko.pl (M. Bernas); kwitos@ath.bielsko.pl (K. Witos) ORCID: 0000-0001-5622-1038 (V. Martsenyuk); 0000-0003-2319-5974 (A. Klos-Witkowska); 0000-0001-8644-0776 (A. Sverstiuk); 0000- 0002-5533-3561 (O. Bahrii-Zaiats); 0000-0002-0099-1647 (M. Bernas); 0000-0001-7959-7329 (K. Witos) © 2021 Copyright for this paper by its authors. Use permitted under Creative Commons License Attribution 4.0 International (CC BY 4.0). CEUR Workshop Proceedings (CEUR-WS.org) number of sensors available on the market that help people monitor their daily fitness, blood glucose levels, and many sensible applications of primary diagnostics at home. In [8] the problems of safe and confidential use of the most modern biosensors for the Internet of Medical Things are considered. The architecture of intelligent monitoring of electronic health of patients with chronic diseases is presented in [9]. Elements of the architecture include portable devices, biosensors, and smartphones for collecting medical and biological indicators [10]. The intelligent system uses scientific machine learning (SciML) as a tool for analytics in the context of complex applications across science, engineering, and medicine [11, 12]. It focuses on Big Data, obtained from various hospitals, and data obtained from the patient, to diagnose and generate warnings of critical and crisis situations. The system of health monitoring using prognostic calculations using Big Data analysis is presented in [13]. The use of Big Data in endoprosthesis is considered in [14]. Google works with US clinics to track searches for the COVID-19 pandemic and to predict in which region the outbreak will occur [15]. In 2013, thanks to an electronic medical card, which stored all the information about the client, scientists established a relationship between brain degradation and diabetes. Scientists from the University of Cape Town, thanks to a correct analysis using Big Data, analyzed the most common types of cancer and found that malignant tumors of the intestine, lung and ovary have clear genetic markers. Medical Big Data helps prevent the development of the disease at an early stage through the analysis of cardiovascular pressure, pulse, respiration and blood sugar [16, 17]. A few years ago, Apple and IBM joined forces to give IBM Watson even more structured information from the iPhone and Apple Watch. One of the most promising vectors of Big Data development in medicine is surgery [18, 19], where the relevant systems analyze all available protocols of operations in order to assess the possible risks of surgical interventions in the form of visualization in different colors. At home, Americans simply connect a stethoscope to the iPhone and send the data to their doctor. At the same time, a large cluster collects information on all patients. It is sorted by criteria and then looks for dependencies, confirm or refute myths. Nowadays, the integration between physical resources and computing leads to the creation of complicated computing systems with distributed parameters. Such systems are controlled or managed by integrated into the Internet computing resources[20, 21]. In recent years, we have seen an interest in cyber-physical systems (CPS), which are tools for controlling and monitoring the studied indicators using computer technology, in which software is closely related to physical objects. CPS for medical applications are used to measure and analyze biochemical parameters in biological fluids, detection of cancer, markers of heart failure, pathogenic bacteria, to determine the level of contamination of food and the environment. In recent years, CPSs have been used effectively, in which the function of selective elements is performed by biosensors. Such systems are called cyber-physical biosensor systems (CPBSS), which have a number of advantages: high selectivity, stability, performance, affordable cost, the ability to use a wide range of users. As a result of the use of CPBSS, Big Data is obtained, which requires processing and interrelation with the measured medical and biological parameters. The analytical analysis of the obtained data in the form of useful knowledge for use by software services that control the studied systems is carried out. For Big Data analysis, SciML is used as a data-driven method for studying computer algorithms that are automatically improved through experience [22, 23]. In CPS, “experience” is usually presented in the form of data, and SciML allows you to create mathematical models based on sample data for forecasting and decision making. SciML has made significant strides in a variety of fields, including computer vision, language recognition, and control systems, as well as related scientific fields in physics, chemistry, and biology [24]. Processing data ingested from CPS of medical and biological processes is one of the important solutions which fits the “5 V’s of Big Data” characterization of Big Data [25]. The volume of such type of data grows drastically. So, the amount of healthcare data gathered from 2013 to 2020 has grown 15 times to 2314 exabytes [26]. The functioning of most CPS in medical applications is in real- type mode, corresponding to the requirements of velocity for BigData. Data sources used for medical CPS varies from numerical and literal signals (ECG, EEG, EMG, etc.) to medical images (CT, MRT, USG, etc.), which fits a variety of data sources. The veracity of medical CPS data is assured by the quality of standardized measurements used. Finally, the Value of Big Data solutions within medical data-driven-CPS is supported by new results on diagnostics, treatment, and prophylaxis, which can be exceptionally obtained using SciML. Algorithms and methods for processing measurement parameters in the CPS of medical and biological processes using Big Data are created on the basis of their mathematical models. Mathematical models of biosensors for CPBSS on rectangular lattice using lattice differential equations with delay have been developed [27, 28]. The mathematical models of biosensors for CPBSS on rectangular lattice using lattice difference equations with time delay are developed in the article [29]. A detailed description of the considered model is presented in [30, 31]. When developing the modern CPS with the desired qualitative characteristics we need to process the medical information of the different kinds, which corresponds to the notion of Big Data. Namely, the medical data fit the main features of Big Data such as variety, value, velocity, and veracity. Medical data in CPBSS come from different sources and types, so in their operation it is necessary to ensure independence and operational stability. Therefore, there is a need to develop a software package for the intelligent system BigData in the study of the sustainability of medical CPS. These studies are particularly relevant for autonomy, which is defined as the improvement of enzyme activity, protein when maintained under certain conditions, as well as operational stability to ensure the performance of CPS when measuring biomedical parameters. The usefulness of the studied systems is determined by the stability of the sensitive element, which is located in the receptor, and the functioning of the corresponding components of the matrix of the system during its use. The purpose of the study: to develop a software package for the intelligent system BigData in the study of the stability of medical CPS. 1. Mathematical models of biosensors of the intelligent Big Data system for medical aplications 1.1. Mathematical model of biosensor on rectangular lattice using differential equations with time delay In models of lattice differential equations, the spatial structure has a discrete character and lattice dynamics is widely used in the examples considered in [30, 31] and the mathematical model of the immunosensor, which characterizes the change in the concentration of macrophages M i , j (t ) and monoclonal antibodies Ai , j (t ) was substantiated: dM i , j (t )  (  - Ai , j (t -  ) -   M i , j (t -  )) M i , j (t )  S  M i , j , ∧ dt (1) dAi , j (t )  (-  f   M i , j (t -  ) -  f Ai , j (t )) Ai , j (t ), t  0, dt where  - the speed of detection of macrophages;  - birth rate constant of macrophages;  - time delays in the formation of the macrophage-monoclonal antibody complex;  M - the average rate of decrease in the concentration of macrophages;  A - the average rate of decrease in the concentration of monoclonal antibodies;  - probabilistic rate of formation of the macrophage-monoclonal antibody complex;  - distance between immunopixels; D - macrophage diffusion coefficient; ∧ { } S M i , j - operator of diffusion processes. The model of the immunosensor (1) using lattice differential equations most fully takes into account all parameters of immunosensory systems, namely: to take into account when measuring the concentration of macrophages M i , j (t ) and monoclonal antibodies Ai , j (t ) and their binding properties, constant fertility of macrophages, time delay formation of a macrophage-monoclonal antibody complex; the average rate of decrease in the concentration of macrophages and monoclonal antibodies; probabilistic rate of formation of the macrophage-monoclonal antibody complex; the distance between the immunopixels; macrophage diffusion coefficient; diffusion process operator. These parameters are used by immunosensory systems. Model (1) is given by initial functions (2) M i , j (t )  M i0, j (t ) ≥ 0, Ai , j (t )  Ai0, j (t ) ≥ 0, (2) t ∈[- ,0), M i , j (0), Ai , j (0)  0. For a square array N × N , the discrete diffusion is used for the spatial operator, taking into account the imbalance constant ndsbn    D-2 M 1, 2  M 2,1  M i , j -1 - 2ndsbn M 1,1 , i, j  1  -2    D M 2, j  M 1, j -1  M 1, j 1  M i , j 1 - 3ndsbn M i , j , i  1, j ∈2, N - 1  -2  D M 1, N -1  M 2, N - 2ndsbn M 1, N , i, j ∈2, N - 1  -2 (3)  D M i -1, N  M i 1, N  M i , N -1 - 3ndsbn M i , N , i ∈2, N - 1, j  N    S M i , j   D- 2 M N -1, N  M N , N -1 - 2ndsbn M N , N , i  N , j  N ∧  -2    D M 00 N -1, j  M N , j -1  M N , j 1  M i , j 1 - 3ndsbn M N , j , i  N , j ∈2, N - 1  D M N -1,1  M N , 2 - 2ndsbn M N ,1 , i  N , j  1  -2  D- 2 M  M i 1,1  M i , 2 - 3ndsbn M i ,1 , i ∈2, N - 1, j  1  i -1,1    D- 2 M i -1, j  M i 1, j  M i , j -1  M i , j 1 - 4ndsbn M i , j , i, j ∈2, N - 1.  Each pixel is exposed to macrophages from four adjacent pixels that are separated by equal distances  . The boundary condition M i , j  0 for the edges of the array i, j  0, N  1 is used. Figure 1: A rectangular grid that connects four adjacent pixels in the biosensor model using a Cartesian coordinate system 1.2. Mathematical model of biosensor on hexagonal lattice using lattice difference equations with time delay For the intelligent Big Data system of medical and biological processes dynamics, we use the mathematical description with the help of nonlinear difference equations with delay. The model of the biosensor on the basis of a hexagonal lattice is considered. In this case, for the numbering of biopixels (i, j , k ) , i, j, k   N , N , i  j  k  0 the cubic coordinate system is used [32, 33]. Let M i , j , k (t ) is the concentration of macrophages, Ai , j , k (t ) is the concentration of monoclonal antibodies in the biopixel (i, j , k ) ; i, j, k   N , N , i  j  k  0 . In the model it is used next parameters:   0 - the speed of detection of macrophages;  - birth rate constant of macrophages;  - time delays in the formation of the macrophage-monoclonal antibody complex;  M - the average rate of decrease in the concentration of macrophages;  A - the average rate of decrease in the concentration of monoclonal antibodies;  - probabilistic rate of formation of the macrophage-monoclonal antibody complex;   0 - distance between immunopixels; D - macrophage diffusion coefficient; ∧ { } S M i , j - operator of diffusion processes. When colonies of monoclonal antibodies are absent, colonies of macrophages are regulated by a logistic equation with time delay: M i , j ,k ( n  1)  (1     v M i , j ,k (n  r ))M i , j ,k (n) , (4) where  and  – positive numbers, and r  0 mean latency of the negative responce of the macrophages’ colonies. The diffusion of macrophages from six adjacent pixels is considered (i  1, j , k  1) , (i  1, j  1, k ) , (i, j  1, k  1) , (i  1, j, k  1) , (i  1, j  1, k ) і (i, j  1, k  1) (Fig. 2) with diffusion speed D 2 . Fig. 2. Hexagonal lattice, which binds six neighboring pixels in the model of the biopixels using the cubic coordinates: 1, 3, 5, 8, 9, 11 -  D M (t )  ; 2 -  2 M i 1, j , k 1 (t )  ;  D  i , j ,k   2   4 -  D M (t )  ; 6 -  2 M i , j 1, k 1 (t )  ; 7 -  2 M i 1, j , k 1 (t )  ; 10 -  2 M i 1, j 1, k (t )  ; 12 -  2 M i , j 1, k 1 (t )  .  D  D  D  D  i 1, j 1, k      2      The mathematical model of late- macrophage - monoclonal antibody interaction for a hexagonal array of biopixels using based on the well-known Marchuk model [34-36] and uses the spatial  operator S proposed in [37]. M i , j , k ( n  1)  M i , j , k ( n) exp  -  Ai , j , k n - r  -   M i , j ,k n - r  Sˆ M i , j ,k ( n), (5) Ai , j ,k ( n  1)  Ai , j ,k ( n) exp-  f   M i , j ,k ( n - r ) -  f Ai , j , k ( n), n  0 where S M i , j , k  is a discrete diffusion for a spatial operator S . ∧  M i , j ,k (n)  M i0, j ,k (n) ≥ 0, Ai , j ,k (n)  Ai0, j ,k (n) ≥ 0, (6) n ∈[-r ,0), M i , j ,k (0), Ai , j ,k (0)  0. 2. Dynamic logical simulation of intelligent Big Data system for investigation the stability of CPS of medical and biological processes. In order to simulate the dynamic logic of an intelligent Big Data system for investigation the stability of CPS of medical and biological processes, the syntax which is proposed by A. Platzer for the general CPS [38, 39]it is used. The CPS uses the Hybrid Programming Language (HP), which has more features than differential equations. Consider the dynamic logical simulation of intelligent Big Data system for investigation the stability of CPS for medical applications on the example of the mathematical model of biosensor on hexagonal lattice using lattice difference equations with time delay. Dynamic program is the first level of HP and it is defined by the following grammar:     a :: M i , j ,k (n  1)  M i , j ,k (n) exp   Ai , j ,k (n  r )   x Ai , j ,k (n  r )  Sˆ M i , j , k (n) , (7)   Ai , j ,k (n  1)  Ai , j ,k (n) exp   f  M i , j ,k (n  r )   y Ai , j ,k (n) &  t . where  t is an evolutionary domain constraint in the form of a formula for the logic of the first order of real arithmetic def  t  M min  Vi , j ,k (n)  M max (8)  Amin  Ai , j ,k (n)  Amax  i, j, k   N , N  n  0, i  j  k  0 The functioning of the biopixel (i, j , k ) is determined by two states, with respect to fluorescence. Namely, s f l is a state of fluorescence and snon f l is one of the non-fluorescence states. The use of the first order of semantics of logic and the satisfaction ratio s  L for the first-order formula L of real arithmetic and state s can be determined for some pixels i, j, k ; i, j, k   N , N , i  j  k  0 states s f l and snon f l as s ft  k fl M i , j ,k (n) Ai , j ,k (n) ≥  fl , (9) snonfl  k fl M i , j ,k (n) Ai , j ,k (n)   fl Discrete changes occur in computer programs when they accept new values for variables. This situation occurs when a fluorescence phenomenon occurs in a pixel i, j , k  ; i, j , k   N , N , i  j  k  0 . The state s f l ,i , j ,k : 1 is assigned a value of 1 to the variable s f l ,i , j ,k . This leads to a discrete, jump-like change, as the value s f l ,i , j ,k does not change smoothly, but rapidly when it suddenly changes from 1 to s f l ,i , j ,k , causing a discrete jump of values s f l ,i , j ,k . In this way, we obtain a discrete model of change s f l ,i , j ,k : 1 , except for the model of change (9). 3. Software complex of the intelligent Big Data system for investigation the stability of CPS of medical and biological processes The software package for the intelligent Big Data system for investigation the stability of CPS of medical and biological processes consists of the following main software modules: unit of ingesting and processing Big Data used to identify input parameters of the model of CPS, the software module of research of dynamic behavior of CPS, the software module of research of dynamic logic of CPS, the block of decision-making on stability of CPS, the block of visualization, lattice images of macrophage/ monoclonal antibody data source, a data source that receives lattice images of macrophage binding to monoclonal antibodies and a database that receives images of fluorescent biopixels. A software package for the study of CPS phase diagrams using the R package (http://www.r- project.org/) has been developed [40]. The block diagram of the software package of the intelligent Big Data system for investigation the stability of CPS of medical and biological processes is shown in Fig. 4. Unit for simulation Unit for obtaining Unit of ingesting and of lattice images of lattice images of processing Big Data antigens / antigen-antibody used to identify antibodies bonds input parameters of the model of CPS Software module for studying the Software module for dynamic behavior studying the discrete of CPBSS dynamics of CPBSS Unit for obtaining Unit for obtaining Unit for receiving bifurcation images of an electrical signal diagrams fluorescent from the converter biopixels Phase diagrams The decision- study unit Visualization unit making unit on the stability of CPBSS Data source of Data source of Data source of lattice images of fluorescent lattice images of antigens / biopixel images antigen-antibody antibodies binding Figure 4: Block diagram of the software package of the intelligent Big Data system for investigation the stability of CPS of medical and biological processes The software module for studying the dynamic behavior of CPS consists of blocks for obtaining bifurcation and phase diagrams. Using software to study the stability of intelligent Big Data system for investigation the stability of CPS of medical and biological processes, obtained the results of computer simulation in the form of bifurcation and phase diagrams of macrophage populations relative to monoclonal antibodies, lattice images of macrophages, monoclonal antibodies, probabilities of macrophage binding to monoclonal antibodies and monoclonal antibodies, the signal from the converter. The software module for studying the dynamic logic of CPS consists of blocks for modeling lattice images of macrophages/ monoclonal antibodies, block for obtaining lattice images of macrophage binding to monoclonal antibodies, obtaining images of fluorescent biopixels and the electrical signal from the converter. The study of CPS stability of medical and biological processes based on the Big Data intelligent system is an important data source: lattice images of macrophage/ monoclonal antibody data source, a data source that receives lattice images of macrophage binding to monoclonal antibodies, and a database that receives images of fluorescent biopixels. 4. Computer simulation of mathematical models of biosensors using an intelligent Big Data system for investigation the stability of CPS for medical applications 4.1. Parameters of mathematical models of biosensors for investigation the stability of CPS for medical applications To study the occurrence of bifurcation and deterministic chaos in compartmental mathematical models of lattice type on a rectangular lattice using differential equations, models (1), (4) at N  4 ml and values of parameters   2 min 1 ,   2 ,  f  1 min 1 ,   0 .8 /  , minꞏmkg 2    0 .5 ml ,  f  0.5 ml , D  0.2 nm ,   0.3nm are considered. The concentrations min mkg min mkg min of macrophages populations Vi, j (t ) , Vi , j , k (t ) and monoclonal antibodies populations Fi, j (t ) , Fi, j , k (t ) are measured in mkg . ml 4.2. Results of computer simulation of mathematical model of biosensor on rectangular lattice using lattice differential equations with delay The computer simulation was implemented for different values  . The long-term behavior of the model (1) – (3) for   0.05 ,  = 0.22 ,  = 0.2865 with a set of parameter values, which are presented above is analyzed. Qualitative changes in the behavior of biopixels and the intelligent Big Data system model for the investigation of the stability of the CPS of medical and biological processes are observed. Figures 5 and 6 show the result of computer simulation of the discrete dynamics of the intelligent Big Data system for the investigation of the stability of CPS of medical and biological processes in the form of lattice images of macrophages and monoclonal antibodies in pixels of the studied system. Figures 5 (a) and 6 (a) show the results of computer simulation of lattice images of macrophages and monoclonal antibodies in the pixels of the system (1) at  = 0.05 , which corresponds to a stable focus. For  = 0.22 there is a less pronounced (Fig. 5 (b) and Fig. 6 (b)), and for  = 0.2865 more pronounced traveling wave of monoclonal antibodies, which is presented in Figures 5 (c) and 6 (c). а) b) c) Figure 5: Lattice images of macrophages in pixels of system (1) at  = 0.05 (a),  = 0.22 (b),  = 0.2865 (c) а) b) c) Figure 6: Lattice images of monoclonal antibodies in pixels of system (1) at   0.05 (a),  = 0.22 (b),  = 0.2865 (c) In the second stage of computer modeling of intelligent Big Data system for investigation of the stability of CPS of medical and biological processes lattice graphs are used on which for each pixel the probability of contact of macrophages with monoclonal antibodies, as Vi , j × Fi , j at  = 0.05 ,  = 0.22 ,  = 0.287 is presented. They are shown in Figures 7 (a - c). а) b) c) Figure 7: Lattice images of connections of macrophages with monoclonal antibodies in pixels of system (1) at  = 0.05 (а),  = 0.22 (b),  = 0.2865 (c) Graphs of fluorescent pixels are shown in Figures 8 (a - c). а) b) c) Figure 8: Fluorescence image of system (1) as a result of numerical simulation at  = 0.05 (а),  = 0.22 (b),  = 0.2865 (c) Figure 8 (a) shows the result of numerical simulation of system (1) at  = 0.05 , which corresponds to a stable focus. For  = 0.22 there is a less pronounced (Fig. 8 (b)), and for  = 0.2865 it is more pronounced traveling wave of fluorescent pixels, which is shown in Figure 8 (c). Figure 8 (c) shows the result of numerical simulation of system (1), which corresponds to the approach to the limit cycle (there is a traveling wave of fluorescent pixels). In the case  = 0.2865, chaotic behavior is observed, which begins with wave-like changes in fluorescent pixels (Fig. 8 (c)) and quickly progresses to chaotic changes. Figure 8 (c) shows the result of numerical simulation of system (1), in which a chaotic wave of fluorescent pixels is observed. To control the measurement process in the intelligent Big Data system for investigation of the stability of the CPS for medical applications, the result of computer simulation of the electrical signal from the converter (Fig. 9 (a, b)), which characterizes the number of fluorescent pixels at values of  = 0.05 ,  = 0.2865 plays an extremely important role. Analyzing the type of electrical signal in Figure 9, we see that when the value  changes qualitatively changes the behavior of the pixels and the whole intelligent Big Data system to investigate the stability of the CPS for medical and biological applications. In Figure 9 (a) there is a steady state at  = 0.05 . In Figure 9 (b) (  = 0.2865) a traveling wave of non-fluorescent pixels is clearly visible. а) b) Figure 9: Electrical signal from the converter, which characterizes the number of fluorescent pixels at  = 0.05 (a),  = 0.2865 (b) As shown by numerical analysis, fluorescent states in biopixels change according to the laws of discrete dynamics. The threshold value for fluorescence is  fl = 1,5 . Taking into account the continuous dynamics of the immunological response, each biopixel is considered as a CPS. 4.3. Results of computer simulation of mathematical model of biosensor on hexagonal lattice using lattice difference equations with delay The long - term behavior of model (4) at r = 5 , r = 17 , r = 22 with a set of parameter values, which are presented above (Fig. 10 - 14), is analyzed. We observe qualitative changes in the behavior of biopixels and the biosensor model on hexagonal lattice using lattice difference equations with delay in general. Figures 10 and 11 show the first stage of computer modeling of discrete dynamics of an intelligent Big Data system for investigation of the stability of CPS of medical and biological processes in the form of lattice images of macrophages and monoclonal antibodies in pixels of the studied system. Figure 10 (a) shows the result of numerical simulation of system (4) at r = 5 , which corresponds to a stable focus. For r = 17 there is less pronounced (Fig. 10 (b)), and for r = 22 more pronounced wavy changes in the lattice images of macrophages and monoclonal antibodies in the pixels of the system (4), as shown in Figures 10 (c). а) b) c) Figure 10: Lattice images of macrophages in pixels of system (4) at r = 5 (a), r = 17 (b), r = 22 (c) а) b) c) Figure 11: Lattice images of monoclonal antibodies in pixels of the system (4) at r = 5 (a), r = 17 (b), r = 22 (c) In the second stage of computer modeling of the intelligent Big Data system, lattice graphs were used to investigate the stability of the CPS of medical and biological processes. Firstly, the corresponding graphs are constructed, on which for each pixel the image of probability of contact of macrophages with monoclonal antibodies, as Vi , j , k × Fi , j , k at r = 5 , r = 17 , r = 22 is presented in Figure 12 (a - c). а) b) c) Figure 12: Lattice images of macrophage‐ monoclonal antibody binding in pixels of system (4) at r = 5 (a), r = 17 (b), r = 22 (c) Figure 12 (a) shows the result of numerical simulation of system (4) at r = 5 , which corresponds to a stable focus. For r = 17 there is less pronounced (Fig. 12 (b)), and for r = 22 more pronounced wavy changes in the images of the probability of contact of macrophages with monoclonal antibodies of system (4), as shown in Figure 12 (c). In the third stage of computer modeling of the intelligent Big Data system for investigation of the stability of CPS of medical and biological processes, lattice graphs of fluorescent pixels were obtained based on the fulfillment of condition (6), which are shown in Figures 13 (a - c). Figure 13 (a) shows the result of numerical simulation of system (4) at r = 5 , which corresponds to a stable focus. For r = 17 there is a less pronounced (Fig. 13 (b)), and for r = 22 - it is more pronounced traveling wave of fluorescent pixels of the system (4), which is shown in Figure 13 (c). а) b) c) Figure 13: Fluorescence image of the system (4) as a result of numerical simulation at r  5 (a), r  17 (b), r  22 (c) Hopf bifurcation is observed with increasing time delay [41]. Figure 13 (c) shows the result of numerical simulation of the system (4) at r = 17 , which corresponds to the limit cycle (there is a traveling wave of fluorescent pixels). In the case r = 22 , chaotic behavior is observed, which begins with wave-like changes in fluorescent pixels (Fig. 13 (c)) and quickly turns into chaotic changes. Figure 14 shows the electric signal from the converter (Fig. 14), which characterizes the number of fluorescent pixels at r = 17 , r = 22 . а) b) Figure 14: Electrical signal from the converter, which characterizes the number of fluorescent pixels at r = 17 (a), r = 22 (b) To control the measurement process using an electrical signal to an intelligent big data system for investigation of the stability of the CPS of medical and biological processes is extremely important the result of computer simulation. Analyzing the type of electrical signal in Figure 14, we can conclude that changing the value qualitatively changes the behavior of the pixels and the entire biosensor model based on lattice difference equations. In Figure 14 (a) ( r = 17 ) occurs, and in Figure 14 (b) ( r = 22 ) there is a traveling wave of non-fluorescent pixels. The threshold value for fluorescence is  fl = 1,5 . A two-dimensional array of biopixels was used to develop the dynamic logic of an intelligent Big Data system to investigate the stability of the CPS of medical and biological processes on a hexagonal lattice using delay equations. As shown by the results of numerical analysis, fluorescent states in biopixels change according to the laws of discrete dynamics. 5. Conclusion In the work, the general scheme of the cyber-physical sensor system proposed in [42] was extended to the usage of Big Data. Peculiarities of immunosensors functioning are taken into account in CPS mathematical models. According to the laws of discrete dynamics, lattice images in immunopixels are modified. The proposed mathematical models take into account the interaction of immunopixels through the diffusion of macrophages. The intelligent Big Data system uses the discrete dynamics of macrophages and antibodies, which is based on dynamic logic, to study the stability of CPS in medicine. The classes of rectangular lattice using differential equations or hexagonal lattice using difference equations are used in the work as modeling of the interaction of macrophages and monoclonal antibodies in immunopixels. This also takes into account the presence of macrophage and monoclonal antibodies colonies, that are localized in pixels, as well as the diffusion of macrophage colonies between pixels using the R package. According to a series of experiments of intelligent Big Data system for investigation of the stability of CPS of medical and biological processes using lattice differential equations with delay it is established that in the case of using a rectangular lattice at constant time delay  ∈[0, 0.22] the solution of the studied system is a stable focus. When   0.23 (in the case of a rectangular lattice) there is a Hopf bifurcation and all subsequent trajectories correspond to stable boundary cycles for all imunopixels. With a further increase in the delay constant  , chaotic behavior of the CPS of medical processes occurs on the basis of lattice differential equations with time delay. According to the results of a series of experiments of intelligent Big Data system for investigation of the stability of CPS of medical processes using lattice difference equations with delay it is established that in case of using hexagonal lattice at r ∈[0, 16] the solution of the investigated system tend to endemic states, which are a stable focus. For r  17 (in the case of a hexagonal lattice) there is a Hopf bifurcation and all subsequent trajectories correspond to stable boundary cycles for all imunopixels. With a further increase in the delay constant r , chaotic behavior of the CPS occurs on the basis of lattice difference equations. The developed computer programs for the study stability of intelligent Big Data system for investigation the stability of CPS of medical processes in the form of phase diagrams, lattice images of macrophages / monoclonal antibodies, lattice images of macrophage binding to monoclonal antibodies, images of fluorescent imunopixels and the electrical signal from the converter should be used in research, design organizations, medical and laboratory centers in the development and testing of cyber-physical systems of medical processes. In the following researches is necessary to offer mathematical models of intelligent systems of Big Data on the basis of other geometrical structures and to carry out research of their stability. 6. Aknowledgements The work was co-funded by the European Union's Erasmus + Programme for Education under KA2 grant (project no. 2020-1-PL01-KA203-082197 “Innovations for Big Data in a Real World”) 7. References [1] S. M. Khan, X. Liu, S. Nath, E. Korot, L. Faes, S. K. Wagner, P. A. Keane, N. J. Sebire, M. J. Burton, A. K. Denniston, A global review of publicly available datasets for ophthalmological imaging: barriers to access, usability, and generalisability, in: The Lancet Digital Health, Volume 3, Issue 1, 2021, pp. e51-e66. doi:10.1016/S2589-7500(20)30240-5. [2] V. Sujatha, S. Prasanna Devi, S. Vinu Kiran, S. Manivannan, Big Data Analytics on Diabetic Retinopathy Study (DRS) on Real-time Data Set Identifying Survival Time and Length of Stay, in: Procedia Computer Science, Volume 87, 2016, pp. 227-232. doi: 10.1016/j.procs.2016.05.153. [3] T. Inomata, J. Sung, M. Nakamura, K. Fujisawa, K. Muto, N. Ebihara, M. Iwagami, M. Nakamura, K. Fujio, Y. Okumura, M. Okano, A. Murakami, New medical big data for P4 medicine on allergic conjunctivitis, in: Allergology International, Volume 69, Issue 4, 2020, pp. 510-518. doi: 10.1016/j.alit.2020.06.001. [4] D. C. Elbers, N. R. Fillmore, F.-Chi Sung, S. S. Ganas, A. Prokhorenkov, Ch. Meyer, R. B. Hall, S. J. Ajjarapu, D. C. Chen, F. Meng, R. L. Grossman, M. T. Brophy, N. V. Do, The Veterans Affairs Precision Oncology Data Repository, a Clinical, Genomic, and Imaging Research Database, in: Patterns, Volume 1, Issue 6, 2020, p. 100083. doi: 10.1016/j.patter.2020.100083. [5] S. M. Willems, S. Abeln, K. A. Feenstra, Remco de Bree, E. F. van der Poel, R. J. Baatenburg de Jong, J. Heringa, M. W.M. van den Brekel, The potential use of big data in oncology, in: Oral Oncology, Volume 98, 2019, pp. 8-12. doi: 10.1016/j.oraloncology.2019.09.003. [6] R. Tripathi, R. S. Mazmudar, K. D. Knusel, J. S. Bordeaux, J. F. Scott, Big data in dermatology: Publicly available health care databases for population health research, in: Journal of the American Academy of Dermatology, Volume 83, Issue 5, 2020, pp. 1546-1556. doi: 10.1016/j.jaad.2020.04.145. [7] A. L. P. Ribeiro, G. M. M. Paixão, P. R. Gomes, M. H. Ribeiro, A. H. Ribeiro, J. A. Canazart, D. M. Oliveira, M. P. Ferreira, E. M. Lima, J. Lopes de Moraes, N. Castro, L. B. Ribeiro, P. W. Macfarlane, Tele-electrocardiography and Big Data: The CODE (Clinical Outcomes in Digital Electrocardiography) study, in: Journal of Electrocardiology, Volume 57, Supplement, 2019, pp. S75- S78. doi: 10.1016/j.jelectrocard.2019.09.008. [8] P. P. Ray, D. Dash, N. Kumar, Sensors for internet of medical things: State-of-the-art, security and privacy issues, challenges and future directions, in: Computer Communications, Volume 160, 2020, pp. 111-131. doi: 10.1016/j.comcom.2020.05.029. [9] J. Lloret, L. Parra, M. Taha, J. Tomás, An architecture and protocol for smart continuous eHealth monitoring using 5G, in: Computer Networks, Volume 129, Part 2, 2017, pp. 340-351. doi: 10.1016/j.comnet.2017.05.018. [10] I. Jayatilleka, M. N. Halgamuge, Chapter 1 - Internet of Things in healthcare: Smart devices, sensors, and systems related to diseases and health conditions, Editor(s): Himansu Das, Nilanjan Dey, Valentina Emilia Balas, In Advances in Ubiquitous Sensing Applications for Healthcare, Real-Time Data Analytics for Large Scale Sensor Data, Academic Press, Volume 6, 2020, pp. 1-35. [11] Center for Scientific Machine Learning. URL: https://www.oden.utexas.edu/research/centers- groups/scientific-machine-learning/. [12] B. Yuan, C. Shen, A. Luna, A. Korkut, D. S. Marks, J. Ingraham, Ch. Sander, CellBox: Interpretable Machine Learning for Perturbation Biology with Application to the Design of Cancer Combination Therapy, in: Cell Systems, Volume 12, Issue 2, 2021, pp. 128-140. doi: 10.1016/j.cels.2020.11.013. [13] M. Srivathsan, K. Yogesh Arjun, Health Monitoring System by Prognotive Computing Using Big Data Analytics, in: Procedia Computer Science, Volume 50, 2015, pp. 602-609. doi:10.1016/j.procs.2015.04.092. [14] P. N. Ramkumar, H. S. Haeberle, M. R. Bloomfield, J. L. Schaffer, A. F. Kamath, B. M. Patterson, V. E. Krebs, Artificial Intelligence and Arthroplasty at a Single Institution: Real-World Applications of Machine Learning to Big Data, Value-Based Care, Mobile Health, and Remote Patient Monitoring, in: The Journal of Arthroplasty, Volume 34, Issue 10, 2019, pp. 2204-2209. doi: 10.1016/j.arth.2019.06.018. [15] S. J. Kurian, A. R. Bhatti, M. A. Alvi, H. H. Ting, C. Storlie, P. M. Wilson, N. D. Shah, H. Liu, M. Bydon, Correlations Between COVID-19 Cases and Google Trends Data in the United States: A State-by-State Analysis, in: Mayo Clinic Proceedings, Volume 95, Issue 11, 2020, pp. 2370-2381. doi: 10.1016/j.mayocp.2020.08.022. [16] C. Baladrón, J. J. Gómez de Diego, I. J. Amat-Santos, Big data and new information technology: what cardiologists need to know, in: Revista Española de Cardiología (English Edition), Volume 74, Issue 1, 2021, pp. 81-89. doi: 10.1016/j.rec.2020.06.036. [17] A. Zamanifar, 7 - Data analytics in IOT-based health care, Editor(s): K. Ch. Lee, S. S. Roy, P. Samui, V. Kumar, in: Data Analytics in Biomedical Engineering and Healthcare, Academic Press, 2021, pp. 119-130. doi: 10.1016/B978-0-12-819314-3.00007-0. [18] A. Perry, P. Kerezoudis, Ch. S. Graffeo, L. P. Carlstrom, M. Peris-Celda, F. B. Meyer, M. Bydon, M. J. Link, Little Insights from Big Data: Cerebrospinal Fluid Leak After Skull Base Surgery and the Limitations of Database Research, in: World Neurosurgery, Volume 127, 2019, pp. e561- e569. doi: 10.1016/j.wneu.2019.03.207. [19] J. Wall, T. Krummel, The digital surgeon: How big data, automation, and artificial intelligence will change surgical practice, in: Journal of Pediatric Surgery, Volume 55, Supplement, 2020, pp. 47-50. doi: 10.1016/j.jpedsurg.2019.09.008. [20] Y. Krainyk, Y. Davydenko and V. Starchenko, Message-level Decoding of Error Patterns for Turbo-Product Codes, in: Proceedings of the 39th International Conference on Electronics and Nanotechnology (ELNANO), Kyiv, Ukraine, 2019, pp. 660–663. doi: 10.1109/ELNANO.2019.8783849. [21] Y. Krainyk, Y. Davydenko and V. Tomas, Configurable Control Node for Wireless Sensor Network, in: Proceedings of the 3rd International Conference on Advanced Information and Communications Technologies (AICT), Lviv, Ukraine, 2019, pp. 258–262. doi: 10.1109/AIACT.2019.8847732. [22] D. Michie, D.J. Spiegelhalter, C. Taylor Machine Learning. Neural and Statistical Classification (1994), pp. 1-298. [23] M. Rezaei-Ravari, M. Eftekhari, F. Saberi-Movahed, Regularizing extreme learning machine by dual locally linear embedding manifold learning for training multi-label neural network classifiers, in: Engineering Applications of Artificial Intelligence, Volume 97, 2021, p. 104062. doi: 10.1016/j.engappai.2020.104062. [24] L. Huang, H. Leng, X. Li, K. Ren, J. Song, D. Wang, A Data-Driven Method for Hybrid Data Assimilation with Multilayer Perceptron, in: Big Data Research, Volume 23, 2021, p. 100179. doi: 10.1016/j.bdr.2020.100179. [25] The seven Vs of Big Data. URL: https://www.bbva.com/en/five-vs-big-data/. [26] Total amount of global healthcare data generated in 2013 and a projection for 2020. URL: https://www.statista.com/statistics/1037970/global-healthcare-data-volume/. [27] V.P.Martsenyuk, I.Ye. Andrushchak., P.M. Zinko, A.S. Sverstiuk, On Application of Latticed Differential Equations with a Delay for Immunosensor Modeling, in: Journal of Automation and Information Sciences, Volume 50, Issue 6, 2018, pp. 55–65. doi: 10.1615/JAutomatInfScien.v50.i6.50. [28] V.P. Martsenyuk, A.S. Sverstiuk, I.Ye. Andrushchak, Approach to the Study of Global Asymptotic Stability of Lattice Differential Equations with Delay for Modeling of Immunosensors, in: Journal of Automation and Information Science, Volume 48, Issue 8, 2019, pp. 58–71. doi: 10.1615/JAutomatInfScien.v51.i2.70. [29] V.P. Martsenyuk, A. Klos-Witkowska, A.S.Sverstiuk, Stability investigation of biosensor model based on lattice difference equations. Difference Equations and Discrete Dynamical Systems with Applications, in: Springer Proceedings in Mathematics and Statistics. 24th ICDEA, Dresden, Germany. Volume 312, 2020, pp. 297–322. doi: 10.1007/978-3-030-35502-9_13. [30] V.P. Martsenyuk, A. Klos-Witkowska, A.S. Sverstiuk, Stability, bifurcation and transition to chaos in a model of immunosensor based on lattice differential equations with delay, in: Electronic Journal of Qualitative Theory of Differential Equations: No. 27 (2018), pp. 1–31. doi: 10.14232/ejqtde.2018.1.27. [31] V.P. Martsenyuk, A.S. Sverstiuk, A. Klos-Witkowska, N.V. Kozodii, O.A. Bagriy-Zayats, I.R. Zubenko, Numerical Analysis of Results Simulation of Cyber-physical Biosensor Systems, in: 1st International Workshop Informatio-Communication Techologies&Embedded Systems, 14–15 November, Mykolaiv, Volume 1, 2019, pp. 149–164. [32] A.Patel, Hexagonal Grids. URL: https://www.redblobgames.com/grids/hexagons/. [33] L. Čomić, B. Nagy, A combinatorial coordinate system for the body-centered cubic grid, in: Graphical Models, Volume 87, 2016, pp. 11-22. doi: 10.1016/j.gmod.2016.08.001. [34] U. Forys, Marchuk's model of immune system dynamics with application to tumour growth, in: Journal of Theoretical Medicine, Volume 4, no. 1, 2002, pp. 85-93. doi: 10.1080/10273660290052151. [35] U. Forys, Stability and bifurcations for the chronic state in Marchuk's model of an immune system, in: Journal of Mathematical Analysis and Applications, Volume 352, Issue 2, 2009, pp. 922- 942. doi: 10.1016/j.jmaa.2008.11.055. [36] Z. Huang, Qigui Yang, Junfei Cao, Stochastic stability and bifurcation for the chronic state in Marchuk’s model with noise, in: Applied Mathematical Modelling, Volume 35, Issue 12, 2011, pp. 5842-5855. doi: 10.1016/j.apm.2011.05.027. [37] A. Prindle, P. Samayoa, I. Razinkov, T. Danino, LS Tsimring, and J. Hasty, A sensing array of radically coupled genetic 'biopixels', in: Nature, Volume 481, no. 7379, Dec. 2011, pp. 39-44. [38] A. Platzer, Logics of Dynamical Systems. 27th Annual IEEE Symposium on Logic in Computer Science, Dubrovnik, 2012, pp. 13–24. doi: 10.1038/nature10722. [39] A. Platzer, The complete proof theory of hybrid systems, in: 27th Annual ACM/IEEE Symposium on Logic in Computer Science, 2012, pp. 541–550. doi: 10.1109/lics.2012.64. [40] V.P. Martsenyuk, A.S. Sverstiuk, O.A. Bahrii-Zaiats, Yu.A. Rudyak, B. Shelestovsky, Software complex in the study of the mathematical model of cyber-physical systems, in: ICT&ES 2nd International Workshop Information-Communication Technologies & Embedded Systems, 12 November, 2020 Mykolaiv, Ukraine. Volume 1, 2020, pp. 87–97. [41] D. Duan, B. Niu, J. Wei, Hopf-Hopf bifurcation and chaotic attractors in a delayed diffusive predator-prey model with fear effect, in: Chaos, Solitons & Fractals, Volume 123, 2019, pp. 206-216. doi: 10.1016/j.chaos.2019.04.012. [42] C. Berger, A.Hees, S.Braunreuther, G.Reinhart, Characterization of Cyber-Physical Sensor Systems, in: Procedia CIRP, Volume 41, 2016, pp. 638–643. doi: 10.1016/j.procir.2015.12.019.