; MONOLIX project generated with IQRtools

; ==PROJECT HEADER START===================================================
; COMMENT             = 'PKPD model - MONOLIX version'
; TOOL                = 'MONOLIX'
; TOOLVERSION         = 'MLX2023R1'
; DATA                = '../../Ex7/data.csv'
; DOSINGTYPES         = 'INFUSION'
; TK0NAMES            = 'NA'
; COVNAMES            = 'WT,HT'
; CATNAMES            = ''
; CATCATEGORIES       = ''
; REGRESSIONNAMES     = 'Vp,Vc,CL,Q'
; OUTPUTS             = 'A_circ'
; ERRORMODELS         = 'rel'
; PARAMNAMES          = 'CIRC0,MTT,SLOPU,GAM'
; PARAMTRANS          = 'exp(phi),exp(phi),exp(phi),(phi)'
; PARAMINVTRANS       = 'log(psi),log(psi),log(psi),(psi)'
; COVARIATENAMES      = 'WT,HT'
; COVARIATESUSED      = ''
; BETACOVNAMES        = ''
; BETACOVTRANS        = ''
; BETACATNAMES        = ''
; BETACATREFERENCE    = ''
; NROBSERVATIONS      = '176'
; ==PROJECT HEADER END=====================================================

; =============================================
<DATAFILE>
; =============================================

[FILEINFO]
file = '../../Ex7/data.csv'
delimiter = comma
header = {IXGDF, USUBJID, ID, TIME, TIMEPOS, TAD, NAME, DV, UNIT, YTYPE, AMT, ADM, EVID, MDV, CENS, TINF, RATE, Vp, Vc, CL, Q, WT, HT}

[CONTENT]
ID = {use=identifier}
TIME = {use=time}
DV = {use=observation, name={y}, ytype={1}, type={continuous}}
YTYPE = {use=observationtype}
AMT = {use=amount}
ADM = {use=administration}
EVID = {use=eventidentifier}
MDV = {use=missingdependentvariable}
CENS = {use=censored}
TINF = {use=infusiontime}
Vp = {use=regressor}
Vc = {use=regressor}
CL = {use=regressor}
Q = {use=regressor}
WT = {use=covariate, type=continuous}
HT = {use=covariate, type=continuous}

; =============================================
<MODEL>
; =============================================

[COVARIATE]
input = {WT,HT}

EQUATION:
tWT = log(WT/75)
tHT = log(HT/175)

[INDIVIDUAL]
input = {CIRC0_pop, MTT_pop, SLOPU_pop, GAM_pop, omega_CIRC0, omega_MTT, omega_SLOPU, omega_GAM}

DEFINITION:
CIRC0 = {distribution=logNormal, typical=CIRC0_pop, sd=omega_CIRC0}
MTT = {distribution=logNormal, typical=MTT_pop, sd=omega_MTT}
SLOPU = {distribution=logNormal, typical=SLOPU_pop, sd=omega_SLOPU}
GAM = {distribution=normal, typical=GAM_pop, no-variability}

[LONGITUDINAL]
input = {b}

file = './project_model.txt'

DEFINITION:
y = {distribution=normal, prediction=A_circ ,errorModel=proportional(b)}

; =============================================
<FIT>
; =============================================
data = {y}
model = {y}

; =============================================
<PARAMETER>
; =============================================
CIRC0_pop = {value=7.21, method=MLE}
MTT_pop = {value=124, method=MLE}
SLOPU_pop = {value=28.9, method=MLE}
GAM_pop = {value=0.239, method=FIXED}
omega_CIRC0 = {value=0.2, method=MLE}
omega_MTT = {value=0.4, method=MLE}
omega_SLOPU = {value=0.3, method=FIXED}
b = {value=0.3, method=MLE}

; =============================================
<MONOLIX>
; =============================================

[TASKS]
populationParameters()
individualParameters(method = conditionalMode)
fim(method = Linearization)
logLikelihood(method = Linearization)
plotResult()

[SETTINGS]
GLOBAL:
seed = 123456
exportpath = 'RESULTSORIG'
nbchains = 2
autochains = no

POPULATION:
exploratoryautostop = no
smoothingautostop = no
smoothingiterations = 20
exploratoryiterations = 50
variability = FirstStage