; MONOLIX project generated with IQRtools

; ==PROJECT HEADER START===================================================
; COMMENT             = 'First covariate model'
; TOOL                = 'MONOLIX'
; TOOLVERSION         = 'MLX2023R1'
; DATA                = '../../Ex12/data.csv'
; DOSINGTYPES         = 'BOLUS,INFUSION'
; TK0NAMES            = 'NA,NA'
; COVNAMES            = ''
; CATNAMES            = ''
; CATCATEGORIES       = ''
; REGRESSIONNAMES     = ''
; OUTPUTS             = 'Cc'
; ERRORMODELS         = 'absrel'
; PARAMNAMES          = 'Fabs1,ka,CL,Vc,Q,Vp'
; PARAMTRANS          = 'exp(phi),exp(phi),exp(phi),exp(phi),exp(phi),exp(phi)'
; PARAMINVTRANS       = 'log(psi),log(psi),log(psi),log(psi),log(psi),log(psi)'
; COVARIATENAMES      = ''
; COVARIATESUSED      = ''
; BETACOVNAMES        = ''
; BETACOVTRANS        = ''
; BETACATNAMES        = ''
; BETACATREFERENCE    = ''
; NROBSERVATIONS      = '1991'
; ==PROJECT HEADER END=====================================================

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

[FILEINFO]
file = '../../Ex12/data.csv'
delimiter = comma
header = {IXGDF, USUBJID, ID, TIME, TIMEPOS, TAD, TIMEUNIT, YTYPE, NAME, DV, UNIT, CENS, MDV, EVID, AMT, ADM, ROUTE, TINF, RATE}

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

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

[INDIVIDUAL]
input = {Fabs1_pop, ka_pop, CL_pop, Vc_pop, Q_pop, Vp_pop, omega_Fabs1, omega_ka, omega_CL, omega_Vc, omega_Q, omega_Vp}

DEFINITION:
Fabs1 = {distribution=logNormal, typical=Fabs1_pop, no-variability}
ka = {distribution=logNormal, typical=ka_pop, sd=omega_ka}
CL = {distribution=logNormal, typical=CL_pop, sd=omega_CL}
Vc = {distribution=logNormal, typical=Vc_pop, sd=omega_Vc}
Q = {distribution=logNormal, typical=Q_pop, no-variability}
Vp = {distribution=logNormal, typical=Vp_pop, sd=omega_Vp}

[LONGITUDINAL]
input = {a, b}

file = './project_model.txt'

DEFINITION:
y = {distribution=normal, prediction=Cc ,errorModel=combined1(a, b)}

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

; =============================================
<PARAMETER>
; =============================================
Fabs1_pop = {value=0.5, method=MLE}
ka_pop = {value=0.2, method=MLE}
CL_pop = {value=0.3, method=MLE}
Vc_pop = {value=2.5, method=MLE}
Q_pop = {value=0.8, method=MLE}
Vp_pop = {value=4, method=MLE}
omega_ka = {value=0.5, method=MLE}
omega_CL = {value=0.5, method=MLE}
omega_Vc = {value=0.5, method=MLE}
omega_Vp = {value=0.5, method=MLE}
a = {value=1, method=MLE}
b = {value=1, method=MLE}

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

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

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

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