MU referencing

Utile pour les méthode EM (inutile pour FOCE).
– MU_x et ETA_x doivent avoir le même index.
– un THETA référencé par un MU, ne peut pas être réutilisé plus bas dans le code. Pour contourner cette limitation, il faut « sortir »la cov du référencement (voir plus bas).
– un paramètre sans ETA ne peut pas être référencé MU. Dans ce cas, il est conseillé de fixer ETA à une valeur très faible.
– les covariables qui varient avec le temps, les doses ou les observations ne peuvent pas être référencées (paramètre bio).
– tous les paramètres ne doivent pas forcement être référencé. Dans ce cas, il faut préciser lesquelles avec l’option MUM de $EST.
ex. $EST… MUM=MMMN seule le 4ème paramètre n’est pas référencé (N).
– avec les modèles MIXTURE, l’option MUM doit être utilisé.
– pas de MU_ dans le block $ERROR
– pas d’utilisation de IF THEN avec MU_.
– dès le début du $PK (avant code verbatim)

Exemple:
CL = THETA(1) * EXP(ETA(1)) se code:
1ère possib.
MU_1 = THETA(1)
CL = EXP(MU_1 + ETA(1))
NB: avec ce code THETA(1) doit être back-transformé pour être interprétable.
2ème possib. avec même valeur initiale de theta
MU_1 = LOG(THETA(1))
CL = EXP(MU_1 + ETA(1))

Exemple
CL = THETA(1)*(WT/80)**THETA(2) * EXP(ETA(1)) se code:
1ère possib.
MU_1 = LOG(THETA(1)*(WT/80)**THETA(2))
CL = EXP(MU_1 + ETA(1)
2ème possib. avec même valeur ini. de theta
TVCL = THETA(1)*(WT/80)**THETA(2)
MU_1 = LOG(TVCL)
CL = EXP(MU_1 + ETA(1))
3ème possib. plus rapide
MU_1 = LOG(THETA(1)) + THETA(2)*LOG(WT/80))
CL = EXP(MU_1 + ETA(1))
4ème possib. en « sortant » la cov.
MU_2 = LOG(THETA(2))
V = (WT/80) * EXP(MU_2 + ETA(2))

Exemple
CL = THETA(4) + ETA(2) se code:
MU_2=THETA(4)
CL=MU_2 + ETA(2)