KR20160107508A - Method for Assessinging Residual Life Using Biological Age - Google Patents

Method for Assessinging Residual Life Using Biological Age Download PDF

Info

Publication number
KR20160107508A
KR20160107508A KR1020150030331A KR20150030331A KR20160107508A KR 20160107508 A KR20160107508 A KR 20160107508A KR 1020150030331 A KR1020150030331 A KR 1020150030331A KR 20150030331 A KR20150030331 A KR 20150030331A KR 20160107508 A KR20160107508 A KR 20160107508A
Authority
KR
South Korea
Prior art keywords
age
examination
value
subject
survival probability
Prior art date
Application number
KR1020150030331A
Other languages
Korean (ko)
Other versions
KR101669526B1 (en
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
Application filed by 주식회사 바이오에이지 filed Critical 주식회사 바이오에이지
Priority to KR1020150030331A priority Critical patent/KR101669526B1/en
Priority to JP2017546614A priority patent/JP6640868B2/en
Priority to PCT/KR2016/000353 priority patent/WO2016140432A2/en
Priority to CN201680020966.9A priority patent/CN108475542A/en
Publication of KR20160107508A publication Critical patent/KR20160107508A/en
Application granted granted Critical
Publication of KR101669526B1 publication Critical patent/KR101669526B1/en
Priority to HK18115784.8A priority patent/HK1256728A1/en

Links

Images

Classifications

    • G06F19/00
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16ZINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00Subject matter not provided for in other main groups of this subclass
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H10/00ICT specially adapted for the handling or processing of patient-related medical or healthcare data
    • G06F19/32

Abstract

The present invention relates to a method for measuring a remaining life of a person by using a biological age and a survival probability and, more specifically, to a method for measuring a survival probability by using a biological age measured through a biological age measurement algorithm to which a health medical examination result is reflected and a method for measuring a remaining life by using the biological age and the survival probability measured by the method. According to the present invention, a remaining life can be predicted more accurately than a remaining life method which does not consider an existing biological age.

Description

생체나이를 이용한 잔여 수명 예측방법{Method for Assessinging Residual Life Using Biological Age}[0001] The present invention relates to a method for predicting residual lifetime using biological age,

본 발명은 생체나이와 생존확률을 이용하여, 사람의 잔여수명을 측정하는 방법에 관한 것으로, 더욱 자세하게는 건강검진 결과를 반영하는 생체나이 측정알고리즘을 통해 측정된 생체나이를 이용한 생존확률 예측방법 및 상기 방법으로 예측된 생존확률을 이용하여 잔여수명을 예측하는 방법에 관한 것이다.
More particularly, the present invention relates to a method for predicting survival probability using a biological age measured by a biological age measurement algorithm that reflects a result of a health examination, and more particularly, And a method for estimating the remaining service life using the predicted survival probability.

최근 사람의 질병 원인에 대한 많은 연구가 이루어지고 있음에도 불구하고 대부분이 흥미위주의 일회성 정보로 취급되고 있어서, 질병예방에는 별다른 효용이 없는 것이 사실이며, 이에 따른 대표적인 종래의 기술은 각 의료분야별로 질병을 예측하고 치료하는 기술이 있는 바, 특허출원번호 제2003-0067652호인 간암조기진단을 위한 간암 예측시스템 및 그 제어방법에 관한 기술로써, 개인별로 간암발생률의 추정과 간암발생의 상대적 위험(Relative risk)을 구하고, 각 개인별로 구한 간암발생률과 간암발생의 상대적 위험을 통해 간암발생위험군에 관한 계층화의 분류(Classify)를 수행하는 간암조기진단을 위한 간암 예측시스템 및 그 제어방법에 관한 기술이고, 상기 기술은 환자의 일반적인 정보, 초음파 시행에 따른 정보, 환자 최초등록시 소견과, 간암진단 시 소견에 관한 정보를 포함하는 임상정보 및 위험군 정보를 데이터베이스에 저장하고, 저장된 임상정보와 위험군 정보를 기반으로 각 위험요인에 해당하는 기여도인 회귀계수를 산출하여 간암발생 위험도를 측정함으로써, 간암발생 예측에 따른 개인별 간암발생을 예방할 수 있으며, 각 개인별로 구한 간암발생률과 간암발생의 상대적 위험을 통해 간암발생위험군에 관한 계층화의 분류를 수행함으로써, 맞춤형 간암발생예측모형의 기반을 구축하는 등의 효과가 있으며, 환자의 주치의는 간암 예측 결과를 이동통신단말을 통해 단문메시지 형태로 수신하거나 전자메일을 통해 수신함으로써 환자의 위험률을 지속적으로 모니터링할 수 있으며, 위험상황 발생시 즉각적으로 조치할 수 있는 효과가 있는 기술에 관한 것이 있었다.Recently, despite the fact that many researches on the causes of human diseases have been made, most of them are treated as one-time information focused on interest, and it is true that there is no utility for disease prevention. The technique of predicting and treating hepatocellular carcinoma is disclosed in Korean Patent Application No. 2003-0067652, which is related to the liver cancer prediction system for early diagnosis of liver cancer and the control method thereof. The present invention relates to a system for predicting liver cancer for early diagnosis of liver cancer and a control method thereof for classifying the stratification of liver cancer risk group through the relative risk of liver cancer incidence and the risk of occurrence of liver cancer, The technology includes general information of the patient, information of the ultrasound implementation, findings at the time of initial registration of the patient, By storing the clinical information and the risk group information including the information on the time of day observation in the database and calculating the regression coefficient which is the contribution factor corresponding to each risk factor based on the stored clinical information and the risk group information, It is possible to prevent the occurrence of individual hepatocellular carcinoma according to the prediction and to establish the basis of the customized hepatocarcinogenesis prediction model by classifying the risk of hepatocarcinoma risk group by the relative risk of hepatocellular carcinoma incidence and hepatocellular carcinoma incidence And the patient's primary care physician can continuously monitor the risk rate of the patient by receiving the liver cancer prediction result through a mobile communication terminal in the form of a short message or receiving it through the e-mail, and can immediately take measures in case of a dangerous situation There was something about technology.

상기와 같이 종래의 예측기술은 의학적으로 특정질병에 대하여 예측하고, 예측된 정보를 기반으로 치료하는 기술들이 대부분으로, 현재 건강상태에 따라 인간의 수명이 어떻게 변화되는지를 측정할 수는 없었다. As described above, the conventional prediction techniques are predominantly technologies for medically predicting a specific disease and treating it based on predicted information, and it has not been possible to measure how the lifespan of a human is changed according to current health conditions.

이에, 본 발명자들은 검강검진에서 얻어지는 데이터를 통하여 측정한 생체나이와 한국인의 수명을 조사한 대규모 데이터베이스를 이용하여 잔여수명 예측을 수행하는 경우, 신뢰성 있게 잔여수명을 예측할 수 있다는 것을 확인하고, 본 발명을 완성하게 되었다.
The present inventors confirmed that the remaining life can be reliably predicted when the remaining life prediction is performed using a large-scale database that examines the age of the living body and the life span of Koreans measured through the data obtained from the examination. It was completed.

본 발명의 목적은 건강검진에서 얻어지는 데이터를 기반으로 측정된 생체나이를 이용한 생존확률의 측정방법을 제공하는데 있다. It is an object of the present invention to provide a method of measuring a survival probability using a biological age measured based on data obtained from a health examination.

본 발명의 또 다른 목적은 상기 방법으로 측정된 생체나이와 생존확률을 이용하여, 잔여수명을 측정하는 방법을 제공하는데 있다.
It is still another object of the present invention to provide a method of measuring the remaining life using the biological age and the survival probability measured by the above method.

상기 목적을 달성하기 위하여, 본 발명은 검강검진 결과를 이용하여 측정된 생체나이와 하기 식(1)를 이용한 생존확률의 예측방법을 제공한다:
In order to achieve the above object, the present invention provides a method for predicting the survival probability using the biological age measured using the result of the examination of urination and the following equation (1)

Figure pat00001
.....(1)
Figure pat00001
.....(One)

상기 식 (1)에서, S0(t)는 통계청 자료 상의 현재 생존해 있는 시점에서부터 t년 후의 생존확률이고, D는 생체나이-출생나이이며, 는 D의 평균, a는 D가 생존확률에 미치는 영향력, (t)는 생체나이를 이용하여 보정계산된 t년 후까지 생존확률을 나타낸다.(1) where S 0 (t) is the survival probability after t years from the presently alive date on the statistical data, D is the age of the organism - birth date, D is the average of D, a is the survival probability (T) represents the probability of survival until t years after calibration using biometric age.

본 발명은 또한, x 축을 기대여명(ER), y 축을 T년 후까지 생존확률로 하여, 작성한 그래프에서, y축의 피검자의 생존확률 값에 대응하는 x 축의 값을 피검자의 검강검진을 받은 날로부터 남은 잔여수명으로 하는 것을 특징으로 하는 잔여수명 예측방법을 제공한다.
The present invention also relates to a method for estimating the survival probability of a subject on the y-axis from the date of receipt of the examination of the subject's chewing gum by using the x-axis as the survival survival rate (ER) and the y- And the residual service life remaining is set as the remaining service life.

본 발명에 따르면, 기존의 생체나이를 고려하지 않은 잔여수명 예측 방법에 비해 훨씬 정확하게 잔여수명을 예측할 수 있다.
According to the present invention, it is possible to predict the remaining service life much more accurately than the existing service life prediction method that does not consider the age of the living body.

도 1은 본 발명에 따른 잔여수명 예측방법을 도식화하여 나타낸 것이다.
도 2는 2013년 통계청 자료에 따른 출생나이 기준 T년 후까지 평균 생존확률을 나타낸 것이다.
도 3은 2013년 통계청 자료에 따른 평균 한국인의 연령별 기대 여명까지 생존확률 그래프이다.
도 4는 본 발명에 따른 잔여수명 측정방법으로 측정된 50세 남성(D=생체나이-출생나이=0)의 생존확률그래프이다.
도 5는 본 발명에 따른 잔여수명 측정방법으로 측정된 50세 남성(D=생체나이-출생나이=5)의 생존확률그래프이다.
도 6은 본 발명에 따른 잔여수명 측정방법으로 측정된 50세 남성(D=생체나이-출생나이=-5)의 생존확률그래프이다.
도 7은 10,000명의 한국인 남자를 대상으로 본 발명에 따라 추정된 기대 수명 분포를 나타낸 그래프이다.
도 8은 10,000명의 한국인 여자를 대상으로 본 발명에 따라 추정된 기대 수명 분포를 나타낸 그래프이다.
1 is a schematic representation of a residual life predicting method according to the present invention.
Fig. 2 shows the average survival probability after T year based on the 2013 National Statistical Office data.
FIG. 3 is a graph of the survival probability by an average Korean population's life expectancy by the 2013 National Statistical Office.
FIG. 4 is a graph of survival probability of a 50 year-old male (D = living organism age-birth age = 0) measured by the residual life measuring method according to the present invention.
FIG. 5 is a graph of survival probability of a 50 year-old male (D = living organism age-birth age = 5) measured by the residual life measuring method according to the present invention.
FIG. 6 is a graph of survival probability of a 50-year-old male (D = living organism age-birth age = -5) measured by the residual life measurement method according to the present invention.
7 is a graph showing an estimated life expectancy distribution according to the present invention for 10,000 Korean men.
8 is a graph showing an estimated life expectancy distribution according to the present invention for 10,000 Korean women.

본 발명에서 "출생나이(CA)"는 태어난 날을 기준으로 해서 건강 검진을 받은 날 까지 생존한 년수를 말하며, "생체나이(BA)"는 건강 검진 항목을 입력 값으로 하여 통계 알고리즘을 적용 후 계산한 개인의 건강 상태를 반영한 나이를 말한다.In the present invention, the term " birth age (CA) "refers to the number of years of survival until the day of receiving the health examination based on the day of birth, and" It refers to the age that reflects the state of health of the person who has been calculated.

본 발명에서, "T년 생존확률(SP)"은 건강 검진을 받은 날을 기준으로 하여 향후 년 까지 생존할 확률을 말하고, "기대여명 (ER)"은 건강 검진을 받은 날을 기준으로 하여 향후 생존할 수 있을 년 수를 말하며, "기대수명 (EL)"은 건강 검진을 받았을 때의 출생 나이(CA)+ 기대여명(ER)을 말한다.In the present invention, the term "T year survival probability (SP)" refers to the probability of surviving until the next year based on the day of the health examination, and "Life expectancy (ER) (EL) refers to the age at birth (CA) + life expectancy (ER) at the time of a health checkup.

본 발명에서, "기대 여명까지 생존 확률 그래프 (SP2ER)"는 x축을 기대 여명(ER), y축을 T년 후의 생존 확률(SP)로 설정한 그래프(성별, 연령별로 구분됨, 통계청 발표자료)를 말한다. In the present invention, the graph "Survival probability graph up to life expectancy (SP2ER)" is a graph in which the x axis is the life expectancy (ER) and the y axis is the survival probability (SP) after T year (classified by sex and age, It says.

본 발명에서, "생체나이를 기반으로 한 T년 후의 생존확률 (SPBA)"은 건강 검진을 받은 날 계산된 출생나이(CA)와 생체나이(BA)와 T년 생존확률(SP)을 이용하여 계산된 생존 확률생체 생존확률을 말한다. In the present invention, the "survival probability (TBA) after T year based on the age of the living body" is calculated by using the calculated birth age (CA), the biological age (BA) and the T year survival probability Calculated survival probability refers to the bio-survival probability.

본 발명에서, "생체나이(BA)를 기반으로 한 기대 여명까지 생존 확률 그래프 (SP2ERBA)"는 x축을 기대여명(ER), y축을 생체나이(BA)를 기반으로 한 T년 생존확률(SPBA)로 설정한 그래프(성별, 연령별로구별됨)를 말한다.In the present invention, the survival probability graph (SP2ERBA) up to the life expectancy based on the biological age (BA) is calculated by multiplying the x axis by the expected life expectancy (ER), the y axis by the T year survival probability ), Which is classified by sex and age.

본 발명에서, "잔여수명(RL)"은 기대 여명까지 생존 확률 그래프(SP2ER)과 생체나이(BA)를 기반으로 한 기대 여명까지 생존 확률 그래프(SP2ERBA) 를 이용하여 계산된 건강검진을 받은 날로부터 남은 생존년수를 말한다.In the present invention, the term "remaining life (RL)" is the day when the health checkup calculated using the survival probability graph (SP2ERBA) to the life expectancy based on the survival probability graph (SP2ER) The number of surviving years.

일관점에서, 본 발명은 검강검진 결과를 이용하여 측정된 생체나이와 하기 식(1)를 이용한 생존확률의 측정방법에 관한 것이다:In one aspect, the present invention relates to a method of measuring the survival probability using the bio-age measured using the result of the examination of urination and the following formula (1)

Figure pat00002
.....(1)
Figure pat00002
.....(One)

상기 식 (1)에서, S0(t)는 통계청 자료 상의 현재 생존해 있는 시점에서부터 t년 후의 생존확률이고, D는 생체나이-출생나이이며, 는 D의 평균, a는 D가 생존확률에 미치는 영향력, (t)는 생체나이를 이용하여 보정계산된 t년 후까지 생존확률을 나타낸다.(1) where S 0 (t) is the survival probability after t years from the presently alive date on the statistical data, D is the age of the organism - birth date, D is the average of D, a is the survival probability (T) represents the probability of survival until t years after calibration using biometric age.

상기 식(1) 에서 S0(t) 는 통계청에서 각 연도마다 남녀별로 제공하는 생존확률 값으로, 통계청 홈페이지를 통해 확인할 수 있다(http://kosis.kr/statHtml/statHtml.do?orgId=101&tblId=DT_1B42&vw_cd=MT_ZTITLE&list_ id=A5&seqNo=&lang_mode=ko&language=kor&obj_var_id=&itm_id=&conn_path=E1#). In the above equation (1), S 0 (t) is the survival probability value provided for each male and female in each year in the National Statistical Office, and can be confirmed through the National Statistical Office's homepage (http://kosis.kr/statHtml/statHtml.do?orgId= 101 & tblId = DT_1B42 & vw_cd = MT_ZTITLE & list_id = A5 & seqNo = & lang_mode = en & language = kor & obj_var_id = & itm_id = & conn_path = E1 #).

식(1)에서 a는 생존자료를 입력값으로 이용하여 콕스 비례 위험 모형(cox proportional hazards model) 알고리즘을 이용하여 자동으로 산출되어 나오는 값으로, SPSS 나 SAS, R package 등과 같은 통상적인 통계 분석을 위한 상용 또는 무료 사 분석 소프트웨어를 실행시켜 구현할 수 있다. 여기서 생존 자료라 함은 건강검진을 받은 피검자 개개인을 대상으로 하여 이후 특정한 시점까지 관측했을 때 사망 또는 생존했는지에 대한 정보가 있는 자료를 의미한다.
In the equation (1), a is the value automatically calculated using the cox proportional hazards model algorithm using the survival data as the input value. It is a statistical analysis such as SPSS, SAS, and R package And can be implemented by running commercial or free analysis software. Here, the term "survival data" refers to data on individuals who have undergone medical examination and have information on whether they have died or survived to a certain point in time.

바람직한 일 양태로서, 본 발명에 따른 식 (1)에서의 a 값 및

Figure pat00003
값은 연령에 따라 표 1과 같은 값을 가질 수 있지만, 이에 한정되는 것은 아니다.
In a preferred embodiment, the value of a in formula (1) according to the present invention and
Figure pat00003
The value may have a value as shown in Table 1 according to age, but is not limited thereto.

Figure pat00004
Figure pat00004

본 발명에 있어서, 상기 생체나이는 하기 식 (2)~(3)에 의하여 측정되는 것을 특징으로 할 수 있다:
In the present invention, the biological age may be measured by the following formulas (2) to (3):

Figure pat00005
...(2)
Figure pat00005
...(2)

Figure pat00006
Figure pat00006

....(3).... (3)

여기서, preBA는 식 (3)으로 보정하기 전의 가생체나이이고, BA는 가생체나이(preBA)를 기반으로 보정한 후의 생체나이이며, xj는 피검진자의 j번째 검진 항목의 검진 수치,

Figure pat00007
는 j번째 검진 항목의 샘플 내 평균값, sd(xj)는 j번째 검진 항목의 샘플 내 표준 편차, βij는 주성분 분석 결과 j번째 검진 항목의 i번째 요인값이 각 샘플의 실제 나이를 설명할 수 있는 정도, m은 주성분 분석의 수행결과 얻어진 요인 개수, P i는 i번째 요인 값이 preBA 계산에 주는 가중치, y는 각 개체의 출생나이(CA)이며,
Figure pat00008
는 생체나이(BA) 모델 구축에 사용되는 모든 개체의 출생나이(CA)의 평균, sd(y)는 y의 표준 편차이고,
Figure pat00009
는 각 샘플에 대하여 실제 나이 및 preBA를 회귀 분석한 결과 생성된 회귀식
Figure pat00010
를 구성하는 계수이다.Here, preBA is the age of the subject before correction by equation (3), BA is the age of the subject after correction based on the age of the subject (preBA), xj is the examination value of the subject,
Figure pat00007
Is the mean value in the sample of the jth item, sd (x j ) is the standard deviation of the sample in the jth item, and β ij is the i-th factor of the item in the j th item of the principal component analysis. ( P ) is the weight given to the preBA calculation, and y is the birth age (CA) of each individual,
Figure pat00008
Is the mean of the birth years (CA), sd (y) is the standard deviation of y for all individuals used in constructing the BA model,
Figure pat00009
Is the regression of the actual age and preBA for each sample,
Figure pat00010
.

본 발명에 있어서, 상기 Bij는 검강검진을 통해 측정된 2개 이상의 변수 (Xj)를 SPSS, SAS 및 R pack age 프로그램을 이용하여 주성분분석 통계알고리즘에 입력하여 산출되는 값인 것을 특징으로 할 수 있으며, m은 주성분 요인수를 나타내는 데, 주성분 분석 알고리즘 실행 시 고유치(eigen value)를 1로 설정했을 때 자동으로 결정되어 나오는 값이다.In the present invention, B ij may be a value calculated by inputting two or more variables (X j) measured through a screening examination into a principal component analysis statistical algorithm using SPSS, SAS and R package programs , m represents the number of principal component factors, which is automatically determined when the eigenvalue is set to 1 when the principal component analysis algorithm is executed.

본 발명의 일 양태에서, 생체나이 측정은 대한민국 특허공개 2014-0126229에 기재된 방법을 사용할 수 있다.
In one embodiment of the present invention, the method described in Korean Patent Publication No. 2014-0126229 can be used for measurement of living body age.

본 발명에 있어서, 피검자가 남자일 경우, x는 허리둘레(WC), 1초간 노력성 기량(FEV1), 간효소(G-GTP), 요소질소(BUN), 고지단백(HDL), 저지단백(LDL), 중성지방(TG), 공복혈당(FBS), 체지방률(FBR), 근육률(BMR), 알부민-글로불린 비율(AGR) 및 수축기혈압(SBP)으로 구성된 군에서 선택되는 건강검진 항목의 결과수치이거나, 허리둘레(WC), 수축기혈압(SBP), 간효소(G-GTP), 고지단백(HDL), 중성지방(TG), 공복혈당(FBS), 헤모글로빈(Hemoglobin), 저지단백(LDL), 체질량지수(BMI) 및 맥압(PP)으로 구성된 군에서 선택되는 건강검진 항목의 결과수치인 것을 특징으로 할 수 있다. In the present invention, when the subject is a male, x is expressed by the following equation: WC, FEV 1, G-GTP, urea nitrogen, BN, (S) selected from the group consisting of LDL, triglyceride (TG), fasting blood glucose (FBS), body fat percentage (FBR), muscle rate (BMR), albumin-globulin ratio (AGR) and systolic blood pressure (HDL), triglyceride (TG), fasting blood sugar (FBS), hemoglobin, and low-density lipoprotein (HDL) LDL), body mass index (BMI), and pulse pressure (PP).

본 발명에 있어서, 피검자가 여자일 경우, x는 허리둘레(WC), 1초간 노력성 기량(FEV1), 감마-GTP)G-GTP), 요소질소(BUN), 고지단백(HDL), 저지단백(LDL), 중성지방(TG), 공복혈당(FBS), 체지방률(FBR), 근육률(BMR), 알부민-글로불린 비율(AGR), 이완기혈압(DBP), 적혈구침강속도(ESR) 및 체질량지소(BMI)로 구성된 군에서 선택되는 건강검진 항목의 결과수치이거나, 허리둘레(WC), 이완기혈압(DBP), 간효소(AST), 고지단백(HDL), 중성지방(TG), 공복혈당(FBS), 헤모글로빈(Hemoglobin), 저지단백(LDL), 체질량지수(BMI) 및 맥압(PP)으로 구성된 군에서 선택되는 건강검진 항목의 결과수치인 것을 특징으로 할 수 있다. In the present invention, when the subject is an exciter, x is expressed by the following equation: WC (waist circumference), FEV1 for 1 second, GTP-GTP), urea nitrogen (BUN) (LDL), triglyceride (TG), fasting blood sugar (FBS), body fat percentage (FBR), muscle rate (BMR), albumin-globulin ratio (AGR), diastolic blood pressure (DBP), erythrocyte sedimentation rate (DBP), liver enzymes (AST), high fat protein (HDL), triglyceride (TG), fasting blood glucose (BMI) And the result of the health screening item selected from the group consisting of FBS, hemoglobin, LDL, body mass index (BMI), and pulse pressure (PP).

