;model definition
<MODEL>

file='./arthritisModel_mlxt.txt'



<DESIGN>
[ADMINISTRATION]
trt_0 = {time = {1, 8, 15, 22, 29}, amount = 0}
trt_1 = {time = {1, 8, 15, 22, 29}, amount = 70}
trt_10 = {time = {1, 8, 15, 22, 29}, amount = 700}

;parameter initial values
<PARAMETER>
; Initialization parameters
a = 1
b = 0.5
I0 = 2.52

; PK model 
alpha = 0.02327
beta = .045
CL = 2.5
V1 = 15
V2 = 25

; Effect model 
sigma1 = 0.154
sigma2 = 0.065
sigma3 = 0.003

; Arthritis model parameters
k1 = 0.183
k2 = 0.092
k3 = 5
k4 = 0.064
k5 = 0.016
tau_pop = 11.2
omega_tau = 3

;prediction outputs and grid
<OUTPUT>
list={TAS,AKS, G, I, D,Cc}
grid=0:.01:40

<RESULTS>
[GRAPHICS]
pTAS = {y={TAS},  ylabel = 'Total Arthritic Score I(t)+D(t)', xlabel = 'time (days)'}
pAKS = {y={AKS},  ylabel = 'Ankylosis Score D(t)', xlabel = 'time (days)'}
pG = {y={G},  ylabel = 'GM-CSF G(t)', xlabel = 'time (days)'}
pI = {y={I},  ylabel = 'Inflammation I(t)', xlabel = 'time (days)'}
gridarrange(pTAS, pAKS, pG, pI, 2,2)

<SETTINGS>
[GRAPHICS]
nb_simulations = 200