Chronic disease co-disease mode identification method and device based on relaxation constraint symmetric low-rank representation, electronic equipment and storage medium

Abstract

The invention discloses a chronic disease co-disease mode recognition method and device based on relaxation constraint symmetric low-rank characterization, electronic equipment and a storage medium, and the method comprises the steps: obtaining medical data resources, constructing a chronic disease co-disease network through the medical data resources, and enabling the chronic disease co-disease network to be used for describing the correlation and influence between diseases, the method comprises the following steps: performing low-rank representation learning on a chronic disease co-disease network, mining and outputting a community structure therein, identifying the community structure in the co-disease network by using a community division mechanism, mining a co-disease mode, introducing relaxation symmetric constraint into a learning target, and well perceiving the inherent symmetric structure characteristics of the network. And a graph regularization technology for capturing local topological characteristics is introduced, so that the inherent geometric structure characteristics of the network are effectively maintained.

Description

Technical Field

[0001] The present invention mainly relates to the field of pattern recognition, and in particular to a method, device, electronic device and storage medium for chronic disease comorbidity pattern recognition based on relaxed constrained symmetric low-rank representation. Background Art

[0002] Chronic non-communicable diseases (referred to as "chronic diseases") have become a major threat to global public health and human health, bringing heavy economic burdens to patients and society. Due to the characteristics of complex etiology, hidden onset, long course of disease and prolonged illness, some patients often suffer from multiple chronic diseases at the same time. The World Health Organization calls this phenomenon "comorbidity". Under the trend of global population aging, the problem of chronic disease comorbidity has gradually become prominent and has become a major challenge in the field of global health. The complex impact of diseases poses a great challenge to disease diagnosis and treatment. Understanding the specific factors and processes of comorbidity, the interaction of diseases and possible synergies is very important for promoting diagnosis, improving the quality of life of patients, improving prevention and reducing the cost of the health care system. Therefore, the use of new technologies such as artificial intelligence and big data to explore the complex comorbidity patterns between chronic diseases provides important decision-making basis for predicting disease risks, disease prevention, treatment and management, which has practical significance for the sustainable development of human health.

[0003] Diseases do not cluster randomly, and comorbidity patterns reveal statistically significant associations and pathophysiological relationships between diseases. Mining comorbidity patterns and exploring their distribution patterns will help classify diseases, reduce the complexity of chronic disease prevention and treatment, and make preventive measures more accurate and efficient. A systematic review of multi-disease epidemiological studies showed that different studies found different properties and patterns of comorbidity. At present, there is no unified definition of methods for identifying comorbidity patterns. Based on different data analysis principles and technical means, common methods can be summarized into the following four categories:

[0004] (1) Odds ratio: It is used to measure the probability of something happening versus not happening. It is often used to analyze the relative risk of two diseases in comorbidity pattern analysis. A study in the Netherlands showed that the relative risk of depression and anxiety, coronary heart disease and heart failure, and chronic obstructive pulmonary disease and heart failure were the highest. The odds ratio algorithm is simple and easy to use and has become one of the commonly used methods in comorbidity pattern research.

[0005] (2) Factor analysis: aims to identify the complex relationships between multiple diseases, combine multiple related factors into a few comprehensive indicators, simplify the data structure and reveal the association between variables. A study of Japanese adults identified five groups of comorbidity patterns for 17 chronic diseases through factor analysis. Factor analysis does not rely on traditional grouping assumptions. It can more effectively reveal the relationship between variables by reducing the dimension and simplifying the data to extract important information.

[0006] (3) Cluster analysis: A data mining technique used to identify data clusters. Common methods include hierarchical clustering and k-means clustering. The comorbidity patterns of COPD patients may include comorbidities with cardiovascular disease, allergic diseases, and other representative diseases. Cluster analysis can explore the commonalities of the underlying factors of these comorbidity patterns and help understand the comorbidity phenomenon of chronic diseases.

[0007] (4) Association rules: used to analyze the direct associations between multiple chronic diseases and reveal the strong associations of comorbidity patterns. A British study evaluated the comorbidity patterns of middle-aged and elderly people through cluster analysis and association rule mining, and found three clusters and 30 disease patterns, revealing the complex relationships and central roles between different diseases, such as diabetes, hypertension, and asthma. There are complex interactions and associations between these chronic diseases, which require comprehensive evaluation and treatment to achieve better health outcomes.

