Data and code for "A kinetic model of continuous radiation damage to populations of cells: Comparison to the LQ model and application to molecular radiotherapy"
Description
Replication data for our study "A kinetic model of continuous radiation damage to populations of cells: Comparison to the LQ model and application to molecular radiotherapy". Introduction The folder Code_and_Data contains the following two routines (functions) and two datasets. The code is written in Matlab (.m functions), while the datasets are saved in Matlab workspace files (.mat). If any reader wants data in any other format, please contact us and we will try to accommodate her/his request. The folder contains the following items: annealing_model_fit.m model_dif_radionuclide_stochastic.m data_unkel.mat data_iodine.mat model_dif_radionuclide_stochastic.m Function containing the implementation of our code. Detailed info about the inputs and outputs of this function can be found in it. annealing_model_fit.m Function containing the implementation of an annealing method to obtain best fits of our model to experimental data. In particular, in this case it is used to fit experimental data taken from Unkel et al of surviving fractions of different cell lines. Detailed info about the inputs and outputs of this function can be found in it. data_unkel.mat Workspace containing the experimental results taken from Unkel et al of surviving fractions of different cell lines, together to best-fits of such data to the LQ model and to our model. Experimental data (doses and log(SF)) for each of the 9 cell lines is save in variables C2, DG, FP, HP, L3, MP, P1, PT, and S2, and altogether are also reported in the structure EXPERIMENTAL_DATA. Variables named par_LQ_”cell_line” contain best-fitting LQ parameters for each cell-line, and variables SF_LQ_”cell_line” contain best-fitting LQ surviving fractions. Variables named par_opt_”cell_line” contain best-fitting parameters of our model for each cell-line, and variables SF_opt_”cell_line” contain best-fitting LQ surviving fractions. data_iodine.mat This workspace contains the sets of parameters used to create figures 2 and 3 (par_1, par_2, par_3, par_4). The variable D contains the simulated dose rate of a radioiodine treatment.