본 발명에서 사용되는 생체나이 예측 모델에는 기존의 검진 데이터인 하나 이상의 샘플이 입력될 수 있다. 상기 생체나이 연산모델에는 생체나이 연산 모델에 포함된 변수(검진 항목)의 검진 수치를 포함하는 샘플인지 여부와 무관하게, 모든 샘플이 입력될 수도 있으나, 생체나이 연산 모델에 포함된 변수의 검진 수치를 모두 포함하는 샘플들만 입력되는 것이 바람직하다. 예를 들어, 상기 제1 생체나이 연산 모델에는 BMI 검진 수치, 폐활량 검진 수치 및 일초량 검진 수치를 모두 가지는 샘플들만 입력되고, 상기 제2 생체나이 연산 모델에는 BMI 검진 수치, 폐활량 검진 수치 및 PSA 검진 수치를 모두 가지는 샘플들만 입력될 수 있다.In the biological age prediction model used in the present invention, one or more samples that are conventional examination data can be input. In the biological age calculation model, all the samples may be input regardless of whether or not the samples include the examination values of the variables (examination items) included in the biological age calculation model, but the examination values of the variables included in the biological age calculation model It is preferable that only the samples including both of them are inputted. For example, only the samples having both the BMI examination value, the spirometric examination value, and the one-minute examination value are input to the first living body age calculation model, and the second living body age calculation model includes the BMI examination value, the spirometry examination value, Only samples with all numerical values can be entered.