[0008] Due to the complex interactions and influences between diseases, traditional methods are difficult to explore the deep relationships and complex connections between diseases, and it is difficult to obtain accurate comorbidity patterns. With the development of complex science, complex network analysis describes complex relationships from a holistic and systematic perspective, which is conducive to comprehensively revealing the complex influence patterns between diseases from a holistic perspective and has become a mainstream tool for exploring chronic disease comorbidity patterns. Community detection algorithms can discover closely connected clusters in the network, which can be used to identify disease comorbidity patterns.

[0009] The present invention designs a new relaxed constrained symmetric low-rank representation method, which strikes a balance between symmetry representation and representation learning ability, and then constructs a community detection model for chronic disease comorbidity networks, and explores comorbidity patterns among chronic diseases from the perspective of complex network representation learning. Summary of the invention

[0010] Based on this, it is necessary to provide a chronic disease comorbidity pattern recognition method, device, electronic device and storage medium based on relaxed constrained symmetric low-rank representation to address existing problems.

[0011] In a first aspect, an embodiment of the present application provides a method for recognizing chronic disease comorbidity patterns based on relaxed constrained symmetric low-rank representation, characterized in that it includes:

[0012] Acquire medical data resources, wherein the medical data resources include data records of chronic disease conditions sorted out from tests, diagnoses, treatments, prescriptions, electronic medical records, and clinical data resources;

[0013] A chronic disease comorbidity network is constructed using medical data resources, wherein the chronic disease comorbidity network is used to describe the correlation and influence between diseases, and can be constructed by calculating the correlation between diseases or patients based on multi-source data, wherein the multi-source data includes basic patient information and test, diagnosis, treatment, prescription data, electronic medical records, clinical data, and omics big data, wherein the chronic disease comorbidity network is represented by a graph G = (V, E), V = {v i |i∈{1,...,n}} represents a set of n disease nodes, E={e ij |i, j∈{1,...,m}} represents a set of m edges; where the adjacency matrix A = [a ij ] is used to represent and store the target comorbidity network G, a ij Represents the disease node v i and v j The interaction between them is equal to I(s i ,s j );

[0014] Performing low-rank representation learning and mining on the chronic disease comorbidity network to output the community structure therein, wherein the low-rank representation learning adopts a symmetric non-negative matrix decomposition algorithm with relaxed constraints to obtain a low-rank representation matrix of the target network;

[0015] Use the community partitioning mechanism to identify the community structure in the comorbidity network and discover the comorbidity pattern.

[0016] Preferably, the method for constructing the chronic disease comorbidity network includes:

[0017] Sort out chronic disease records, and sort out data records from tests, diagnoses, treatments, prescriptions, electronic medical records, and clinical data;

[0018] Count the frequency of disease occurrence, and count the frequency of each disease occurrence Q(s) from the above data records. i ) and the frequency of disease combinations occurring together Q(s i ,s j ), where a disease i If it appears in a disease record, it is counted once. i and j If they co-occur in one disease record, they are counted once, where i, j = 1, 2, ..., n represents the number of diseases of interest;

[0019] Calculate the mutual information value between two diseases as follows:

[0020]

[0021] Among them, p(s i)=Q(s i ) / F,p(s j )=Q(s j ) / F respectively calculate disease s i and j The observation frequency, where F represents the total observation frequency, p(s i ,s j )=Q(s i ,s j ) / F calculates the number of diseases observed at the same time i and j frequency;

[0022] Constructing a comorbidity network to describe the complex impact relationships of chronic diseases.

[0023] Preferably, the method of using a symmetric low-rank characterization method with relaxed constraints to characterize the comorbidity network includes:

[0024] Design the objective function, which is as follows:

[0025]

[0026] stW=P, H=Q, P=Q; P, Q≥0

[0027] Wherein, A is the adjacency matrix corresponding to the comorbidity network; decision parameters W and H represent learnable feature matrices, and P and Q are auxiliary parameter matrices, whose sizes are the same as the feature matrices W and H, respectively; represents the F-norm of the matrix; λ is a positive constant representing the coefficient of the graph regularization term; Tr(·) is used to calculate the trace of the matrix; L represents the Laplacian matrix of the adjacency matrix A. The purpose of introducing auxiliary matrices P and Q is to separate the training and constraints of decision parameters W and H, and to achieve symmetric characterization under the relaxed constraints "W=P, H=Q, P=Q", as well as numerical characteristic characterization under the non-negative constraints "P, Q≥0";

[0028] Construct the augmented Lagrangian function and use the ADMM alternating direction multiplier method to construct the augmented Lagrangian function:

[0029]

