Tip: Comments are indicated in blue. For additional details, click on the keywords highlighted in green. Reasonable default values are specified for all model inputs.
/*
TLK Model Calibration Example: Use this sample file to begin analyzing delayed-plating cell survival data (no repopulation effects; GF=0.0). For small datasets (e.g., only low or high dose rate survival data are available), it may be prudent to perform the data analysis using the LPL or RMR model instead of the TLK. Nevertheless, the TLK model can still provide some useful insights into the radiobiology -- even for small datasets. For example, the TLK might be used to investigate the impact on parameter estimates of biphasic DSB rejoining kinetics. Recommendations for TLK model fitting are provided below. For large datasets (multiple dose rates and/or split-dose survival data), the TLK can provide better fits to the data than the LPL or RMR model.Small datasets (low or high dose rate data only): Use the default parameters indicated below and treat the ETA and PHI parameters as adjustable (2 or 3 parameter fit). Large datasets (multiple dose rates): Use the default parameters indicated below. An adequate fit to most datasets can be obtained by treating 3 to 5 of the TLK model inputs as adjustable. Adjust ETA, PHI (both kinds of DSB), and RHT for simple DSBs to perform a 3 parameter fit. For a 4 parameter fit, adjust ETA, PHI, and the RHT for both simple and complex DSBs. And for a 5 parameter fit, adjust ETA, PHI, the RHT for simple and complex DSBs, and the initial yield of complex DSBs ( DSB). Start with a 3 or 4 parameter fit and increase the number of adjustable parameters until you obtain a good fit to the data. */ MODEL: DRM=TLK CKM=QECK !Model Specification CELL: DNA=5.667D+09 NC=46 IRAD: DSB=25 FCB=0.2 !Radiation damage to DNA CLRP: RHT=HT4R ETA=2.75E-04 PHI=0.01 A0=PMR GAM=0.25 CCKM: TCC=10 TPOT=20 GF=0.0 N0=1 KAP=1.0D+38 VOL=1 RBM: KX=6.0E+06 FMAX=1 !Analysis of PFGE data
! TLK model inputs (simple and complex double strand breaks) HT4R = 0.25 7 / !Half-times for DSB repair PMR = 0.75 0.25 / !Probability DSB is rejoined correctly ! Parameters that impact on solution accuracy SOLVER: DBLV=0 ETOL=1.0E-09 MXSS=999999 HMAX=10 ! Parameters that control the time steps and stopping criteria SIMCON: FSDX=1 TSAX=12 TCUT=96 FRDL=1.0D+00 RCUT=1.0D-09 BCONF: MITR=10000 MXFN=100000 SGT=1.0D-04 SST=1.0D-09 FCT=1.0D-09 MCAT: NPAR=4 !npar=2 or 3 for small datasets, npar=3-5 for large datasets AMP1: CID=CLRP PNAM=ETA PLB=0.0 PUB=1.0D-02 PIG=2.0E-04 AMP2: CID=CLRP PNAM=PHI PLB=0.0 PUB=0.25 PIG=0.01 AMP3: CID=CLRP PNAM=RHT PLB=0.01 PUB=12 PIG=0.15 LI=1 !Simple DSB AMP4: CID=CLRP PNAM=RHT PLB=0.01 PUB=12 PIG=8 LI=2 !Complex DSB AMP5: CID=IRAD PNAM=DSB PLB=20 PUB=60.0 PIG=25 !Total DSB
/* Figure of Merit (FOM) Weight 1: cell survival data 2: DSB rejoining data (FAR) 3: DSB rejoining data (Avg. DSB/cell) 4: Transformed cells per survivor 5: Transformed cells per irradiated cell */ FOM: UW=LFWT LFWT = 1 0 0 0 0 /
|