Hi Ruben, I think I misread Tingjies original posting as taking ABS(ETA), whereas his initial attempt was actually ABS(1+ETA), which is less problematic. The latter would not bias simulations much if IIV is e.g. 30% CV, agreed.
However, as Tingjies is mainly interested in estimation, I believe that without the ABS-correction, no subject will have the EBE at ETA <= -1 for a parameter that could not be <=0. Unless possibly in a subject which is a) uninformative on that parameter and b) where the eta is also part of an omega-block - a scenario which seems unlikely to me, but may occur in theory. Implementing the ABS-korrection ETA=-1.2 would give the same solution (parameter value) as ETA=-0.8, but at a higher OFV for that subject. It seems to me, if negative parameter values are only a problem in the eta search for the EBE, whereas the EBE for individual parameters are always positive, then it should be more straightforward to use FOCE, with the addition e.g.: IF(PARA.LT.0.001) PARA=0.001 Probably, no subject will have such a low individual parameter value, when looking into the table output? If there are any such subjects I would look for errors in the data set and nonmem code (as outlined in my initial reply). The above concerns estimation. In simulation (unless %CV is low), we may get a fraction of subject with PARA=0.001, which may be an unreasonably low parameter value. Whether that is acceptable or not depends on the objectives and in this case there was no need for simulations even for model evaluation (?), so I will not elaborate further. Cheers Jakob