상기 생체나이 연산 모델에 하나 이상의 샘플이 입력된 후에는, 입력된 데이터에 대한 주성분 분석(Principal Component Analysis)이 수행될 수 있다. 상기 주성분 분석의 결과, 상기 생체나이 연산 모델에 입력된 데이터를 설명할 수 있는 하나 이상의 요인(factor)을 얻을 수 있다. 상기 얻어진 요인의 개수는 상기 생체나이 연산 모델의 변수 개수 보다 적을 수 있다. 이하, 상기 요인의 개수는 m으로, 상기 생체나이 연산 모델의 변수 개수는 n으로 기재한다.  After one or more samples are input to the biological age computation model, Principal Component Analysis on the input data may be performed. As a result of the principal component analysis, one or more factors that can explain the data input to the biological age calculation model can be obtained. The number of the obtained factors may be smaller than the number of variables of the biological age calculation model. Hereinafter, the number of the factors is m, and the number of variables of the biological age calculation model is described as n.

본 발명의 일 실시예에 따르면, 상기 주성분 분석을 통하여, 상기 생체나이 연산 모델에 입력된 각 샘플의 j번째 변수값들에 포함되어 있는 i번째 요인 값들이 실제 나이를 설명할 수 있는 정도를 의미하는 가중치(Bij)가 연산될 수 있다. Bij값은 일반적인 통계분석에서 알려진 주성분 분석(principal component analysis) 방법을 통해 산출되는 값이다. 즉, Bij는 주성분 분석 프로세스를 수행하면 자동으로 계산 되는 각 변수별 요인 적재행렬(factor loading matrix)이다.According to an embodiment of the present invention, through the principal component analysis, the i-th factor values included in the j-th variable values of each sample input to the biological age calculation model indicate the degree to which the actual age can be explained (B ij ) can be calculated. The value of B ij is the value calculated by the principal component analysis method known from general statistical analysis. That is, B ij is a factor loading matrix for each variable that is automatically calculated when the principal component analysis process is performed.

본 발명의 일 실시예에 따르면, 상기 주성분 분석을 통하여, i번째 요인 값이 생체나이 계산에 주는 가중치(P i )도 연산될 수 있다. 이하, 가중치(Pi)를 연산하는 방법을 설명한다.According to an embodiment of the present invention, the weight ( P i ) of the i-th factor value and the bio-age calculation may be calculated through the principal component analysis. Hereinafter, a method of calculating the weight P i will be described.

가중치(P i)는 주성분 분석 결과 계산된 각 샘플의 요인(factor)별 요인점수(factor score)와 실제 나이 사이의 결정 계수(coefficient of determination, R2) 를 이용하여 연산될 수 있다. 상기 결정 계수를 이용한 가중치(Pi )의 연산 방법은 아래와 같다.The weight ( P i ) can be calculated using a coefficient of determination (R 2 ) between the factor score and the actual age of each sample factor calculated from the principal component analysis. A calculation method of the weight value (P i ) using the decision coefficient is as follows.