[0030] Wherein, A is the adjacency matrix corresponding to the comorbidity network; decision parameters W and H represent feature matrices, and P and Q are auxiliary parameter matrices that are the same as the feature matrices W and H; represents the F norm of the matrix; λ is a positive constant, representing the coefficient of the graph regularization term; Tr(·) calculates the trace of the matrix; L represents the Laplace matrix of the adjacency matrix A; ○ represents the Hadamard product of the matrix; α, β and θ correspond to the weights of the three penalty terms, respectively, controlling their proportional sizes. The non-negativity of P and Q can be achieved by mapping their value ranges to non-negative space;

[0031] The update rules for parameters are derived. In the ADMM optimization framework, let t and t+1 represent the current and updated iteration states of the parameters respectively. All parameters are iteratively updated using the update rules until the convergence condition is reached. The convergence condition is: (1) the absolute value of the difference between the objective function values ​​after two consecutive iterative trainings is less than 0.1; (2) or the number of iterations reaches 20. The update rules are as follows:

[0032]

[0033] in represents the gradient, η is the learning rate in gradient ascent, outputs feature parameters W, H, P and Q, and Q is used as the community partition indicator matrix.

[0034] Preferably, the method of using a community partitioning mechanism to identify the community structure in the comorbidity network to discover the comorbidity pattern includes:

[0035] After training to obtain the characteristic parameters W, H, P and Q, Q is used as the indicator matrix for node-disease comorbidity module affiliation division, and the nodes in the comorbidity network are judged one by one and assigned to the comorbidity module. The assignment rules are as follows:

[0036]

[0037] Where V represents the node set in the comorbidity network, v j represents any node; C s is the comorbidity combination s, s∈{1,2,…,K}, K is the total number of comorbidity combinations; q jk is an element in the low-rank matrix Q, representing node v j The probability of being classified into the comorbidity combination k, q js is the maximum value among them, indicating that node v j The probability of being classified into the comorbidity combination s is the highest. After the comorbidity combinations of all nodes in the target comorbidity network are divided, all possible comorbidity combinations identified can be obtained. C = {C1, C2, ..., C K}.

[0038] In a second aspect, an embodiment of the present application provides a chronic disease comorbidity pattern recognition device based on relaxed constrained symmetric low-rank representation, characterized in that it includes:

[0039] An acquisition unit, used to acquire medical data resources, wherein the medical data resources include data records of chronic disease conditions sorted out from inspections, diagnoses, treatments, prescriptions, electronic medical records, and clinical data resources;

[0040] A construction unit constructs a chronic disease comorbidity network using the medical data resources, wherein the chronic disease comorbidity network is used to describe the correlation and influence between diseases, and can be constructed by calculating the correlation between diseases or patients according to multi-source data, wherein the multi-source data includes basic patient information and test, diagnosis, treatment, prescription data, electronic medical records, clinical data, and omics big data, and the chronic disease comorbidity network is represented by a graph G=(V, E), wherein V={v i |i∈{1,...,n}} represents a set of n disease nodes, E={e ij |i, j∈{1,...,m}} represents a set of m edges; the adjacency matrix A = [a ij ] is used to represent and store the target comorbidity network G, where a ij Represents the disease node v i and v j The interaction between them is equal to I(s i ,s j );

[0041] A representation unit, used for performing low-rank representation learning on the chronic disease comorbidity network to mine and output the community structure therein, wherein the low-rank representation learning adopts a symmetric non-negative matrix decomposition algorithm with relaxed constraints to obtain a low-rank representation matrix of the target network;

[0042] The partitioning unit is used to identify the community structure in the comorbidity network and discover the comorbidity pattern using the community partitioning mechanism.

[0043] In a third aspect, an embodiment of the present application provides an electronic device, characterized in that it includes:

[0044] processor;

[0045] a memory for storing instructions executable by the processor;

[0046] The processor is used to read the executable instructions from the memory and execute the executable instructions to implement the above method steps.

[0047] In a fourth aspect, an embodiment of the present application provides a computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program, and the computer program is used to execute the above method.

[0048] Compared with the prior art, the present invention has the following advantages: (1) Compared with the standard symmetric non-negative matrix factorization (SNMF) algorithm, the technical solution RSN of the present invention has an expanded feature space. Multiple feature matrices are used to represent the target undirected network, thereby maintaining its representation learning ability; (2) Relaxed constraints are used to achieve symmetry characterization. Compared with the standard non-negative matrix factorization (NMF) algorithm, by introducing relaxed symmetry constraints in the learning objectives, the designed model can well perceive the inherent symmetric structural characteristics of the network; (3) The spatial geometric structure constraints are effectively maintained, and the graph regularization technology that captures local topological features is introduced to effectively maintain the inherent geometric structure characteristics of the network. BRIEF DESCRIPTION OF THE DRAWINGS

