WO2023013164A1 - ＨｂＡ１ｃリスク推定装置、ＨｂＡ１ｃリスク推定方法及びプログラム - Google Patents

ＨｂＡ１ｃリスク推定装置、ＨｂＡ１ｃリスク推定方法及びプログラム Download PDF

Info

Publication number
WO2023013164A1
WO2023013164A1 PCT/JP2022/015095 JP2022015095W WO2023013164A1 WO 2023013164 A1 WO2023013164 A1 WO 2023013164A1 JP 2022015095 W JP2022015095 W JP 2022015095W WO 2023013164 A1 WO2023013164 A1 WO 2023013164A1
Authority
WO
WIPO (PCT)
Prior art keywords
hba1c
risk
risk estimation
biological information
information
Prior art date
Application number
PCT/JP2022/015095
Other languages
English (en)
French (fr)
Inventor
徳隆 安藤
尚久 正箱
Original Assignee
日清食品ホールディングス株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Priority claimed from JP2022040636A external-priority patent/JP2023024256A/ja
Application filed by 日清食品ホールディングス株式会社 filed Critical 日清食品ホールディングス株式会社
Priority to CN202280006142.1A priority Critical patent/CN117693313A/zh
Priority to US18/041,939 priority patent/US20230402185A1/en
Publication of WO2023013164A1 publication Critical patent/WO2023013164A1/ja

Links

Images

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H40/00—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices
    • G16H40/60—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices
    • G16H40/63—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices for local operation
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems

