RE: current sigma
Try the following examples, according to your need, in $PRED or $PK:
SIGVAL11=SIGMA(1,1)
OVAL22=OMEGA(2,2)
IF(NEWIND.eq.0) WRITE(*,*) SIGVAL11,OVAL22
...
$TABLE SIGVAL11 OVAL22
Robert J. Bauer, Ph.D.
Senior Director
Pharmacometrics R&D
ICON Early Phase
820 W. Diamond Avenue
Suite 100
Gaithersburg, MD 20878
Office: (215) 616-6428
Mobile: (925) 286-0769
[email protected]
www.iconplc.com
Quoted reply history
-----Original Message-----
From: [email protected] [mailto:[email protected]] On
Behalf Of [email protected]
Sent: Sunday, November 04, 2018 7:03 PM
To: [email protected]
Subject: [NMusers] current sigma
Hello,
I would like to know how to know current sigma values in a customized PRED
routine. In noncom.for file, it is documented that initial and final theta,
omega and sigma values are set in ROCM6, and current omega values are set in
ROCM22. However, I cannot find current sigma values.
We developed a new program which can restore numerous time-course data by
expanding the analysis of nonlinear mixed effect model. Detail of our
analysis, named SReFT, is described in recent Clin Pharm Ther
(PMID:29951994). Currently, SReFT is implemented in my own program, Napp,
but I am considering its implementation in NONMEM. For the implementation of
SReFT, it is necessary to know current sigma value in the PRED routine.
Akihiro Hisaka,
Chiba University, Japan.
email: [email protected]