먼저, 주성분 분석을 통해 분석에 사용된 각 샘플의 요인 별 요인점수를 연산한다. 예를 들어, 샘플이 100명이고 요인이 m개 라고 한다면 m개 요인 별로 각각 100개의 요인 점수가 계산될 것이다.First, we calculate factor scores for each sample used in the analysis by principal component analysis. For example, if you have 100 samples and m factors, then 100 factor scores will be calculated for each of m factors.

다음으로, 각 요인별 요인점수와 각 샘플의 실제 나이를 회귀 분석 하여 결정 계수를 m개(

Figure pat00011
) 계산 한다. 이 m개의 결정 계수의 합을 S 라고 하면, i번째 요인이 생체나이 계산에 주는 가중치(Pi )는 아래의 수학식 1을 이용하여 연산될 수 있다.
Next, the factor scores of each factor and the actual age of each sample are regression analyzed to find the coefficient of determination m (
Figure pat00011
). Assuming that the sum of the m decision coefficients is S, the weight P i of the i-th factor to the biometric age calculation can be calculated using the following equation (1).

[수학식 1][Equation 1]

Figure pat00012

Figure pat00012

가중치(Pi)는 주성분 분석 결과 연산된 각 요인 별 고유치(eigen value)를 이용하여 연산될 수도 있다. 이 때, 가중치(Pi)는 아래의 수학식 2를 이용하여 연산될 수 있다.
The weight P i may be calculated using the eigenvalue of each factor calculated from the principal component analysis. At this time, the weight P i can be calculated using the following equation (2).

[수학식 2]&Quot; (2) "

Figure pat00013

Figure pat00013

상기 수학식 2에서 ei는 요인 i의 고유치(eigen value)이고, m은 요인의 개수이다.In Equation (2), e i is the eigenvalue of the factor i, and m is the number of factors.

상기 주성분 분석이 완료되면, 상기 요인을 이용하여 상기 생체나이 연산 모델에 따른 생체나이 산출식이 생성될 수 있다.When the principal component analysis is completed, a bio-age calculation expression according to the bio-age calculation model can be generated using the factors.

일 실시예에 따르면, 상기 생체나이는 식(2)에 기반하여 수학식 3에 따라 산출될 수 있다.
According to one embodiment, the biological age can be calculated according to equation (3) based on equation (2).

Figure pat00014
Figure pat00014

...식(2)
... (2)

(여기서, preBA는 가생체나이, xj는 피검진자의 j번째 검진 항목의 검진 수치,

Figure pat00015
는 j번째 검진 항목의 샘플 내 평균값, sd(xj)는 j번째 검진 항목의 샘 플 내 표준 편차, Bij는 주성분 분석 결과 j번째 검진 항목의 i번째 요인값이 각 샘플의 실제 나이를 설명할 수 있는 정도, m은 주성분 분석의 수행결과 얻어진 요인 개수, Pi는 i번째 요인 값이 BA 계산에 주는 가중치)
(Where preBA is the age of the living organism, x j is the examination value of the j th examination item of the examinee,
Figure pat00015
Is the mean value in the sample of the jth item, sd (x j ) is the sample standard deviation of the jth item, and B ij is the actual age of each sample The number of factors obtained as a result of the principal component analysis, and P i is the weight given to BA calculation of the i-th factor value)

수학식 3Equation 3

Figure pat00016
Figure pat00016

(상기 수학식 3에서 BA는 생체나이, preBA 및 preBA'는 가생체나이, xj는 피검진자의 j번째 검진 항목의 검진 수치,

Figure pat00017
는 j번째 검진 항목의 샘플 내 평균값, sd(xj)는 j번째 검진 항목의 샘플 내 표준 편차, Bij는 주성분 분석 결과 j번째 검진 항목의 i번째 요인값이 각 샘플의 실제 나이를 설명할 수 있는 정도, m은 주성분 분석의 수행결과 얻어진 요인 개수, Pi는 i번째 요인 값이 preBA 계산에 주는 가중치, y는 각 개체의 출생나이(CA)이며,
Figure pat00018
는 생체나이(BA) 모델 구축에 사용되는 모든 개체의 출생나이(CA)의 평균, sd(y)는 y의 표준 편차,
Figure pat00019
Figure pat00020
는 각 샘플에 대하여 실제 나이 및 preBA를 회귀 분석한 결과 생성된 회귀식
Figure pat00021
를 구성하는 계수)In the above equation (3), BA is the age of the living organism, preBA and preBA 'are the age of the living organism, x j is the examination value of the j th examination item of the examinee,
Figure pat00017
Is the average value in the sample of the jth item, sd (x j ) is the standard deviation of the sample in the jth item, and B ij is the i-th factor of the item in the j th item of the principal component analysis. (P ) is the weight given to the preBA calculation, and y is the birth age (CA) of each individual,
Figure pat00018
Is the mean of the birth years (CA) of all individuals used in constructing the BA model, sd (y) is the standard deviation of y,
Figure pat00019
And
Figure pat00020
Is the regression of the actual age and preBA for each sample,
Figure pat00021
≪ / RTI >

상기 수학식 3에서, 일 실시예에 따르면, 상기 y평균은 모든 샘플의 실제 나이에 대한 평균이고, 상기 sd(y)는 모든 샘플의 실제 나이에 대한 표준 편차일 수 있으나, 다른 실시예에 따르면, 상기 y평균은 상기 생체나이 연산 모델에 포함된 모든 검진 항목의 수치 데이터를 유효한 값으로 가지고 있는 모든 샘플의 실제 나이에 대한 평균이고, 상기 sd(y)는 상기 생체나이 연산 모델에 포함된 모든 검진 항목의 수치 데이터를 유효한 값으로 가지고 있는 모든 샘플의 실제 나이에 대한 표준 편차일 수도 있다.In Equation (3), according to one embodiment, the y average is an average of the actual age of all samples, and sd (y) may be a standard deviation of the actual age of all samples. (Y) is an average of all of the samples having valid values of all the items of the examination items included in the biological age calculation model, and sd (y) is an average of all the samples included in the biological age calculation model It may also be the standard deviation of the actual age of all samples with valid values for the numerical data of the examination item.

정리하면, 상기 생체나이 연산 모델에 샘플 데이터를 입력하고, 입력된 데이터에 대한 주성분 분석을 거치면 하나 이상의 요인이 얻어지고, 상기 주성분 분석 결과 및 상기 요인을 이용하여 생체나이 산출식을 얻을 수 있다. 상기 생체나이 산출식의 예시는 식 (2) 또는 수학식 3을 참조할 수 있다.In summary, one or more factors are obtained by inputting the sample data into the biological age calculation model and subjecting the input data to principal component analysis, and the biological age calculating equation can be obtained using the principal component analysis result and the factors. An example of the bio-age calculating equation may refer to equation (2) or equation (3).

상기 생체나이 연산 모델에, 이번의 피검진자 검진 결과 데이터에 포함된 검진 항목을 변수로서 하나씩 추가해나가면서, 최적의 결과를 도출하는 생체나이 연산 모델을 구성하는 변수 집합을 찾는 과정으로 이해할 수 있다. 즉, 처음에는 비어 있는 상기 생체나이 연산 모델에 상기 원본 집합의 검진 항목 중 최선의 하나를 선정하여, 선정된 검진 항목을 변수로서 추가하는 것을 반복하는, 소위 단계별 전진 선택 방법을 이용한다. It can be understood as a process of finding a set of variables constituting a biological age calculation model for deriving an optimum result while adding the examination items included in the result data of the examinee's examination results to the living body age calculation model one by one as a variable. That is, a so-called step-by-step forward selection method is used in which the best one of the items of the original set of the original set is selected for the above-described living body age calculation model that is empty and the selected items are added as variables.

이하, 상기 단계별 전진 선택 방법에 대하여 설명한다. 원본 집합에 BMI, 폐활량, 일초량, 알부민(Albumin), PSA의 총 5가지 검진 항목이 원소로서 포함된 것을 전제로 한다. Hereinafter, the step-by-step forward selection method will be described. It is assumed that the original set includes five items of examination, BMI, lung capacity, daily dose, albumin, and PSA.

첫번째 단계로, 하나의 변수를 포함하는 테스트 모델이 생성된다. 다음으로, 생성된 테스트 모델에 해당 변수의 값, 즉 해당 검진 항목의 검진 [0075] 수치를 포함하는 샘플들이 입력되며, 입력된 샘플들에 대하여 주성분 분석을 수행하여 요인을 얻고, 얻어진 요인 및 주성분 분석 결과를 이용하여 생체나이 산출식을 생성한다. 이하, 상기 테스트 모델에 새로 추가되는 변수를 테스트 검진 항목으로 기재한다. 다음으로, 상기 입력된 샘플들 각각을 상기 생체나이 산출식에 입력하여 각 샘플의 생체나이를 연산한다.  As a first step, a test model containing a single variable is created. Next, samples including the value of the corresponding variable, that is, the examination value of the examination item, are input to the generated test model, and the principal components analysis is performed on the input samples to obtain the factors, Using the results of the analysis, a bioassay equation is generated. Hereinafter, variables newly added to the test model are described as test screening items. Next, each of the input samples is input to the bio-age calculating equation to calculate bio-age of each sample.