[0049] By referring to the following drawings, the exemplary embodiments of the present invention can be more completely understood. The drawings are used to provide a further understanding of the embodiments of the present application, and constitute a part of the specification. Together with the embodiments of the present application, they are used to explain the present invention and do not constitute a limitation of the present invention. In the drawings, the same reference numerals generally represent the same parts or steps.

[0050] Figure 1 A flowchart of a method for identifying chronic disease comorbidity patterns based on relaxed constrained symmetric low-rank representation according to an exemplary embodiment of the present application;

[0051] Figure 2 This is a schematic diagram of a chronic disease comorbidity pattern recognition device based on relaxed constrained symmetric low-rank representation according to another exemplary embodiment of the present application. DETAILED DESCRIPTION

[0052] The exemplary embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. Although the exemplary embodiments of the present disclosure are shown in the accompanying drawings, it should be understood that the present disclosure can be implemented in various forms and should not be limited by the embodiments described herein. On the contrary, these embodiments are provided in order to enable a more thorough understanding of the present disclosure and to fully convey the scope of the present disclosure to those skilled in the art.

[0053] In the description of the present invention, it should be noted that the terms "center", "upper", "lower", "left", "right", "vertical", "horizontal", "inner", "outer", etc., indicating the orientation or positional relationship, are based on the orientation or positional relationship shown in the drawings, and are only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore cannot be understood as limiting the present invention. In addition, the terms "first", "second", and "third" are used for descriptive purposes only, and cannot be understood as indicating or implying relative importance.

[0054] In the description of the present invention, it should be noted that, unless otherwise clearly specified and limited, the terms "installed", "connected", and "connected" should be understood in a broad sense, for example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be a direct connection, or it can be indirectly connected through an intermediate medium, or it can be the internal communication of two components. For ordinary technicians in this field, the specific meanings of the above terms in the present invention can be understood according to specific circumstances.

[0055] In addition, the technical features involved in the different embodiments of the present invention described below can be combined with each other as long as they do not conflict with each other.

[0056] An embodiment of the present application provides a method for recognizing chronic disease comorbidity patterns based on relaxed constrained symmetric low-rank representation, which is described below in conjunction with the accompanying drawings.

[0057] Reference Figure 1 , which shows a chronic disease comorbidity pattern recognition method based on relaxed constrained symmetric low-rank representation provided by some embodiments of the present application, such as Figure 1 As shown, the method may include the following:

[0058] S101: Acquire medical data resources, wherein the medical data resources include data records of chronic disease conditions sorted out from tests, diagnoses, treatments, prescriptions, electronic medical records, and clinical data resources;

[0059] S102: Use medical data resources to build a chronic disease comorbidity network, where the chronic disease comorbidity network is used to describe the correlation and influence between diseases. It can be built by calculating the correlation between diseases or patients based on multi-source data. The multi-source data includes basic patient information as well as test, diagnosis, treatment, prescription data, electronic medical records, clinical data, and omics big data;

[0060] Specifically, the chronic disease comorbidity network is represented by a graph G = (V, E), where V = {v i |i∈{1,...,n}} represents a set of n disease nodes, E={e ij |i, j∈{1,...,m}} represents a set of m edges; where the adjacency matrix A = [a ij ] is used to represent and store the target comorbidity network G, a ij Represents the disease node v i and v j The interaction between them is equal to I(s i ,s j );

[0061] S103: Perform low-rank representation learning on the chronic disease comorbidity network to mine and output the community structure therein, wherein the low-rank representation learning uses a symmetric non-negative matrix factorization algorithm with relaxed constraints to obtain a low-rank representation matrix of the target network;

[0062] S104: Use the community division mechanism to identify the community structure in the comorbidity network and discover the comorbidity pattern. Specifically, the method for constructing the chronic disease comorbidity network includes the following steps:

[0063] (1) Sort out chronic disease records and extract data records from tests, diagnoses, treatments, prescriptions, electronic medical records, and clinical data;

[0064] (2) Count the frequency of disease occurrence. Calculate the frequency of each disease occurrence Q(s) from the data records. i ) and the frequency of disease combinations occurring together Q(s i ,s j ), where a disease i If it appears in a disease record, it is counted once. i and j If they co-occur in one disease record, they are counted once, where i, j = 1, 2, ..., n represents the number of diseases of interest;

