Intraluminal thrombus effect on the progression of abdominal aortic aneurysms by using a multistate continuous-time Markov chain model ====================================================================================================================================== * Liangliang Zhang * Byron A. Zambrano * Jongeun Choi * Whal Lee * Seungik Baek * Chae Young Lim ## Abstract **Objective** Recent studies reported the intraluminal thrombus (ILT) among others strongly affects abdominal aortic aneurysms (AAA) expansion rates. Thus, we investigate characteristics of ILT with AAA expansion. **Methods** We applied homogeneous multistate continuous-time Markov chain models to longitudinal data of 26 Korean AAA patients as a retrospective clinical study. We considered four states of AAA and maximal thickness of ILT (maxILT), fraction of wall area covered by ILT (areafrac) and fraction of ILT volume (volfrac) as possible covariates. **Results** Based on likelihood-ratio statistics, areafrac is the most significant biomarker and maxILT is the second most significant. Besides, within AAAs that developed an ILT layer, we found that an AAA expands relatively fast at an early stage but the rate goes slower once AAA reaches in a large size. **Conclusion** Results recommends surgical intervention when any patient with areafrac more than 60% or maxILT more than 30mm. Although this recommendation should be considered with caution given the limited sample size, one can use the proposed model as a tool to find such recommendations with their own data. Keywords * Abdominal aortic aneurysms * Intraluminal thrombus * Markov chain model * survival analysis ## Introduction Abdominal aortic aneurysm (AAA), the dilatation of an aorta at an abdominal level is a common threatening disease that affects 9.5% of the elderly population (>65 years)1. This dilatation of the abdominal aorta can cause death when it ruptures. The rupture occurs when the stress on the AAA wall overcomes the wall strength. Since the maximum diameter is being positively associated with the wall stress2 and with AAA expansion3, it has been used as a risk factor for AAA growth and subsequently for rupture. Current medical treatment suggests that once a patient is diagnosed with AAA, s/he should be under surveillance until the aneurysm reaches 5.5 cm in diameter4 in the U.S. Some studies, however, challenge the 5.5 cm threshold criterion since small AAAs still rupture. For example, 10–24% of ruptured AAAs were less than 5 cm in diameter5, although more cases succumb to rupture prior to surgical intervention as the diameter increases. In fact, the screening criterion varies by countries6. Therefore, we take advantage of this association and investigate the statistical distribution of the maximum diameters of AAAs by quantitatively dichotomizing the patients into four different prognosis groups: early, mild, severe and fatal. In this study, patients were observed over time and morphological characteristics of AAA were collected as well. With the time-to-event data, researchers and practitioners usually would like to make inferences on the progression path (e.g. deterioration or recovery) of the disease. Survival analysis is a commonly used statistical tool to achieve such goal. However, usual survival analysis can only deal with binary events (alive or death). Instead, the multistate model7 is able to model longitudinal studies where individuals may experience several events, and conduct forward or backward conversions between events. Therefore, we adopted a multistate continuous-time Markov Chain model to estimate the transition from one state to another for AAA. In this model, the transition intensities provide the hazards8 and the survival probabilities for AAA progression. We can also calculate the mean sojourn time9 in a given state. In addition, we incorporate the intraluminal thrombus (ILT) information as covariates into the model through transition intensities as ILT is strongly associated with AAA expansion rates. Choosing the appropriate covariate to predict AAA growth rate and time-to-event is an important task, because an accurate prediction can allow personalized clinical management and proper timing of surgery10. Patient-specific prediction of AAAs have been performed via Bayesian methods combining sequential CT images11 and biomedical and computational models12 without taking into account the intraluminal thrombus (ILT). There has been, however, substantial heterogeneity of AAA expansion rates and an accurate growth prediction is still remaining a top priority to improve the prediction capability. Numerous variables have been suggested as potential predictors for growth and its rupture13 and, among them, the ILT is being recognized as a potential metric of AAA growth14. For examples, Stevens *et al*.15 obtained a set of follow-up CT images from four patients and estimated AAA growth along with geometrical characteristics of ILT. They found that the intraluminal ILT volume and maximum ILT thickness were correlated with AAA volume growth. With 14 patients, Zambrano *et al*.16 classified scan images of AAAs into two groups of AAAs by the ILT areal fractions and found that AAAs grew significantly faster in the group with larger ILT compared to the group with lesser ILT. Nevertheless, those studies of ILT on AAA growth were found only from on a small set of AAAs and few studied on which geometrical characteristics of ILT were more linked to the prediction of AAA expansion. Therefore, using a larger dataset of follow–up AAAs, this study calls to investigate which geometrical characteristics of ILT can enhance the prediction of AAA expansion and their transition probabilities. Particularly, we study the association between AAA growth and ILT by considering a homogeneous multistate continuous-time Markov chain model with covariates related to ILT. Generally, the progression of AAA growth can be related to other factors such as age, gender, smoking, and vascular diseases17. Our retrospective data, on the other hand, were originally obtained for characterizing relevant morphological parameters so that the clinical and biochemical data were not available. As our focus is on influence of ILT information for AAA progression, we extracted ILT thicknesses, percentage of the luminal area covered during the AAA progression and the volume of ILT per each scan and applied a statistical model to predict AAA expansion based on these geometrical characteristics of ILT. By incorporating contributing factors such as ILT’s volume or fraction, a fitted model can be used to predict the progression of AAA and its rate from one state to another state. ## Materials & Methods This study uses longitudinal computed tomography (CT) images from 26 de-identified Korean patients with each patient having up to 7 surveillance scans. Thus, this study is the retrospective review of clinical data. The study was reviewed and exempted by Internal Review Board at Michigan State University. For each scan, a maximum spherical diameter *D* was measured and used to categorize AAAs into 4 stages18: 1) ‘early’ (*D* < 40mm), 2) ‘mild’ (40mm ≤ *D* <47mm), ‘severe’ (47mm ≤ *D* < 58 mm), 4) ‘fatal’ (*D* ≥ 58mm). Gharahi *et al*.18 proposed a diameter measurement called the spherical diameter that the largest sphere fits within the aorta and found that the spherical diameter measurement has the minimum uncertainty compared to others, which has an advantage for the growth prediction, although the spherical diameters are slightly smaller than the traditional orthogonal diameters (average 4 mm with the maximum value of 9 mm). Note that there are variations in AAA screening among countries. However, the typical orthogonal diameter criterion for surgical recommendation is 50mm or 55mm4,6. Given such variations in criterion and use of spherical diameter in our study, it is reasonable to assume that a patient in severe state is recommended for surgical intervention and a patient in fatal state should be guided for a prompt surgery. Also, for each scan, maximal thickness of ILT, fraction of wall area covered by ILT and fraction of ILT volume were measured by the approach introduced in Zambrano *et al*. (2016)16 and denoted by maxILT, areafrac and volfrac, respectively. AAA maximum spherical diameter information of all patients at different scan times are shown in Figure 1. There are cases where a state of AAA growth is moved to another state (e.g. move from state ‘early’ to state ‘mild’ or from state ‘mild’ to state ‘severe’.). Note that one patient (P14) experienced a decrease of the maximum diameter before an expansion. There are also cases that scan times are relatively short. One patent (P21) has much smaller diameter compared to the other patients. The measurements of this patient could be an outlier or due to the characteristics of the spherical diameter (smaller than that of the traditional orthogonal diameter18). As P21 stayed in the “early” state according to our criterion, the results would not change much with/without this patient so we decide to keep the patient in the analysis. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/02/2020.08.31.20185330/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2020/09/02/2020.08.31.20185330/F1) Figure 1. Plot of AAA spherical diameter vs. age for each patient. Each line corresponds to a patient. Patient ID is provided on the lower left corner of each line. Dots on each line indicate when the CT image was taken. A multistate continuous-time Markov chain model is used to model the progression of diseases (see, e.g. Jackson *et al*., 200319) which can provide information of disease progression by transition probabilities between different states and mean sojourn time at one state. Specifically, we used a homogeneous model by assuming a transition intensity (the rate of a transition probability) is independent of time. This model allows a transition probability that changes over time but the rate remains constant. In addition, we incorporate ILT information into the model by considering the logarithm of an intensity as a linear function of explanatory variables from ILT. To investigate which model fits the data better, we consider likelihood ratio (LR) statistics given the proposed statistical model. The LR statistics is the ratio of two likelihoods (or -2 log likelihood ratio), which is used to perform a statistical hypothesis test for model comparison. The likelihood in the numerator is from the null model and the likelihood in the denominator is from the alternative model. Thus, a larger value of the likelihood ratio statistics (or a smaller value of -2 log likelihood ratio) indicates the data support the null model. The estimation algorithm is converged produced the stable results. ## Results Recall that we consider the maximal thickness of ILT (maxILT), the fraction of wall area covered by ILT (areafrac) and the volume fraction of ILT (volfrac) into the model as covariates. These variables are continuous variables. Note that these three variables are highly correlated with each other. In particular, the correlation between the fraction of wall area covered by ILT and the volume fraction of ILT is 0.89. High correlation between covariates can lead unstable estimates if they are included in the model together. Thus, we consider a model with one covariate each. Although three variables are highly correlated, as shown in Figure 2, they have different effects on the transition probability curves at different times. We want to investigate which variable is a significant biomarker for the progression of AAA. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/02/2020.08.31.20185330/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2020/09/02/2020.08.31.20185330/F2) Figure 2. Fitted transition probability curves at different years. In each graph, the red circle line (p12) denotes transition probabilities of moving from ‘early’ state to ‘mild’ state; the blue triangle line (p23) denotes transition probabilities of moving from ‘mild’ state to ‘severe’ state; the purple plus line (p32) denotes transition probabilities of moving from ‘severe’ state to ‘mild’ state: the green cross line (p34) denotes transition probabilities of moving from ‘severe’ state to ‘fatal’ state. The x-axis is years. The y-axis is the transition probability. ### Model comparison Table 1 shows likelihood ratio statistics (−2 log likelihood ratio) and the corresponding p-value for each model when the null model is the model without covariates. This result shows that the model with one of such covariates is overall statistically significant. Although we cannot directly compare which covariate fits the data better than the others by a hypothesis test as they are not nested models, the larger LR statistics could imply a better fit as the LR statistics of each model is calculated with the same null model. View this table: [Table 1.](http://medrxiv.org/content/early/2020/09/02/2020.08.31.20185330/T1) Table 1. Likelihood ratio statistics of different models. Another evidence that these three covariates contribute much to the progress of the AAA enlargement is the comparison of transition probability curves under different models (see Figures 2a, 2b, 2c and 2d). The fitted transition probability curve without covariates has relatively lower peaks than all the others, especially when the transition is from the ‘early’ state to the ’mild’ state. In addition, by finding the peak of the curves *p*12 and *p*23, we conclude that patients in ‘early’ state are most likely to move to ‘mild’ state between 1 and 2 years, while patients in ‘mild’ state are most likely to move to ‘severe’ state around 4 years. The model with areafrac has the highest transition probability from the ‘mild’ state to the ‘severe’ state among the three models given the severe criterion within 3 years (Figures 2b, 2c and 2d). This may indicate the progression from the ‘mild’ state to the ‘severe’ state is sensitive to the changes of areafrac. In other words, when a patient already has a mild AAA, areafrac can be a useful biomarker to predict his/her state after a certain time. A prevalence plot could also show goodness of fit for a multistate continuous Markov Chain model and Pearson tests can provide goodness of fit for the models as well. These results indicate that the Markov model with maxILT or areafrac as a covariate passed the test of a goodness of fit. We provide detailed results in the supplementary document. Among the models with ILT characteristics as a covariate, we investigate the model using wall area covered by ILT (areafrac) with more details. Results of the other two models (maxILT and volfrac) are provided in the supplementary document. ### Fraction of wall area covered by ILT The value of areafrac ranges from 0 to 0.74, and about 75% of cases are less than 0.4. So the value 0.4 (75% quantile) could be used as a criterion for the severity of ILT. As an extreme case, we also consider 0.6 which is 95% quantile. With these values for areafrac, we provide estimated transition probabilities for the duration of 3 years (Table 2a, 2b). Also we provide mean sojourn time in Table 2c. According to the estimated transition probability ((1,1)-th entry of Table 2a), we can see that there is no chance that a patient with areafrac being 0.4 is in ‘early’ state. Besides, such a patient will instantly transit to or already stay in either ‘mild’ or ‘severe’ state within 3 years. By comparing transition probabilities when areafrac is 0.6 (Table 2b) with those when areafrac is 0.4 (Table 2a), we also conclude that disease state of AAA progresses significantly with the increase of areafrac. Also, mean sojourn time for ‘early’ state is significantly small compared to those for other states. Although the range of AAA size for each state is rather different, this still implies that a patient stays ‘early’ for a shorter time with 0.4 of areafrac. View this table: [Table 2.](http://medrxiv.org/content/early/2020/09/02/2020.08.31.20185330/T2) Table 2. Results from multistate continuous-time Markov chain models with the covariate areafrac. (a) and (b) denote estimated transition probabilities (ETP) at the different values of the covariate and time. (c) denote the estimated mean sojourn time (EMST). To see the effects of areafrac on the survival probability, we plot fitted survival probability curves with Kaplan-Meier (KM) curves, the empirical estimate of the survival probability20. By comparing Figure 3a with Figure 3b, we find that the model with larger areafrac as a covariate has a lower survival probability which implies that a patient with larger areafrac is exposed to a higher risk of entering into the last stage (‘fatal’ state). This result is consistent with what we have found earlier. One can compare two survival curves, for example, the one with areafrac = 0.4 and the other with areafrac = 0.6 by two sample Kolmogorov-Smirnov test21, which indicates the two curves are statistically different. The detailed results of the test are provided in the supplementary document. ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/02/2020.08.31.20185330/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2020/09/02/2020.08.31.20185330/F3) Figure 3. Comparison of empirical and fitted survival probability for two multistate continuous-time Markov chain models with areafrac. The red solid line denotes the fitted survival curve. The blue dashed line denotes the empirical survival. The red dotted line denotes the 95% confidence interval of the fitted survival. A bigger set of patient data should be tested in order to reach a stronger conclusion. However, despite the fact that a small data set was used, a trend was clearly seen. If we believe that a surgical intervention should be recommended when survival probability lower than 40% at 4 years, Figure 3b suggests that a patient with 0.6 of areafrac is being considered for such recommendation. For maxILT, 30mm was the value we obtain for such recommendation. Note that depending on the practitioner’s criterion on the survival probability (e.g. 40%) and the time period (e.g. 4 years), the recommendation will change. However, the model we used can provide such information adaptively. ## Discussion A homogeneous multistate continuous-time Markov chain model is considered for the analysis of AAA progression to investigate transition of progression from one state to the another state in AAA growth. The model allows to estimate mean sojourn time for a patient being in each state. Given the estimated mean sojourn times for ’mild’ and ’severe’ states from the model with areafrac, we may say that an approximated time for 1mm increase in the ’mild’ state is shorter than that for the ’severe’ state, which implies that AAA expands relatively fast at an early stage but the rate gets slower once AAA is enlarged enough. The model with ILT area fraction has the highest transition probabilities among the three models within 3 years at the value of the third quartile (severe condition). This finding may imply that the progression from ‘mild’ state to ‘severe’ state is sensitive to the changes of ILT area fraction. In other words, when someone already has a mild AAA, ILT area fraction may be a useful biomarker to predict his/her state after a certain time. Likelihood ratio statistics supports that the model with one of three covariates fits the data better than the model without ILT information. This may suggest that when we consider AAA growth, we should not ignore ILT information such as maximal ILT thickness, ILT area fraction, ILT volume fraction information of the patients. More specifically, if we believe that the patient with the fitted survival probability lower than 40% at 4 year should be recommended for surgical intervention, we may suggest any patients with ILT area fraction more than 60% is highly recommended for surgical intervention, considering the low survival probabilities in these two progression stages. This recommendation could be changed depending on the threshold of survival probability and the time. Also, the results would be different than the reported due to the dependency of the chosen stages threshold criteria. More detailed analysis on the appropriate threshold should be conducted with a larger cohort before reaching stronger conclusion. Despite this limitation, our results are in agreement with other previous studies such as Domonkos *et al*. (2019)22 and Metaxa *et al*. (2015)23. The thrombus accumulation is an ongoing process presented in 75% of detected AAA and the prevalence of finding such accumulation increases as the aneurysm size increases24. It is, therefore, expected to observe AAAs without ILT at the early stage (early) of this disease. This study found that most of our AAAs were remained without ILT accumulation during the early stage and while they transitioned to the following stage, a significant ILT accumulation was observed. Over this early stage, a faster AAA expansion in comparison with other stages was found. This could suggest that the faster growth may promote ILT accumulation. Additionally, this observation and the fact that ILT accumulation rate was found to be the same to the AAA expansion rate16 would explain the proposed idea of using sudden increase of thrombus as a potential predictor of rupture. Once the initial ILT accumulation occurs (usually in the transition between ‘early’ to ‘mild’ state), our results have shown that the next transition stage would be sensitive to the amount of lumen area covered by the thrombus (areafrac). The intraluminal thrombus is a complex biological entity composed also by many inflammatory cells, including macrophages and neutrophils25 that not just interfere with the direct wall shear stress and strain relationship between lumen wall and blood flow but that also interact biochemically (promoting wall thinning, cell inflammation, degradation of the extracellular matrix26) and biomechanically (modifying wall shear stress27) with the arterial wall. It is then understandable that a larger area covered by thrombus would have a greater impact on the AAA growth process and thus in our transition time. Furthermore, our results may indicate that the risk of a patient to transition to the ‘fatal’ last stage is outweighed as a larger ILT thickness (maxILT) is observed. These results agree with signs of lower pO2 levels and signs of hypoxia in aneurysm with ILT thickness greater than 4 mm28. The data were originally obtained for characterizing relevant morphological parameters so that the clinical and biochemical data such as smoking habit (current, ex-smoking and how long) and blood pressure were not diligently collected. Also, we have only three female patients so that we could not investigate the effect by gender. Thus, we were unable to add more patient-specific information. If we have enough observations from the control group (individual with no AAA) and patients whose AAA ruptured or received surgical treatment, we would have more reliable estimates of transition probabilities as well as mean sojourn times. Also, thorough investigation on the effects of ILT information as well as other patient-specific information to the AAA growth can be done with a larger size of data. Therefore, we would like to point out that our findings bring an interesting hypothesis which should be verified by further studies with a larger group of AAA patients. On the other hand, to our knowledge, there were only a few studies, available for the morphological data, obtained from follow-up patients. In fact, this study was one of large data sets, in terms of characterizing morphological parameters in the follow-up study. The reason of the scarcity of such data should be due to requiring manual time consuming segmentation or development of automatic process of quantifying the morphological parameters. We think that, even with the small size of samples, we was able to demonstrate the new statistical tool in the role of thrombus effect and predicting the progression of AAAs. For a future study, we wish to get clinical data with various risk factors and we plan to develop automatic morphological quantification so that we can extract ILT information easily, which helps to make a larger data for the analysis. A growth rate/expansion rate could be also an informative criterion to determine the stages of AAA as well. However, we decided to select a maximum diameter given that it is the clinically important metric to assess the onset, progression, and risk of rupture in AAAs. The use of this widely known relationship to stratify AAAs severity would make it easier for the reader to relate our results to other patients’ cohort or to patients in clinical practice. ## Data Availability The data that support the findings of this study are available on request from the corresponding author. ## Declaration of conflicting interests The authors declare that there is no conflict of interest. ## Funding The authors gratefully acknowledge the support, in part, by the National Heart, Lung, and Blood Institute of the National Institutes of Health (R01HL115185 and R21HL113857) and National Science Foundation (CMMI-1150376). Choi and Lim were supported in part by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (2018R1A4A1025986 and 2019R1A2C1002213). ## Acknowledgements All the codes are written and performed in R (r-project.org version 3.4.4.). The R package that we mainly used is “msm” and usage is introduced in Jackson *et al*. (2011) 29. * Received August 31, 2020. * Revision received August 31, 2020. * Accepted September 2, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NoDerivs 4.0 International), CC BY-ND 4.0, as described at [http://creativecommons.org/licenses/by-nd/4.0/](http://creativecommons.org/licenses/by-nd/4.0/) ## References 1. 1.Alcorn HG, Wolfson SJ, Sutton-Tyrrell K et al. Risk factors for abdominal aortic aneurysms in older adults enrolled in the cardiovascular health study. Arter Thromb Vasc Biol 1996; 16(8): 963–970. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYXR2YmFoYSI7czo1OiJyZXNpZCI7czo4OiIxNi84Lzk2MyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA5LzAyLzIwMjAuMDguMzEuMjAxODUzMzAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 2. 2.Long A, Rouet L, Lindholt JS et al. Measuring the maximum diameter of native abdominal aortic aneurysms: review and critical analysis. Eur J Vasc Endovasc Surg 2012; 43(5): 515–524. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ejvs.2012.01.018&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22336051&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F02%2F2020.08.31.20185330.atom) 3. 3.Brady AR, Thompson SG, Fowkes FG et al. Abdominal aortic aneurysm expansion: risk factors and time intervals for surveillance. Circulation 2004; 110(1): 16–21. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTQ6ImNpcmN1bGF0aW9uYWhhIjtzOjU6InJlc2lkIjtzOjg6IjExMC8xLzE2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDkvMDIvMjAyMC4wOC4zMS4yMDE4NTMzMC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 4. 4.Brewster DC, Cronenwett JL, Hallett JJ et al. Guidelines for the treatment of abdominal aortic aneurysms. Report of a subcommittee of the joint council of the american association for vascular surgery and society for vascular surgery. J VASC SURG 2003; 37(5): 1106–1117. 5. 5.Nicholls SC, Gardner JB, Meissner MH et al. Rupture in small abdominal aortic aneurysms. Journal of vascular surgery 1998; 28(5): 884–888. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0741-5214(98)70065-5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9808857&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F02%2F2020.08.31.20185330.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000076878800016&link_type=ISI) 6. 6.Stather P, Dattani N, Bown M et al. International variations in aaa screening. European Journal of Vascular and Endovascular Surgery 2013; 45: 231–234. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ejvs.2012.12.013&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23332307&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F02%2F2020.08.31.20185330.atom) 7. 7.Meira-MachadoL, deUña-Álvarez J, Cadarso-Suvrez C et al. Multi-state models for the analysis of time-to-event data. Statistical methods in medical research 2009; 18(2): 195–222. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/0962280208092301&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18562394&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F02%2F2020.08.31.20185330.atom) 8. 8.Spotswood LS, Julia ER, Michael G et al. Hazard ratio in clinical trials. Antimicrob Agents Chemother 2004; 48(3): 2787–2796. [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWFjIjtzOjU6InJlc2lkIjtzOjk6IjQ4LzgvMjc4NyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA5LzAyLzIwMjAuMDguMzEuMjAxODUzMzAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 9. 9.Stephen WD. Screening, Sojourn Time. USA: John Wiley & Sons, Ltd, 2005. 10. 10.Lee R, Jarchi D, Perera R et al. Applied Machine Learning for the Prediction of Growth of Abdominal Aortic Aneurysm in Humans. EJVES Short Reports 2018; 39: 24—-28. 11. 11.Do HN, Ijaz A, Gharahi H et al. Prediction of abdominal aortic aneurysm growth using dynamical gaussian process implicit surface. IEEE T BIO-MED ENG 2018; 66(3): 609–622. 12. 12.Zhang L, Jiang Z, Choi J et al. Patient-specific prediction of abdominal aortic aneurysm expansion using bayesian calibration. IEEE J BIOMED HEALTH 2019;. 13. 13.Shum J, Martufi G, Di Martino E et al. Quantitative assessment of abdominal aortic aneurysm geometry. Ann Biomed Eng 2010; 39(1): 277–286. 14. 14.Speelman L, Schurink GW, Bosboom EM et al. The mechanical role of thrombus on the growth rate of an abdominal aortic aneurysm. J VASC SURG 2010; 51(1): 19–26. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jvs.2009.08.075&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19944551&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F02%2F2020.08.31.20185330.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000273708500005&link_type=ISI) 15. 15.Stevens RR, Grytsan A, Biasetti J et al. Biomechanical changes during abdominal aortic aneurysm growth. PLoS ONE 2017; 12(11): e0187421. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F02%2F2020.08.31.20185330.atom) 16. 16.Zambrano BA, Gharahi H, Lim C et al. Association of intraluminal thrombus, hemodynamic forces, and abdominal aortic aneurysm expansion using longitudinal CT images. ANN BIOMED ENG 2016; 44: 1502–1514. 17. 17.Vega de Céniga M, Gómez R, Estallo L et al. Growth rate and associated factors in small abdominal aortic aneurysms. EUR J VASC ENDOVASC 2006; 31(3): 231–236. 18. 18.Gharahi H, Zambrano B, Lim C et al. On growth measurements of abdominal aortic aneurysms using maximally inscribed sphere. MED ENG PHYS 2015; 37: 683–691. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.medengphy.2015.04.011&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26004506&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F02%2F2020.08.31.20185330.atom) 19. 19.Jackson CH, Sharples LD, Thompson SG et al. Multistate markov models for disease progression with classification error. J ROY STAT SOC D-STA 2003; 52(2): 193–209. 20. 20.Manish KG, Pardeep K and Jugal K. Understanding survival analysis: Kaplan–Meier estimate. Int J Ayurveda Res 2010; 1(4): 274–278. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4103/0974-7788.76794&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21455458&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F02%2F2020.08.31.20185330.atom) 21. 21.Stephens MA. Edf statistics for goodness of fit and some comparisons. Journal of the American statistical Association 1974; 69(347): 730–737. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2307/2286009&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1974U523600021&link_type=ISI) 22. 22.Domonkos A, Staffa R and Kubíč Effect of intraluminal thrombus on growth rate of cek L. abdominal aortic aneurysms. International Angiology 2019; 38: 39–45. 23. 23.Metaxa E, Kontopodis N, Tzirakis K et al. Effect of intraluminal thrombus asymmetrical deposition on abdominal aortic aneurysm growth rate. Journal of Endovascular Therapy 2015; 22: 406–412. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/1526602815584018&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25900725&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F02%2F2020.08.31.20185330.atom) 24. 24.Shindo S, Matsumoto H, Kubota K et al. Is the size of an abdominal aortic aneurysm associated with coagulopathy? World J Surg 2005; 29(7): 925–929. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00268-005-7718-z&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15951938&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F02%2F2020.08.31.20185330.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000230819400024&link_type=ISI) 25. 25.Quintana RA and Taylor WR. Cellular Mechanisms of Aortic Aneurysm Formation. CIRC RES 2019; 124: 607—-618. 26. 26.Virag L, Wilson JS, Humphrey JD et al. Potential biomechanical roles of risk factors in the evolution of thrombus-laden abdominal aortic aneurysms. INT J NUMER METH BIO 2017; 33(12): e2893. 27. 27.Georgakarakos E, Ioannou CV, Volanis S et al. The influence of intraluminal thrombus on abdominal aortic aneurysm wall stress. J VASC SURG 2009; 28(4): 325–333. 28. 28.Vorp DA, Lee PC, Wang DH et al. Association of intraluminal thrombus in abdominal aortic aneurysm with local hypoxia and wall weakening. J VASC SURG 2001; 34(2): 291–299. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1067/mva.2001.114813&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11496282&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F02%2F2020.08.31.20185330.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000170473900017&link_type=ISI) 29. 29.Jackson CH et al. Multi-state models for panel data: the msm package for r. Journal of statistical software 2011; 38(8): 1–29.