Thyrosim.jl is an updated version of THYROSIM. It heavily relies on the amazing DifferentialEquations.jl package to solve systems of ODEs.
This tool is developed to
- Better optimize replacement LT4 and LT4+LT3 dosing for treating hypothyroid patients, based on their individual hormone levels, heights, weights and gender
- Aid in understanding more precisely how gender and BMI impact thyroid hormone transient and steady state dynamical regulation over time in these patients.
Download and install Julia. This package supports Julia v1.6+ for Mac, Linux, and window machines. Within Julia, copy and paste the following:
using Pkg
pkg"add https://github.com/biona001/Thyrosim.jl"For running the examples below, the following packages are also necessary
pkg"add Plots"Load packages
using Thyrosim, PlotsIn plots below, orange lines indicate normal range.
sex = false # true = male, false = female
h = 1.63 # unit = meters
w = 60 # unit = KG
sol = simulate(h, w, sex) # solve ODE
plt = output_plot(sol, title="Female, $h meters, $w KG") # make plot
sex = false # true = male, false = female
h = 1.63 # unit = meters
w = 90 # unit = KG
T4dose = 120.0 # unit = mcg
dial = [0.0; 0.88; 0.0; 0.88] # no T4 secretion, normal T4 excretion, no T3 secretion, normal T3 excretion
sol = simulate(h, w, sex, dial=dial, T4dose=T4dose) # solve ODE
plt = output_plot(sol, title="T4 monotherapy on thyroidectomized patient") # make plot
sex = false # true = male, false = female
h = 1.63 # unit = meters
w = 45 # unit = KG
T4dose = 75.0 # unit = mcg
T3dose = 5.0 # unit = mcg
RTF = 0.10 # 10% residual thyroid function
dial = [RTF; 0.88; RTF; 0.88] # T4 secretion, normal T4 excretion, T3 secretion, normal T3 excretion
sol = simulate(h, w, sex, dial=dial, T4dose=T4dose, T3dose=T3dose) # solve ODE
plt = output_plot(sol, title="T4/T3 combination therapy on hypothyroid patient") # make plot
sex = true # true = male, false = female
h = 1.78 # unit = meters
w = 70 # unit = KG
T4dose = 600 # unit = mcg
interval = 24*5 # unit = hours
dial = [0.0; 0.88; 0.0; 0.88] # no T4 secretion, normal T4 excretion, no T3 secretion, normal T3 excretion
sol = simulate(h, w, sex, dial=dial, warmup=false, T4dose=T4dose, dosing_interval=interval) # solve ODE
plt = output_plot(sol, title="Absurdly large T4 dose in 5 day intervals") # make plot
If you use Thyrosim.jl for academic purposes, please cite the following paper
If you encounter a bug or need user support, please open a new issue on Github. Please provide as much detail as possible for bug reports, ideally a sequence of reproducible code that lead to the error.
PRs and feature requests are welcomed!