[0065] Among them, the mutual information value between two diseases is calculated as follows:

[0066]

[0067] Among them, p(s i )=Q(s i ) / F,p(s j )=Q(s j ) / F respectively calculate disease s i and j The observation frequency, where F represents the total observation frequency, p(s i ,s j )=Q(s i ,s j ) / F calculates the number of diseases observed at the same time i and j frequency;

[0068] (3) Construct a comorbidity network to describe the complex impact relationships of chronic diseases.

[0069] Specifically, the symmetric low-rank representation method with relaxed constraints to characterize the comorbidity network includes the following steps:

[0070] (1) Design objective function;

[0071] Specifically, the objective function is as follows:

[0072]

[0073] stW=P, H=Q, P=Q; P, Q≥0

[0074] Among them, A is the adjacency matrix corresponding to the comorbidity network; the decision parameters W and H represent the learnable feature matrices, and P and Q are auxiliary parameter matrices, whose sizes are the same as the feature matrices W and H respectively; represents the F-norm of the matrix; λ is a positive constant representing the coefficient of the graph regularization term; Tr(·) is used to calculate the trace of the matrix; L represents the Laplacian matrix of the adjacency matrix A. The purpose of introducing auxiliary matrices P and Q is to separate the training and constraints of decision parameters W and H, and to achieve symmetric characterization under the relaxed constraints "W=P, H=Q, P=Q", as well as numerical characteristic characterization under the non-negative constraints "P, Q≥0";

[0075] (2) Constructing the augmented Lagrangian function;

[0076] Specifically, the ADMM alternating direction multiplier method is used to construct the augmented Lagrangian function:

[0077]

[0078] Among them, A is the adjacency matrix corresponding to the comorbidity network; the decision parameters W and H represent the feature matrix, and P and Q are the auxiliary parameter matrices that are the same as the feature matrices W and H; represents the F norm of the matrix; λ is a positive constant, representing the coefficient of the graph regularization term; Tr(·) calculates the trace of the matrix; L represents the Laplace matrix of the adjacency matrix A; ○ represents the Hadamard product of the matrix; α, β and θ correspond to the weights of the three penalty terms, respectively, controlling their proportional sizes. The non-negativity of P and Q can be achieved by mapping their value ranges to non-negative space;

[0079] The update rules for parameters are derived. In the ADMM optimization framework, let t and t+1 represent the current and updated iteration states of the parameters respectively. All parameters are iteratively updated using the update rules until the convergence condition is reached. The convergence condition is: (1) the absolute value of the difference between the objective function values ​​after two consecutive iterative trainings is less than 0.1; (2) or the number of iterations reaches 20. The update rules are as follows:

[0080]

[0081] in, represents the gradient, η is the learning rate in gradient ascent, outputs feature parameters W, H, P and Q, and Q is used as the community partition indicator matrix.

[0082] Specifically, the community partitioning mechanism is used to identify the community structure in the comorbidity network to discover the comorbidity pattern, including:

[0083] After training to obtain the characteristic parameters W, H, P and Q, Q is used as the indicator matrix for node-disease comorbidity module affiliation division, and the nodes in the comorbidity network are judged one by one and assigned to the comorbidity module. The assignment rules are as follows:

[0084]

[0085] Among them, V represents the node set in the comorbidity network, v j represents any node; C s is the comorbidity combination s, s∈{1,2,…,K}, K is the total number of comorbidity combinations; q jk is an element in the low-rank matrix Q, representing node v j The probability of being classified into the comorbidity combination k, q js is the maximum value among them, indicating that node v j The probability of being classified into the comorbidity combination s is the highest. After the comorbidity combinations of all nodes in the target comorbidity network are divided, all possible comorbidity combinations identified can be obtained. C = {C1, C2, …, C K}.

[0086] Compared with the prior art, the present invention has the following advantages: (1) Compared with the standard symmetric non-negative matrix factorization (SNMF) algorithm, the technical solution RSN of the present invention has an expanded feature space. Multiple feature matrices are used to represent the target undirected network, thereby maintaining its representation learning ability; (2) Relaxed constraints are used to achieve symmetry characterization. Compared with the standard non-negative matrix factorization (NMF) algorithm, by introducing relaxed symmetry constraints in the learning objectives, the designed model can well perceive the inherent symmetric structural characteristics of the network; (3) The spatial geometric structure constraints are effectively maintained, and the graph regularization technology that captures local topological features is introduced to effectively maintain the inherent geometric structure characteristics of the network.

