Academic Journal of Science and Technology ISSN: 2771-3032 | Vol. 9, No. 1, 2024 241 Prediction of Medical Trauma Rehabilitation by GM (1,1) Model Xiaoyu Yang1, * 1School of Surveying and Land Information Engineering, Henan Polytechnic University, Jiaozuo, Henan 454000, China *Corresponding author: YANG Xiaoyu (E-mail: yxy-yxy0915@163.com) Abstract: With the continuous development of science and technology, scholars at home and abroad have adopted different methods to predict the recovery of trauma patients. GM(1,1) gray model was constructed to predict the trauma recovery of patients with anal fistula.For the anal fistula trauma data of small samples, the model constructed in this paper can still meet the prediction demand.The constructed prediction model can realize the multi-stage prediction of the wound condition of the patient, and the accuracy of the model is high, which is helpful for the medical staff to adopt the appropriate treatment strategy according to the rehabilitation condition of the patient. Keywords: GM(1,1); Prediction; Trauma; Test of model. 1. Introduction Trauma is mainly a type of damage caused to organs and tissues of the human body by various external forces or mechanical factors, including cuts, punctures, contusions, etc[1].Trauma not only exists during times of war, but also frequently occurs during times of peace. Trauma not only has a high incidence, but also varies greatly in severity, making it a prominent global issue.Traumatic patients often have complex injuries, rapid disease progression, and a risk of death[2].According to incomplete statistics, approximately one million people worldwide die from trauma every year, and in China, only a few hundred thousand people die from trauma each year[3].Trauma has become the leading cause of death for young adults, attracting widespread attention from countries around the world and various sectors of society, and becoming a focus of discussion among many scholars[4]. Scholars have established dynamic trauma prediction models based on machine learning algorithms such as random forests and support vector machines.There are also scholars who use texture synthesis technology to simulate wound healing image prediction algorithms for myofibroblasts.Guo Chengyu and others developed a trauma bleeding volume grading prediction model based on deep learning methods to assist in predicting trauma bleeding volume[5]. The GM (1,1) model is suitable for predicting small sample data and meets the needs of predicting trauma rehabilitation, so this prediction model was adopted in the study.Predicting the trend of changes in patient trauma area and circumference over time can help medical staff evaluate the degree of recovery from trauma, further implement predictive treatment for patients, and reduce medical safety risks.At the same time, by predicting and reducing trauma mortality and disability rates, we can reduce the harm that trauma brings to humanity[6]. 2. Experimental Data The entire experimental data was collected from 4 patients with 5 stages of anal fistula medical trauma data, as shown in Table 1. Table 1. Medical trauma data of patients with anal fistula Patient Number of periods Area of trauma/cm2 Perimeter of trauma/cm 1 1 2 3 4 5 0.7235 0.6497 0.5581 0.4873 0.4024 4.6694 4.2183 3.7236 3.2941 2.8493 2 1 2 3 4 5 1.0510 1.0034 0.9525 0.9049 0.8550 5.7021 5.4276 5.1318 4.8589 4.5472 3 1 2 3 4 5 0.8574 0.7998 0.7369 0.6731 0.6233 4.6941 4.3983 4.0816 3.7514 3.4456 4 1 2 3 4 5 0.5948 0.5322 0.4487 0.3806 0.3071 4.0537 3.6614 3.2018 2.8135 2.3631 3. Basic Theories and Methods 3.1. Basic Principles In the late 1970s and early 1980s, Professor Deng Julong, a renowned scholar in China, proposed the Grey System Theory, which is an original discipline in China and an important discipline in the fields of deterministic theory and systems science[7]. Usually, a system where some information is known but some information is unclear is called a grey system. Grey system theory is used to solve problems of data scarcity and uncertainty. The grey system theory model, also known as the grey model or grey dynamic model, abbreviated as the GM model. Grey prediction is a method of establishing a GM model that extends from the past to the future based on known or uncertain information from the past, in order to understand the development trend of events and provide a basis for planning decisions and state evaluation. Grey prediction is 242 mainly built on the basis of modeling. In addition to predicting future values such as population and food consumption, it can also predict the time when outliers occur, that is, disaster prediction[8]. 3.2. Research Method The grey system theory solves the problems of data scarcity and uncertainty. The GM (1,1) model is the most widely used grey dynamic prediction model in grey system theory. Due to its advantages of requiring less raw data and high accuracy, it has been widely applied in the field of medical research in recent years [9].In the GM (1,1) model, G represents Grey and M represents Model. The first 1 in parentheses represents the first order of the differential equation, while the second 1 represents a variable in the differential equation. Its basic method is to accumulate the initial data to generate a data sequence, and then modify the accumulated sequence to a nearest neighbor generated sequence. The corresponding curve of the new data sequence can be approximated infinitely by a specific curve. The approximation curve is used as the basic model, and the predicted values are rolled and accumulated several times to predict the development trend. This study used the GM (1,1) model to construct a trauma prediction model, predicting the trauma of patients, and used the small residual probability P and a posterior difference ratio C value to evaluate the accuracy of the model. 3.3. Construction of GM (1,1) Model 3.3.1. Feasibility Analysis of Modeling The rank ratio formula is nk kx kx k ,...,3,2, )( )1( )( )0( )0( )0(    , if all )()0( k are within the interval ),( 1 2 1 2   nn ee (e is the natural logarithm, n is the sample size analyzed), then through the rank ratio test, the original time series can be predicted using the model GM (1,1) for grey prediction. If the level ratio is outside the interval, it does not meet the grey prediction condition and needs to be processed. 3.3.2. Data Preprocessing (1) Record non negative original time series ;,...2,1)},(),...,2(),1({)( )0()0()0()0( nknxxxkX  (2) Record the one-time accumulation of the original time series to generate a sequence )}(),...,2(),1({)( )1()1()1()1( nxxxkX  , ;,...,2,1 nk  (3) Revise a sequence generated by accumulation to a sequence generated by the nearest mean )}(),...,2(),1({)( )1()1()1()1( nzzzkz  , .,...,3,2),1(5.0)(5.0)( )1()1()1( nkkxkxkz  3.3.3. Modelling Principles The GM (1,1) model constructs differential equations on preprocessed data sequences, obtains the time response function of the differential equations, and then uses the cumulative method to perform reverse operations to obtain predicted values. The specific operations are as follows: (1) The grey differential equation of the GM (1,1) model is bkazkx  )()( )1()0( .The parameter in the formula is the development grey number, and parameter b is the endogenous control grey number, )()0( kx is grey derivative, )()1( kz is whitening background value; (2) Introducing a data matrix                   1)( ...... 1)3( 1)2( )1( )1( )1( nz z z B ,data vector                )( ... )3( )2( 0 0 0 nx x x Yn , parameter vector        b a u ; (3) Estimating the Parameter Vector of GM (1,1) Model Using Least Squares Method n TT YBBB b a u 1)(ˆ ˆ ˆ        ; (4) Substitute the estimated values of parameters a and b into the whitening equation btax dt tdx  )( )( )1( )1( ,the solution to the equation, also known as the time response function, is obtained as a b e a b xtx at   ))0(()( )1()1( , )1()0( )0()1( xx  ; (5) The time response sequence of GM (1,1) grey differential equation is: nk a b e a b xkx ak ,...,2,1,))1(()1(ˆ )0()1(   ; (6) The original data sequence is nkkxkxkx ,...,2,1),(ˆ)1(ˆ)1(ˆ )1()1()0(  . 3.3.4. Accuracy Inspection The GM (1,1) model may have certain errors in practical application, so after selecting the model, the accuracy of the model needs to be verified before the problem can be predicted and analyzed. The GM (1,1) model generally adopts three testing methods, namely residual test, correlation test, and posterior test. Only when all three tests are passed, can the constructed model meet the requirements and can be used for problem prediction and analysis. (1) Residual Test Set absolute residuals )(ˆ)()( )0()0( kxkxk  ,relative residual )( )( )0( kx k k   , Average relative residual nk n n k kk ,...,2,1, 1 1     ,when the average relative residual value is less than 0.2, the model is called a qualified residual model. (2) Correlation Test Correlation coefficient nk kk kk k ,...,2,1, )}(max{)( )}(max{)}(min{ )(       ,  is resolution coefficient, the smaller the value of  , the greater the difference between the correlation coefficients, 243 and the stronger the discriminative ability, 5.0 ,correlation    n k k n r 1 )( 1  ,when the correlation is greater than 0.6, the correlation test is passed. (3) Posterior Error Test Posteriori difference test is a test of the statistical characteristics of residual distribution, recording the mean square deviation of the original time series   2 1 1 2)0()0( 1 1 )(                  n xkx S n k , Mean squared deviation of residuals   nk n kk S n k ,...,2,1, 1 )()( 2 1 1 2 2                   , )0(x is the average of the original data, )(k is the mean of absolute residuals. By calculating the posterior difference ratio 1 2 S S C  ,Small residual probability  )()(6745.0 1 kkSpP  determine the accuracy of the constructed model,the smaller the C value, the higher the credibility of the constructed model. The larger the P value, the better the accuracy of the model fit. The reference table for model accuracy level discrimination is shown in Table 2: Table 2. Reference Table for Precision Level Discrimination of GM (1,1) Model Model Accuracy level P C Excellent P≥0.95 C≤0.35 Qualified 0.80≤P<0.95 0.350.65 4. Construction of Predictive Models After calculation, the maximum level ratio of the original time series was 1.2393, and the minimum level ratio was 1.0474, both falling within (0.7165, 1.3956). The level ratio passed the test.Build a differential equation model and calculate the development grey number a, endogenous control grey number b, correlation degree r, and posterior error ratio C for each prediction model, as shown in Table 3.The predictive equations of the GM (1,1) model for the recovery of medical trauma in four patients are: .8310.317773.27)1( ;8661.32713.3)1( ;1460.614519.56)1( ;7843.109269.9)1( ;6547.1019526.95)1( ;4876.204366.19)1( ;6241.399547.34)1( ;2642.55407.4)1( 1418.0)1( 4 1780.0)1( 4 0813.0)1( 3 0840.0)1( 3 0583.0)1( 2 0530.0)1( 2 1288.0)1( 1 1547.0)1( 1                 k p k a k p k a k p k a k p k a ekx ekx ekx ekx ekx ekx ekx ekx )1()1( 1 kx a represents the predictive equation for the trauma area of the first patient; )1()1( 1 kx p represents the predictive equation for the circumference of the first patient's trauma and so on. 5. Discussion and Conclusion Research has shown that the area and circumference of anal fistula wounds in four patients gradually decrease over time, and the wounds heal steadily.At present, scholars mostly predict trauma complications, trauma outcomes, and trauma stress disorders. Lee Na Hyeon, Kim Seon Hee et al. used chest computed tomography to quantify lung contusion volume and predict respiratory complications in patients with chest trauma [10].Scholars have used injury severity score risk assessment tools to predict trauma outcomes [11]. In order to scientifically grasp the trend of trauma development, this article establishes a trauma GM (1,1) prediction model. The establishment of a trauma prediction model can improve the effectiveness of trauma treatment and increase the success rate of trauma treatment. Therefore, the study of trauma prediction models has important practical significance. Although this study has made certain progress and achievements, there are still some shortcomings. The GM (1,1) model has a strong dependence on historical data, and the established model does not consider the relationship between various factors, resulting in a large error; Based on the research results of this article, the next step will be to optimize the initial conditions of the GM (1,1) prediction model and improve its prediction accuracy. References [1] Bai Xiangjun,Gao Wei,and Li Zhanfei.Promoting the construction and graded treatment of trauma centers to improve the level of trauma treatment [J].Chinese Journal of Emergency Medicine,2013,22 (6):567-569. [2] Wang Zhengguo.Disease Trauma in Modern Society [J].Journal of the Chinese Academy of Medical Sciences,2002 24 (1):1-3. [3] Yao Yuan,Development of Trauma Database and Input and Analysis of Trauma Data [D].Third Military Medical University,2006:121. [4] Borse,N.N.and A.A.Hyder,Call for more research on injury from the developing world:results of a bibliometric analysis.The Indian journal of medical research,2009.129(3). [5] Farrow,L.,et al.,Interpretation and reporting of predictive or diagnostic machine-learning research in Trauma & Orthopaedics.Bone Joint J,2021.103-B(12):p.1754-1758. 244 [6] Guo Chengyu.Development and Verification of a Dynamic Prediction Model for Traumatic Fatal Massive Hemorrhage based on Vital Sign Time Series Data [J]. Journal of PLA Medicine, 2023:1-12. [7] Deng Julong.Introduction to Grey System Theory. Inner Mongolia Electric Power [J],1993 (03):51-52. [8] Xiaoyun,Z.,et al.,Prediction of China’s Grain Consumption from the Perspective of Sustainable Development—Based on GM(1,1) Model.Sustainability,2022.14(17). [9] Liang Binghua,Huang Lifeng.Prediction and analysis of medical and health resources in China based on GM (1,1) grey prediction model [J].Modern Preventive Medicine,2021,48(20): 3655-3659. [10] Lee,N.H.,et al.,Prediction of respiratory complications by quantifying lung contusion volume using chest computed tomography in patients with chest trauma. Scientific Reports, 2023.13(1). [11] Santos,J. et al.,Comparison of National Surgical Quality Improvement Program Surgical Risk Calculator and Trauma and Injury Severity Score Risk Assessment Tools in Predicting Outcomes in High-Risk Operative Trauma Patients.The American surgeon,2023.