다음으로, 샘플들의 상기 분포에 대한 회귀 분석을 통하여 단순 회귀식 BA = a + by를 생성한다. 그 결과, 샘플들의 회귀선(

Figure pat00022
)과 실제 나이(BA = y) 와의 결정 계수(R2)을 구할 수 있다. 첫번째 단계의 결과 얻어진 각 검진항목 별 결정 계수 값이 표시되어 있다. 첫번째 단계의 결과, 가장 큰 값을 가지는 폐활량이 선택될 것이다. 이는, 폐활량 검진 수치는 생체나이 연산 모델에 포함되는 것을 의미한다. 폐활량 검진 항목은 상기 생체나이 연산 모델에 추가되었으므로, 상기 원본 집합에서 제거된다. 두번째로 상기 생체나이 연산 모델에 포함될 두번째 검진항목이 선택된다. 두번째 단계의 테스트 모델은 상기 폐활량 검진 항목 및 그 밖의 검진 항목 중 하나를 포함하는, 총 2개의 변수를 가지는 모델이다. 따라서, 두번째 단계에서는 (폐활량, BMI), (폐활량, 일초량), (폐활량, 알부민), (폐활량, PSA) 총 4번의 테스트가 진행될 것이다.Next, a simple regression equation BA = a + by is created through regression analysis of the distribution of the samples. As a result, the regression line of the samples
Figure pat00022
) And it can be obtained for the coefficient of determination (R 2) between the actual age (BA = y). The result of the first step shows the determination coefficient values for each examination item. As a result of the first step, the lung capacity with the largest value will be selected. This means that the spirometric examination value is included in the biological age calculation model. The spirometry check item is removed from the original set because it has been added to the biological age calculation model. Second, a second examination item to be included in the biological age calculation model is selected. The second-stage test model is a model having a total of two variables, including one of the above-mentioned spirometry examination items and other examination items. Therefore, in the second stage (spirometry, BMI), (lung capacity, one second), (lung capacity, albumin), (lung capacity, PSA)

다음으로, 생성된 생체나이 연산 모델에 해당 변수의 값, 즉 해당 검진 항목의 검진 수치를 포함하는 샘플들이 입력되며, 입력된 샘플들에 대하여 주성분 분석을 수행하여 요인을 얻고, 얻어진 요인 및 주성분 분석 결과를 이용하여 생체나이 산출식을 생성한다.Next, samples including the value of the corresponding variable, that is, the examination value of the corresponding examination item, are input to the generated biological age calculation model, and the principal components analysis is performed on the input samples to obtain the factors, The result is used to generate a bio-age expression.

다음으로, 상기 입력된 샘플들 각각을 상기 생체나이 산출식에 입력하여 각 샘플의 생체나이를 연산하고, 샘플들의 상기 분포에 대한 회귀 분석을 통하여 단순 회귀식

Figure pat00023
를 생성한다. 그 결과, 샘플들의 회귀선 (
Figure pat00024
)과 실제 나이(BA = y) 와의 결정 계수(R2)를 구할 수 있다.Next, each of the input samples is input to the bio-age calculating equation to calculate bio-age of each sample, and a simple regression equation
Figure pat00023
. As a result, the regression line of the samples
Figure pat00024
) And it can be calculated the coefficient of determination (R 2) between the actual age (BA = y).

두번째 단계의 결과, 기존의 결정 계수인 0.4보다 모델의 결정 계수 값을 증가 시키는 것 중, 증가량이 가장 큰 검진 항목인, BMI가 선정될 것이다. 이는, BMI 수치가 폐활량 수치와 함께 생체나이 연산 모델에 포함되는 것을 의미한다. BMI 검진항목은 상기 생체나이 연산 모델에 추가되었으므로, 상기 원본 집합에서 제거된다. 세번째로 상기 생체나이 연산 모델에 포함될 세번째 검진항목이 선택된다. 세번째 단계의 테스트 모델은 상기 폐활량 검진 항목, BMI 및 그 밖의 검진 항목 중 하나를 포함하는, 총 3개의 변수를 가지는 모델이다. 따라서, 두번째 단계에서는 (폐활량, BMI, 일초량), (폐활량, BMI, 알부민), (폐활량, BMI, PSA) 총 3번의 테스트가 진행될 것이다. 세번째 단계 역시 두번째 단계와 유사항 반식으로 진행된다. 일초량의 결정 계수 값이 기존의 결정 계수 값인 0.5에 비하여 유일하게 증가하는 값이므로, 일초량이 상기 생체나이 연산 모델에 포함되는 3번째 검진 항목으로 선정된다. 일초량은 상기 원본 집합에서 제거된다. As a result of the second step, BMI, which is the largest increase in the number of determinants of the model, will be selected. This implies that the BMI value is included in the biological age calculation model together with the spirometric value. The BMI check item is removed from the original set because it has been added to the biological age calculation model. Thirdly, a third examination item to be included in the biological age calculation model is selected. The third-stage test model is a model having a total of three variables, including one of the above-mentioned spirometric examination item, BMI, and other examination items. Therefore, in the second phase, three tests will be conducted (lung capacity, BMI, one second), (lung capacity, BMI, albumin), (lung capacity, BMI, PSA). The third step is also the second step and the benefit process. The one-second amount of determination coefficient is selected as the third examination item included in the biological age calculation model since the determination value of one-minute amount is uniquely increased compared to the conventional determination coefficient value of 0.5. One second is removed from the original set.

네번째로 상기 생체나이 연산 모델에 포함될 세번째 검진항목이 선택된다. Fourth, a third examination item to be included in the biological age calculation model is selected.

남은 두개의 검진 항목 모두, 기존의 결정 계수 값인 0.6 보다 낮은 결정 계수를 가지는 생체나이 연산을 만들고 있다. 네번째 단계까지 진행한 결과, 원본 집합에서 더 이상 추가될 검진 항목이 존재하지 않는 상황이다. 따라서, 네번째 단계에서 상기 단계별 전진 진행 단계는 종료된다. 그 결과, 이번 피검진자에 대한 생체나이 연산 모델에는 5개의 검진 항목 중, BMI, 폐활량, 일초량 총 3개만 변수로 포함되는 것이 바람직하다. 상기 단계별 전진 진행 단계가 종료되는 시점에서의 상기 생체나이 연산 모델이 상기 수신된 검진 데이터에 대하여 적용될 생체나이 연산 모델로 확정된다. Both of the remaining examinations are making biometric calculations with a coefficient of determination lower than 0.6, which is the conventional coefficient of determination. As a result of proceeding to the fourth step, there is no further examination item to be added in the original set. Therefore, in the fourth step, the stepwise advancing step is terminated. As a result, it is preferable that only three of the five examination items, BMI, lung capacity, and one second, are included in the biological age calculation model for the examinee. The biological age calculation model at the end of the stepwise progression step is determined as a biological age calculation model to be applied to the received examination data.

다음으로, 상기 수신된 검진 데이터를 상기 확정된 생체나이 연산 모델에 따른 생체나이 산출식에 입력하여 피검진자의 생체나이를 연산할 수 있다. Next, the received examination data can be inputted into a living body age calculating formula according to the determined living body age calculating model to calculate the living body age of the examinee.

일 실시예에 따르면, 상기 검진 데이터를 상기 샘플 데이터에 추가함으로써, 다음 피검진자를 위한 샘플로 활용할 수 있다.According to one embodiment, by adding the examination data to the sample data, it can be utilized as a sample for the next examinee.

본 발명에 따른 단계별 전진 선택 단계는 피검진자의 검진 데이터에 포함된 검진 항목 중에서 생체나이 연산 모델에 포함될 적어도 일부의 검진 항목을 하나씩 추가하는 과정이다. 데이터에 따라, 단계별 전진 선택 단계를 진행할수록 상기 결정 계수 값이 계속 증가하는 경우도 존재할 수 있으므로, 피검진자의 검진 데이터에 포함된 검진 항목 모두가 상기 생체나이 연산 모델에 포함되는 경우도 발생 할 수 있음은 물론이다.The step-by-step forward selection step according to the present invention is a step of adding at least a part of the examination items to be included in the biological age calculation model among the examination items included in the examination data of the examinee. There may be a case where the decision coefficient value continues to increase as the stepwise advancing selection step is progressed according to the data. Therefore, it may happen that all the examination items included in the examination data of the examinee are included in the biometric age calculation model Of course it is.

한편, 피검진자의 검진 데이터와 유사한 패턴을 보일 수 있는 샘플을 이용하는 것이 보다 정확한 생체나이 연산 모델을 생성하는데 유리할 것이다. 따라서, 상기 단계별 전진 선택 단계에서 상기 테스트 모델에 입력되는 샘플은 상기 수신된 검진 데이터의 피검진자와 유사한 신상 정보를 가지는 샘플으로 제한될 수 있다. 예를 들어, 상기 테스트 모델에 입력되는 샘플은 상기 피검진자의 성별와 동일하고, 상기 피검진자의 실제 나이를 기준으로 기 지정된 범위의 실제 나이를 가지는 것으로 제한될 수 있다. On the other hand, using a sample that can show a pattern similar to that of the examinee's test data would be advantageous for generating a more accurate biometric age calculation model. Therefore, in the step-by-step forward selection step, the sample input to the test model may be limited to the sample having the personal information similar to the examinee of the received examination data. For example, the sample input to the test model is the same as the sex of the examinee, and may be limited to having the actual age of the predetermined range based on the actual age of the examinee.

