Dear Members of Certara Forum
I have an issue during the development of popPK modeling. I tried to fix the parameters that have large ETA shrinkage to zero. However when I fixed it to zero, the model excution always failed with following error message : out of range final -LL value : 1.0E+70
In NONMEM, I can fix the ETA like this example code :
$OMEGA
0.04
0.04
0 FIX
How can I handle this issue?
And, is there the way to check ETA bar and its p-value in NLME? (I have checked it in manual)
I have an issue during the development of popPK modeling. I tried to fix the parateres that have large ETA shrinkage to zero. However when I fixed it to zero, the model excution always failed with following error message like this : out of range final -LL value : 1.0E+70
In NONMEM, I can fix the ETA like this example code :
$OMEGA
0.04
0.04
0 FIX
Dear Go-Wun, Instead of trying to fix to zero, I wold first try to remove random effects from the parameter with the worst shrinkage, as I have done below for CL2
alternatively if I understand you correctly and you still want to freeze a parameter it’s quite simpe either in the UI or the code; below I froze (fixed) V2
fixef(tvV = c(, 5.35391, ))
fixef(tvV2(freeze) = c(, 10.9549, )) # commenting that v2 is frozen
fixef(tvCl = c(, 2.249, ))
fixef(tvCl2 = c(, 2.56732, ))
Simon.
Dear Go Wun Choi
Can you provide a reproducible example on our end that reproduces this -eLL error ?
There is ways to fix a random effect to zero or enable disable it.
In Phoenix PML each parameter has a name
so you can just comment out the parameter that you don’t need without afecting the other part of the model there is no sequential numbering of OMEGA so you will not have to edit your other parts of code
CL =tvCL #*exp (nCL)
#ranef(ncl=c(0.1))
some alrogrithms like SAEM or QRPEM do not like 0 variance fixing you might need to fix to a small value instead
if you have all of your n on one block or line you might need to edit the ranef
ranef( nCL,nV)
to
ranef(nv)
Samer