Definitions

  • the present invention relates to an HbA1c risk estimation device, an HbA1c risk estimation method and a program.
  • Hemoglobin is a type of protein in red blood cells and plays a role in transporting oxygen throughout the body. Hemoglobin binds to glucose in the blood and irreversibly changes to glycated hemoglobin. Since hemoglobin circulates in the body for about 120 days, measuring the glycated hemoglobin ratio (HbA1c) is thought to reflect long-term glucose metabolism and blood glucose levels that cannot be determined from fasting blood glucose levels. According to the blood test criteria established by the Society of Ningen Dock, if HbA1c is 5.5% or less, there is no abnormality, and if it is 5.6% or more, it is determined that there is a risk.
  • measuring HbA1c required collecting the subject's blood and conducting a biochemical analysis.
  • Patent Document 1 discloses a method for noninvasively estimating HbA1c by collecting the urine of a subject and using the urinary zinc level as an index.
  • the present invention has been made in view of the above circumstances, and an object of the present invention is to determine the health risk of HbA1c from non-invasive biological information with extremely high accuracy and without blood or urine collection.
  • An object of the present invention is to provide an HbA1c risk estimating device, an HbA1c risk estimating method, and a program that enable rapid estimation.
  • An HbA1c risk estimation device includes an information acquisition unit that acquires predetermined user attribute information and non-invasive biological information, an estimation model storage unit that stores an HbA1c risk estimation model, and an HbA1c risk estimation model. and an estimation processing unit that calculates an HbA1c risk estimate value of a given user based on the given user's attribute information and/or non-invasive biological information.
  • the HbA1c risk estimation device is characterized by further comprising a training data storage unit that stores a training data set, and a learning processing unit that generates an HbA1c risk estimation model by machine learning based on the training data set.
  • Attribute information includes age and sex, or a combination thereof
  • non-invasive biometric information includes BMI, circulating blood volume, blood pressure, pulse wave data, electrocardiogram data, bioimpedance, or a combination thereof. characterized by
  • the training dataset is characterized by including subject's attribute information, non-invasive biological information, and HbA1c measurement value measured from blood.
  • the noninvasive biological information is characterized by further including oxygen saturation (SpO2).
  • the estimation accuracy of the HbA1c risk estimation value is characterized by the accuracy with which the presence or absence of risk can be classified with an ROC_AUC of 0.7 or more.
  • the learning processing unit adds labels indicating the presence or absence of HbA1c risk based on HbA1c measurement values measured from blood, and in the labels, the difference between the number with HbA1c risk and the number without HbA1c risk is greater than or equal to a predetermined value, the sample data in the training data set is increased so as to reduce the difference.
  • the learning processing unit generates a first HbA1c risk estimation model and a second HbA1c risk estimation model by machine learning based on different types of training data sets, and the estimation processing unit generates the first HbA1c risk estimation model and the second HbA1c risk estimation model. is used to calculate the HbA1c risk estimate for a given user.
  • a biological information estimator for estimating at least one biological information selected from BMI, blood pressure, pulse wave data, electrocardiogram data, biological impedance, and oxygen saturation included in the biological information;
  • the biometric information estimated by the estimation unit is obtained as the biometric information of the predetermined user.
  • HbA1c risk estimation device In addition to the HbA1c risk estimation device, it is characterized by further comprising a biological information measurement device that measures non-invasive biological information to constitute a noninvasive HbA1c risk estimation system.
  • the HbA1c risk estimation method of the present invention includes steps of storing a training data set including subject attribute information, non-invasive biological information, and HbA1c measured values measured from blood, and performing machine learning based on the training data set. generating an HbA1c risk estimation model by and using the HbA1c risk estimation model to calculate an HbA1c risk estimate for a given user based on the given user's attribute information and/or non-invasive biometric information. characterized by comprising
  • the program in the present invention comprises steps of storing a training data set including subject attribute information, non-invasive biological information, and HbA1c measured values measured from blood; generating an estimation model; and using the HbA1c risk estimation model to calculate an HbA1c risk estimate for a given user based on a given user's attribute information and/or non-invasive biometric information. It is characterized by executing
  • an HbA1c risk estimation device capable of extremely accurately estimating the health risk of HbA1c by machine learning using non-invasive biological information. be able to.
  • FIG. 1 is a block diagram showing a schematic configuration of an HbA1c risk estimation system; FIG. It is a figure for demonstrating ESG (Electroscangram).
  • FIG. 2 is a hardware configuration diagram of an HbA1c risk estimation device; 4 is a flow chart showing a procedure for generating and executing an HbA1c risk estimation model by machine learning; It is the hierarchical structure of the neural network (NN) used for HbA1c risk estimation.
  • FIG. 10 is a flow chart showing an execution procedure of HbA1c risk estimation processing;
  • FIG. 4 is an ROC_AUC curve of estimation results in Example 1.
  • FIG. 10 is an ROC_AUC curve of estimation results in Example 2.
  • FIG. 1 is a block diagram showing a schematic configuration of an HbA1c risk estimation system 1 according to this embodiment.
  • the HbA1c risk estimation system 1 includes a terminal device 10 , a biological information measurement device 20 , an HbA1c risk estimation device 30 and a display device 39 .
  • "user” refers to a person who uses the HbA1c risk estimation system in order to non-invasively obtain an HbA1c risk estimate.
  • subjects are used in the HbA1c risk estimation system with attribute information such as age and gender, non-invasive biological information, and HbA1c measured values measured from blood after obtaining prescribed procedures and consent.
  • attribute information such as age and gender, non-invasive biological information, and HbA1c measured values measured from blood after obtaining prescribed procedures and consent.
  • the terminal device 10 is an information terminal capable of inputting user attribute information (name, ID, age, gender, etc.) and outputting the input information to the HbA1c risk estimation device 30 via a wired or wireless communication network.
  • Any information terminal may be used as long as it is. Examples include tablet terminals, smart phones, mobile terminals including wearable terminals, and PCs (Personal Computers).
  • the height, weight, and the like may be measured by the biological information measuring device 20, which will be described later.
  • the biological information measuring device 20 measures the user's non-invasive biological information.
  • non-invasive biometric information is biometric information obtained by a method that does not require the insertion of instruments into the skin or into body orifices.
  • Non-invasive biological information can be measured using commercially available stature scales, weight scales, sphygmomanometers, pulse oximeters, pulse wave meters, electrocardiographs, impedance measuring devices, galvanic skin measuring devices, and the like.
  • ESTEC BC-3 which can simultaneously measure pulse wave data, electrocardiogram data, bioimpedance, and oxygen saturation (SpO2), can also be used. These devices can noninvasively measure biometric data without imposing a psychological or physical burden on the user.
  • the noninvasive biological information includes any one or a combination of BMI (body mass index), circulating blood volume, blood pressure, pulse wave data, electrocardiogram data, and bioimpedance, and oxygen It may also include saturation (SpO2).
  • BMI body mass index
  • circulating blood volume blood pressure
  • pulse wave data pulse wave data
  • electrocardiogram data electrocardiogram data
  • bioimpedance oxygen It may also include saturation (SpO2).
  • Circulating blood volume [L] 0.3669 x h 3 + 0.03129 x w + 0.6041
  • Circulating blood volume [L] 0.3561 x h 3 + 0.03308 x w + 0.1833
  • Pulse wave data is obtained by irradiating red light ( ⁇ 660 nm) from a Red LED and near-infrared light ( ⁇ 905 nm) from an IR LED to a protruding part of the body, such as a finger, using a pulse wave meter or pulse oximeter. is measured by measuring the phototransistor.
  • Pulse wave data include pulse, elastic index, peripheral vascular resistance, accelerogram, b/a, e/a, -d/a, Takazawa accelerometer aging index, ejection fraction, LVET, DEI elastic index), or a combination thereof.
  • the elastic index is a numerical value obtained by dividing the height by the time from the systolic peak to the diastolic peak detected in the finger plethysmogram.
  • Peripheral vascular resistance is determined by mean arterial pressure/cardiac output ⁇ 80.
  • DEI depleted elastic index
  • 0.3 to 0.7 is normal, 0.3 or less suggests hypertension or arteriosclerosis, and 0.7 or more suggests acute anxiety.
  • the acceleration pulse wave is a second derivative of photoplethysmogram (SDPTG) of finger plethysmogram (photoplethysmogram: PTG).
  • Acceleration pulse wave includes early positive wave (a-wave), early negative wave (b-wave), mid-systolic re-rising wave (c-wave), late-systolic re-descending wave (d-wave), and diastolic early positive wave (e-wave).
  • the above b/a, e/a, and -d/a are calculated from the ratio of each wave height.
  • vascular aging is evaluated using the Takazawa acceleration pulse wave aging index (bcde)/a. can do.
  • Ejection fraction is the fraction of blood ejected from the ventricle with each heartbeat and is proportional to the Accelerometer Pulsatility Aging Index.
  • LVET is Left Ventricular Ejection Time, the time for blood in the left ventricle to be ejected into the aorta after the aortic valve opens.
  • Electrocardiogram data can be measured by electrode electrocardiography (ECG) or photoplethysmography (PPG).
  • ECG data includes respiratory rate, heart rate, RR interval, standard deviation of RR interval, MxDMn ratio, power spectrum of low frequency band, power spectrum of high frequency band, heart rate variability index LF/HF, total power, or combinations thereof.
  • the RR interval is the interval from one QRS wave to the next QRS wave in an electrocardiogram.
  • the MxDMn ratio is the ratio of the longest RR interval to the shortest RR interval in time and is an index of irregular heartbeats.
  • the total power is the calculated value of the total power of the power spectrum of frequencies 0 to 0.4 Hz (VLF, LF, HF) measured for 2 minutes.
  • This value reflects the total autonomic nervous system activity dominated by sympathetic activity.
  • the high frequency power spectral ratio (0.1875 to 0.50 Hz: HF)
  • the low frequency power spectral ratio (0.05 to 0.1875 Hz: LF)
  • the LF/HF ratio the power spectrum ratio of the very low frequency band (0-0.05 Hz: VLF)
  • the impedance (conductance) of a living body can be measured, for example, by passing weak electricity between two of the six electrodes on both legs, both hands, and left and right forehead.
  • ⁇ S anode/cathode conductance
  • ⁇ S cathode/anode conductance
  • ⁇ S/m electrical conductivity
  • muscle mass, body fat mass, total water content, phase angle, and resistance value can be measured at the same time.
  • the dielectric constant ( ⁇ Si) when the current is applied between the right hand and the left hand and between the right forehead and the left forehead is also possible. It is preferable to measure the impedance (conductance) of the living body using 22 patterns of conductivity from six electrodes.
  • Bioimpedance includes body fat mass (kg), body fat mass (%), lean mass, lean mass, muscle mass, total water content (kg), total water content (%), intracellular water content (% ), cardiac output, 1 forehead left-2 left hand/SCR A, 1 forehead left-2 left hand/delta SCR C-SCR A, 3 forehead right-4 right hand/SCRC, 5 left hand-6 left foot/delta SCR C- SCR A, 9 left forehead - 10 right forehead / SCRC, 13 left foot - 14 right foot / SCR C, 13 left foot - 14 right foot / SCR A, 15 right hand - 16 forehead left / Delta SCR C - SCR A, 15 right hand - 16 forehead Left SCR C, 19 Right foot -20 Left hand/Delta SCR C-SCR A , ESG2+4+15+17 ( ⁇ S/m), ESG6+13+19 (%), ESG 6+8+19+21 (%), ESG6+8+19+21 ( ⁇ S/m), ESG
  • ESG is an abbreviation for electroscangram.
  • the "+" in ESG2+4+15+17 means which electrode attached to the body was used to measure.
  • ESG 2 + 4 + 15 + 17, as shown in Fig. 2 is the left hand when electricity is applied from the left hand to the left forehead, the right hand when electricity is applied from the right hand to the right forehead, the right hand when electricity is applied from the right hand to the left forehead, and the left hand to the right forehead. It means the average value of the conductance measured with the left hand when the current is applied.
  • Non-Patent Document 1 has a detailed description of these conductances.
  • Non-Patent Documents 2 and 3 The details of the ESG (electroscangram) measurement method are described in Non-Patent Documents 2 and 3.
  • “1 left side -2 left hand/SCR A” is the conductance (or conductivity) of the path measured when electricity flows so that "1 left side” is the cathode and "2 left hand” is the anode direction.
  • “5 left hand - 6 left foot / delta SCR C-SCR A” is the difference in conductance measured when current is applied between "5 left hand” and "6 left foot” so that anode-cathode, cathode-anode. be.
  • BMI and blood pressure can be measured using a height weight scale and a blood pressure monitor.
  • the non-invasive biological information may also include oxygen delivery calculated from SpO2 and cardiac output.
  • the measured noninvasive biological information is output to the HbA1c risk estimation device 30 through a wired or wireless communication network.
  • the biological information measuring device 20 may be a built-in measuring device or a portable measuring device such as a wearable terminal.
  • the HbA1c risk estimation device 30 includes a first acquisition unit 31, a second acquisition unit 32, a user data storage unit 33, a training data storage unit 34, a learning processing unit 35, an estimation model storage unit 36, and an estimation process.
  • a section 37 and an estimated data storage section 38 are provided.
  • the first acquisition unit 31 acquires user attribute information from the terminal device 10 .
  • the second acquisition unit 32 acquires non-invasive biological information of the user from the biological information measuring device 20 .
  • the user data storage unit 33 stores user attribute information and non-invasive biological information acquired from the first acquisition unit 31 and the second acquisition unit 32 .
  • the training data storage unit 34 consists of attribute information of a plurality of subjects acquired in advance, non-invasive biological information, and sample test information such as HbA1c obtained by blood test as a training data set for machine learning. A plurality of training datasets are stored.
  • the specimen test information may further include test information obtained from blood, urine, stool, and the like.
  • the learning processing unit 35 acquires the training data set stored in the training data storage unit 34 and creates an HbA1c risk estimation model using the training data set. Specifically, when estimating HbA1c risk, the acquired training data set is subjected to gradient boosting such as XGBoost, neural network, and logistic regression, or machine learning by ensemble learning of these learning results, attribute information and non-invasive learn the relationship between biological information and HbA1c risk.
  • the estimated model storage unit 36 stores the HbA1c risk estimated model generated by the learning processing unit 35 .
  • Non-invasive biological information includes BMI, circulating blood volume, blood pressure, pulse wave data, electrocardiogram data, and biological impedance, or a combination thereof. Additionally, oxygen saturation (SpO2) is included as needed.
  • the estimation processing unit 37 uses the estimation model generated by the learning processing unit 35 to estimate the user's HbA1c risk based on predetermined user attribute information and/or non-invasive biological information. Then, the estimated HbA1c risk value is stored in the estimated data storage unit 38 .
  • the display device 39 can display the HbA1c risk estimate along with the user's attribute information and non-invasive biological information. These data may be displayed on the terminal device 10 owned by the user.
  • FIG. 3 is a hardware configuration diagram of the HbA1c risk estimation device 30.
  • the HbA1c risk estimation device 30 is composed of a computer 300 having one or more processors 301 , memory 302 , storage 303 , input/output port 304 and communication port 305 .
  • the processor 301 performs processing related to HbA1c estimation according to this embodiment by executing a program.
  • the memory 302 temporarily stores programs and computation results of the programs.
  • Storage 303 stores a program for executing processing by HbA1c risk estimation device 30 .
  • the storage 303 may be of any kind as long as it is readable by a computer, for example, various recording media such as a recording medium (magnetic disk, optical disk, etc.), random access memory, flash type memory, read-only memory, etc. can be used.
  • the input/output port 304 inputs information from the terminal device 10 and the biological information measurement device 20 and outputs the HbA1c estimated value to the display device 39 .
  • a communication port 305 transmits and receives data to and from an information terminal such as another computer (not shown). As a communication method, wireless communication or wired communication can be used.
  • the HbA1c risk estimator 30 can be implemented by a commercially available desktop PC or notebook PC, and it takes several seconds to calculate the HbA1c risk estimate using the estimation model.
  • the processor 301 of the HbA1c risk estimation device 30 the above-described first acquisition unit 31, second acquisition unit 32, learning processing unit 35, estimation processing unit 37, etc. function during operation.
  • HbA1c risk estimation without blood test by generating HbA1c risk estimation model by machine learning based on non-invasive biological data including BMI (body mass index), blood pressure, pulse wave data, electrocardiogram data, bioimpedance, etc. It is possible to do Further, as described below, it is possible to determine whether or not HbA1c is normal even when the number of each data included in the noninvasive biological data is limited.
  • BMI body mass index
  • blood pressure blood pressure
  • pulse wave data electrocardiogram data
  • bioimpedance bioimpedance
  • FIG. 4 is a flow chart showing a procedure for generating an HbA1c risk estimation model by machine learning.
  • the learning processing section 35 preprocesses the input data. Specifically, the learning processing unit 35 converts the HbA1c obtained by the blood test to 0 (no risk) when the HbA1c is 5.5% or less, and 1 (risk) when the HbA1c is 5.6% or more.
  • the learning processing unit 35 applies SMOTE (Chawla, NV. et al. 2002) to the learning data, Training samples may be generated artificially.
  • a label indicating the presence or absence of HbA1c risk (e.g., 0 or 1 above) is added based on HbA1c obtained by blood test, and in the label, there is HbA1c risk (label : If the difference between the number of 1) and the number of no HbA1c risk (label: 0) is a predetermined value or more, even if the sample data in the training data set is increased (generated) so as to reduce this difference good.
  • the learning processing unit 35 performs machine learning using logistic regression, neural network, and gradient boosting decision tree
  • step ST103 learning result of machine learning using logistic regression, neural network, and gradient boosting decision tree Ensemble learning (voting) is further performed for
  • FIG. 5 shows the structure of a neural network (NN) used as an HbA1c risk estimation model. Rectangles indicate layers for data conversion, and rounded rectangles indicate input/output data. D is the number of inspection items.
  • the NN structure transforms data in the order of D dimension, 96 dimension, 96 dimension, 96 dimension, and 1 dimension via four layer groups (Compile A1, Compile A2, Compile A3, Compile B1).
  • Compile A1, A2, and A3 include a fully connected layer “Linear” that performs fully connected processing, a "Kernel Regularizer” that performs regularization, and a "ReLU” layer that performs ReLU processing.
  • the hidden layer includes a dropout layer that suppresses some input values to 0 to prevent overfitting.
  • Logistic Regression provided by Scikit-learn, an open source machine learning library for Python, can be used. Also, if necessary, the number of dimensions may be compressed by principal component analysis.
  • the HbA1c risk obtained from the blood test was compared with the HbA1c risk estimate estimated by machine learning, and each parameter of Logistic Regression (C, regularization method, max_iter, solver) was adjusted to maximize the f1 score.
  • C is a trade-off parameter that determines the strength of regularization, and the larger the value, the weaker the strength of regularization.
  • Regularization method means L1 regularization or L2 regularization, which is selected.
  • max_iter is the maximum number of iterations to learn.
  • the solver chooses a convergence method (eg L-BFGS method, Newton CG method, liblinear, sag and saga) that minimizes the cross-entropy error.
  • a convergence method eg L-BFGS method, Newton CG method, liblinear, sag and saga.
  • Machine learning using gradient boosted regression trees can use software libraries such as XGBoost, CatBoost, and LightBGM.
  • the risk value (0: no risk, 1: with risk) obtained by comparing the HbA1c value obtained by blood test with a predetermined threshold is compared with the HbA1c risk value estimated by machine learning, and the f1 score is
  • Each parameter of XGBoost max_depth, subsample, colsample_bytree, learning_rate) was adjusted to maximize.
  • max_depth is the depth of the decision tree
  • subsample is the proportion of randomly sampled samples in each tree
  • colsample_bytree is the proportion of randomly sampled columns in each tree
  • learning_rate indicates the learning rate.
  • the learning processing unit 35 stores the HbA1c risk estimation model generated by the learning process described above in the estimation model storage unit 36 . It should be noted that the machine learning algorithms described above are merely examples, and the present invention is not limited to these.
  • step ST201 the first acquisition unit 31 of the HbA1c risk estimation device 30 acquires user attribute information from the terminal device 10.
  • the second acquisition unit 32 of the HbA1c risk estimation device 30 acquires non-invasive biological information of the user. Then, the user's attribute information and non-invasive biological information are stored in the user data storage unit 33 .
  • step ST203 using the HbA1c risk estimation model stored in the estimation model storage unit 36, the estimation processing unit 37 determines the probability of belonging to class 0 (no risk) or class 1 (with risk), that is, HbA1c risk estimation. Calculate the value.
  • step ST204 the calculated HbA1c risk estimate is stored in the estimate data storage unit 38, and in step ST205, the HbA1c risk estimate is output to an external terminal such as the display device 39 for display.
  • HbA1c risk estimation An example of HbA1c risk estimation is described below.
  • Attribute information includes ID, name, age, sex, or a combination of these.
  • Non-invasive biometric information includes BMI, circulating blood volume, blood pressure, pulse wave data, electrocardiogram data, bioimpedance, and oxygen saturation (SpO2), or combinations thereof.
  • Height and weight which are the criteria for calculating BMI, were measured using a height scale and weight scale, respectively, and blood pressure was measured using a sphygmomanometer.
  • Pulse wave data, electrocardiogram data, bioimpedance, and oxygen saturation (SpO2) were measured with Estech BC-3 (Ryobi Systems). Instead of Estech BC-3, a commercially available plethysmograph, electrocardiograph, impedance measuring device, and pulse oximeter may be used in combination. Alternatively, the non-invasive biometric information may be acquired using a predetermined wearable terminal. Bioimpedance (conductance) was measured by passing a weak electric current between two of six electrodes on both feet, both hands, and left and right forehead. The voltage and current were 1.28 V, 200 ⁇ A and the conductance was measured for 32 milliseconds per second.
  • Estec BC-3 measured pulse wave data, electrocardiogram data, bioimpedance, and oxygen saturation (SpO2) for 2 minutes for each subject.
  • a device equipped with electrocardiogram, pulse wave monitor, and pulse oximeter functions was attached to the subject's left index finger, and two electrodes were attached to the forehead. placed on top.
  • ⁇ Learning model 1> In learning model 1, machine learning was performed using only the logistic regression in FIG. Since only logistic regression is used, ensemble learning in FIG. 4 is not performed. At this time, the data shown below were selected and used as non-invasive biological data. Note that this model does not use attribute information.
  • BMI blood pressure: diastolic blood pressure, pulse pressure, pulse wave data: elastic index, pulse, electrocardiogram data: heart rate, LF/HF, bioimpedance: 5 left hand - 6 left leg / delta SCRC-SCRA, left leg conductance, R
  • the noninvasive biometric data further includes pulse pressure/pulse data, cardiac output included in bioimpedance, and oxygen delivery estimated from oxygen saturation (SpO2).
  • ⁇ Learning model 2> In learning model 2, as shown in Fig. 4, a blood HbA1c risk estimation model is created by performing machine learning using logistic regression, neural network, and gradient boosting, and performing ensemble learning by voting on these learning results. bottom. At this time, the data shown below were selected and used as non-invasive biological data.
  • Non-invasive biological data Body weight, BMI, circulating blood volume
  • Blood pressure systolic blood pressure, diastolic blood pressure, mean arterial pressure
  • Pulse wave data elastic index ⁇ Electrocardiogram data ⁇ Heart rate, LF/HF ⁇ Bioimpedance ⁇ Muscle mass (%), total water content (kg), total water content (%), cardiac output, 3 right side of forehead - 4 right hand/SCRC , left forehead - right forehead /SCRA, impedance index for right forehead (sum of right forehead -4 right forehead /SCRC and left forehead for forehead - right forehead for 10/SCRC), conductance for left foot, R ( ⁇ ), once Cardiac output
  • Non-invasive biometric data further includes cardiac output included in bioimpedance and oxygen delivery estimated from oxygen saturation (SpO2).
  • Example 1 In Example 1, (1) attribute information of 321 subjects, (2) height and weight scale, sphygmomanometer, and noninvasive biological information measured by Estec BC-3, and (3) noninvasive biological information measurement and Using the HbA1c training data set obtained from the blood test performed on the same day, the HbA1c risk estimation model was generated through the machine learning of learning model 1 above. Then, the estimated accuracy of the HbA1c risk estimation model was evaluated by the ROC_AUC curve. As a result, ROC_AUC showed 0.78, which exceeded 0.7 indicating that good classification was achieved.
  • FIG. 7 shows the ROC_AUC curve of the estimation result of Example 1. In FIG.
  • Example 2 In Example 2, (1) attribute information of a total of 712 subjects, (2) height and weight scale, blood pressure monitor, and noninvasive biological information measured by Estech BC-3, and (3) noninvasive biological information An HbA1c risk estimation model was generated through machine learning of learning model 2 above using a training data set of HbA1c obtained from a blood test performed on the same day as the measurement. Then, the estimated accuracy of the HbA1c risk estimation model was evaluated by the ROC_AUC curve. As a result, ROC_AUC showed 0.80, which is the same value as 0.8, which indicates extremely good classification.
  • FIG. 8 shows the ROC_AUC curve of the estimation result of Example 2. As shown in FIG.
  • HbA1c risk may be estimated using a plurality of learning models.
  • HbA1c risk can be estimated with higher accuracy than using one learning model.
  • BMI is generally not obtained by wearable devices, etc., but is obtained from the height and weight entered by the user.
  • the above-mentioned HbA1c risk is obtained by obtaining only biological information. Therefore, convenience for the user is improved.
  • the method of estimating BMI is not particularly limited, for example, BMI is known to be correlated with the inclination of the user's abdomen (at a predetermined position).
  • a predetermined acceleration sensor is provided on the abdomen of the user (or a wristband-type wearable terminal equipped with an acceleration sensor is applied to the abdomen), and the inclination of the abdomen is calculated based on the data output from the acceleration sensor. may be obtained and estimated BMI.
  • pulse wave data and oxygen saturation may also be estimated using a wristband-type wearable terminal equipped with a pulse wave sensor or medium oxygen concentration sensor, as described above.
  • blood pressure may also be estimated using a wristband-type wearable terminal or the like in the same manner as described above. Since it is known that there is a correlation between the speed of the pulse wave sent through the artery by the heartbeat and the blood pressure, using a predetermined sensor that measures the speed of the pulse wave sent through the artery by the heartbeat, Blood pressure may be estimated.
  • electrocardiogram data may also be estimated using a wristband-type wearable terminal or the like in the same manner as described above.
  • the electrocardiogram data can be estimated based on the data obtained from the electrodes provided on the surface opposite to the display surface of the wristband wearable terminal and the electrodes provided on the display surface side.
  • the wrist of a hand for example, left hand wearing a wristband-type wearable terminal is in contact with the electrode provided on the opposite side, and the hand attached to the electrode provided on the display surface side
  • Electrocardiogram data may be estimated from data obtained by contacting the fingertips of the opposite hand (eg, right hand).
  • bioimpedance may also be estimated using a wristband-type wearable terminal or the like equipped with various electrodes.
  • bioimpedance may be estimated based on biometric information obtained from the chest and wrist using a wristband-type wearable terminal or the like.
  • biometric information that can be acquired by the wearable terminal and the biometric information to be estimated (BMI, blood pressure, pulse wave data, electrocardiogram data, and biometric information of at least one or more of bioimpedance and oxygen saturation)
  • BMI blood pressure
  • pulse wave data pulse wave data
  • electrocardiogram data electrocardiogram data
  • biometric information of at least one or more of bioimpedance and oxygen saturation As teacher data, classifiers generated using various machine learning algorithms may be used to estimate the biometric information of the estimation target described above.
  • the above-described second acquisition unit may acquire the estimated biological information.
  • the series of processes described above can be executed by hardware or by software.
  • the functional configuration described above is merely an example and is not particularly limited. In other words, it is sufficient for the information processing system to have a function capable of executing the series of processes described above as a whole, and what kind of functional block is used to realize this function is not particularly limited to the example described above.
  • the locations of the functional blocks are not particularly limited to those shown in FIG. 1, and may be arbitrary.
  • the functional blocks of the server may be transferred to another terminal, device, or the like.
  • functional blocks of other terminals or devices may be transferred to a server or the like.
  • one functional block may be composed of hardware alone, software alone, or a combination thereof.
  • the computer may be a computer built into dedicated hardware. Also, the computer may be a computer capable of executing various functions by installing various programs, such as a server, a general-purpose smart phone, or a personal computer.
  • a recording medium containing such a program not only consists of a removable medium (not shown) that is distributed separately from the device main body in order to provide the program to the user, etc., but is also preinstalled in the device main body and stored in the user's memory. It is composed of a recording medium etc. provided for Since the program can be distributed via a network, the recording medium may be installed in or accessible to a computer connected or connectable to the network.
  • the steps of writing a program recorded on a recording medium are not necessarily processed chronologically, but may be executed in parallel or individually. It also includes the processing to be performed.
  • the term "system” means an overall device composed of a plurality of devices, a plurality of means, or the like.
  • HbA1c risk estimation system 10 terminal device 20 biological information measurement device 30 HbA1c risk estimation device 31 first acquisition unit 32 second acquisition unit 33 user data storage unit 34 training data storage unit 35 learning processing unit 36 estimation model storage unit 37 estimation processing Unit 38 Estimated data storage unit 39

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Public Health (AREA)
  • Medical Informatics (AREA)
  • Biomedical Technology (AREA)
  • General Health & Medical Sciences (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Pathology (AREA)
  • Databases & Information Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Surgery (AREA)
  • Optics & Photonics (AREA)
  • Veterinary Medicine (AREA)
  • Animal Behavior & Ethology (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Business, Economics & Management (AREA)
  • General Business, Economics & Management (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

【課題】従来、HbA1cの測定には、被験者の血液を採取し、生化学的に分析を行う必要があった。しかしながら、この方法では、侵襲的に被験者の皮膚に針等を刺す必要があり、被験者にとって心理的または肉体的な負担を伴うという課題があった。 【解決手段】本発明によれば、事前に取得した複数の被験者の属性情報、非侵襲な生体情報、および血液検査の検査データに基づき機械学習によりHbA1cリスク推定モデルを生成することで、所定のユーザの属性情報および／または非侵襲な生体情報から非侵襲にHbA1cリスクを推定することが可能となる。

Description

ＨｂＡ１ｃリスク推定装置、ＨｂＡ１ｃリスク推定方法及びプログラム
　本発明は、HbA1cリスク推定装置、HbA1cリスク推定方法及びプログラムに関する。
　ヘモグロビンは赤血球内のたんぱく質の１種であり、全身に酸素を運搬する役割を担っている。ヘモグロビンは血中のブドウ糖と結合し、不可逆的に糖化ヘモグロビンへと変化する。ヘモグロビンは120日間ほど体内を循環することから、糖化ヘモグロビンの割合（HbA1c）を測定することによって、空腹時血糖値ではわからない長期的な糖代謝、血糖値を反映すると考えられている。人間ドック学会が制定する血液検査での判定基準によると、HbA1cが５．５％以下であれば異常なし、５．６％以上であればリスク有と判定される。
　従来、HbA1cの測定には、被験者の血液を採取し、生化学的に分析を行う必要があった。しかしながら、この方法では、侵襲的に被験者の皮膚に針等を刺す必要があり、被験者にとって心理的または肉体的な負担を伴うという課題があった。
　一方、特許文献１には、被験者の尿を採取し、尿中亜鉛値を指標として非侵襲にHbA1cを推定する方法が開示されている。
特開2021-047056
Psychology Research and Behavior Management 2011:4 81-86、Summary of the clinical investigations E.S.Teck Complex March, 20, 2010要旨 R. N. Chua, Y. W. Hau, C. M. Tiew and W. L. Hau, "Investigation of Attention Deficit/Hyperactivity Disorder Assessment Using Electro　Interstitial Scan Based on Chronoamperometry Technique," in IEEE　Access, vol. 7, pp. 144679-144690, 2019, doi:　10.1109/ACCESS.2019.2938095. Maarek A. Electro interstitial scan system: assessment of 10 years of research and development. Med Devices (Auckl). 2012;5:23-30. doi:10.2147/MDER.S29319
　しかしながら、特許文献１に記載された方法では、被験者の尿を採取する必要があり、亜鉛やクレアチニン量等の測定には分析装置が必要であるなどの課題があった。
　本発明は上記のような事情に鑑みてなされたものであって、本発明の目的は、採血や採尿を行わずに、非侵襲な生体情報からHbA1cの健康へのリスクを極めて精度良く、かつ迅速に推定することが可能なHbA1cリスク推定装置、HbA1cリスク推定方法及びプログラムを提供することにある。
　本発明に係るHbA1cリスク推定装置は、所定のユーザの属性情報および非侵襲な生体情報を取得する情報取得部と、HbA1cリスク推定モデルを記憶する推定モデル記憶部と、HbA1cリスク推定モデルを用いて、所定のユーザの属性情報および／または非侵襲な生体情報に基づいて、所定のユーザのHbA1cリスク推定値を計算する推定処理部と、を備えることを特徴とする。
　本発明に係るHbA1cリスク推定装置は、訓練データセットを記憶する訓練データ記憶部と、訓練データセットに基づいて、機械学習によりHbA1cリスク推定モデルを生成する学習処理部と、をさらに備えることを特徴とする。
　属性情報は、年齢、および性別のいずれか、またはこれらの組み合わせを含み、非侵襲な生体情報はBMI、循環血液量、血圧、脈波データ、心電図データ、生体インピーダンス、またはこれらの組み合わせを含むことを特徴とする。
　訓練データセットは、被験者の属性情報、非侵襲な生体情報、および血液から測定されたHbA1c測定値を含むことを特徴とする。
　非侵襲な生体情報は、酸素飽和度（ＳｐＯ２）をさらに含むことを特徴とする。
　HbA1cリスク推定値の推定精度は、リスクの有無をROC_AUCが０．７以上で分類可能な精度であることを特徴とする。
　学習処理部は、訓練データセットにおいて、血液から測定されたHbA1c測定値に基づいてHbA1cリスクの有無を示すラベルを付加し、ラベルにおいて、HbA1cリスク有りの数と、HbA1cリスク無しの数との差が所定の値以上の場合は、差を小さくするように訓練データセットにおけるサンプルデータを増加させることを特徴とする。
　学習処理部は、異なる種類の訓練データセットに基づいて、それぞれ機械学習により第１HbA1cリスク推定モデルおよび第２HbA1cリスク推定モデルを生成し、推定処理部は、第１HbA1cリスク推定モデルおよび第２HbA1cリスク推定モデルを用いて、所定のユーザのHbA1cリスク推定値を計算することを特徴とする。
　生体情報に含まれるBMI、血圧、脈波データ、心電図データ、および生体インピーダンス、酸素飽和度のうち少なくとも１つ以上の生体情報を推定する生体情報推定部をさらに備え、情報取得部は、生体情報推定部によって推定された生体情報を所定のユーザの生体情報として取得することを特徴とする。
　HbA1cリスク推定装置に加えて、非侵襲な生体情報を測定する生体情報測定装置をさらに備え、非侵襲HbA1cリスク推定システムを構成することを特徴とする。
　本発明におけるHbA1cリスク推定方法は、被験者の属性情報、非侵襲な生体情報、および血液から測定されたHbA1cの測定値を含む訓練データセットを記憶するステップと、訓練データセットに基づいて、機械学習によりHbA1cリスク推定モデルを生成するステップと、HbA1cリスク推定モデルを用いて、所定のユーザの属性情報および／または非侵襲な生体情報に基づいて、所定のユーザのHbA1cリスク推定値を計算するステップを含むことを特徴とする。
　本発明におけるプログラムは、被験者の属性情報、非侵襲な生体情報、および血液から測定されたHbA1cの測定値を含む訓練データセットを記憶するステップと、訓練データセットに基づいて、機械学習によりHbA1cリスク推定モデルを生成するステップと、HbA1cリスク推定モデルを用いて、所定のユーザの属性情報および／または非侵襲な生体情報に基づいて、所定のユーザのHbA1cリスク推定値を計算するステップと、をコンピュータに実行させることを特徴とする。
　本発明によれば、非侵襲な生体情報を用いて機械学習により極めて精度良くHbA1cの健康へのリスクを極めて精度良く推定することが可能なHbA1cリスク推定装置、HbA1cリスク推定方法及びプログラムを提供することができる。
HbA1cリスク推定システムの概略構成を示すブロック図である。 ESG（エレクトロスキャングラム）を説明するための図である。 HbA1cリスク推定装置のハードウェア構成図である。 機械学習によるHbA1cリスク推定モデルの生成実行手順を示すフローチャートである。 HbA1cリスク推定時に使用するニューラルネットワーク（NN）の階層構造である。 HbA1cリスク推定処理の実行手順を示すフローチャートである。 実施例１における推定結果のROC_AUC曲線である。 実施例２における推定結果のROC_AUC曲線である。
　以下、実施形態について図面を参照しつつ説明する。なお、実施形態は例示であり、本発明は下記の構成に限定されるものではない。
＜装置の機能＞
　図１～図６を参照して、本実施形態に係るHbA1cリスク推定システム１およびHbA1cリスク推定装置３０を説明する。図１は、本実施形態に係るHbA1cリスク推定システム１の概略構成を示すブロック図である。HbA1cリスク推定システム１は、端末装置１０と、生体情報測定装置２０と、HbA1cリスク推定装置３０と、表示装置３９とを備える。
　ここで、「ユーザ」とは、非侵襲にHbA1cリスクの推定値を得るために、HbA1cリスク推定システムを利用する人をいう。また、「被験者」とは、所定の手続きと同意を得た上で、年齢や性別等の属性情報、非侵襲な生体情報、および血液から測定されたHbA1c測定値を、HbA1cリスク推定システムにおいて利用される訓練データセットとして提供する人をいう。
　端末装置１０は、ユーザの属性情報（氏名、ＩＤ、年齢、性別等）が入力可能であり、かつ、入力された情報を有線または無線の通信ネットワークを通じてHbA1cリスク推定装置３０へ出力可能な情報端末であれば、いかなる情報端末であっても構わない。例えば、タブレット端末、スマートフォン、ウェアラブル端末等を含む携帯端末等、またはＰＣ（Ｐｅｒｓｏｎａｌ　Ｃｏｍｐｕｔｅｒ）が挙げられる。なお、身長や体重等は、後述する生体情報測定装置２０により測定してもよい。
　生体情報測定装置２０は、ユーザの非侵襲な生体情報を測定する。ここで、非侵襲な生体情報とは、皮膚内への、または体の開口部への器具の挿入を必要としない方法で取得した生体情報である。非侵襲な生体情報は、市販の身長計、体重計、血圧計、パルスオキシメーター、脈波計、心電図計、インピーダンス測定機、ガルバニック皮膚測定機等を利用して測定することが可能である。また、脈波データ、心電図データ、生体インピーダンス、酸素飽和度（ＳｐＯ２）を同時に測定可能な、エステックBC-3（両備システムズ）を利用することもできる。これらの装置は、ユーザに心理的、または肉体的な負担を与えることなく、非侵襲な生体データを測定することが可能である。
　本発明の実施形態において、非侵襲な生体情報は、ＢＭＩ（ボディマス指数）、循環血液量、血圧、脈波データ、心電図データ、および生体インピーダンスのいずれか、またはこれらの組み合わせを含み、さらに、酸素飽和度（ＳｐＯ２）を含んでもよい。
　ＢＭＩは身長ｈ[m]及び体重ｗ[kg]より、下記の計算式により算出される。
　ＢＭＩ＝w/ h2 [kg/m2]
　循環血液量は、下記のナドラー式によって算出される。
＜男性の場合＞循環血液量 [L] ＝０．３６６９×h 3 + ０．０３１２９×w + ０．６０４１
＜女性の場合＞循環血液量 [L] ＝０．３５６１×h 3 + ０．０３３０８×w  + ０．１８３３
　血圧には、収縮期血圧、拡張期血圧、脈圧、および平均動脈圧のいずれか、またはその組み合わせが含まれる。
　脈圧は下記の式により計算される。
　脈圧＝収縮期血圧―拡張期血圧
　平均動脈圧は下記の式により計算される。
　平均動脈圧＝拡張期血圧＋脈圧×１／３
　脈波データは、脈波計やパルスオキシメーターにより、指など体の突出した部位にRed LEDから赤色光(～660nm)、IR LEDから近赤外光(～905nm) を照射し、その透過光をフォトトランジスタ測定することで測定される。
　脈波データには、脈拍、弾性指数、末梢血管抵抗、加速度脈波、b/a、e/a、-d/a、高沢式加速度脈波加齢指数、駆出率、LVET、DEI（重拍弾性指数）のいずれか、またはこれらの組み合わせが含まれる。
　ここで弾性指数とは、身長を、指尖容積脈波において収縮期のピークから拡張期のピークが検出されるまでの時間で割って得られる数値である。末梢血管抵抗は、平均動脈圧／心拍出量×８０によって求められる。DEI（重拍弾性指数）は、拡張期血管の弾性を示す指標であり、PWV測定装置により測定可能である。0.3～0.7が正常、0.3以下で高血圧または動脈硬化の可能性があり、0.7以上では急性不安神経症の可能性が示唆される。
　加速度脈波は、指尖容積脈波（photoplethysmogram：PTG）の二次微分波（second derivative of photoplethysmogram：SDPTG）である。加速度脈波は、初期陽性波（ａ波）、初期陰性波（ｂ波）、収縮中期再上昇波（ｃ波）、収縮後期再降下波（ｄ波）、および拡張初期陽性波（ｅ波）より構成され、上記b/a、e/a、-d/aは、各波高の割合より計算される。加齢に伴ってb/aの上昇およびc/a、-d/a、e/aの低下が認められるため、高沢式加速度脈波加齢指数（b-c-d-e）/aにより血管の加齢を評価することができる。駆出率は、各心拍で心室から送られる血液の割合であり、加速度脈波加齢指数に比例する。LVETは左心室駆出時間であり、大動脈弁解放後、左室内の血液が大動脈へ駆出される時間である。
　心電図データは、電極による心電図法（ECG）または光電脈波法（PPG）により測定可能である。
　心電図データには、呼吸数、心拍数、RR間隔、RR間隔の標準偏差、MxDMn比率、低周波帯のパワースペクトル、高周波帯のパワースペクトル、心拍数変動指標LF/HF、トータルパワーのいずれか、またはその組み合わせが含まれる。ここで、RR間隔とは、心電図のQRS波から次のQRS波までの間隔である。MxDMn比率とは、時間内の最も長いRR間隔と最も短いRR間隔の比であり、不規則な心拍の指数である。トータルパワーとは、２分間の分間測定における周波数0～0.4Hｚ（VLF,LF,HF）のパワースペクトルのトータルパワーの計算値である。この値は交感神経活動が主に占める自律神経系活動全体を反映する。
　心電図よりパワースペクトル密度を計算することで、高周波のパワースペクトル割合（０．１８７５～０．５０Ｈｚ：HF）、低周波のパワースペクトル割合（０．０５～０．１８７５Ｈｚ：LF）、LF/HF比、超低周波帯のパワースペクトル割合（０～０．０５Ｈｚ：VLF）を計算することができる。
　生体のインピーダンス（コンダクタンス）は、例えば、両脚、両手、左右額の計６か所の電極のうち、２か所の間に微弱な電気を流すことで測定することができる。６か所の電極中で２か所に電流を流すと、（１）アノード/カソードコンダクタンス（μS）、（２）カソード/アノードのコンダクタンス（μS）、（３）上記（１）で測定されたコンダクタンスと上記（２）で測定されたコンダクタンスの差分（デルタSCRA－SCRC）、（４）電気伝導度（μS/m）を測定することが可能である。また、筋量、体脂肪量、総水分量、位相角、抵抗値も同時に測定可能である。また、右手⇔左手間、右額⇔左額間で通電している際の誘電率（μSi）も可能である。６か所の電極から２２パターンの通電性を利用し、生体のインピーダンス（コンダクタンス）を測定することが好ましい。
　生体インピーダンスには、体脂肪量（kg）、体脂肪量（%）、除脂肪体重、除脂肪率、筋量、総水分量（kg）、総水分量（%）、細胞内水分量（%)、心拍出量、1額左側-2左手/SCR A、1額左側-2左手/デルタSCR C-SCR A、3額右側-4右手/SCRC 、5左手-6左足/デルタSCR C-SCR A、9額左側-10額右側/SCRC 、13左足-14右足/SCR C、13左足-14右足/SCR A、15右手-16額左側/デルタSCR C-SCR A、15右手-16額左側SCR C、19右足-20左手/デルタSCR C-SCR A 、ESG2＋4＋15＋17（μS/m）、ESG6＋13＋19（%）、ESG 6＋8＋19＋21 （%）、ESG6＋8＋19＋21（μS/m）、ESG 9＋10（μS/m）、ESG9＋10（%）、左足のコンダクタンス、R（Ω）、位相角、額経路の誘電率、額経路への電気伝導度(9)、片手-片手経路の誘電率手から手への電気伝導度（11,12）、1回拍出量（心拍出量÷心拍数）、右額のインピーダンス指数（3額右側-4右手/SCRCと9額左側-10額右側/SCRCの和）のいずれか一つまたはこれらの組み合わせを含む。
　ここで、SCRとは皮膚コンダクタンス反応（skin conductance response）の略称であり、ESGとはエレクトロスキャングラムの略称である。ESG2＋4＋15＋17の「＋」は、身体に添付されたどの電極で測定されたかを意味する。例えば、ESG2＋4＋15＋17とは、図2に示すように左手から左額へ通電した際の左手、右手から右額へ通電した際の右手、右手から左額へ通電した際の右手、左手から右額へ通電した際の左手で測定されたコンダクタンスの平均値を意味する。これらコンダクタンスに関しては非特許文献１に詳細な記載がある。なお、ESG（エレクトロスキャングラム）の測定方法については、非特許文献２、３に詳細が記載されている。
　「1額左側-2左手/SCR A」は、「1額左側」がカソード、「2左手」がアノード方向となるように電気が流れた際に測定された経路のコンダクタンス（または伝導性）であり、「5左手-6左足/デルタSCR C-SCR A」は、「5左手」、「6左足」間でアノード-カソード、カソード-アノードとなるよう通電した際に測定されたコンダクタンスの差分である。
　ＢＭＩおよび血圧については、身長体重計、血圧計によって測定可能である。また、非侵襲な生体情報には、ＳｐＯ２および心拍出量より計算される酸素運搬量が含まれてもよい。
　測定された非侵襲な生体情報は、有線または無線の通信ネットワークを通じてHbA1cリスク推定装置３０へ出力される。生体情報測定装置２０は、備え付けの測定装置であっても、ウェアラブル端末のように可搬可能な測定装置であってもよい。
　HbA1cリスク推定装置３０は、第1取得部３１と、第2取得部３２と、ユーザデータ記憶部３３と、訓練データ記憶部３４と、学習処理部３５と、推定モデル記憶部３６と、推定処理部３７と、推定データ記憶部３８とを備える。第1取得部３１は、端末装置１０よりユーザの属性情報を取得する。また、第2取得部３２は、生体情報測定装置２０より、ユーザの非侵襲な生体情報を取得する。
　ユーザデータ記憶部３３は、第1取得部３１及び第2取得部３２より取得したユーザの属性情報と非侵襲な生体情報を記憶する。
　訓練データ記憶部３４は、機械学習のための訓練データセットとして、事前に取得した複数の被験者の属性情報と、非侵襲な生体情報と、血液検査により得られたHbA1c等の検体検査情報より構成される訓練データセットが複数記憶されている。なお、検体検査情報には、血液、尿、便等から得られた検査情報をさらに含めても良い。
　学習処理部３５は、訓練データ記憶部３４に記憶されている訓練データセットを取得し、訓練データセットを用いてHbA1cリスクの推定モデルを作成する。具体的には、HbA1cリスクを推定する場合、取得した訓練データセットをXGBoostなどの勾配ブースティング、ニューラルネットワーク、およびロジスティック回帰、またはこれらの学習結果のアンサンブル学習による機械学習により、属性情報および非侵襲な生体情報と、HbA1cリスクの関係を学習する。
　推定モデル記憶部３６は、学習処理部３５により生成したHbA1cリスクの推定モデルを記憶する。
　非侵襲な生体情報には、ＢＭＩ、循環血液量、血圧、脈波データ、心電図データおよび生体インピーダンスのいずれか、またはその組み合わせが含まれる。さらに、必要に応じて酸素飽和度（ＳｐＯ２）が含まれる。
　推定処理部３７は、学習処理部３５により生成した推定モデルを用いて、所定のユーザの属性情報および／または非侵襲な生体情報に基づいて、ユーザのHbA1cリスクを推定する。そして、HbA1cリスク推定値は推定データ記憶部３８に記憶される。
　表示装置３９は、ユーザの属性情報や非侵襲な生体情報と共に、HbA1cリスク推定値を表示することが可能である。なお、これらのデータは、ユーザが保有する端末装置１０に表示されてもよい。
＜装置のハードウェア構成＞
　図３は、HbA1cリスク推定装置３０のハードウェア構成図である。図３に示すように、HbA1cリスク推定装置３０は、一つまたは複数のプロセッサ３０１と、メモリ３０２と、ストレージ３０３と、入出力ポート３０４と、通信ポート３０５とを有するコンピュータ３００により構成される。プロセッサ３０１は、プログラムを実行することにより本実施形態に係るHbA1c推定に関する処理を行う。メモリ３０２は、プログラム及びプログラムの演算結果を一時的に記憶する。ストレージ３０３は、HbA1cリスク推定装置３０による処理を実行するプログラムを格納している。ストレージ３０３は、コンピュータにより読み取り可能であればどのようなものであっても良く、例えば、記録媒体（磁気ディスク、光ディスク等）、ランダムアクセスメモリー、フラシュタイプメモリー、リードオンリーメモリー等の種々の記録媒体を利用することができる。入出力ポート３０４は、端末装置１０および生体情報測定装置２０より情報を入力したり、表示装置３９へHbA1c推定値の出力を行う。通信ポート３０５は、図示しない他のコンピュータ等の情報端末との間でデータの送受信を行う。通信の方法は、無線通信、有線通信を利用することができる。なお、HbA1cリスク推定装置３０は、市販のデスクトップPCまたはノートPCにより実装することが可能であり、推定モデルを用いたHbA1cリスクの推定値の算出に要する時間は数秒である。
　なお、HbA1cリスク推定装置３０のプロセッサ３０１においては、動作する際に、上述の第１取得部３１、第２取得部３２、学習処理部３５、推定処理部３７等が機能する。
　ＢＭＩ（ボディマス指数）、血圧、脈波データ、心電図データ、および生体インピーダンス等を含む非侵襲生体データに基づき機械学習によりHbA1cリスク推定モデルを生成することで、血液検査を行うことなくHbA1cリスクの推定を行うことが可能となる。また、下記に示すように、非侵襲生体データに含まれる各データの数を限定した場合であっても、HbA1cが正常か否かの判定を行うことも可能である。
＜機械学習によるHbA1cリスク推定モデルの生成＞
　図４は、機械学習によるHbA1cリスク推定モデルの生成実行手順を示すフローチャートである。
　ステップST101では、学習処理部３５は、入力データの前処理を行う。具体的には、学習処理部３５は、血液検査により得られたHbA1cについて、HbA1c５．５％以下を0（リスク無し）、５．６％以上を1（リスク有り）に変換する。また、学習処理部３５は、0と分類された人数と１と分類された人数が乖離し不均衡な場合は、学習データに対してSMOTE （Chawla, NV. et al. 2002）を適用し、人工的にトレーニングサンプルを生成してもよい。すなわち、入力データ（訓練データセット）において、血液検査により得られたHbA1cに基づいてHbA1cリスクの有無を示すラベル（例えば、上述の０または１）を付加し、当該ラベルにおいて、HbA1cリスク有り（ラベル：１）の数と、HbA1cリスク無し（ラベル：０）の数との差が所定の値以上の場合は、この差を小さくするように訓練データセットにおけるサンプルデータを増加（生成）させてもよい。なお、性別については、女性を０、男性を１とした。また、年齢、体重、BMI、循環血液量については、それぞれの値をそのまま使用する。
　ステップST102では、学習処理部３５は、ロジスティック回帰、ニューラルネットワーク、および勾配ブースティング決定木による機械学習を行い、ステップST103では、ロジスティック回帰、ニューラルネットワーク、および勾配ブースティング決定木による機械学習による学習結果について、さらにアンサンブル学習（投票）を行う。
　図５は、HbA1cリスク推定モデルとして使用するニューラルネットワーク（NN）の構造である。長方形はデータの変換を行う層群を示し、角丸長方形は入出力データを示している。Dは検査項目の数である。NNの構造は、４つの層群（Compile A1、Compile A2、CompileA3、Compile B1）を介し、D次元、９６次元、９６次元、９６次元、１次元の順にデータを変換している。
　また、Compile A1、A2、A3は、全結合処理を行う全結合層「Linear」、正則化を行う「Kernel Regularizer」、ReLU処理を行う「ReLU 」層を含み、層群Compile B1は、全結合処理を行う全結合層「Linear」、および最適化処理を行う「Adagrad」層を含む。層群Compile A1の全結合層の入力ユニットが入力層、層群Compile B1の出力ユニットが出力層、これらの間にあるユニットが中間層（隠れ層）に相当する。中間層には一部の入力値を０に抑え過学習を防止するDropout層が含まれる。
　ロジスティック回帰による機械学習では、例えばPythonのオープンソース機械学習ライブラリであるScikit-learnにて提供されているLogistic Regressionを用いることができる。また、必要に応じて、主成分分析により次元数を圧縮してもよい。血液検査により得られたHbA1cリスクを機械学習により推定されたHbA1cリスク推定値と比較し、f1スコアが最大となるようにLogistic Regressionの各パラメーター（C、正則化方法、max_iter、solber）を調整した。ここで、Cとは正則化の強度を決定するトレードオフパラメータであり、値が大きいほど正則化強度は弱くなる。正則化方法とは、L1正則化またはL2正則化を意味し、これを選択する。max_iterは、学習を反復する最大回数である。solberでは交差エントロピー誤差を最小化する収束方法（例えば、L-BFGS法、ニュートンCG法、liblinear、sagおよびsaga）を選択する。なお、下記の実施例１、２ではliblinear法を選択した。
　勾配ブースティング回帰木による機械学習では、XGBoost、CatBoost、LightBGM等のソフトウェアライブラリを用いることができる。血液検査により得られたHbA1c値を所定の閾値と比較することにより得られたリスク値（０：リスク無し、１：リスク有り）を機械学習により推定されたHbA1cリスク値と比較し、f1スコアが最大となるようにXGBoostの各パラメータ（max_depth、subsample、colsample_bytree、learning_rate）を調整した。max_depthは決定木の深さであり、subsampleは各木においてランダムに抽出される標本の割合であり、colsample_bytreeは各木においてランダム抽出される列の割合であり、learning_rateは学習率を示す。ただし、max_depthは1～10、subsampleは0.1～1.0、colsample_bytreeは0.3～1.0、learning_rateは0.1～0.7の範囲で調整を行った。
　学習処理部３５は、上述の学習処理によって生成されたHbA1cリスク推定モデルを、推定モデル記憶部３６に記憶する。
　なお、上述の機械学習アルゴリズムは一例であって、これらに限定されるものではない。
＜HbA1cリスク推定モデルを用いたHbA1cリスクの推定＞
　図６に示すように、ステップST201では、HbA1cリスク推定装置３０の第１取得部３１は、端末装置１０よりユーザの属性情報を取得する。ステップST202では、HbA1cリスク推定装置３０の第２取得部３２は、ユーザの非侵襲な生体情報を取得する。そして、ユーザの属性情報と非侵襲な生体情報をユーザデータ記憶部３３に記憶する。そして、ステップST203では、推定モデル記憶部３６に格納されているHbA1cリスク推定モデルを用いて、推定処理部３７によりクラス０（リスク無し）またはクラス１（リスク有り）に属する確率、すなわちHbA1cリスク推定値を算出する。ステップST204では、算出されたHbA1cリスク推定値は推定データ記憶部３８に記憶され、ステップST205では、HbA1cリスク推定値を表示装置３９等の外部端末へ出力し、表示させる。
＜実施例（HbA1cリスク推定）＞
　以下にHbA1cリスク推定の実施例を記載する。但し、本発明におけるHbA1cリスク推定の態様は、以下の実施例に限定されるものではない。
　属性情報には、ＩＤ、氏名、年齢、性別のいずれか、またはこれらの組み合わせが含まれ、非侵襲な生体情報には、BMI、循環血液量、血圧、脈波データ、心電図データ、生体インピーダンス、および酸素飽和度（ＳｐＯ２）のいずれか、またはこれらの組み合わせが含まれる。BMIの算出基準となる身長、体重については、それぞれ身長計および体重計により計測し、血圧については血圧計により測定した。また、脈波データ、心電図データ、生体インピーダンス、酸素飽和度（ＳｐＯ２）については、エステックBC-3（両備システムズ）により測定した。なお、エステックBC-3に替えて、市販の脈波計、心電図計、インピーダンス測定装置、パルオキシメーターを組み合わせて使用してもよい。また、所定のウェアラブル端末を用いて、上述の非侵襲な生体情報を取得してもよい。
　生体インピーダンス（コンダクタンス）の測定は、両足、両手、左右額の計6か所の電極のうち、２か所の間に微弱な電気を流すことで測定した。電圧と電流は1.28V, 200 μAとし、コンダクタンスは1秒につき32ミリ秒間測定した。６か所の電極中で２か所に電流を流し、（１）アノード/カソードコンダクタンス（μS）、（２）カソード/アノードのコンダクタンス（μS）、（３）上記（１）で測定されたコンダクタンスと上記（２）で測定されたコンダクタンスの差分（デルタSCRA－SCRC）、（４）電気伝導度（μS/m）を測定した。
　また、筋量、体脂肪量、総水分量、位相角、および抵抗値を測定し、右手⇔左手間、右額⇔左額間で通電している際の誘電率（μSi）も測定した。
　エステックBC-3により、各被験者について、脈波データ、心電図データ、生体インピーダンス、および酸素飽和度（ＳｐＯ２）の測定を２分間行った。測定に際しては、被験者の左手人差し指には心電図、脈波計、およびパルオキシメーターの機能を備えたデバイスを、額には電極を２個装着し、椅子に座った状態で両手両足を電極板の上に置かせた。
＜学習モデル１＞
　学習モデル１では、図４におけるロジスティック回帰のみを用いて機械学習を行わせた。ロジスティック回帰のみを用いるため、図４におけるアンサンブル学習は行わない。
　この際、非侵襲な生体データとして下記に示すデータを選択し、使用した。なお、このモデルでは属性情報は使用しない。（A）非侵襲な生体データ・BMI・血圧・・・拡張期血圧、脈圧・脈波データ・・・弾性指数、脈拍・心電図データ・・・心拍数、LF/HF・生体インピーダンス・・・5左手-6左足/デルタSCRC-SCRA、左足のコンダクタンス、R（Ω）
　また、非侵襲な生体データには、脈圧／脈拍のデータ、生体インピーダンスに含まれる心拍出量、および酸素飽和度（ＳｐＯ２）より推定される酸素運搬量がさらに含まれる。
＜学習モデル２＞
　学習モデル２では、図４に示すように、ロジスティック回帰、ニューラルネットワークおよび勾配ブースティングによる機械学習を行わせ、これらの学習結果について投票によるアンサンブル学習を行うことで、血中HbA1cリスク推定モデルを作成した。この際、非侵襲な生体データとして下記に示すデータを選択し、使用した。（Ａ）属性情報・年齢、性別（Ｂ）非侵襲な生体データ・体重、BMI、循環血液量・血圧・・・収縮期血圧、拡張期血圧、平均動脈圧・脈波データ・・・弾性指数・心電図データ・・・心拍数、LF/HF・生体インピーダンス・・・筋量（％）、総水分量（kg）、総水分（%）、心拍出量、3額右側-4右手/SCRC、9額左側-10額右側/SCRA、右額のインピーダンス指数（3額右側-4右手/SCRC　と　9額左側-10額右側/SCRCの和）、左足のコンダクタンス、R（Ω）、１回拍出量
　また、非侵襲な生体データには、生体インピーダンスに含まれる心拍出量、および酸素飽和度（ＳｐＯ２）より推定される酸素運搬量がさらに含まれる。
＜実施例１＞
　実施例１では、被験者３２１名の（１）属性情報、（２）身長体重計、血圧計、およびエステックBC-3により測定した非侵襲な生体情報、および（３）非侵襲な生体情報測定と同日に行われた血液検査により得られたHbA1cの訓練データセットを用いて、上記学習モデル１の機械学習を通じて、HbA1cリスク推定モデルを生成した。
　そして、HbA1cリスク推定モデルの推定精度について、推定結果をROC_AUC曲線により評価した。その結果、ROC_AUCは0.78を示し、良好な分類ができていることを示す0.7を超えた。実施例１の推定結果のROC_AUC曲線を図７に示す。
＜実施例２＞
　実施例２では、のべ被験者７１２名の（１）属性情報、（２）身長体重計、血圧計、およびエステックBC-3により測定した非侵襲な生体情報、および（３）非侵襲な生体情報測定と同日に行われた血液検査により得られたHbA1cの訓練データセットを用いて、上記学習モデル２の機械学習を通じて、HbA1cリスク推定モデルを生成した。
　そして、HbA1cリスク推定モデルの推定精度について、推定結果をROC_AUC曲線によって評価した。その結果、ROC_AUCは0.80を示し、極めて良好な分類ができていることを示す0.8と同値となった。実施例２の推定結果のROC_AUC曲線を図８に示す。
（変形例）
　上述の実施形態では、学習モデル１または学習モデル２を用いてHbA1cリスクを推定する例について説明したが、複数の学習モデルを用いて、HbA1cリスクを推定してもよい。
　これにより、１つの学習モデルを用いてHbA1cリスクを推定するよりも、高精度に推定できる。
　また、上述の実施形態および実施例において、学習モデルのデータセットとして、BMIを用いると推定精度が向上する事例が見られたことから、当該BMIを推定する機能部がHbA1cリスク推定装置に設けられてもよい。
　BMIは、一般に、ウェアラブル端末等で取得されず、ユーザの入力した身長や体重から求めているが、当該BMIを推定することで、生体情報のみを取得することにより、上述のHbA1cリスクを取得することができるため、ユーザにとって利便性が向上する。
　BMIの推定方法は特に限定されないが、例えば、BMIは、ユーザの腹部の（所定の位置の）傾きと相関があることが知られている。そこで、例えば、ユーザの腹部に所定の加速度センサを設け（または加速度センサを備えるリストバンド型のウェアラブル端末等を腹部に当てて）、当該加速度センサから出力されるデータに基づいて、腹部の傾きを求め、BMIを推定してもよい。
　また、上述と同様に脈波データ、酸素飽和度についても、脈波センサまたは中酸素濃度センサを備えるリストバンド型のウェアラブル端末等を用いて推定してもよい。
　また、上述と同様に血圧についても、リストバンド型のウェアラブル端末等を用いて推定してもよい。これは、心拍によって動脈を通じて送られる脈波の速度と血圧との相関があることが知られていることから、当該心拍によって動脈を通じて送られる脈波の速度を測定する所定のセンサを用いて、血圧を推定してもよい。
　また、上述と同様に心電図データについても、リストバンド型のウェアラブル端末等を用いて推定してもよい。例えば、リストバンド型のウェアラブル端末の表示面と反対側の面に設けられた電極と、表示面側に設けられた電極から得られるデータに基づいて心電図データを推定することができる。具体的には、上記反対側面に設けられた電極にリストバンド型のウェアラブル端末を装着した手（例えば、左手）の手首が接し、上記表示面側に設けられた電極に当該装着した手とは反対側の手（例えば、右手）の指先が接することで得られるデータから心電図データを推定するとよい。
　また、上述と同様に生体インピーダンスについても、各種電極を備えるリストバンド型のウェアラブル端末等を用いて、生体インピーダンスを推定してもよい。例えば、リストバンド型のウェアラブル端末等を用いて、胸部と手首から得られた生体情報に基づいて、生体インピーダンスを推定するとよい。
　なお、ウェアラブル端末で取得可能な生体情報と、上述の推定対象の生体情報（BMI、血圧、脈波データ、心電図データ、および生体インピーダンス、酸素飽和度のうち少なくとも１つ以上の生体情報）とを教師データとして、種々の機械学習アルゴリズムを用いて生成した分類器を用いて、上述の推定対象の生体情報を推定してもよい。
　なお、この場合、上述の第２取得部は、推定された生体情報を取得してもよい。
（その他）
　また例えば、上述した一連の処理は、ハードウェアにより実行させることもできるし、ソフトウェアにより実行させることもできる。換言すると、上述の機能的構成は例示に過ぎず、特に限定されない。即ち、上述した一連の処理を全体として実行できる機能が情報処理システムに備えられていれば足り、この機能を実現するためにどのような機能ブロックを用いるのかは特に上述の例に限定されない。また、機能ブロックの存在場所も、図１に特に限定されず、任意でよい。例えば、サーバの機能ブロックを他の端末や装置等に移譲させてもよい。逆に他の端末や装置の機能ブロックをサーバ等に移譲させてもよい。また、一つの機能ブロックは、ハードウェア単体で構成してもよいし、ソフトウェア単体で構成してもよいし、それらの組み合わせで構成してもよい。
　一連の処理をソフトウェアにより実行させる場合には、そのソフトウェアを構成するプログラムが、コンピュータ等にネットワークや記録媒体からインストールされる。コンピュータは、専用のハードウェアに組み込まれているコンピュータであってもよい。また、コンピュータは、各種のプログラムをインストールすることで、各種の機能を実行することが可能なコンピュータ、例えばサーバの他汎用のスマートフォンやパーソナルコンピュータであってもよい。
　このようなプログラムを含む記録媒体は、ユーザ等にプログラムを提供するために装置本体とは別に配布される図示せぬリムーバブルメディアにより構成されるだけでなく、装置本体に予め組み込まれた状態でユーザ等に提供される記録媒体等で構成される。プログラムはネットワークを介して配信可能であることから、記録媒体は、ネットワークに接続された、或いは接続可能なコンピュータに搭載、或いはアクセス可能なものであってもよい。
　なお、本明細書において、記録媒体に記録されるプログラムを記述するステップは、その順序に沿って時系列的に行われる処理はもちろん、必ずしも時系列的に処理されなくとも、並列的あるいは個別に実行される処理をも含むものである。また、本明細書において、システムの用語は、複数の装置や複数の手段等より構成される全体的な装置を意味するものとする。
１　　　HbA1cリスク推定システム
１０　　端末装置
２０　　生体情報測定装置
３０　　HbA1cリスク推定装置
３１　　第１取得部
３２　　第２取得部
３３　　ユーザデータ記憶部
３４　　訓練データ記憶部
３５　　学習処理部
３６　　推定モデル記憶部
３７　　推定処理部
３８　　推定データ記憶部
３９　　表示装置
３００　コンピュータ
３０１　プロセッサ
３０２　メモリ
３０３　ストレージ
３０４　入出力ポート
３０５　通信ポート

Claims (12)

  1. 　所定のユーザの属性情報および非侵襲な生体情報を取得する情報取得部と、
    　HbA1cリスク推定モデルを記憶する推定モデル記憶部と、
    　前記HbA1cリスク推定モデルを用いて、前記所定のユーザの属性情報および／または非侵襲な生体情報に基づいて、前記所定のユーザのHbA1cリスク推定値を計算する推定処理部と、を備えることを特徴とするHbA1cリスク推定装置。
  2. 　前記属性情報は、年齢、および性別のいずれか、またはこれらの組み合わせを含み、
    　前記非侵襲な生体情報は、BMI、循環血液量、血圧、脈波データ、心電図データ、および生体インピーダンスのいずれか、またはこれらの組み合わせを含む、ことを特徴とする請求項１に記載のHbA1cリスク推定装置。
  3. 　前記HbA1cリスク推定値の推定精度は、リスクの有無をROC_AUCが０．７以上で分類可能な精度である、ことを特徴とする請求項１または２に記載のHbA1cリスク推定装置。
  4. 　訓練データセットを記憶する訓練データ記憶部と、
    　前記訓練データセットに基づいて、機械学習により前記HbA1cリスク推定モデルを生成する学習処理部と、をさらに備えることを特徴とする請求項１から３のいずれか一項に記載のHbA1cリスク推定装置。
  5. 　前記訓練データセットは、被験者の属性情報、非侵襲な生体情報、および血液から測定されたHbA1c測定値を含む、ことを特徴とする請求項４に記載のHbA1cリスク推定装置。
  6. 　前記非侵襲な生体情報は、酸素飽和度（ＳｐＯ２）をさらに含む、ことを特徴とする請求項５に記載のHbA1cリスク推定装置。
  7. 　前記学習処理部は、
    　　前記訓練データセットにおいて、血液から測定されたHbA1c測定値に基づいて前記HbA1cリスクの有無を示すラベルを付加し、
    　　前記ラベルにおいて、前記HbA1cリスク有りの数と、前記HbA1cリスク無しの数との差が所定の値以上の場合は、前記差を小さくするように前記訓練データセットにおけるサンプルデータを増加させる、ことを特徴とする請求項４から６のいずれか一項に記載のHbA1cリスク推定装置。
  8. 　前記学習処理部は、異なる種類の訓練データセットに基づいて、それぞれ機械学習により第１HbA1cリスク推定モデルおよび第２HbA1cリスク推定モデルを生成し、
    　前記推定処理部は、前記第１HbA1cリスク推定モデルおよび前記第２HbA1cリスク推定モデルを用いて、前記所定のユーザのHbA1cリスク推定値を計算する、ことを特徴とする請求項４から７のいずれか一項に記載のHbA1cリスク推定装置。
  9. 　前記生体情報に含まれるBMI、血圧、脈波データ、心電図データ、および生体インピーダンス、酸素飽和度のうち少なくとも１つ以上の生体情報を推定する生体情報推定部をさらに備え、
    　前記情報取得部は、前記生体情報推定部によって推定された生体情報を前記所定のユーザの生体情報として取得する、ことを特徴とする請求項１から８のいずれか一項に記載のHbA1cリスク推定装置。
  10. 　請求項１から９のいずれか一項に記載の前記HbA1cリスク推定装置と、
    　非侵襲な生体情報を測定する生体情報測定装置と、を備えることを特徴とする非侵襲HbA1cリスク推定システム。
  11. 　被験者の属性情報、非侵襲な生体情報、および血液から測定されたHbA1c測定値を含む訓練データセットを記憶するステップと、
    　前記訓練データセットに基づいて、機械学習によりHbA1cリスク推定モデルを生成するステップと、
    　前記HbA1cリスク推定モデルを用いて、所定のユーザの属性情報および／または非侵襲な生体情報に基づいて、前記所定のユーザのHbA1cリスク推定値を計算するステップと、を含む、HbA1cリスク推定方法。
  12. 　被験者の属性情報、非侵襲な生体情報、および血液から測定されたHbA1cの測定値を含む訓練データセットを記憶するステップと、
    　前記訓練データセットに基づいて、機械学習によりHbA1cリスク推定モデルを生成するステップと、
    　前記HbA1cリスク推定モデルを用いて、所定のユーザの属性情報および／または非侵襲な生体情報に基づいて、前記所定のユーザのHbA1cリスク推定値を計算するステップと、をコンピュータに実行させるプログラム。

     
PCT/JP2022/015095 2021-08-06 2022-03-28 ＨｂＡ１ｃリスク推定装置、ＨｂＡ１ｃリスク推定方法及びプログラム WO2023013164A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202280006142.1A CN117693313A (zh) 2021-08-06 2022-03-28 HbA1c风险推定装置、HbA1c风险推定方法以及程序
US18/041,939 US20230402185A1 (en) 2021-08-06 2022-03-28 HbA1c RISK ESTIMATION DEVICE, HbA1c RISK ESTIMATION METHOD, AND COMPUTER PROGRAM

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
JP2021130045 2021-08-06
JP2021-130045 2021-08-06
JP2022-040636 2022-03-15
JP2022040636A JP2023024256A (ja) 2021-08-06 2022-03-15 ＨｂＡ１ｃリスク推定装置、ＨｂＡ１ｃリスク推定方法及びプログラム

Publications (1)

Publication Number Publication Date
WO2023013164A1 true WO2023013164A1 (ja) 2023-02-09

Family

ID=85155665

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2022/015095 WO2023013164A1 (ja) 2021-08-06 2022-03-28 ＨｂＡ１ｃリスク推定装置、ＨｂＡ１ｃリスク推定方法及びプログラム

Country Status (3)

Country Link
US (1) US20230402185A1 (ja)
TW (1) TW202307430A (ja)
WO (1) WO2023013164A1 (ja)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013508687A (ja) * 2009-10-15 2013-03-07 クレッシェンド バイオサイエンス インコーポレイテッド 炎症性疾患活動性を測定およびモニタリングするためのバイオマーカーおよび方法
WO2020183609A1 (ja) * 2019-03-12 2020-09-17 三菱電機株式会社 移動体制御装置および移動体制御方法
WO2020203728A1 (ja) * 2019-03-29 2020-10-08 株式会社タニタ 健康情報提供システム及び健康情報提供プログラム
JP2021047056A (ja) 2019-09-17 2021-03-25 株式会社ファンケル 非侵襲的なＨｂＡ１ｃ値の推定方法
JP2021072100A (ja) * 2019-10-28 2021-05-06 株式会社デンソー 情報処理装置、統合モデル生成方法、及び統合モデル生成プログラム

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013508687A (ja) * 2009-10-15 2013-03-07 クレッシェンド バイオサイエンス インコーポレイテッド 炎症性疾患活動性を測定およびモニタリングするためのバイオマーカーおよび方法
WO2020183609A1 (ja) * 2019-03-12 2020-09-17 三菱電機株式会社 移動体制御装置および移動体制御方法
WO2020203728A1 (ja) * 2019-03-29 2020-10-08 株式会社タニタ 健康情報提供システム及び健康情報提供プログラム
JP2021047056A (ja) 2019-09-17 2021-03-25 株式会社ファンケル 非侵襲的なＨｂＡ１ｃ値の推定方法
JP2021072100A (ja) * 2019-10-28 2021-05-06 株式会社デンソー 情報処理装置、統合モデル生成方法、及び統合モデル生成プログラム

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
MAAREK A.: "Electro interstitial scan system: assessment of 10 years of research and development", MED DEVICES (AUCKL, vol. 5, 2012, pages 23 - 30
PSYCHOLOGY RESEARCH AND BEHAVIOR MANAGEMENT, vol. 4, 2011, pages 81 - 86
R. N. CHUAY. W. HAUC. M. TIEWW. L. HAU: "Investigation of Attention Deficit,Hyperactivity Disorder Assessment Using Electro Interstitial Scan Based on Chronoamperometry Technique", IEEE ACCESS, vol. 7, 2019, pages 144679 - 144690, XP011750230, DOI: 10.1109/ACCESS.2019.2938095

Also Published As

Publication number Publication date
US20230402185A1 (en) 2023-12-14
TW202307430A (zh) 2023-02-16

Similar Documents

Publication Publication Date Title
JP2022153288A (ja) 血糖値推定装置、血糖値推定方法及びプログラム
WO2023013164A1 (ja) ＨｂＡ１ｃリスク推定装置、ＨｂＡ１ｃリスク推定方法及びプログラム
WO2023013166A1 (ja) 尿酸値推定装置、尿酸値推定方法及びプログラム
WO2023013167A1 (ja) クレアチニンリスク推定装置、クレアチニンリスク推定方法及びプログラム
WO2022210570A1 (ja) 血糖値推定装置、血糖値推定方法及びプログラム
WO2022270098A1 (ja) 血中中性脂肪推定装置、血中中性脂肪推定方法及びプログラム
WO2022270099A1 (ja) コレステロールリスク推定装置、コレステロールリスク推定方法及びプログラム
WO2023013165A1 (ja) γＧＴ推定装置、γＧＴ推定方法及びプログラム
EP4383274A1 (en) Hba1c risk estimation device, hba1c risk estimation method, and program
JP2023024256A (ja) ＨｂＡ１ｃリスク推定装置、ＨｂＡ１ｃリスク推定方法及びプログラム
EP4383264A1 (en) yGT ESTIMATION DEVICE, yGT ESTIMATION METHOD, AND PROGRAM
EP4383275A1 (en) Creatinine risk estimation device, creatinine risk estimation method, and program
EP4383265A1 (en) Uric acid value estimation device, uric acid value estimation method, and program
Shamim et al. Diagnostic accuracy of smartphone-connected electrophysiological biosensors for prediction of blood glucose level in a type-2 diabetic patient using machine learning: A pilot study
JP2023024258A (ja) 尿酸値推定装置、尿酸値推定方法及びプログラム
JP2023024259A (ja) クレアチニンリスク推定装置、クレアチニンリスク推定方法及びプログラム
JP2023002452A (ja) 血中中性脂肪推定装置、血中中性脂肪推定方法及びプログラム
JP2023024257A (ja) γＧＴ推定装置、γＧＴ推定方法及びプログラム
JP2023002453A (ja) コレステロールリスク推定装置、コレステロールリスク推定方法及びプログラム
Vakulenkoa et al. New application of blood pressure monitor with software environment Oranta-AO based on arterial oscillography methods
CN117693313A (zh) HbA1c风险推定装置、HbA1c风险推定方法以及程序
CN117677849A (zh) 肌酐风险推定装置、肌酐风险推定方法以及程序
Manamperi et al. A robust neural network-based method to estimate arterial blood pressure using photoplethysmography
CN117678033A (zh) 尿酸值推定装置、尿酸值推定方法以及程序
CN117678032A (zh) γGT推定装置、γGT推定方法以及程序

Legal Events

Date Code Title Description
WWE Wipo information: entry into national phase

Ref document number: 202280006142.1

Country of ref document: CN

121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 22852593

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

ENP Entry into the national phase

Ref document number: 2022852593

Country of ref document: EP

Effective date: 20240306