The goal of this report is to develop code in Python that implements and evaluates a multivariate linear error propagation for commonly used DCE-MRI models. Two pharmoacokinetic models (PK-models) were tested; modified Kety model and two-compartment exchange model. The analysis compared a Monte carlo simulated signal and a simple analytic model for uncertainty. The purpose of comparing was to have the Monte carlo method and the linear error estimation method reasonably close in terms of estimation of the error of the PK-parameters (parameters are specific for the PK-model). The analysis included an inspection of the resulting plots of the Coefficent of Variation (CV) of the estimation of the PK-parameters, when increasing the ratio of the true value and the added error for each of the input parameters. Generally the CV increased, I.e. the precision decreased, as the noise level increased. For both PK-models, error in injection time of contrast agent were the input parameter with the least tolerance of noise. By comparing the results from both methods, the conclusion is to implement the modified Kety model as the PK-model used in MICE Toolkit (NONPI Medical AB, Umeå, Sweden) and keep the noise to signal level below 15\% in order to obtain precision of the results to less than 33\% error. We have verified the uncertainty estimation method works well for the Modified Kety model but not so well for the two-compartment exchange model. The scope of this study only included two PK-models and one signal model, therefore it would be beneficial to test the linear estimation method used here on other PK-models and signal models to find what models it works well for. / Målet med denna rapport är att utveckla kod i Python som implementerar och utvärderar multivariat linjär felpropagering för ofta använda DCE-MRI modeller. Två farmakokinetiska modeller (PK-modeller) testades; den modifierade Kety-modellen och två-utrymme utbytesmodellen. Analysen jämförde en Monte carlo-simulerad signal och en enkel analytisk modell för osäkerhet. Syftet med jämförelsen var att ha Monte carlo-metoden och den linjära feluppskattningsmetoden någorlunda nära när det gäller uppskattning av PK-parametrarnas fel (parametrarna är specifika för PK-modellen). Analysen innebar en inspektion av de resulterande kurvorna för variationskoefficienten (CV) för uppskattningen av PK-parametrarna, när förhållandet mellan det sanna värdet och det pålagda felet för var och en av ingångsparametrarna ökades. Generellt ökade CV:n, det vill säga precisionen minskade när bruset/felet ökade. För båda PK-modellerna var det fel i injektionstid för kontrastmedel den ingångsparameter som hade minst tolerans för brus. Genom att jämföra resultaten från båda metoderna är slutsatsen att använda den modifierade Kety-modellen som PK-modell i MICE Toolkit (NONPI Medical AB, Umeå, Sverige) och hålla brus-till-signal nivån under 15\% för att få precision av resultaten till mindre än 33\% fel. Vi har verifierat att metoden för osäkerhetsuppskattning fungerar bra för den modifierade Kety-modellen men inte så bra för två-utrymme utbytesmodellen. Omfattningen av denna studie inkluderade endast två PK-modeller och en signalmodell, därför skulle det vara fördelaktigt att testa den linjära uppskattningsmetoden som används här på andra PK-modeller och signalmodeller för att hitta vilka modeller den fungerar bra för.
| Identifer | oai:union.ndltd.org:UPSALLA1/oai:DiVA.org:umu-201567 |
| Date | January 2022 |
| Creators | Höglund, Jonas |
| Publisher | Umeå universitet, Institutionen för fysik |
| Source Sets | DiVA Archive at Upsalla University |
| Language | English |
| Detected Language | Swedish |
| Type | Student thesis, info:eu-repo/semantics/bachelorThesis, text |
| Format | application/pdf |
| Rights | info:eu-repo/semantics/openAccess |
Page generated in 0.003 seconds