BYOM package GUTS walkthrough
- Author: Tjalling Jager
- Date: April 2017
- Web support: http://www.debtox.info/byom.html
Step-by-step walk through the code of the GUTS package for the BYOM platform. This walk through is made with the 'publish' option in Matlab, which might also be very convenient to keep track of your work (as a modeller's log book).
The two example files in the main directory of this package are a good place to start. They use the reduced GUTS model (toxicokinetics and damage dynamics combined in a single one-compartment model), and only the two simplest cases (pure SD or pure IT). These examples apply the analytical solution to the model, which is very fast and allows you to try out the various options without very long waiting times. More complete GUTS models, including options for working with pulsed exposure, can be found in the directories 'reduced' and 'full'.
This walk through consists of the following files:
- byom_guts_start.m: an example script to demonstrate GUTS on a straightforward survival data set. This demo shows the profile likelihood for individual parameters and calculation of LC50s.
- byom_guts_extra.m: an example script to demonstrate more elabrate options (especially for confidence intervals). This demo shows the likelihood-region method to provide the joint-confidence region of the parameters, and how it can be used to make intervals on various types of plots.
- derivatives.m: the actual model in the form of a system of ordinary differential equations (ODEs).
- simplefun.m: the actual model in the form of a system of explicit functions (e.g., an analytical solution).
- call_deri.m: calls derivatives.m or simplefun.m to calculate the model output. Here, the model output can be transformed to match the data, if needed.
- pathdefine.m: a piece of code that searches for the engine directory, and adds it to the Matlab path. No need to make changes here, just make sure it is in every directory from which you run BYOM scripts.