Information processing device, control method for information processing device, program, calculation device, and calculation method
Patent Information
- Application Number
- JP2019208332
- Authority / Receiving Office
- JP · JP
- Patent Type
- Patents
- Current Assignee / Owner
- Priority Date
- 2018-12-26
- Filing Date
- 2019-11-18
- Publication Date
- 2025-03-04
- Estimated Expiration
- 2039-11-18
AI Technical Summary
【0013】 本発明に係る情報処理装置及び算出装置によれば、夾雑物を分離するための前処理や、ピーク分割法などの演算処理をしなくても、スペクトル情報から被検物質について、精度の高い情報を得ることができる。
Smart Images

Figure 0007638620000003 
Figure 0007638620000004 
Figure 0007638620000005
Abstract
Description
[Technical field]
[0001] The present invention relates to an information processing device, a control method for an information processing device, a program, a calculation device, and a calculation method. [Background technology]
[0002] Spectral analysis is widely used as a method to determine the concentration and amount of specific components (hereinafter referred to as test substances) contained in various samples. In spectral analysis, the response to a certain stimulus is detected, and information about the components that make up the sample (spectral information) can be obtained based on the obtained signal. Spectral information includes the intensity of electromagnetic waves, including light, which characterize the stimulus and response, as well as temperature, mass, and the number of counts of fragments with a specific mass. Spectral analysis also includes using electron collisions as a stimulus to record the amount of fragments produced by decomposition relative to their mass, and obtaining information such as structure.
[0003] In some spectral analysis, the components are first separated using the differences in their three-dimensional size, charge, and hydrophilic / hydrophobic properties, and then irradiated with electromagnetic waves. This is called separation analysis. For example, in liquid chromatography (HPLC), the test substance is separated from other substances (hereafter referred to as impurities) by optimizing analytical conditions such as the column type, mobile phase type, temperature, and flow rate. The concentration and amount of the separated test substance can be determined by measuring the spectrum. In addition, if it is difficult to separate the test substance from the impurities, pretreatment may be performed to remove some of the impurities in advance, or the separation conditions may be optimized. If separation from the impurities is not possible even after optimizing the pretreatment and separation conditions, peak division by calculation is attempted.
[0004] Conventional peak division methods include a method of setting a baseline (FIG. 5), a method of dividing the peaks vertically using the minimum value between the peaks (FIG. 6), and a method of fitting an appropriate function such as the Gaussian function described in Patent Documents 1 and 2 using the least squares method to divide the peaks.
[0005] Here, HPLC is often used to analyze samples derived from living organisms. However, samples derived from living organisms, such as urine and blood, contain many impurities, and may contain unknown impurities derived from ingested substances, so operators who are skilled in the consideration of separation conditions, pretreatment, peak division methods, etc., are required to separate the test substance from the impurities.
[0006] In addition, there are many cases where samples contain a large number of impurities, such as in the analysis of pesticide residues in food and in the analysis of impurities in the environment. For this reason, there has been a strong demand for a method that allows even beginners to easily and accurately analyze the target substance in a sample of impurities without the need for pretreatment. [Prior art documents] [Patent documents]
[0007] [Patent Document 1] Japanese Patent Application Publication No. 6-324029 [Patent Document 2] JP 2006-177980 A [Patent Document 3] JP 2018-152000 A Summary of the Invention [Problem to be solved by the invention]
[0008] As described above, conventionally, in order to obtain information on a test substance from spectral information, pretreatment for separating impurities and calculation processing such as peak division methods were required. [Means for solving the problem]
[0009] In order to solve the above problem, an information acquisition device according to one embodiment of the present invention has an information acquisition means for acquiring quantitative information of a test substance estimated by inputting spectral information of a sample containing the test substance and impurities into a learning model, and the learning model is generated based on learning spectral information generated based on the spectral information of the test substance.
[0010] Furthermore, a control method for an information processing device according to another aspect of the present invention includes an information acquisition step of acquiring quantitative information of a test substance estimated by inputting spectral information of a sample containing the test substance and impurities into a learning model, and the learning model is generated based on learning spectral information generated based on the spectral information of the test substance.
[0011] In order to solve the above problem, a calculation device according to another aspect of the present invention is characterized in having a receiving unit that receives spectral information of a sample containing a test substance, an acquisition unit that acquires a learning model related to the spectral information of the test substance, and a calculation unit that calculates quantitative information of the test substance based on the spectral information of the sample and the learning model.
[0012] A calculation method according to still another aspect of the present invention includes a receiving step of receiving spectral information of a sample including a test substance, and an acquiring step of acquiring a learning model related to the spectral information of the test substance. The method further comprises a calculation step of calculating quantitative information of the test substance based on the spectral information of the sample and the learning model. Effect of the Invention
[0013] According to the information processing device and calculation device of the present invention, highly accurate information about a test substance can be obtained from spectral information without preprocessing for separating impurities or calculation processing such as peak division. [Brief description of the drawings]
[0014] [Figure 1] 1 is a schematic block diagram of an information processing device according to a first embodiment of the present invention. [Diagram 2] FIG. 2 is a diagram showing an example of a flowchart of a processing procedure for generating a learning model in the first embodiment of the present invention. [Diagram 3] FIG. 1 is a diagram showing an example of a flowchart of a processing procedure for acquiring quantitative information of a test substance in the first embodiment of the present invention. [Figure 4] FIG. 1 is a diagram showing an example of a flowchart of a process for calculating the amount of a test substance in a sample in the first embodiment of the present invention. [Diagram 5] FIG. 1 is a diagram showing an example of a conventional peak division method for establishing a baseline. [Figure 6] FIG. 1 is a diagram showing an example of a conventional peak splitting method in which a peak is split vertically by utilizing a minimum value between peaks. [Figure 7] FIG. 11 is a schematic block diagram of a calculation device according to a second embodiment of the present invention. [Figure 8] FIG. 11 is a diagram showing an example of a flowchart of a process for calculating the amount of a test substance in a sample in embodiment 2 of the present invention. [Figure 9] FIG. 1 is a diagram showing an example of a chromatogram including three test substances in Example 1 of the present invention. [Figure 10] FIG. 1 is a diagram showing an example of a chromatogram in which a plurality of different components are added in addition to three types of test substances in Example 1 of the present invention. [Figure 11] FIG. 1 is a diagram showing the results of predicting peak heights of test substances using a learning model in Example 1 of the present invention. [Figure 12] FIG. 11 is a diagram showing a method for dividing spectrum information in Example 2 of the present invention. [Figure 13] FIG. 11 is a diagram showing the results of predicting peak heights of test substances using a learning model in Example 2 of the present invention. [Figure 14] FIG. 11 is a diagram showing the results of predicting peak heights of test substances using a learning model in Example 3 of the present invention. [Figure 15] FIG. 13 is a diagram showing an example of a learning model for predicting related information from a peak height in Example 4 of the present invention. [Figure 16] FIG. 13 is a diagram showing an example of a result of predicting related information from a peak height in Example 4 of the present invention. [Figure 17] FIG. 11 is a diagram showing the prediction results of peak heights of test substances using a learning model in Example 5 of the present invention. [Figure 18]FIG. 6 shows the predicted peak heights of test substances using the conventional method (FIG. 5). [Figure 19] FIG. 1 shows a chromatogram obtained by HPLC analysis in Example 6 of the present invention. [Figure 20] FIG. 13 is a diagram showing the prediction results of peak heights of test substances using a learning model in Example 6 of the present invention. [Figure 21] FIG. 1 shows a chromatogram obtained in Example 7 of the present invention. [Figure 22] FIG. 13 is a diagram showing the prediction results of peak heights of test substances using a learning model in Example 7 of the present invention. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS
[0015] First Embodiment Hereinafter, embodiments of the present invention will be described with reference to the drawings. However, the scope of the present invention is not limited to the embodiments described below.
[0016] First, in describing this embodiment, terminology will be explained.
[0017] (sample) The sample in this embodiment is a mixture containing a plurality of types of compounds. In this embodiment, the sample is a mixture containing a test substance and other substances (impurities), and is not particularly limited as long as it is a mixture. In addition, the components of the mixture do not need to be specified, and unknown components may be contained. For example, the sample may be a mixture of blood, urine, saliva, etc. derived from a living body, or may be food or drink. Analysis of a sample derived from a living body contains clues to know the nutrition and health condition of the sample donor, so the analysis is valuable both medically and nutritionally. For example, urinary vitamin B3 is involved in the metabolism of carbohydrates, lipids, and proteins, and energy production, so measurement of its urinary metabolite, N1-methyl-2-pyridone-5-carboxamide, is useful for nutritional guidance for maintaining health.
[0018] (Test substance) The test substance in this embodiment is one or more known components contained in a sample, such as at least one selected from the group consisting of proteins, DNA, viruses, fungi, water-soluble vitamins, fat-soluble vitamins, organic acids, fatty acids, amino acids, sugars, pesticides, and environmental hormones.
[0019] For example, if one wishes to know the amount of nutrients, test substances that can be used include water-soluble vitamins such as thiamine (vitamin B1), riboflavin (vitamin B2), vitamin B3 metabolites N1-methylnicotinamide and N1-methyl-2-pyridone-5-carboxamide, vitamin B6 metabolites 4-pyridoxic acid and N1-methyl-4-pyridone-3-carboxamide, pantothenic acid (vitamin B5), pyridoxine (vitamin B6), biotin (vitamin B7), pteroylmonoglutamic acid (vitamin B9), cyanocobalamin (vitamin B12), and ascorbic acid (vitamin C); amino acids such as L-tryptophan, lysine, methionine, phenylalanine, threonine, valine, leucine, isoleucine, and L-histidine; and minerals such as sodium, potassium, calcium, magnesium, and phosphorus.
[0020] (Quantitative information) In this embodiment, the quantitative information may be the amount of the test substance contained in the sample, the concentration of the test substance contained in the sample, the presence or absence of the test substance in the sample, etc. In addition, other quantitative information of the test substance may be the concentration or ratio of the amount of the test substance contained in the sample to a reference amount of the test substance, or the ratio of the amount or concentration of the test substance contained in the sample, etc.
[0021] (Spectral information) The spectral information in this embodiment is at least one type selected from the group consisting of a chromatogram, a photoelectron spectrum, an infrared absorption spectrum (IR spectrum), a nuclear magnetic resonance spectrum (NMR spectrum), a fluorescence spectrum, an X-ray fluorescence spectrum, an ultraviolet / visible absorption spectrum (UV / Vis spectrum), a Raman spectrum, an atomic absorption spectrum, a flame emission spectrum, an optical emission spectrum, an X-ray absorption spectrum, an X-ray diffraction spectrum, a paramagnetic resonance absorption spectrum, an electron spin resonance spectrum, a mass spectrum, and a thermal analysis spectrum.
[0022] (Information processing system, information processing device) Next, an information processing system according to this embodiment will be described with reference to Fig. 1. Fig. 1 is a diagram showing the overall configuration of an information processing system including an information processing device according to this embodiment.
[0023] The information processing system includes an information processing device 10, a database 22, and an analysis device 23. The information processing device 10 and the database 22 are connected to each other so as to be able to communicate with each other via a communication means. In this embodiment, the communication means is configured by a LAN (Local Area Network) 21. The information processing device 10 and the analysis device 23 are connected by a communication means conforming to a standard such as USB (Universal Serial Bus). The LAN may be a wired LAN or a wireless LAN, or may be a WAN (Wide Area Network). The USB may be a LAN.
[0024] The database 22 manages the spectral information acquired by the analysis performed by the analysis device 23. The database 22 also manages a learning model (trained model) generated by a learning model generation unit 42 (described later). The information processing device 10 acquires the spectral information and the learning model managed in the database 22 via the LAN 21.
[0025] (Learning model) The learning model in this embodiment is a regression learning model, and one generated by machine learning such as deep learning can be used. A learning model is a model constructed by learning a machine learning algorithm using teacher data and making appropriate predictions. There are various types of machine learning algorithms used in the learning model. For example, deep learning using a neural network can be used. A neural network is composed of an input layer, an output layer, and multiple hidden layers, and each layer is connected by a formula called an activation function. When using teacher data with labels (output corresponding to input), the coefficients of the activation function are determined so that the relationship between the input and the output is established. By determining the coefficients using multiple teacher data, a learning model that can predict the output for the input with high accuracy can be generated.
[0026] (Analyzer) The analysis device 23 is a device for analyzing samples, test substances, etc. The analysis device 23 corresponds to an example of an analysis means. As described above, in this embodiment, the information processing device 10 and the analysis device 23 are connected so as to be able to communicate with each other. However, the analysis device 23 may be provided inside the information processing device 10, or the information processing device 10 may be provided inside the analysis device 23. Furthermore, the analysis result (spectral information) may be transferred from the analysis device 23 to the information processing device 10 via a recording medium such as a non-volatile memory.
[0027] The analytical device 23 in this embodiment is not limited as long as it can acquire spectral information, and a device using a chemical analysis method or a physical analysis method can be used. In this embodiment, the device using a chemical analysis method uses at least one method selected from the group consisting of, for example, chromatography such as liquid chromatography and gas chromatography, and capillary electrophoresis. In this embodiment, the device using a physical analysis method uses at least one method selected from the group consisting of, for example, photoelectron spectroscopy, infrared absorption spectroscopy, nuclear magnetic resonance spectroscopy, fluorescence spectroscopy, X-ray fluorescence spectroscopy, visible / ultraviolet absorption spectroscopy, Raman spectroscopy, atomic absorption spectroscopy, flame emission spectroscopy, emission spectroscopy, X-ray absorption spectroscopy, X-ray diffraction, electron spin resonance spectroscopy using paramagnetic resonance absorption, mass spectrometry, and thermal analysis.
[0028] For example, liquid chromatography is equipped with a mobile phase container, a liquid delivery pump, a sample injection unit, a column, a detector, and an A / D converter. Detectors that can be used include electromagnetic wave detectors that use ultraviolet light, visible light, infrared light, etc., as well as electrochemical detectors and ion detectors. In this case, the spectral information obtained is the output intensity from the detector versus time.
[0029] The information processing device 10 includes, as its functional configuration, a communication IF 31, a ROM 32, a RAM 33, a storage unit 34, an operation unit 35, a display unit 36, and a control unit 37.
[0030] The communication IF (Interface) 31 is realized by, for example, a LAN card and a USB interface card. The communication IF 31 manages communication between the information processing device 10 and external devices (for example, the database 22 and the analysis device 23) via the LAN 21 and the USB. The ROM (Read Only Memory) 32 is realized by a non-volatile memory or the like, and stores various programs and the like. The RAM (Random Access Memory) 33 is realized by a volatile memory or the like, and temporarily stores various information. The storage unit 34 is realized by, for example, a HDD (Hard Disk Drive) or the like, and stores various information. The operation unit 35 is realized by, for example, a keyboard, a mouse, or the like, and inputs instructions from the user into the device. The display unit 36 is realized by, for example, a display, or the like, and displays various information to the user. The operation unit 35 and the display unit 36 provide functions as a GUI (Graphical User Interface) under the control of the control unit 37.
[0031] (Control unit) The control unit 37 is realized by, for example, at least one CPU (Central Processing Unit) or the like, and controls the processing in the information processing device 10. The control unit 37 includes, as its functional configuration, a spectrum information acquisition unit 41, a learning model generation unit 42, a learning model acquisition unit 43, an estimation unit 44, an information acquisition unit 45, and a display control unit 46.
[0032] (Spectral information acquisition unit) The spectral information acquisition unit 41 acquires the analysis results of a sample containing at least a test substance and impurities, specifically, the spectral information of the sample, from the analysis device 23. Note that the spectral information of the sample may be acquired from the database 22 in which the analysis results are stored in advance. Similarly, the spectral information of the test substance is acquired. This spectral information of the test substance is the spectral information when a single test substance is present. The spectral information acquisition unit 41 then outputs the acquired spectral information of the sample to the estimation unit 44. The spectral information of the test substance is also output to the learning model generation unit 42. In the following, a case where the material contains multiple types of test substances will be described, but there may be only one type of test substance.
[0033] (divided part) The spectral information acquisition unit 41 according to this embodiment may further include a division unit (not shown) that divides the spectral information of the sample. The division unit can divide the spectral information of the sample containing the test substance and impurities into spectral information for each test substance.
[0034] One example of a method for dividing the spectral information of a sample into multiple parts is a method in which the user specifies the division points of the spectral information. For example, there is a method for specifying the range of the spectrum for each test substance, or a method for specifying the center of the spectrum for each test substance and extracting a certain range around the center. In the method for specifying the center, there is a method for setting the range to be extracted in advance, and a method for automatically determining the range to be extracted depending on the position of the spectral information (retention time if the spectral information is a chromatogram). For example, when obtaining a chromatogram by liquid chromatography, the longer the retention time, the more likely the chromatogram is to become broad, so it is better to widen the range to be extracted. In other words, the range to be extracted from the chromatogram can be changed according to the retention time of the chromatogram. For example, when dealing with a test substance whose peak appears at a position where the retention time of the chromatogram is long, it is better to widen the range to be extracted from the chromatogram.
[0035] Another example of a method for dividing the spectral information of a sample into multiple parts is a method in which the division points of the spectral information are automatically determined. For example, a threshold is set for the spectral information, and if the threshold is exceeded, it can be determined that the spectral information of the test substance is present there. For example, if the spectral information has a Gaussian distribution, the midpoint of the two points that exceed the threshold becomes the center of the peak of the spectral information. From there, the standard deviation (σ) can be obtained from the peak height and full width at half maximum. In this case, it is advisable to set the extraction range from the peak center based on σ. The extraction range is preferably about ±2σ to ±4σ.
[0036] Concerning the range for extracting spectral information, it is highly likely that the optimal range will vary depending on the type of test substance. Therefore, after determining the extraction range, a learning model can be created in the generation unit described later, and the extraction range can be corrected after checking the accuracy. Specifically, after determining the extraction range, spectral information is extracted within that range, and a learning model is created in the generation unit described later. The extracted spectral information (A) plus noise (B) is input to the created learning model, and the correlation coefficient between output C and the information in A (peak height, peak area, etc.) is calculated. The extraction range can be changed, the correlation coefficient confirmed using the same method, and it can be determined whether the extraction range should be widened or narrowed based on the change in the correlation coefficient.
[0037] Furthermore, machine learning may be used to divide the spectral information. By learning division points for various pieces of spectral information, it is possible to divide the spectral information for each test substance.
[0038] (Learning model generation unit) The learning model generation unit 42 generates teacher data using the spectral information of the test substance acquired by the spectral information acquisition unit 41. Then, the learning model generation unit 42 executes deep learning using the teacher data to generate a learning model. A detailed explanation of the generation of teacher data and the generation of a learning model will be given later. Then, the learning model generation unit 42 outputs the generated learning model to the learning model acquisition unit 43. Note that the learning model generation unit 42 may output the generated learning model to the database 22.
[0039] Here, an example of a method for generating a learning model will be described. First, where information on each test substance is included is extracted from the analysis results including the test substance. The method of extraction is the same as the means of division in the division unit described above. Therefore, if the spectral information has been divided in advance in the division unit, it is not necessary to extract it again. If it has not been divided, a spectrum for each test substance is extracted from the spectral information by the same means as the division unit. Data is generated by increasing or decreasing the amount of information on each test substance based on the extracted information (peak position) of each test substance. For example, if there are three peaks of test substances A, B, and C, a peak shape is generated in which the peak height of test substance A is 0 times, 0.2 times, 0.4 times, 0.6 times, 0.8 times, 1.2 times, 1.4 times, 1.6 times, and 1.8 times the original peak height. Similarly, peak shapes are generated for test substances B and C. In this example, 10 types of peak shapes are generated for each test substance. It is preferable that the generated peak shape is in a range that includes the expected amount of the test substance. Furthermore, the more peak shapes that are generated, the more accurate the learning model will be.
[0040] From here, there are two main ways to generate a learning model.
[0041] One method is to use training data that combines the peaks of each test substance that have been generated. In the above example, there are three test substances A, B, and C, so there are a total of 10 x 10 x 10, or 1,000 possible combinations. An arbitrary waveform generated using random numbers is added to this and used as input data. The peak heights of each test substance used in the input data are prepared as output data, and these are used as training data. In this case, only one learning model is generated.
[0042] The other method is to generate a learning model for each peak of each test substance. The output is the height of the peak prepared for each test substance, and the input data is prepared by adding an arbitrary waveform generated by random numbers to the peak shape prepared for each test substance. These are used as the training data. In this case, one learning model is generated for each test substance, so in the above example, three learning models are generated.
[0043] In either method, the arbitrary waveform generated by the random numbers is based on the information obtained by each analysis method. For example, in liquid chromatography, the spectrum information obtained is often Gaussian distributed. In that case, it is preferable to add up multiple Gaussian curves in which the peak height, median, and standard deviation are determined by random numbers, and further add the analysis result of the test substance, and use the result as input and the analysis result of the test substance as output to obtain the teaching data. The output may be the information to be obtained, for example, only the peak height. More specifically, the number of additions is preferably a number that cannot be separated on the chromatogram and may result in overlapping peaks, but usually about 2 to 8 are sufficient for one test substance. If the number exceeds 8, it becomes difficult to predict the shape of the test substance peak, and the quantitative accuracy decreases. If the number is 3 or less, accurate quantitative analysis cannot be performed for a chromatogram with insufficient separation. Preferably, 3 to 6, and even better, 4 to 5. The shape of the arbitrary waveform is a Gaussian function shown in Equation 1.
[0044]
number
[0045] Here, a is determined by random numbers within a range of 0 to x% of the expected peak height of the test substance, and b is determined by random numbers within a range of x% of the trimmed range. For example, when the range of ±3σ is trimmed from the center of the peak of the test substance, b can be any value within a range of -8σ×x% to +8σ×x%. x is 100 to 300, preferably 100 to 250, and even better 100 to 200. c is preferably determined by random numbers within a range of 0.1 to 10 times, more preferably 0.2 to 8 times, and even more preferably 0.5 to 5 times the standard deviation of the test substance peak. After adding multiple arbitrary waveforms for each test substance, the waveforms are combined with the desired value (for example, the peak heights of multiple test substances) to form a set of training data. Approximately 1000 to 2000 sets of training data are prepared, and similar training data is also created for peaks with different contents in the same combination of test substances. It is advisable to prepare training data of about 5,000 to 20,000 samples for about 5 to 10 test substances with different concentrations.
[0046] Furthermore, the arbitrary waveform generated by random numbers may be prepared to match the position of the spectrum information to be added. For example, in liquid chromatography, peaks with later retention times tend to be broader. If the standard deviation of the arbitrary waveform generated by random numbers is made a function of time, and the standard deviation becomes larger as the time becomes later, a more accurate learning model can be created.
[0047] (Learning model acquisition part) The learning model acquisition unit 43 acquires the learning model generated by the learning model generation unit 42. When the learning model is stored in the database 22, the learning model acquisition unit 43 acquires the learning model from the database 22. Then, the learning model acquisition unit 43 outputs the acquired learning model to the estimation unit 44.
[0048] (Estimation Department) The estimation unit 44 inputs the spectral information of the sample acquired by the spectral information acquisition unit 41 into the learning model acquired by the learning model acquisition unit 43, thereby causing the learning model to estimate quantitative information of the test substance contained in the sample. Then, the estimation unit 44 outputs the estimated quantitative information to the information acquisition unit 45. The estimation unit 44 corresponds to an example of an estimation means that estimates quantitative information of the test substance by inputting the spectral information of the sample into the learning model.
[0049] The information acquiring unit 45 acquires quantitative information estimated by the learning model. That is, the information acquiring unit 45 corresponds to an example of an information acquiring means for acquiring quantitative information of the test substance estimated by inputting the spectral information of a sample including the test substance and impurities into the learning model. Then, the information acquiring unit 45 outputs the acquired quantitative information to the display control unit 46.
[0050] (Display control unit) The display control unit 46 causes the display unit 36 to display the quantitative information acquired by the information acquisition unit 45. The display control unit 46 corresponds to an example of a display control means.
[0051] At least some of the units included in the control unit 37 may be realized as independent devices. Also, each unit may be realized as software that realizes the function. In this case, the software that realizes the function may run on a server via a network such as the cloud. In this embodiment, each unit is realized by software in a local environment.
[0052] 1 is merely an example. For example, the storage unit 34 of the information processing device 10 may have the functions of the database 22, and the storage unit 34 may hold various types of information.
[0053] Next, the processing procedure in this embodiment will be described with reference to Fig. 2. Fig. 2 is a flowchart of the processing procedure related to generation of a learning model.
[0054] (S201) (Analyze the test substance alone) In step S201, the analysis device 23 analyzes the test substance alone and acquires the spectral information of the test substance. The analysis conditions may be appropriately selected from the viewpoints of sensitivity, analysis time, and the like. In this case, the analysis device 23 performs the analysis by changing the concentration of the test substance in several ways. The number of concentrations required varies depending on the properties of the substance, and the like, but it is generally desirable to change the concentration at three or more points. The analysis of the test substance may be performed by analyzing each test substance individually or simultaneously. Then, the analysis device 23 outputs the acquired spectral information to the information processing device 10. The information processing device 10 receives the spectral information from the analysis device 23 and stores it in the RAM 33 or the storage unit 34. The spectral information acquisition unit 41 acquires the spectral information thus stored. Note that, as described above, the spectral information that is the analysis result may be stored in the database 22. In this case, the spectral information acquisition unit 41 acquires the spectral information from the database 22. The timing of the analysis of the test substance by the analysis device 23 may be any timing as long as it is performed before the generation of the teacher data in step S202.
[0055] (S202) (Generate training data) In step S202, the learning model generating unit 42 generates a plurality of teacher data using the spectral information of the test substance acquired by the spectral information acquiring unit 41. A method of generating teacher data will be specifically described. The teacher data is generated by adding an arbitrary waveform generated by random numbers to the spectral information of the test substance. For example, in liquid chromatography, the waveform indicated by the spectral information (chromatogram) is often a Gaussian distribution. Therefore, the learning model generating unit 42 generates a plurality of random noises by adding together a plurality of Gaussian curves (Gaussian functions) whose peak heights, medians, and standard deviations are determined by random numbers. Then, the learning model generating unit 42 generates a plurality of waveforms by adding together each of the plurality of random noises and the waveform indicated by the spectral information of the test substance. The plurality of waveforms thus generated are used as the spectral information (learning spectral information) of a virtual sample including the test substance and impurities. In other words, the plurality of generated spectral information is determined as input data constituting the teacher data. Furthermore, the learning model generation unit 42 determines the peak height (quantitative information) identified from the spectral information of the test substance on which the generated spectral information is based as the correct answer data constituting the teacher data. In this way, the learning model generation unit 42 generates a plurality of teacher data which are pairs of input data and correct answer data. Then, in step S201, since the learning model generation unit 42 acquires spectral information according to the concentration of the test substance, it generates a plurality of teacher data for each concentration. Note that, in the waveform of a chromatogram, the width of the peak tends to increase as the retention time increases, and therefore the learning model generation unit 42 may widen the width of the waveform to be generated.
[0056] Patent Document 3 discloses a method of performing machine learning by linking mass spectrum data of a specimen with the presence or absence of cancer. However, a large amount of training data is required to improve the accuracy of machine learning. Patent Document 3 provides 90,000 types of data as training data. In other words, machine learning can accurately analyze complex analysis results, but the drawback is that a large amount of training data must be prepared. In this embodiment, there is no need to prepare a large amount of training data, which is a drawback of machine learning, and therefore the burden on the user can be reduced.
[0057] Although the teacher data is generated in this manner, the teacher data may be generated by analyzing a plurality of samples with the analysis device 23 to obtain spectral information of the samples for learning together with quantitative information of the test substance. Also, the spectral information of a virtual sample may be generated by a method different from the above-mentioned method.
[0058] (S203) (Generate learning model) In step S203, the learning model generating unit 42 generates a learning model by performing machine learning according to a predetermined algorithm using the multiple teacher data generated for each concentration in step S202. In this embodiment, a neural network is used as the predetermined algorithm. The learning model generating unit 42 generates a learning model that estimates quantitative information of the test substance contained in the sample based on the input of the spectrum information of the sample by having the neural network learn using the multiple teacher data. Note that since the learning method of the neural network is a well-known technology, detailed description will be omitted in this embodiment. In addition, for example, SVM (support vector machine) and DNN (deep neural network) may be used as the predetermined algorithm. As another algorithm, CNN (convolutional neural network) or the like may be used. When there are multiple types of test substances, a learning model is constructed for each substance. Then, the learning model generating unit 42 stores the generated learning model in the RAM 33, the storage unit 34, or the database 22.
[0059] In this manner, a learning model is generated that estimates quantitative information about a test substance contained in a sample based on the spectral information of the sample.
[0060] (S301) (Analyze the sample) In step S301, the analysis device 23 analyzes a target sample and acquires spectral information of the sample. The analysis conditions are the same as those in step S201 described above. Then, the analysis device 23 outputs the acquired spectral information to the information processing device 10. The information processing device 10 receives the spectral information from the analysis device 23 and stores it in the RAM 33 or the storage unit 34. The spectral information acquisition unit 41 acquires the spectral information thus stored. As described above, the spectral information that is the analysis result may be stored in the database 22. In this case, the spectral information acquisition unit 41 acquires the spectral information from the database 22. The timing at which the analysis device 23 analyzes the sample may be any timing as long as it is performed before the estimation of quantitative information in step S302.
[0061] (S302) (Estimating quantitative information) In step S302, the learning model acquisition unit 43 acquires a learning model stored in the RAM 33, the storage unit 34, or the database 22. Then, the estimation unit 44 inputs the spectrum information of the sample acquired in step S301 into the acquired learning model, thereby estimating quantitative information of the test substance contained in the sample. Furthermore, as necessary, the estimation unit 44 converts the estimated quantitative information into a format to be displayed on the display unit 36. The format to be displayed on the display unit 36 may be a concentration such as g / L or mol / L, or a ratio to a reference amount (standard amount). If the value estimated by the learning model is in these display formats, conversion is not necessary. Then, the information acquisition unit 45 acquires the estimated quantitative information from the estimation unit 44 and stores it in the RAM 33 or the storage unit 34.
[0062] In this way, by using the learning model obtained through machine learning, even if the peaks of the test substance and impurities cannot be completely separated, it is possible to obtain accurate quantitative information about the test substance without complex and advanced knowledge of analysis.
[0063] As a result, even an unskilled person can easily perform quantitative analysis of a test substance with high accuracy.
[0064] (S303) (Display quantitative information) In step S304, the display control unit 46 causes the display unit 36 to display quantitative information on the test substance contained in the sample, estimated by the learning model in step S302. At this time, the information may be displayed in a graph or table format (data analysis method).
[0065] Next, the processing in the calculation device will be described in more detail with reference to Fig. 4. Fig. 4 is a flowchart of the processing for calculating the amount of a test substance in a sample.
[0066] First, the analysis unit analyzes the test substance (step S401). The analysis conditions may be appropriately selected from the viewpoints of sensitivity, measurement time, etc. At this time, the concentration of the test substance may be changed in several ways to perform the measurement.
[0067] Next, the analysis unit analyzes the target sample containing the test substance (step S402). The measurement conditions are set to the same conditions as those in step S401.
[0068] Next, the generating unit generates teacher data using the measurement results of the test substance obtained in step S401. Using the generated teacher data, machine learning according to a predetermined algorithm is performed to construct a learning model (step S403). As a specific learning method, for example, a neural network or a support vector machine may be used as a general machine learning method. In addition, a deep learning method with multiple hidden layers such as a DNN (deep neural network) or a CNN (convolutional neural network) may be used.
[0069] Next, the calculation unit applies the learning model generated in step S403 to the measurement data obtained in step S402 to calculate the amount of the test substance (step S404). At this time, the amount is converted into a format to be displayed on the display unit. The format to be displayed on the display unit may be a concentration such as g / L or mol / L, or a ratio to a standard amount.
[0070] Next, the display unit displays the amount of the test substance calculated in step S404 above and presents it to the user (S405). At that time, the amount may be displayed in a graph or table format.
[0071] The present invention can also be realized by a process in which a program for implementing one or more of the functions of the above-described embodiments is supplied to a system or device via a network or a storage medium, and one or more processors in a computer of the system or device read and execute the program. The present invention can also be realized by a circuit (e.g., ASIC) for implementing one or more of the functions.
[0072] The present invention will be described in more detail below with reference to examples and comparative examples, but the present invention is not limited to the following examples.
[0073] Second Embodiment FIG. 7 is a schematic block diagram of a calculation device according to the second embodiment.
[0074] (Configuration of the calculation device) The calculation device 700 according to this embodiment includes a receiving unit 703 that receives an analysis result of a test substance alone for a sample 702 that contains a test substance 701, and an acquiring unit 704 that acquires a learning model related to the analysis result of the test substance. The calculation device 700 further includes a calculation unit 705 that calculates quantitative information about the test substance using the analysis result of the sample and the learning model.
[0075] In recent years, analytical methods using machine learning have been developed for complex analytical results. Patent Document 3 discloses a method of performing machine learning by linking mass spectrum data of a specimen with the presence or absence of cancer. However, a large amount of training data is required to improve the accuracy of machine learning. Patent Document 3 provides 90,000 types of data as training data. In other words, machine learning can accurately analyze complex analytical results, but the drawback is that a large amount of training data must be prepared.
[0076] In this embodiment, even if the peaks of the test substance and the peaks of the impurities are not completely separated, accurate peak information of the test substance can be obtained without complex and advanced knowledge of analysis by using a learning model obtained by machine learning. In addition, there is no need to prepare a large amount of training data, which is a drawback of machine learning.
[0077] Furthermore, in the past, when the peak of the test substance and the peak of the impurities could not be completely separated, work requiring complex and advanced knowledge of analysis was required, such as pre-processing and optimization of analytical conditions, etc. According to the calculation device of this embodiment, it is possible to obtain information on the test substance from the spectrum information without pre-processing for separating the impurities or calculation processing such as peak division method.
[0078] As a result, even an unskilled person can easily perform quantitative analysis of a test substance with high accuracy.
[0079] The calculation device according to this embodiment may further include an analysis unit 706 that analyzes samples. Even if the calculation device does not include an analysis unit, the analysis results of a separately prepared analysis unit 707 may be acquired, or the results may be acquired from a database (DB1) 708. In addition, the learning model may be generated by a generation unit 709, or the learning model may be acquired from a database (DB2) 710.
[0080] The calculation device according to this embodiment will be described in detail below.
[0081] (Reception Department) The reception unit receives information on the analysis results of the sample, specifically, the spectral information of the sample. In order to obtain the analysis results, if the calculation device according to this embodiment has an analysis unit described later, the spectral information of the sample can be acquired from the analysis unit. Furthermore, if the calculation device according to this embodiment does not have an analysis unit, a database (DB1 in the figure) in which the analysis results are stored in advance may be prepared, and the analysis results, i.e., the spectral information of the sample, may be acquired from the database. The database may be present in a storage device built into the calculation device, in a storage device external to the calculation device, or on a cloud via a network.
[0082] (Acquisition Department) The acquisition unit acquires a learning model based on the analysis results of the sample obtained by the reception unit. The learning model may be generated by a generation unit described later. Alternatively, a database (DB2 in the figure) that stores learning models of test substances in advance may be prepared, and the learning model may be obtained from the database. The database may be present in a storage device built into the calculation device, in a storage device external to the calculation device, or on a cloud via a network.
[0083] (Learning model) In this embodiment, the learning model can be one generated by machine learning such as deep learning.
[0084] A machine learning algorithm is constructed to learn from training data and make appropriate predictions, and is called a learning model here. There are various types of machine learning algorithms used in learning models. For example, deep learning using neural networks can be used. A neural network consists of an input layer, an output layer, and multiple hidden layers, and each layer is connected by a formula called an activation function. When using training data with labels (output corresponding to input), the coefficients of the activation function are determined so that the relationship between input and output is established. By determining the coefficients using multiple training data, it is possible to create a learning model that can predict the output for an input with high accuracy.
[0085] As the learning model in this embodiment, for example, information including spectrum information (spectrum information of the test substance) when a single test substance is present and random noise can be used. The random noise is not particularly limited, but a waveform combining multiple Gaussian functions can be used. In addition, the width of the waveform of the chromatogram may be widened in consideration of the tendency for the width of the peak to increase as the retention time increases.
[0086] (Spectral information) In this embodiment, the spectral information is at least one selected from the group consisting of a chromatogram, a photoelectron spectrum, an infrared absorption spectrum, a nuclear magnetic resonance spectrum, an X-ray fluorescence spectrum, a visible / ultraviolet absorption spectrum, and a Raman spectrum, or at least one selected from the group consisting of an atomic absorption spectrum, a flame emission spectrum, an emission spectroscopy spectrum, an X-ray absorption spectrum, an X-ray diffraction spectrum, a paramagnetic resonance absorption spectrum, a mass spectrum, and a thermal analysis spectrum.
[0087] (Calculation section) The calculation unit applies the learning model obtained in the acquisition unit to the analysis results of the sample, and calculates quantitative information about the test substance in the sample.
[0088] (Analysis Department) The analysis unit in this embodiment is various analyzers for obtaining the analysis results of the sample and the test substance. The analysis unit may be provided in the same computer as at least one of the other parts of the calculation device, or may be connected via a wireless or wired Internet network, or may transmit and receive the analysis results via a recording medium such as a non-volatile memory.
[0089] The analysis unit in this embodiment uses at least one technique selected from the group consisting of chromatography, capillary electrophoresis, photoelectron spectroscopy, infrared absorption spectroscopy, nuclear magnetic resonance spectroscopy, X-ray fluorescence spectroscopy, visible / ultraviolet absorption spectroscopy, and Raman spectroscopy, and also uses at least one technique selected from the group consisting of atomic absorption spectrometry, flame emission spectroscopy, optical emission spectroscopy, X-ray absorption spectroscopy, X-ray diffraction, paramagnetic resonance absorption, mass spectrometry, and thermal analysis.
[0090] There are various instruments used for analysis, including visible and ultraviolet absorption spectroscopy (UV / Vis spectroscopy), infrared absorption spectroscopy (IR spectroscopy), nuclear magnetic resonance spectroscopy (NMR spectroscopy), and Raman spectroscopy. Other instruments include fluorescence spectroscopy, atomic absorption spectroscopy, flame analysis, flame fluorescence spectroscopy, emission spectroscopy, emission X-ray analysis, X-ray diffraction, and fluorescence X-ray diffraction. In addition, there are fluorescence X-ray spectroscopy, paramagnetic resonance absorption spectroscopy, mass spectroscopy, thermal analysis, gas chromatography, liquid chromatography, and photoelectron spectroscopy.
[0091] For example, liquid chromatography is equipped with a mobile phase container, a liquid delivery pump, a sample injection unit, a column, a detector, and an A / D converter. Detectors include electromagnetic wave detectors that use ultraviolet, visible light, infrared, etc., as well as electrochemical detectors and ion detectors. In this case, the spectral information obtained is the output intensity from the detector versus time.
[0092] (Generation part) The generation unit creates teacher data based on the analysis results of the test substance, and generates a learning model by performing deep learning. The teacher data is generated by adding an arbitrary waveform generated by random numbers to the analysis results of the test substance. For example, in liquid chromatography, the spectrum information obtained is often Gaussian distributed. In that case, it is preferable to add multiple Gaussian curves in which the peak height, median, and standard deviation are determined by random numbers, and further add the analysis results of the test substance, and use the result as input to generate teacher data with the analysis results of the test substance as output. The output may be the information to be obtained, for example, only the peak height. More specifically, data is prepared by trimming the peak of the test substance at the center. The wider the trimming range, the higher the accuracy of quantification in the subsequent calculation unit, but the number of teacher data required to increase the accuracy increases. The trimming range is about 6 to 30 times the standard deviation (σ) of the peak of the test substance. It is preferable that it is 10 to 20 times the standard deviation (σ) of the peak of the test substance, and even better that it is 14 to 18 times. Next, an arbitrary waveform is added to the trimmed data. The number of added peaks is preferably a number that is not separable on the chromatogram and may result in overlapping peaks, but usually 2 to 8 is sufficient. If the number of added peaks exceeds 8, it becomes difficult to predict the shape of the peak of the test substance, and the quantitative accuracy may decrease. Furthermore, if the number of added peaks is 3 or less, accurate quantification may not be possible for insufficiently separated chromatograms. Preferably, 3 to 6 peaks are added, and even better, 4 to 5 peaks are added. The shape of the arbitrary waveform can be a Gaussian function shown in Equation 1.
[0093]
number
[0094] Here, a is determined by a random number within a range of 0 to x% of the expected peak height of the test substance, and b is determined by a random number within a range of x% of the trimmed range. For example, when the range of ±3σ is trimmed from the center of the peak of the test substance, b can be any value within a range of -8σ×x% to +8σ×x%. x is determined by a random number within the range. x is 100 to 300, preferably 100 to 250, and more preferably 100 to 200. c is preferably determined by a random number within a range of 0.1 to 10 times, more preferably 0.2 to 8 times, and even more preferably 0.5 to 5 times the standard deviation of the test substance peak. After adding a plurality of arbitrary waveforms, it can be set with the value to be obtained (for example, the peak height of the test substance) to form a set of teacher data. About 1000 to 2000 pieces of teacher data are prepared for one test substance peak, and similar teacher data is also created for peaks of the same test substance with different contents. It is advisable to prepare training data of about 5,000 to 20,000 test substances with different contents of about 5 to 10 types.
[0095] (Storage part) The calculation device according to the present embodiment may further include a storage unit in which the learning model is stored. In this case, the acquisition unit can acquire the learning model from the storage unit.
[0096] (Display control unit) The calculation device according to this embodiment may further include a display control unit that controls the display of the quantitative information.
[0097] (sample) The sample is not particularly limited as long as it is a mixture containing a plurality of types of compounds. In addition, the components of the mixture do not need to be specified, and unknown components may be included. For example, it may be a mixture of blood, urine, saliva, etc. derived from a living body, or it may be food or drink. Analysis of a sample derived from a living body contains clues to know the nutrition and health condition of the sample donor, so the analysis is valuable both medically and nutritionally. For example, urinary vitamin B3 is involved in the metabolism of carbohydrates, lipids, and proteins, and energy production, so the measurement of its urinary metabolite, N1-methyl-2-pyridone-5-carboxamide, is useful for nutritional guidance to maintain health.
[0098] (Test substance) The test substance in this embodiment is one or more known components contained in the sample, for example, at least one selected from the group consisting of proteins, DNA, viruses, fungi, water-soluble vitamins, fat-soluble vitamins, organic acids, fatty acids, amino acids, sugars, pesticides, and environmental hormones.
[0099] For example, if you want to know the amount of nutrients, the test substances include thiamine (vitamin B1), riboflavin (vitamin B2), and vitamin B3 metabolites N1-methylnicotinamide and N1-methyl-2-pyridone-5-carboxamide. Other water-soluble vitamins include N1-methyl-4-pyridone-3-carboxamide, pantothenic acid (B5), 4-pyridoxic acid (B6), biotin (B7), pteroylmonoglutamic acid (B9), and ascorbic acid (C). Other amino acids include L-tryptophan, lysine, methionine, phenylalanine, threonine, valine, leucine, isoleucine, and L-histidine. Other minerals include sodium, potassium, calcium, magnesium, and phosphorus, as well as folic acid, cyanocobalamin, and ascorbic acid.
[0100] (Quantitative information) It is at least one selected from the group consisting of the amount of the test substance contained in the sample, the concentration of the test substance contained in the sample, and the presence or absence of the test substance in the sample. Also, it is at least one selected from the group consisting of the concentration or amount ratio of the test substance contained in the sample to a reference amount, and the amount or concentration ratio of the test substance contained in the sample.
[0101] (Calculation method) The calculation method according to this embodiment includes at least the following steps. (1) A receiving step of receiving spectral information of a sample containing a test substance. (2) An acquisition step of acquiring a learning model regarding the spectral information of the test substance. (3) A calculation step of calculating quantitative information of the test substance based on the spectral information of the sample and the learning model.
[0102] Here, the spectrum information is exemplified by the information described above.
[0103] Furthermore, an analysis step may be included in which an analysis is performed to obtain spectral information of the sample, and the above-mentioned analysis is exemplified.
[0104] The method may further include a generating step of generating a learning model, and the acquiring step may include a step of acquiring the learning model generated in the generating step. The learning model may be information including spectrum information of the test substance and random noise. Here, the random noise may be a waveform obtained by combining a plurality of Gaussian functions.
[0105] Examples of the test substance include the substances described above. Examples of the quantitative information include the information described above.
[0106] The learning model may be generated by machine learning such as deep learning.
[0107] The calculation method may further include a display control step of controlling the display of the quantitative information.
[0108] An example of a calculation method according to an embodiment of the present invention will be described in detail below.
[0109] (Data analysis method) Next, the processing in the calculation device will be described in more detail with reference to Fig. 8. Fig. 8 is a flowchart of the processing for calculating the amount of a test substance in a sample.
[0110] (Step S801) First, the analysis unit analyzes the test substance alone (step S801). The analysis conditions may be appropriately selected from the viewpoints of sensitivity, measurement time, etc. In this case, the concentration of the test substance is changed in several ways and measured. The number of changes required varies depending on the properties of the substance, etc., but it is generally desirable to change the concentration at three or more points. When there are multiple types of test substances, it is desirable to measure each of them, but if the signals of the test substances are sufficiently separated, they may be measured simultaneously.
[0111] (Step S802) Next, the analysis unit analyzes the target sample containing the test substance (step S802). The measurement conditions are set to the same conditions as those in step S801.
[0112] Next, the generation unit generates multiple pieces of teacher data using the measurement results of the test substance alone obtained in step S801. Using the generated teacher data, machine learning according to a predetermined algorithm is performed to construct a learning model (step S803). As a specific learning method, for example, a neural network or a support vector machine may be used as a general machine learning method. In addition, a deep learning method with multiple hidden layers such as a deep neural network (DNN) or a convolutional neural network (CNN) may be used. When there are multiple types of test substances, a learning model is constructed for each substance.
[0113] (Step S803) Next, the calculation unit applies the learning model generated in step S3 to the measurement data obtained in step S802 to calculate the amount of the test substance (step S804). At this time, the amount is converted into a format to be displayed on the display unit. The format to be displayed on the display unit may be a concentration such as g / L or mol / L, or a ratio to a reference amount (standard amount).
[0114] (Step S804) Next, the display unit displays the amount of the test substance calculated in step S804 above and presents it to the user. At that time, the amount may be displayed in a graph or table format.
[0115] The present invention can also be realized by a process in which a program for implementing one or more of the functions of the above-described embodiments is supplied to a system or device via a network or a storage medium, and one or more processors in a computer of the system or device read and execute the program. The present invention can also be realized by a circuit (e.g., ASIC) for implementing one or more of the functions.
[0116] The present invention will be described in more detail below with reference to examples and comparative examples, but the present invention is not limited to the following examples. EXAMPLES
[0117] Example 1 Fig. 9 is a chromatogram including three test substances (501, 502, 503). A threshold value 504 of 0.1 was set, and the peak position of each test substance was extracted from the intersection of the threshold value and the chromatogram. For example, for peak 501, the median value 506 of the intersection point 505 of the threshold value and the chromatogram was determined to be the position of peak 501. The peak positions of peaks 502 and 503 were identified in the same manner. The respective peak positions were 100, 250, and 400.
[0118] For each peak, the range of the peak position ±75 was multiplied by a value from 0 to 1.8 (in increments of 0.2), and 10 types of waveforms with different peak heights were obtained for each test substance. Waveforms with different peak heights for test substances 501, 502, and 503 were combined to generate 10 x 10 x 10, or 1,000 waveforms. Hereinafter, this will be referred to as test substance data. For each test substance data, five normal distribution waveforms with median, standard deviation, and peak height set by random numbers were added to create sample data. The added normal distribution waveforms were determined by random numbers within the range of peak height (a in formula 1) 0.0 to 1.5, median (b in formula 1) 10 to 490, and standard deviation (c in formula 1) 10 to 100. 100 types of sample data were prepared for one test substance data. Each sample data and the peak height of the test substance data contained therein were paired to create 100,000 training data, and a learning model was generated by performing machine learning. The machine learning method used was a fully connected neural network, with the relu function and linear function as activation functions. The loss function was the mean squared error, and the optimization algorithm was Adam. Repeated calculations of about 20 epochs were required to obtain sufficient quantitative accuracy.
[0119] Figure 10 is an example of a chromatogram in which multiple different components are added in addition to the three test substances in Figure 9. The amounts of the three test substances have been adjusted to arbitrary amounts. This chromatogram in which multiple different components have been added in addition to the three test substances is applied to the learning model, and Figure 11 shows a comparison of the calculated peak heights of the test substances (predicted values) and the arbitrarily adjusted peak heights of the test substances (correct values). Predictions could be made with a correlation coefficient of about 0.7 to 0.9 against the correct values.
[0120] Example 2 As in Example 1, the chromatogram in FIG. 9 was divided into three types of test substances, and each was multiplied by a value from 0 to 1.8 (in increments of 0.2) to obtain test substance data with 10 different peak heights for each test substance. Sample data was obtained by adding four normal distribution waveforms in which the median, standard deviation, and peak height were set by random numbers to each test substance data. 1000 types of sample data were prepared for one test substance data. Each sample data and the peak height of the test substance data contained therein were paired to obtain 10,000 pieces of teacher data, and a learning model was generated by performing machine learning for each test substance. A fully connected neural network was used as the machine learning method, and a relu function and a linear function were used as the activation function. The mean square error was used as the loss function, and Adam was used as the optimization algorithm. Repeated calculations of about 100 epochs were required to obtain sufficient quantitative accuracy.
[0121] As in Example 1, a chromatogram was prepared by adding a number of different components in addition to three kinds of test substances containing arbitrary amounts. As shown in FIG. 12, the chromatogram was divided into the same regions as the test substance data (501, 502, 503). Then, the chromatograms (601, 602, 603) in the regions corresponding to each test substance were applied to the learning model created for each test substance, and the peak height of each test substance was calculated (predicted value). FIG. 13 shows a comparison of the peak heights of the test substances adjusted arbitrarily as the correct value. It was possible to predict with a correlation coefficient of about 0.8 to 0.9 with respect to the correct value.
[0122] Example 3 The standard deviation of the normal distribution waveform added to the test substance data in Example 2 was set as follows based on the trend of the chromatogram. Standard deviation = (random number 0 to 5) + time x 0.2 + 5 Other than that, the procedure was carried out as in Example 2, and the correct values were compared with the predicted values (FIG. 14). Predictions were possible with a correlation coefficient of about 0.9.
[0123] Example 4 In the first to third embodiments, the peak height was obtained as the final output, but here, a means for predicting related information from the peak height is implemented. It is known that the manufacturer of the test substances 501 to 503 can be estimated based on the ratio of their amounts. A learning model is prepared by machine learning using the amounts of the test substances 501 to 503 in the sample of each manufacturer and the manufacturer (five types A to E in this embodiment) as training data. FIG. 15 shows an example of a learning model using a fully connected neural network. 1101 is an input layer, which inputs information corresponding to the amounts of the test substances 501 to 503. 1102 is an intermediate layer. In this embodiment, there is one layer, but multiple layers may be prepared depending on the number of test substances and classifications. The number of nodes may also be changed as appropriate. 1103 is an output layer, and classifications A to E represent each manufacturer. The activation function connecting the intermediate layer to the output layer is a Softmax function, and the output is a probability. By applying the amounts of each test substance (calculated from the peak heights) obtained in Examples 1 to 3 to the learning model of this example, it was possible to display which manufacturer's sample had a high probability of containing the test substance. An example of the results is shown in FIG.
[0124] Example 5 First, an example will be described in which the above-mentioned data processing method is applied to simulation data in order to evaluate the effect of the method.
[0125] As the test substance data, 11 types of data with normal distribution waveforms were prepared, with a median of 250, a standard deviation of 20, and peak heights of 0.0 to 1.0 in increments of 0.1. The data range (trimming range) was set to 0 to 500.
[0126] For each test substance data, four normal distribution waveforms with the median, standard deviation, and peak height set by random numbers were added to obtain the sample data. The peak height (a in Equation 1) was determined by random numbers within the range of 0.0 to 1.5, the median (b in Equation 1) was determined by random numbers within the range of 10 to 490 (±240 relative to the peak median of the test substance), and the standard deviation (c in Equation 1) was determined by random numbers within the range of 10 to 100. 1000 types of sample data were prepared for one test substance data. Each sample data and the peak height of the test substance data contained therein were paired to obtain 11000 pieces of teacher data, and a learning model was generated by machine learning. A fully connected neural network was used as the machine learning method, and the relu function and linear function were used as the activation function. The mean square error was used as the loss function, and Adam was used as the optimization algorithm. Repeated calculations for about 100 epochs were required to obtain sufficient quantitative accuracy.
[0127] Data created using the same method as the sample data was applied to the obtained learning model to determine the peak heights of the test substances contained in the sample data. In Figure 17, the horizontal axis shows the peak heights of the test substances used when creating the sample data (correct values), and the vertical axis shows the peak heights of the test substances obtained using the learning model (predicted values). The correlation coefficient between the correct values and the predicted values was 0.99.
[0128] In addition, a comparative example in which the conventional peak division method (method of establishing a baseline (FIG. 5)) was applied to data created using the same method as the sample data is shown in FIG. 18. The correlation coefficient between the correct value and the predicted value was 0.68.
[0129] Example 6 In order to evaluate the effect of the above-mentioned data processing method, an example in which the method was applied to the quantification of 1,6-dihydro-1-methyl-6-oxo-3-pyridinecarboxamide (hereinafter referred to as 2py), a type of metabolite of vitamin B3 contained in urine, will be described. Measured waveforms of 2py (0.0 μg / mL, 2.5 μg / mL, 5.0 μg / mL, 7.5 μg / mL, 10 μg / mL) were prepared by measuring them by HPLC. The retention time was trimmed in the range from 100 to 400, and four normal distribution waveforms with median, standard deviation, and peak height set by random numbers were added to the waveform of 2py actually measured by HPLC to create sample data. The peak median of 2py was 250, the standard deviation was about 10, and the expected maximum peak height was 500. Based on this, the peak height (a in formula 1) was determined by random numbers within the range of 0 to 750, the median (b in formula 1) was determined by random numbers within the range of 130 to 370, and the standard deviation (c in formula 1) was determined by random numbers within the range of 5 to 50. 1000 types of sample data were prepared for one 2py data. Each sample data and the peak height of the 2py data contained therein were paired to create 5000 pieces of training data, and a learning model was generated by performing machine learning. The machine learning method was the same as in Example 5.
[0130] Next, samples containing 2py (0.0 μg / mL, 2.5 μg / mL, 5.0 μg / mL, 7.5 μg / mL, 10 μg / mL) and other impurities were analyzed by HPLC, and the chromatograms shown in Figure 19 were obtained. Other impurities include (a: uric acid, b: trans-urocanic acid, c: N1-methyl-4-pyridone-3-carboxamide, d: adenosine). When each chromatogram was applied to the generated learning model to predict the height of the peak corresponding to 2py, the correlation coefficient between the correct value and the predicted value was 0.99 (Figure 20). Although 2py was adjacent to the peak of the above impurity (c), it was predicted without any problems.
[0131] Example 7 As an example with more impurities, an example will be described in which the method is applied to the quantification of 2py added to Kirin Co., Ltd.'s "100% Rich Vegetable Juice with 48 Kinds of Additive-Free Vegetables (hereinafter referred to as vegetable juice)". As for 2py, the waveform was measured in advance by HPLC as in Example 6, and a learning model was generated. The retention time was trimmed in the range of 0 to 200. The peak height (a in Equation 1) was determined by random numbers in the range of 0.0 to 1.5, the median (b in Equation 1) was determined by random numbers in the range of 0 to 290, and the standard deviation (c in Equation 1) was determined by random numbers in the range of 10 to 100, and the four normal distribution waveforms were added to create sample data. The above vegetable juice was previously filtered with a 0.8 μm filter to remove coarse particles, and pretreated with Monospin C18 from GL Science Co., Ltd. The pretreatment was performed by centrifugation at 10,000 rpm. 2py (0.0 μg / mL, 2.5 μg / mL, 5.0 μg / mL, 7.5 μg / mL, 10 μg / mL) was added to the pretreated vegetable juice and analyzed by HPLC to obtain the chromatogram in FIG. 21. The arrow in the figure indicates the retention time of 2py. As in Example 6, each chromatogram was applied to the generated learning model to predict the height of the peak corresponding to 2py, and the correlation coefficient between the correct value and the predicted value was 0.99 (FIG. 22). It was found that the peak height of the test substance 2py could be predicted without any problems even in a system with a large amount of impurities and a baseline that does not become 0. [Explanation of symbols]
[0132] 10. Information processing device 21 LAN 22 Database 23 Analyzer 31 Communication Interface 32 ROM 33 RAM 34 Storage section 35 Control section 36 Display section 37 Control Unit 41 Spectral information acquisition unit 42 Learning model generation unit 43 Learning model acquisition unit 44 Estimation part 45 Information Acquisition Department 46 Display control section 700 Calculation Device 701 Test substance 702 Samples 703 Reception 704 Acquisition Department 705 Calculation Department 706,707 Analysis Department 708 Database 709 Generation part 710 Database
Claims
1. A learning model acquisition unit for acquiring a learning model; an information acquisition unit that acquires quantitative information of a test substance estimated by inputting spectral information of a sample containing a plurality of types of test substances and impurities into the learning model, the information acquisition unit has a division unit that divides the spectral information of the sample into spectral information for each of the test substances; The learning model is generated using training data that is a combination of correct answer data and input data, The information processing device is characterized in that the ground truth data is quantitative information of the test substance identified based on spectral information of the test substance, and the input data is spectral information of a virtual sample generated using spectral information of a sample that does not contain the impurities and contains the test substance and random noise obtained by combining a plurality of Gaussian functions.
2. The information processing apparatus according to claim 1 , wherein the dividing unit divides the spectral information at a division point designated by a user.
3. 3 . The information processing apparatus according to claim 1 , wherein the dividing section extracts and divides a position designated by a user from the spectrum information and a certain range including the position.
4. The information processing apparatus according to claim 1 , wherein the dividing unit extracts and divides a range exceeding a preset threshold value from the spectrum information.
5. 5. The information processing apparatus according to claim 1, further comprising a display control means for displaying the acquired quantitative information on a display unit.
6. 6. The information processing apparatus according to claim 5, wherein said display control means selects and displays information designated by a user from the quantitative information on the plurality of types of test substances.
7. 7. The information processing apparatus according to claim 1, wherein each of the plurality of Gaussian functions is obtained by determining a peak height, a median, and a standard deviation using random numbers.
8. 8. The information processing apparatus according to claim 1, further comprising an estimation means for estimating quantitative information of the test substance by inputting spectral information of the sample to the learning model.
9. 9. The information processing device according to claim 1, wherein the spectral information is at least one selected from the group consisting of a chromatogram, a photoelectron spectrum, an infrared absorption spectrum, a nuclear magnetic resonance spectrum, a fluorescence spectrum, an X-ray fluorescence spectrum, an ultraviolet / visible absorption spectrum, a Raman spectrum, an atomic absorption spectrum, a flame emission spectrum, an emission spectroscopy spectrum, an X-ray absorption spectrum, an X-ray diffraction spectrum, a paramagnetic resonance absorption spectrum, an electron spin resonance spectrum, a mass spectrum, and a thermal analysis spectrum.
10. 10. The information processing apparatus according to claim 1, further comprising an analysis means for performing an analysis to obtain spectral information of the sample.
11. 11. The information processing device according to claim 10, wherein the analysis means uses at least one technique selected from the group consisting of chromatography, capillary electrophoresis, photoelectron spectroscopy, infrared absorption spectroscopy, nuclear magnetic resonance spectroscopy, fluorescence spectroscopy, X-ray fluorescence spectroscopy, visible / ultraviolet absorption spectroscopy, Raman spectroscopy, atomic absorption spectrometry, flame emission spectroscopy, optical emission spectroscopy, X-ray absorption spectroscopy, X-ray diffraction, electron spin resonance spectroscopy, mass spectrometry, and thermal analysis.
12. 12. The information processing device according to claim 1, wherein the test substance is at least one selected from the group consisting of proteins, DNA, viruses, fungi, water-soluble vitamins, fat-soluble vitamins, organic acids, fatty acids, amino acids, sugars, pesticides, and environmental hormones.
13. The information processing device according to any one of claims 1 to 12, characterized in that the test substance is at least one selected from the group consisting of thiamine, riboflavin, N1-methylnicotinamide, N1-methyl-2-pyridone-5-carboxamide, N1-methyl-4-pyridone-3-carboxamide, pantothenic acid, pyridoxine, 4-pyridoxic acid, biotin, cyanocobalamin, and ascorbic acid.
14. The information processing device according to any one of claims 1 to 13, characterized in that the quantitative information is at least one type selected from the group consisting of the amount of the test substance contained in the sample, the concentration of the test substance contained in the sample, the presence or absence of the test substance in the sample, and the ratio of the amount of the test substance contained in the sample to a reference amount of the test substance.
15. A learning model acquisition unit for acquiring a learning model; an information acquisition unit that acquires quantitative information of a test substance estimated by inputting spectral information of a sample including the test substance and impurities into the learning model, The learning model is generated using training data that is a combination of correct answer data and input data, the ground truth data is quantitative information of the test substance identified based on the spectral information of the test substance, and the input data is spectral information of a hypothetical sample generated using the spectral information of a sample containing the test substance but not the impurities and random noise obtained by combining a plurality of Gaussian functions; 11. An information processing apparatus comprising: a processor for processing a plurality of Gaussian functions, each of which is obtained by determining a peak height, a median, and a standard deviation using random numbers.
16. The method includes an information acquisition step of acquiring quantitative information of a test substance estimated by inputting spectral information of a sample containing a plurality of types of test substances and impurities into a learning model, the learning model being generated using teacher data constituted by a combination of ground truth data and input data, The information acquiring step includes a dividing step of dividing the spectral information of the sample into spectral information for each of the test substances, An information processing method characterized in that the ground truth data is quantitative information of the test substance identified based on spectral information of the test substance, and the input data is spectral information of a virtual sample generated using spectral information of a sample that does not contain the impurities and contains the test substance and random noise obtained by combining a plurality of Gaussian functions.
17. 17. The information processing method according to claim 16, wherein the dividing step includes a step of dividing the spectral information at a division point designated by a user in the spectral information.
18. 18. The information processing method according to claim 16, wherein the dividing step includes a step of extracting and dividing a position designated by a user from the spectrum information and a certain range including the position.
19. 19. The information processing method according to claim 16, wherein the dividing step includes a step of extracting and dividing a range exceeding a preset threshold value from the spectrum information.
20. 20. The information processing method according to claim 16, further comprising a display control step of displaying the obtained quantitative information on the test substance on a display unit.
21. 21. The information processing method according to claim 20, wherein said display control step includes a step of selecting and displaying information designated by a user from the quantitative information on the test substance.
22. 22. The information processing method according to claim 16, wherein each of the plurality of Gaussian functions is obtained by determining a peak height, a median, and a standard deviation using random numbers.
23. 23. The information processing method according to claim 16, further comprising an estimation step of estimating quantitative information of the test substance by inputting spectral information of the sample into the learning model.
24. 24. The information processing method according to claim 16, wherein the spectral information is at least one selected from the group consisting of a chromatogram, a photoelectron spectrum, an infrared absorption spectrum (IR spectrum), a nuclear magnetic resonance spectrum (NMR spectrum), a fluorescence spectrum, an X-ray fluorescence spectrum, an ultraviolet / visible absorption spectrum (UV / Vis spectrum), a Raman spectrum, an atomic absorption spectrum, a flame emission spectrum, an emission spectroscopy spectrum, an X-ray absorption spectrum, an X-ray diffraction spectrum, a paramagnetic resonance absorption spectrum, an electron spin resonance spectrum, a mass spectrum, and a thermal analysis spectrum.
25. 25. The information processing method according to claim 16, further comprising an analysis step of performing an analysis to obtain spectral information of the sample.
26. 26. The information processing method according to claim 25, wherein the analysis step comprises performing at least one technique selected from the group consisting of chromatography, capillary electrophoresis, photoelectron spectroscopy, infrared absorption spectroscopy, nuclear magnetic resonance spectroscopy, fluorescence spectroscopy, X-ray fluorescence spectroscopy, visible / ultraviolet absorption spectroscopy, Raman spectroscopy, atomic absorption spectrometry, flame emission spectroscopy, optical emission spectroscopy, X-ray absorption spectroscopy, X-ray diffraction, electron spin resonance spectroscopy, mass spectrometry, and thermal analysis.
27. 27. The information processing method according to claim 16, wherein the test substance is at least one of proteins, DNA, viruses, fungi, water-soluble vitamins, fat-soluble vitamins, organic acids, fatty acids, amino acids, sugars, pesticides, and environmental hormones.
28. The information processing method according to any one of claims 16 to 27, characterized in that the test substance is at least one of thiamine, riboflavin, N1-methylnicotinamide, N1-methyl-2-pyridone-5-carboxamide, N1-methyl-4-pyridone-3-carboxamide, pantothenic acid, pyridoxine, 4-pyridoxic acid, biotin, cyanocobalamin, and ascorbic acid.
29. 29. The information processing method according to any one of claims 16 to 28, characterized in that the quantitative information is at least one type selected from the group consisting of the amount of the test substance contained in the sample, the concentration of the test substance contained in the sample, the presence or absence of the test substance in the sample, and the ratio of the amount of the test substance contained in the sample to a reference amount of the test substance.
30. The method includes an information acquisition step of acquiring quantitative information of a test substance estimated by inputting spectral information of a sample containing the test substance and impurities into a learning model, the learning model being generated using teacher data composed of a combination of ground truth data and input data, the ground truth data is quantitative information of the test substance identified based on the spectral information of the test substance, and the input data is spectral information of a hypothetical sample generated using the spectral information of a sample containing the test substance but not the impurities and random noise obtained by combining a plurality of Gaussian functions; An information processing method, wherein each of the plurality of Gaussian functions is obtained by determining a peak height, a median, and a standard deviation using random numbers.
Citation Information
Patent Citations
Method for analyzing chromatogram and chromatographic device
JP1994094696A
Method and apparatus of analyzing and displaying chromatogram
JP1994324029A
A neural network system using spectral analysis to identify complex organic structures.
JP1994505815A
Chromatographic data processing device, chromatographic data processing method, and chromatographic analyzer
JP2006177980A
Ensemble method and apparatus for classifying materials and quantifying the components of mixtures
JP2010520471A