[0087] In the above embodiment, a method is provided, and correspondingly, the present application also provides a device. The device provided in the embodiment of the present application can implement the above method, and the device can be implemented by software, hardware, or a combination of software and hardware. For example, the device may include integrated or separate functional modules or units to perform the corresponding steps in the above methods.

[0088] In some implementation modes of the embodiments of the present application, the device provided by the embodiments of the present application is based on the same inventive concept as the method provided by the aforementioned embodiments of the present application and has the same beneficial effects.

[0089] like Figure 2 As shown, the device 20 may include:

[0090] An acquisition unit 201 is used to acquire medical data resources, wherein the medical data resources include data records of chronic disease conditions sorted out from inspections, diagnoses, treatments, prescriptions, electronic medical records, and clinical data resources;

[0091] A construction unit 202 constructs a chronic disease comorbidity network using medical data resources, wherein the chronic disease comorbidity network is used to describe the correlation and influence between diseases, and can be constructed by calculating the correlation between diseases or patients based on multi-source data, wherein the multi-source data includes basic patient information and test, diagnosis, treatment, prescription data, electronic medical records, clinical data, and omics big data;

[0092] Specifically, the chronic disease comorbidity network is represented by a graph G = (V, E), where V = {v i |i∈{1,…,n}} represents a set of n disease nodes, E={e ij |i,j∈{1,…,m}} represents a set of m edges; the adjacency matrix A=[a ij ] is used to represent and store the target comorbidity network G, where a ij Represents the disease node v i and v j The interaction between them is equal to I(s i ,s j );

[0093] A representation unit 203 is used to perform low-rank representation learning on the chronic disease comorbidity network to mine and output the community structure therein, wherein the low-rank representation learning adopts a symmetric non-negative matrix decomposition algorithm with relaxed constraints to obtain a low-rank representation matrix of the target network;

[0094] The partitioning unit 204 is used to identify the community structure in the comorbidity network and discover the comorbidity pattern by using the community partitioning mechanism.

[0095] An embodiment of the present application also provides an electronic device corresponding to the method provided in the above embodiment. The device may be an electronic device used for a server, such as a server, including an independent server and a distributed server cluster, etc., to execute the above method; the electronic device may also be an electronic device used for a client to execute the above method.

[0096] The electronic device includes: a processor, a memory, a bus and a communication interface, wherein the processor, the communication interface and the memory are connected via a bus; the memory stores a computer program that can be run on the processor, and the processor executes the aforementioned method of the present application when running the computer program.

[0097] The memory may include high-speed random access memory (RAM), and may also include non-volatile memory, such as at least one disk memory. The communication connection between the system network element and at least one other network element is realized through at least one communication interface (which may be wired or wireless), and the Internet, wide area network, local area network, metropolitan area network, etc. may be used.

[0098] The bus may be an ISA bus, a PCI bus or an EISA bus, etc. The bus may be divided into an address bus, a data bus, a control bus, etc. The memory is used to store a program, and the processor executes the program after receiving an execution instruction. The method disclosed in any implementation of the above-mentioned embodiment of the present application may be applied to the processor or implemented by the processor.

[0099] The processor may be an integrated circuit chip with signal processing capabilities. In the implementation process, each step of the above method can be completed by an integrated logic circuit of hardware in the processor or an instruction in the form of software. The above processor may be a general-purpose processor, including a central processing unit (CPU), a network processor (NP), etc.; it may also be a digital signal processor (DSP), an application-specific integrated circuit (ASIC), a readily available programmable gate array (FPGA) or other programmable logic devices, discrete gates or transistor logic devices, discrete hardware components. The methods, steps and logic block diagrams disclosed in the embodiments of the present application can be implemented or executed. The general-purpose processor may be a microprocessor or the processor may also be any conventional processor, etc. The steps of the method disclosed in the embodiments of the present application can be directly embodied as being executed by a hardware decoding processor, or executed by a combination of hardware and software modules in a decoding processor. The software module may be located in a mature storage medium in the art such as a random access memory, a flash memory, a read-only memory, a programmable read-only memory or an electrically erasable programmable memory, a register, etc. The electronic device provided in the embodiment of the present application and the method provided in the embodiment of the present application are based on the same inventive concept and have the same beneficial effects as the method adopted, run or implemented.

[0100] The embodiments of the present application also provide a computer-readable medium corresponding to the method provided in the aforementioned embodiments, on which a computer program (ie, a program product) is stored. When the computer program is executed by a processor, the aforementioned method is executed.

[0101] It should be noted that examples of the computer-readable storage medium may also include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other optical or magnetic storage media, which are not listed here one by one.