상기 단계별 전진 선택 단계에서, 상기 결정 계수를 기존의 테스트 모델의 그것에 비하여 최대한 증가시키는 지 여부뿐만 아니라, 추가적인 선택 기준이 사용될 수 있다.In the step-by-step forward selection step, an additional selection criterion can be used as well as whether the decision coefficient is increased as much as possible compared to that of the existing test model.

일 실시예에 따르면, 상기 테스트 검진 항목의 수치 데이터를 가지는 샘플의 수가 기지정된 한계치 이하인 경우 상기 테스트 검진 항목에서 제외 될 수 있다. 다른 실시예에 따르면, 상기 테스트 모델에 포함된 검진 항목의 검진 수치를 모두 가지는 샘플의 수가 기지정된 한계치 이하인 경우 상기 테스트 검진 항목에서 제외 될 수 있다. 샘플의 개수가 상기 한계치 미만인 경우, 통계 데이터로서의 신뢰도를 부여하기 어렵기 때문이다.According to an embodiment, when the number of samples having the numerical data of the test examination item is less than the predetermined threshold value, it may be excluded from the test examination item. According to another embodiment, if the number of samples having all the examination values of the examination items included in the test model is less than the predetermined threshold value, it can be excluded from the test examination item. When the number of samples is less than the threshold value, it is difficult to give reliability as statistical data.

또한, 일 실시예에 따르면, 상기 테스트 모델에 포함된 검진 항목 중 적어도 하나와 상관 계수(R)가 기 지정된 한계치를 초과하는 상기 테스트 검진 항목은 상기 테스트 검진 항목에서 제외 될 수 있다. 이는, 테스트 모델에 이미 포함된 검진 항목과 상관 관계가 높은 검진 항목이 새로이 테스트 모델에 포함되게 되면, 기존의 검진 항목이 테스트 모델에서 차지하는 영향력이 감소하는 역효과를 방지하기 위함이다. 또한, 테스트 모델에 상관 관계가 높은 검진 항목들이 포함되는 경우, 주성분 분석 과정에서 수행하는 행렬 연산 과정 중 singularity 현상이 발생되어 고유치(eigen value)가 제대로 연산되지 않는 문제를 방지하기 위함이기도 하다.Also, according to an embodiment, the test screening item for which the correlation coefficient R with at least one of the screening items included in the test model exceeds a predetermined threshold may be excluded from the test screening item. This is to prevent the adverse effect that the influence of the existing test items in the test model is reduced when the test items that are highly correlated with the test items already included in the test model are newly included in the test models. In addition, when the test model includes high correlation items, the problem is that the singularity phenomenon occurs in the matrix operation process performed in the principal component analysis process, and the eigen value is not calculated properly.

상기 상관 계수 연산의 정확성을 위하여, 상기 상관 계수는 상기 테스트 모델에 포함된 전체 검진 항목의 검진 수치를 가지는 샘플 데이터를 이용하여 산출될 수 있다.For the accuracy of the correlation coefficient calculation, the correlation coefficient may be calculated using sample data having the examination values of all the examination items included in the test model.

본 발명의 일 양태에서, 상기 현재 생존해 있는 시점에서부터 t년 후까지 생존확률(통계청 자료)는 예를 들어 표 2의 기대여명표 및 도 2를 사용할 수 있다. In one aspect of the present invention, the survival probability (statistical office data) from the presently-living point to the t-year point later can use, for example, the expected life table of Table 2 and FIG.

또한 본 발명의 일양태에서의 연령별 생존확률을 표 3에 나타내었다. Table 3 shows the probability of survival according to one embodiment of the present invention.

Figure pat00025
Figure pat00025

Figure pat00026
Figure pat00026

Figure pat00027
Figure pat00027

Figure pat00028
Figure pat00028

Figure pat00029
Figure pat00029

Figure pat00030
Figure pat00030

Figure pat00031
Figure pat00031

Figure pat00032

Figure pat00032

Figure pat00033
Figure pat00033

Figure pat00034

Figure pat00034

Figure pat00035
Figure pat00035

Figure pat00036
Figure pat00036

Figure pat00037
Figure pat00037

Figure pat00038
Figure pat00038

Figure pat00039
Figure pat00039

Figure pat00040
Figure pat00040

Figure pat00041

Figure pat00041

Figure pat00042
Figure pat00042

Figure pat00043
Figure pat00043

Figure pat00044
Figure pat00044

Figure pat00045
Figure pat00045

Figure pat00046
Figure pat00046

Figure pat00047
Figure pat00047

Figure pat00048
Figure pat00048

Figure pat00049
Figure pat00049

Figure pat00050
Figure pat00050

Figure pat00051
Figure pat00051

Figure pat00052

Figure pat00052

Figure pat00053
Figure pat00053

Figure pat00054
Figure pat00054

Figure pat00055
Figure pat00055

다른 관점에서, 본 발명은 x 축을 기대여명(ER), y 축을 상기 검강검진 결과를 이용하여 측정된 생체나이와 식(1)를 이용한 생존확률의 측정방법으로 측정된 T년 후까지 생존확률로 하여, 작성한 그래프에서, y축의 생존확률 값에 대응하는 x 축의 값을 검강검진을 받은 날로부터 남은 잔여수명으로 하는 것을 특징으로 하는 잔여수명 측정방법에 관한 것이다.In another aspect, the present invention provides a method for estimating the survival probability (T) of a living body by T year measured by a method of measuring survival probability using the bio-age and expression (1) And the value of the x-axis corresponding to the survival probability value of the y-axis in the created graph is set as the remaining remaining life from the date of the examination.

