2 replies to this topic

[Go-Wun Choi]

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

[Simon Davis]

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.

Attached Files

•  noRanEff_Cl2.jpg   105.1KB   1 downloads
•  freeze.jpg   51.38KB   0 downloads

[smouksassi1]

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