[0102] It should be noted that the flowcharts and block diagrams in the accompanying drawings show the possible architecture, functions and operations of the systems, methods and computer program products according to multiple embodiments of the present application. In this regard, each box in the flowchart or block diagram can represent a module, a program segment or a part of a code, and the module, a program segment or a part of a code contains one or more executable instructions for implementing a specified logical function. It should also be noted that in some alternative implementations, the functions marked in the box can also occur in a different order from the order marked in the accompanying drawings. For example, two consecutive boxes can actually be executed substantially in parallel, and they can sometimes be executed in the opposite order, depending on the functions involved. It should also be noted that each box in the block diagram and / or flowchart, and the combination of the boxes in the block diagram and / or flowchart can be implemented with a dedicated hardware-based system that performs a specified function or action, or can be implemented with a combination of dedicated hardware and computer instructions.

[0103] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working processes of the systems, devices and units described above can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

[0104] In the several embodiments provided in the present application, it should be understood that the disclosed devices and methods can be implemented in other ways. The device embodiments described above are merely schematic. For example, the division of the units is only a logical function division. There may be other division methods in actual implementation. For example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some communication interfaces, and the indirect coupling or communication connection of devices or units can be electrical, mechanical or other forms.

[0105] The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0106] In addition, each functional unit in each embodiment of the present application may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit.

[0107] If the functions are implemented in the form of software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present application, or the part that contributes to the prior art or the part of the technical solution, can be embodied in the form of a software product, which is stored in a storage medium and includes several instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the methods described in each embodiment of the present application. The aforementioned storage medium includes: various media that can store program codes, such as a USB flash drive, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk.

[0108] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present application, rather than to limit it. Although the present application has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or replace some or all of the technical features therein by equivalents. These modifications or replacements do not deviate the essence of the corresponding technical solutions from the scope of the technical solutions of the embodiments of the present application, and they should all be included in the scope of the claims and specification of the present application.

Claims

1. A method for recognizing chronic disease comorbidity patterns based on relaxed constrained symmetric low-rank representation, characterized in that: include: Acquire medical data resources, wherein the medical data resources include data records of chronic disease conditions sorted out from tests, diagnoses, treatments, prescriptions, electronic medical records, and clinical data resources; A chronic disease comorbidity network is constructed using medical data resources, wherein the chronic disease comorbidity network is used to describe the correlation and influence between diseases, and can be constructed by calculating the correlation between diseases or patients based on multi-source data, wherein the multi-source data includes basic patient information and test, diagnosis, treatment, prescription data, electronic medical records, clinical data, and omics big data, wherein the chronic disease comorbidity network is represented by a graph G = (V, E), V = {v i |i∈{1,…,n}} represents a set of n disease nodes, E={e ij |i,j∈{1,…,m}} represents a set of m edges; where the adjacency matrix A=[a ij ] is used to represent and store the target comorbidity network G, a ij Represents the disease node v i and v j The interaction between them is equal to I(s i ,s j ); Performing low-rank representation learning and mining on the chronic disease comorbidity network to output the community structure therein, wherein the low-rank representation learning adopts a symmetric non-negative matrix decomposition algorithm with relaxed constraints to obtain a low-rank representation matrix of the target network; Use the community partitioning mechanism to identify the community structure in the comorbidity network and discover the comorbidity pattern.

2. The method according to claim 1, characterized in that The method for constructing the chronic disease comorbidity network includes: Sort out chronic disease records, and sort out data records from tests, diagnoses, treatments, prescriptions, electronic medical records, and clinical data; Count the frequency of disease occurrence, and count the frequency of each disease occurrence Q(s) from the above data records. i ) and the frequency of disease combination co-occurrence Q(s i ,s j ), where a disease i If it appears in a disease record, it is counted once. i and j If they appear together in one disease record, they are counted once, where i,j = 1, 2, …, n represents the number of diseases of interest; Calculate the mutual information value between two diseases as follows: Among them, p(s i )=Q(s i ) / F,p(s j )=Q(s j ) / F respectively calculate disease s i and j The observation frequency, where F represents the total observation frequency, p(s i ,s j )=Q(s i ,s j ) / F calculates the number of diseases observed at the same time i and j frequency; Constructing a comorbidity network to describe the complex impact relationships of chronic diseases.