본 발명의 잔여수명 예측방법은 도 1에 나타낸 바와 같이, 측정대상의 성별,연령자료와 건강검진 자료와 생체나이 계산 알고리즘(식(2)~(3))을 이용하여, 개인의 생체나이를 계산하고, 상기 생체나이와 식(1)을 이용하여 연령별 t년 후 생존확률을 계산한 후, x 축을 기대여명(ER), y 축을 상기 방법으로 측정한 t년 후 생존확률로 하고, 상기 방법으로 측정된 t년 생존확률을 나타낸 그래프를 이용하여 성별 연령별 잔여수명을 계산하는 것이다. As shown in FIG. 1, the remaining life prediction method of the present invention is a method for estimating the remaining life of an individual using the gender, age data, health examination data, and bio-age calculation algorithm (equations (2) to (3) , And the survival probability after t years by age is calculated using the above-mentioned biological age and equation (1), and the survival probability after t year after the y-axis is measured by the above method is used as the x axis and the y axis, And the remaining life for each age group is calculated using a graph showing the t-year survival probability.

본 발명의 잔여수명 계산 시스템에서는 2013년 통계청 자료를 기반으로한 연령별 생존확률 값을 사용하였지만 (표 2 및 도 4 참조), 이에 한정되는 것은 아니며 매년 또는 일정 기간 마다 통계청에서 새로이 발표되는 값을 이용하여 정확도를 제고할 수 있음은 통상의 기술자에게는 자명한 것이다.
In the remaining life calculation system of the present invention, the age probability of survival based on the data of the National Statistical Office is used (see Table 2 and FIG. 4). However, the present invention is not limited to this, It is obvious to those skilled in the art that the accuracy can be improved.

구체적인 예를 통해 본 발명에 따른 잔여수명 예측방법을 설명한다. 한국인 50세 남성(출생나이 50세)의 평균 잔여수명은 30.57세(표 1 참조)로 나타났는데 이러한 50세 남성의 생존확률은 약 0.5595(표 2 참조)이다. A method for predicting residual service life according to the present invention will be described with reference to specific examples. The average remaining life of a 50-year-old male (age 50) is 30.57 years (see Table 1). The probability of survival for a 50-year-old man is about 0.5595 (see Table 2).

이후 x 축을 기대여명(ER), y 축을 상기 방법으로 측정한 t년 후 생존확률로 하고, 상기 방법으로 측정된 t년 생존확률을 나타낸 그래프에서, y축에서 이 생존확률 값의 위치를 잡고 여기서부터 x 축에 평행한 가상선을 설정하여, 50세 남성의 생존 확률 그래프와 만나는 점의 x 축 값이 바로 본 발명에 따른 잔여수명이 되는 것이다(도 4). In the graph showing the t-year survival probability measured by the above method, the x-axis is the survival probability (ER), the y-axis is the survival probability after t year measured by the above method, The x-axis value of the point meeting the 50-year-old male survival probability graph is the remaining life according to the present invention (Fig. 4).

만약, 한국인 50세 남성 중 생체나이가 출생나이보다 5세 더 많은 (D=생체나이-출생나이=5) 사람의 잔여수명을 생존확률 그래프를 이용하여 측정한다면, 잔여수명이 30.5세보다 현저히 낮은 것을 확인할 수 있었다(도 5).If the survival probability of a 50-year-old man is 5 years older than the age of birth (D = age of the organism - age of birth = 5), the remaining life is significantly lower than 30.5 years (Fig. 5).

또한, 한국인 50세 남성 중 생체나이가 출생나이보다 5세 더 적은 (D=생체나이-출생나이=-5)사람의 잔여수명을 생존확률 그래프를 이용하여 측정한 결과, 잔여수명이 30.5세보다 현저히 높은 것을 확인할 수 있었다(도 6).In addition, the survival probability of a 50-year-old Korean man whose living age was less than 5 years of age (D = body age - birth age = -5) (Fig. 6).

하나의 구체예에서, 본 발명에 따른 잔여수명 예측 방법을 이용하여 산출한 결과, 출생나이 50세인 남자의 기대수명을 확인한 결과, D(생체나이 - 출생나이)가 -4.5 세 일 때, 기대수명은 91.9세이었으며, D=0세 일 때 80.8세 (한국인 평균)이었으며, D=4.5세 일 때 69.3세였다. In one embodiment, the life expectancy of a man with a birth age of 50 years is calculated using the residual life prediction method according to the present invention. As a result, when D (age of the human body - birth age) is -4.5 years, Was 91.9 years old and was 80.8 years old (average Korean) when D = 0 years old and 69.3 years old when D = 4.5 years old.

또 다른 양태에서, 출생 나이 50세인 여자의 기대 수명을 확인한 결과, D = -4.3 세 일 때 94.3세였고, D = 0세 일 때 86.3세 (한국인 평균)이 었으며, D= 4.2세 일 때 79.1세였다. In another aspect, the life expectancy of a woman with a birth age of 50 years was 94.3 years at D = -4.3 years, 86.3 years at Korean D = 0 years, 79.1 years at D = 4.2 years, It was three.

본 발명의 다른 양태에서는 본 발명의 잔여수명 측정방법을 이용하여 한국인 남/여의 10,000명을 연령별로 나누어 잔여수명을 측정하였으며, 그 결과를 각각 도 7 및 도 8에 나타내었다.
In another embodiment of the present invention, the remaining life is measured by dividing 10,000 Korean male / female by age group using the remaining life measuring method of the present invention, and the results are shown in FIG. 7 and FIG. 8, respectively.

이상으로 본 발명 내용의 특정한 부분을 상세히 기술하였는 바, 당업계의 통상의 지식을 가진 자에게 있어서 이러한 구체적 기술은 단지 바람직한 실시 양태일 뿐이며, 이에 의해 본 발명의 범위가 제한되는 것이 아닌 점은 명백할 것이다. 따라서, 본 발명의 실질적인 범위는 첨부된 청구항들과 그것들의 등가물에 의하여 정의된다고 할 것이다.
While the present invention has been particularly shown and described with reference to specific embodiments thereof, those skilled in the art will appreciate that such specific embodiments are merely preferred embodiments and that the scope of the present invention is not limited thereto will be. Accordingly, the actual scope of the present invention will be defined by the appended claims and their equivalents.

Claims (8)

검강검진 결과를 이용하여 측정된 생체나이와 하기 식(1)를 이용한 생존확률의 측정방법;
Figure pat00056
.....(1)

상기 식 (1)에서, S0(t)는 통계청 자료 상의 현재 생존해 있는 시점에서부터 t년 후의 생존확률이고, D는 생체나이-출생나이이며, 는 D의 평균, a는 D가 생존확률에 미치는 영향력, (t)는 생체나이를 이용하여 보정계산된 t년 후까지 생존확률을 나타냄.
A method of measuring the survival probability by using the biological age and the following equation (1) measured using the result of the examining of the urine;
Figure pat00056
.....(One)

(1) where S 0 (t) is the survival probability after t years from the presently alive date on the statistical data, D is the age of the organism - birth date, D is the average of D, a is the survival probability (T) represents the probability of survival until t years after calibration using the age of the living body.
제1항에 있어서, 상기 생체나이는 하기 식 (2)~(3)에 의하여 측정되는 것을 특징으로 하는 방법:
Figure pat00057
...(2)
Figure pat00058

...(3)

여기서, preBA는 식 (3)으로 보정하기전의 가생체나이이고, BA는 가생체나이(preBA)를 식(3)으로 보정한 후 의 생체나이이며, xj는 피검진자의 j번째 검진 항목의 검진 수치,
Figure pat00059
는 j번째 검진 항목의 샘플 내 평균값, sd(xj)는 j번째 검진 항목의 샘 플 내 표준 편차, Bij는 주성분 분석 결과 j번째 검진 항목의 i번째 요인값이 각 샘플의 실제 나이를 설명할 수 있는 정도, m은 주성분 분석의 수행결과 얻어진 요인 개수, Pi는 i번째 요인 값이 preBA 계산에 주는 가중치, y는 각 개체의 출생나이(CA)이며,
Figure pat00060
는 생체나이(BA) 모델 구축에 사용되는 모든 개체의 출생나이(C A)의 평균, sd(y)는 y의 표준 편차이고,
Figure pat00061
는 각 샘플에 대하여 실제 나이 및 preBA를 회귀 분석한 결과 생성된 회귀식
Figure pat00062
를 구성하는 계수임.
2. The method according to claim 1, wherein the biological age is measured by the following formulas (2) to (3):
Figure pat00057
...(2)
Figure pat00058

... (3)

Here, preBA is the age of the subject before correction to equation (3), BA is the age of the subject after adjusting the age of the subject (preBA) to equation (3), and x j is the jth item of examination Examination figures,
Figure pat00059
Is the mean value in the sample of the jth item, sd (x j ) is the sample standard deviation of the jth item, and B ij is the actual age of each sample about, m is the weight, y is the number of factors obtained performance results of principal component analysis, P i is the i-th factor to the value that can be calculated preBA is birth age (CA) of each object,
Figure pat00060
Is the mean of the birth years (CA), sd (y) is the standard deviation of y for all individuals used in constructing the BA model,
Figure pat00061
Is the regression of the actual age and preBA for each sample,
Figure pat00062
.
제2항에 있어서, Bij는 검강검진을 통해 측정된 2개 이상의 변수 (Xj)를 입력값으로 하여 주성분 분석 통계알고리즘을 통해 산출되는 값인 것을 특징으로 하는 방법.
3. The method according to claim 2, wherein B ij is a value calculated through a principal component analysis statistical algorithm with two or more variables (X j) measured through a screening examination as an input value.
제2항에 있어서, 피검자가 남자일 경우, x는 허리둘레(WC), 1초 간 노력성 기량(FEV1), 간효소(G-GTP), 요소질소(BUN), 고지단백(HDL), 저지단백(LDL), 중성지방(TG), 공복 혈당(FBS), 체지방률(FBR), 근육률(BMR), 알부민-글로불린 비율(AGR) 및 수축기혈압(SBP)으로 구성된 군에서 선택되는 건강검진 항목의 결과수치인 것을 특징으로 하는 방법.
3. The method according to claim 2, wherein when the subject is a male, x is at least one of a waist circumference (WC), an FEV1 for one second, G-GTP, urea nitrogen, A physical examination selected from the group consisting of low fat protein (LDL), triglyceride (TG), fasting blood sugar (FBS), body fat percentage (FBR), muscle rate (BMR), albumin-globulin ratio (AGR), and systolic blood pressure Item is a result value of the item.
제2항에 있어서, 피검자가 남자일 경우, x는 허리둘레(WC), 수축기혈압(SBP ), 간효소(G-GTP), 고지단백(HDL), 중성지방(TG), 공복혈당(FBS), 헤모글로빈(Hemoglobin), 저지단백(LDL), 체질량지수(BMI) 및 맥압(PP)으로 구성된 군에서 선택되는 건강검진 항목의 결과수치인 것을 특징으로 하는 방법.
The method according to claim 2, wherein when the subject is a male, x is at least one selected from the group consisting of WC, SBP, G-GTP, HDL, TG, ), Hemoglobin, low-density lipoprotein (LDL), body mass index (BMI), and pulse pressure (PP).
제2항에 있어서, 피검자 가 여자일 경우, x는 허리둘레(WC), 1초간 노력성 기량(FEV1), 감마-GTP)G-GTP), 요소질소(BUN), 고지단백( HDL), 저지단백(LDL), 중성지방(TG), 공복혈당(FBS), 체지방률(FBR), 근육률(BMR), 알부민-글로불린 비율(A GR), 이완기혈압(DBP), 적혈구침강속도(ESR) 및 체질량지소(BMI)로 구성된 군에서 선택되는 건강검진 항목의 결과수치인 것을 특징으로 하는 방법.
3. The method according to claim 2, wherein when the subject is an exciter, x is expressed as a function of the waist circumference (WC), the FEV1 for 1 second, gamma-GTP) G-GTP, urea nitrogen, (LDL), triglyceride (TG), fasting blood sugar (FBS), body fat percentage (FBR), muscle rate (BMR), albumin-globulin ratio (A GR), diastolic blood pressure (DBP), erythrocyte sedimentation rate And body mass index (BMI). ≪ RTI ID = 0.0 > 21. < / RTI >
제2항에 있어서, 피검자가 여자일 경우, x는 허리둘레(WC), 이완기혈압(DBP), 간효소(AST), 고지단백(HDL), 중성지방(TG), 공복혈당(FBS), 헤모글로빈(Hemoglobin), 저지단백(LDL), 체질량지수(BMI) 및 맥압(PP)으로 구성된 군에서 선택되는 건강검진 항목의 결과수치인 것을 특징으로 하는 방법.
The method according to claim 2, wherein when the subject is an female, x is at least one selected from the group consisting of waist circumference (WC), diastolic blood pressure (DBP), liver enzyme (AST), high fat protein (HDL), triglyceride (TG) Wherein the value is a result of a medical examination item selected from the group consisting of hemoglobin, low density lipoprotein (LDL), body mass index (BMI) and pulse pressure (PP).
x 축을 기대여명(ER), y 축을 제1항 내지 제7항 중 어느 한 항의 방법으로 측정된 t년 후 까지 생존확률로 하여 작성한 그래프에서, y축의 피검자의 생존확률 값에 대응하는 x 축의 값을 피검자의 검강검진을 받은 날로부터 남은 잔여수명으로 하는 것을 특징으로 하는 잔여수명 측정방법.

