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.
/*
RMR Model Calibration Example: Use this sample file to begin analyzing delayed-plating cell survival data (no repopulation effects; GF=0.0). If the dataset includes multiple dose rates (low and high dose rates), set the initial DSB yield (DSB parameter) to a value in the range from about 25 to 60 DSB/Gy-cell and adjust the PHI, RHT (or LAM), and ETA parameters. For small datasets, set the DSB parameter to a value in the range from 25 to 60 DSB/Gy-cell, and set RHT to a value in the range from 0.25 to 6 hour. Then adjust the PHI and ETA parameters. For small datasets, it is very worthwhile to perform a few sensitivity studies to ensure that an optimal fit has been identified. For example, performs fits with different values for DSB and RHT.
*/ MODEL: DRM=RMR CKM=QECK !Model Specification
CELL: DNA=5.667D+09 NC=46
IRAD: DSB=25 !Radiation damage to DNA
CLRP: RHT=2.0 ETA=2.75E-04 PHI= 0.01 A0=0.00 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
! 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.0E+00 RCUT=1.0D-09 BCONF: MITR=10000 MXFN=100000 SGT=1.0D-04 SST=1.0D-09 FCT=1.0D-09
MCAT: NPAR=3
AMP1: CID=CLRP PNAM=BMR PLB=0 PUB=1.0E-02 PIG=1.0E-04
AMP2: CID=CLRP PNAM=PHI PLB=0 PUB=1.0 PIG=0.01
AMP3: CID=CLRP PNAM=RHT PLB=0.01 PUB=12.0 PIG=2
/*
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 /
|