3. The method according to claim 2, characterized in that The method of using a symmetric low-rank characterization method with relaxed constraints to characterize the comorbidity network includes: Design the objective function, which is as follows: Wherein, A is the adjacency matrix corresponding to the comorbidity network; decision parameters W and H represent learnable feature matrices, and P and Q are auxiliary parameter matrices, whose sizes are the same as the feature matrices W and H, respectively; represents the F-norm of the matrix; λ is a positive constant representing the coefficient of the graph regularization term; Tr(·) is used to calculate the trace of the matrix; L represents the Laplacian matrix of the adjacency matrix A. The purpose of introducing auxiliary matrices P and Q is to separate the training and constraints of decision parameters W and H, and to achieve symmetric characterization under the relaxed constraints "W=P,H=Q,P=Q", and numerical characteristic characterization under the non-negative constraints "P,Q≥0"; Construct the augmented Lagrangian function and use the ADMM alternating direction multiplier method to construct the augmented Lagrangian function: Wherein, A is the adjacency matrix corresponding to the comorbidity network; decision parameters W and H represent feature matrices, and P and Q are auxiliary parameter matrices that are the same as the feature matrices W and H; represents the F norm of the matrix; λ is a positive constant representing the coefficient of the graph regularization term; Tr(·) calculates the trace of the matrix; L represents the Laplace matrix of the adjacency matrix A; ° represents the Hadamard product of the matrix; α, β and θ correspond to the weights of the three penalty terms, respectively, controlling their proportional sizes. The non-negativity of P and Q can be achieved by mapping their value ranges to non-negative space; The update rules for parameters are derived. In the ADMM optimization framework, let t and t+1 represent the current and updated iteration states of the parameters respectively. All parameters are iteratively updated using the update rules until the convergence condition is reached. The convergence condition is: (1) the absolute value of the difference between the objective function values ​​after two consecutive iterative trainings is less than 0.1; (2) or the number of iterations reaches 20. The update rules are as follows: in, represents the gradient, η is the learning rate in gradient ascent, outputs feature parameters W, H, P and Q, and Q is used as the community partition indicator matrix.

4. The method according to claim 3, characterized in that The method of using the community partitioning mechanism to identify the community structure in the comorbidity network to discover the comorbidity pattern includes: After training to obtain the characteristic parameters W, H, P and Q, Q is used as the indicator matrix for node-disease comorbidity module affiliation division, and the nodes in the comorbidity network are judged one by one and assigned to the comorbidity module. The assignment rules are as follows: Where V represents the node set in the comorbidity network, v j represents any node; C s is the comorbidity combination s, s∈{1,2,…,K}, K is the total number of comorbidity combinations; q jk is an element in the low-rank matrix Q, representing node v j The probability of being classified into the comorbidity combination k, q js is the maximum value among them, indicating that node v j The probability of being classified into the comorbidity combination s is the highest. After the comorbidity combinations of all nodes in the target comorbidity network are divided, all possible comorbidity combinations identified can be obtained. C = {C1, C2, …, C K }.

5. A chronic disease comorbidity pattern recognition device based on relaxed constrained symmetric low-rank representation, characterized in that: include: An acquisition unit, used to acquire medical data resources, wherein the medical data resources include data records of chronic disease conditions sorted out from inspections, diagnoses, treatments, prescriptions, electronic medical records, and clinical data resources; A construction unit is provided for constructing a chronic disease comorbidity network using medical data resources, wherein the chronic disease comorbidity network is used to describe the correlation and influence between diseases, and can be constructed by calculating the correlation between diseases or patients according to multi-source data, wherein the multi-source data includes basic patient information and test, diagnosis, treatment, prescription data, electronic medical records, clinical data, and omics big data, and the chronic disease comorbidity network is represented by a graph G = (V, E), wherein V = {v i |i∈{1,…,n}} represents a set of n disease nodes, E={e ij |i,j∈{1,…,m}} represents a set of m edges; the adjacency matrix A=[a ij ] is used to represent and store the target comorbidity network G, where a ij Represents the disease node v i and v j The interaction between them is equal to I(s i ,s j ); A representation unit, used for performing low-rank representation learning on the chronic disease comorbidity network to mine and output the community structure therein, wherein the low-rank representation learning adopts a symmetric non-negative matrix decomposition algorithm with relaxed constraints to obtain a low-rank representation matrix of the target network; The partitioning unit is used to identify the community structure in the comorbidity network and discover the comorbidity pattern using the community partitioning mechanism.

6. An electronic device, characterized in that: include: processor; a memory for storing instructions executable by the processor; The processor is configured to read the executable instructions from the memory and execute the executable instructions to implement the method according to any one of claims 1 to 4.

7. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and the computer program is used to execute the method according to any one of claims 1 to 4 above.

Images