In the graph prepared by taking the x axis as the survival probability (ER) and the y axis as survival probabilities until t years measured by the method of any one of claims 1 to 7, the value of the x axis corresponding to the survival probability value of the subject on the y axis Is set as the remaining life remaining after the date of receipt of the examination by the examinee.

KR1020150030331A 2015-03-04 2015-03-04 Method for Assessinging Residual Life Using Biological Age KR101669526B1 (en)

Priority Applications (5)

Application Number Priority Date Filing Date Title
KR1020150030331A KR101669526B1 (en) 2015-03-04 2015-03-04 Method for Assessinging Residual Life Using Biological Age
JP2017546614A JP6640868B2 (en) 2015-03-04 2016-01-13 Prediction method of remaining life using living age
PCT/KR2016/000353 WO2016140432A2 (en) 2015-03-04 2016-01-13 Method for predicting residual life using biological age
CN201680020966.9A CN108475542A (en) 2015-03-04 2016-01-13 The method that remaining life is predicted by using biological age
HK18115784.8A HK1256728A1 (en) 2015-03-04 2018-12-10 Method for predicting residual life using biological age

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020150030331A KR101669526B1 (en) 2015-03-04 2015-03-04 Method for Assessinging Residual Life Using Biological Age

Publications (2)

Publication Number Publication Date
KR20160107508A true KR20160107508A (en) 2016-09-19
KR101669526B1 KR101669526B1 (en) 2016-10-26

Family

ID=56848379

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020150030331A KR101669526B1 (en) 2015-03-04 2015-03-04 Method for Assessinging Residual Life Using Biological Age

Country Status (5)

Country Link
JP (1) JP6640868B2 (en)
KR (1) KR101669526B1 (en)
CN (1) CN108475542A (en)
HK (1) HK1256728A1 (en)
WO (1) WO2016140432A2 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20190131860A (en) * 2018-05-17 2019-11-27 재단법인차세대융합기술연구원 Method, system and non-transitory computer-readable recording medium for providing life age
KR102418342B1 (en) * 2021-08-13 2022-07-08 주식회사 유투메드텍 Apparatus for predicting glucose change using dietary information and method thereof

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018137203A1 (en) * 2017-01-25 2018-08-02 深圳华大基因研究院 Method for determining population sample biological indicator set and predicting biological age and use thereof
KR102106428B1 (en) * 2018-02-19 2020-05-06 주식회사 셀바스에이아이 Method for predicting health age
CN108847284B (en) * 2018-05-02 2021-03-23 莱博生物科技股份有限公司 Human body biological age measuring and calculating device and system
CN110277151A (en) * 2019-06-11 2019-09-24 浙江大学 Human Physiology Analysis of age method, system and model based on routine physical examination index
JP7277351B2 (en) * 2019-12-23 2023-05-18 Tvs Regza株式会社 Viewing history analysis device and viewing history analysis device program
CN111243738B (en) * 2019-12-31 2023-04-25 东软集团股份有限公司 Method, device and equipment for constructing survival analysis model and predicting survival rate
JP6901169B1 (en) * 2020-02-25 2021-07-14 日新ビジネス開発株式会社 Age learning device, age estimation device, age learning method and age learning program
KR102371440B1 (en) * 2021-08-28 2022-03-07 유진바이오소프트 주식회사 The method and system for generating model predicting personalized biological age
KR20230114081A (en) 2022-01-24 2023-08-01 주식회사 에비드넷 Health age and health score calculation method changes in health
KR20240012704A (en) 2022-07-21 2024-01-30 주식회사 로그미 An apparatus and a method for predicting biological age
CN115274098B (en) * 2022-08-19 2023-04-07 北京咯咯网络科技有限公司 Intelligent system for predicting height based on bone age and height

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20060111940A (en) * 2005-04-25 2006-10-31 김건형 System for predicting remains life and predicting server controlling method thereof
KR20140126229A (en) * 2013-11-20 2014-10-30 주식회사 바이오에이지 Biological age calculation model generation method and system thereof, biological age calculation method and system thereof

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4062910B2 (en) * 2001-11-29 2008-03-19 株式会社日立製作所 HEALTH MANAGEMENT SUPPORT METHOD AND DEVICE AND HEALTH LIFE LIFE PREDICTION DATA GENERATION METHOD AND DEVICE
KR20060131154A (en) * 2005-06-15 2006-12-20 주식회사 리얼 에이지 System for measuring biological age and service method thereof
JP4619219B2 (en) * 2005-07-19 2011-01-26 株式会社エヌ・ティ・ティ・データ Subject selection device
JP5231835B2 (en) * 2008-02-22 2013-07-10 株式会社東芝 Health guidance support system
US8409104B2 (en) * 2008-05-02 2013-04-02 Conopco Inc. Heart age assessment
JP5185785B2 (en) * 2008-11-19 2013-04-17 オムロンヘルスケア株式会社 Health condition judgment device
EP2781602A1 (en) * 2013-03-21 2014-09-24 Universität Konstanz Method for the determination of biological age in human beings

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20060111940A (en) * 2005-04-25 2006-10-31 김건형 System for predicting remains life and predicting server controlling method thereof
KR20140126229A (en) * 2013-11-20 2014-10-30 주식회사 바이오에이지 Biological age calculation model generation method and system thereof, biological age calculation method and system thereof

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
강중철 외 2인, "장수리스크를 고려한 사망률 추정방법에 관한 연구", 리스크 관리 연구 Vol. 17(1), pp153~178 *
'바이오에이지, 생체나이 측정시스템 중국 진출', 매일경제, 2009.03.16. *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20190131860A (en) * 2018-05-17 2019-11-27 재단법인차세대융합기술연구원 Method, system and non-transitory computer-readable recording medium for providing life age
KR102418342B1 (en) * 2021-08-13 2022-07-08 주식회사 유투메드텍 Apparatus for predicting glucose change using dietary information and method thereof

Also Published As

Publication number Publication date
JP2018512659A (en) 2018-05-17
CN108475542A (en) 2018-08-31
HK1256728A1 (en) 2019-10-04
JP6640868B2 (en) 2020-02-05
KR101669526B1 (en) 2016-10-26
WO2016140432A2 (en) 2016-09-09

Similar Documents

Publication Publication Date Title
KR101669526B1 (en) Method for Assessinging Residual Life Using Biological Age
Kaess et al. Aortic stiffness, blood pressure progression, and incident hypertension
CN110808097A (en) Gestational diabetes prediction system and method
Priya et al. A novel approach to predict diabetes by using Naive Bayes classifier
Bikia et al. Noninvasive estimation of aortic hemodynamics and cardiac contractility using machine learning
Gutiérrez-Tobal et al. Ensemble-learning regression to estimate sleep apnea severity using at-home oximetry in adults
CN113782183A (en) Pressure damage risk prediction device and method based on multi-algorithm fusion
Zheng et al. A novel machine learning-based systolic blood pressure predicting model
Bacchetti et al. Non-Markov multistate modeling using time-varying covariates, with application to progression of liver fibrosis due to hepatitis C following liver transplant
Mamun Cuff-less blood pressure measurement based on hybrid feature selection algorithm and multi-penalty regularized regression technique
Sayed et al. Interchangeability between Pentacam and IOLMaster in phakic intraocular lens calculation
Dessiso et al. Bayesian joint modelling of longitudinal and survival data of HIV/AIDS patients: a case study at bale robe general hospital, Ethiopia
Teitelbaum et al. Limited sampling strategies supporting individualized dose adjustment of intravenous busulfan in children and young adults
KR101398986B1 (en) Medical biological age measurement system and apparatus
Donaldson et al. Continuous glucose monitoring-based composite metrics: a review and assessment of performance in recent-onset and long-duration type 1 diabetes
CN107003315A (en) Insulin secreting ability analytical equipment, the insulin secreting ability analysis system and insulin secreting ability analysis method for possessing the device
Aggarwal A machine learning based approach for the identification of insulin resistance with non-invasive parameters using homa-IR
Mai et al. Non-Laboratory-Based Risk Factors for Automated Heart Disease Detection
Rihana et al. Artificial intelligence framework for COVID19 patients monitoring
Kaypakli et al. The effect of heart rate and pulse pressure on mean arterial pressure: the combined formula for calculation of mean arterial pressure
Pisaruk et al. Impact of COVID-19 on the signs of human ageing
Shruthi et al. Diabetes prediction using machine learning technique
Cooke et al. Methodological considerations for the measurement of arterial stiffness using applanation tonometry
Bouwmans et al. PRM32 Feasibility and Validity of the TiC-P
JP2011257816A (en) Apparatus and program for estimating state of medical examinee

Legal Events

Date Code Title Description
A201 Request for examination
E902 Notification of reason for refusal
E90F Notification of reason for final refusal
E701 Decision to grant or registration of patent right
GRNT Written decision to grant
FPAY Annual fee payment

Payment date: 20190605

Year of fee payment: 4