; MONOLIX project generated with IQRtools

; ==PROJECT HEADER START===================================================
; COMMENT             = 'Covariance model'
; TOOL                = 'MONOLIX'
; TOOLVERSION         = 'MLX2023R1'
; DATA                = '../../Ex9/data.csv'
; DOSINGTYPES         = 'BOLUS'
; TK0NAMES            = 'NA'
; COVNAMES            = ''
; CATNAMES            = ''
; CATCATEGORIES       = ''
; REGRESSIONNAMES     = 'WT0,AGE0'
; OUTPUTS             = 'Cc'
; ERRORMODELS         = 'absrel'
; PARAMNAMES          = 'ka,CL,Vc,Q1,Vp1,beta1,beta2'
; PARAMTRANS          = 'exp(phi),exp(phi),exp(phi),exp(phi),exp(phi),(phi),(phi)'
; PARAMINVTRANS       = 'log(psi),log(psi),log(psi),log(psi),log(psi),(psi),(psi)'
; COVARIATENAMES      = ''
; COVARIATESUSED      = ''
; BETACOVNAMES        = ''
; BETACOVTRANS        = ''
; BETACATNAMES        = ''
; BETACATREFERENCE    = ''
; NROBSERVATIONS      = '1266'
; ==PROJECT HEADER END=====================================================

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

[FILEINFO]
file = '../../Ex9/data.csv'
delimiter = comma
header = {IXGDF, USUBJID, ID, TIME, TIMEPOS, TAD, TIMEUNIT, YTYPE, NAME, DV, UNIT, CENS, MDV, EVID, AMT, ADM, TINF, RATE, DOSE, WT0, AGE0, BMI0, HT0, SEX}

[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}
WT0 = {use=regressor}
AGE0 = {use=regressor}

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

[INDIVIDUAL]
input = {ka_pop, CL_pop, Vc_pop, Q1_pop, Vp1_pop, beta1_pop, beta2_pop, omega_ka, omega_CL, omega_Vc, omega_Q1, omega_Vp1, omega_beta1, omega_beta2}

DEFINITION:
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}
Q1 = {distribution=logNormal, typical=Q1_pop, sd=omega_Q1}
Vp1 = {distribution=logNormal, typical=Vp1_pop, sd=omega_Vp1}
beta1 = {distribution=normal, typical=beta1_pop, sd=omega_beta1}
beta2 = {distribution=normal, typical=beta2_pop, sd=omega_beta2}

[LONGITUDINAL]
input = {a, b}

file = './project_model.txt'

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

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

; =============================================
<PARAMETER>
; =============================================
ka_pop = {value=0.2, method=MLE}
CL_pop = {value=20, method=MLE}
Vc_pop = {value=20, method=MLE}
Q1_pop = {value=20, method=MLE}
Vp1_pop = {value=2000, method=MLE}
beta1_pop = {value=0.1, method=MLE}
beta2_pop = {value=0.1, method=MLE}
omega_ka = {value=0.5, method=MLE}
omega_CL = {value=0.5, method=MLE}
omega_Vc = {value=0.5, method=MLE}
omega_Q1 = {value=0.5, method=MLE}
omega_Vp1 = {value=0.1, method=MLE}
omega_beta1 = {value=0.001, method=FIXED}
omega_beta2 = {value=0.001, method=FIXED}
a = {value=1, method=MLE}
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 = 1
autochains = no

POPULATION:
exploratoryautostop = no
smoothingautostop = no
smoothingiterations = 200
exploratoryiterations = 500
variability = FirstStage