M3 method for LLOQ

$ERROR
; Dans la bdd, les valeurs de DV en-deçà de la LOQ doivent être à « 0 »
; BQL: colonne pour identifier les lignes avec une DV en-deça de la LOQ (Below Quantification Limit)
; BQL=1 si DV inférieur à la BQL, sinon BQL=0
; MDV = 0 pour les lignes BQL

; Faire un run initial avec MDV = 1 pour les lignes BQL
; les lignes seront ignorées à l’estimation et cela permettra de
; vérifier le code avec FOCEI INTERACTION.
; les IPRED des points correspondants à des BQL doivent être bas.

IPRED=F
LOQ= xxx ; concentration de la LOQ
P_RUV=F*THETA(1)
A_RUV=THETA(2)
SD=SQRT(P_RUV*P_RUV + A_RUV*A_RUV)
IF (BLQ.EQ.0) THEN ; FLG=0, MDV=0, non-LLOQ values
F_FLAG=0
Y = IPRED + SD*ERR(1)
ENDIF
IF (BLQ.EQ.1) THEN ; FLG=1, MDV=0, LLOQ values
; PRED sera la proba que la concentration soit inférieur à la LOQ
; IPRED sera la concentration prédite individuelle
F_FLAG=1
DUM = (LOQ-IPRED)/SD
CUMD=PHI(DUM)
Y = CUMD
;MDVRES = 1 ; to be tested
ENDIF

$SIGMA 1 FIXED

$EST METHOD=1 INTER LAPLACIAN NSIG=4 SIGL=12 MAXVAL=9999 NOBT NOTBT NOABORT

$COV UNCONDITIONAL MATRIX=S PRINT=E