Contact pressure explains half of the abdominal aortic aneurysms wall thickness inter-study variability

This is a code used to generate distribution functions of parameters used in the study Contact pressure explains half of the abdominal aortic aneurysms wall thickness inter-study variability. Jupyter Notebook with Python code (anizo_merge.ipynb)

Steps to reproduce:

Run anizo_merge.ipynb in Jupyter Notebook.

Package versions used:

numpy : 1.24.3
ndsplines : 0.1.5
pymc : 5.1.2
pytensor : 2.10.1
scipy : 1.10.1
arviz : 0.15.1
pandas : 2.0.0
matplotlib: 3.7.1
seaborn : 0.12.2

MCMC trace with 200000 samples is stored in tr_An01_StHeDimaMartTong_FT.nc file.

Results of computationally intensive steps can be loaded from the attached files.

本代码用于生成题为《接触压力可解释半数腹主动脉瘤壁厚度的研究间变异》（Contact pressure explains half of the abdominal aortic aneurysms wall thickness inter-study variability）的研究中所用参数的分布函数。附带Python代码的Jupyter Notebook文件为anizo_merge.ipynb。

复现步骤：
在Jupyter Notebook中运行anizo_merge.ipynb。

所用软件包版本如下：
numpy : 1.24.3
ndsplines : 0.1.5
pymc : 5.1.2
pytensor : 2.10.1
scipy : 1.10.1
arviz : 0.15.1
pandas : 2.0.0
matplotlib: 3.7.1
seaborn : 0.12.2

含200000个采样点的马尔可夫链蒙特卡洛（Markov Chain Monte Carlo, MCMC）轨迹存储于tr_An01_StHeDimaMartTong_FT.nc文件中。

计算密集型步骤的结果可从附带文件中加载。