NONMEM Users Network Archive

Hosted by Cognigen

RE: [NMusers] Problem of STS in NONMEM

From: Mills, Richard <Richard.Mills_at_iconplc.com>
Date: Tue, 22 Dec 2015 07:54:18 +0000

Hi Anyue,

You need to amend MAXEVAL=0 in $EST (I suggest MAXEVAL=9999) in order to allow estimation.

Kind regards,
Richard
Richard Mills PhD
Senior Scientist, PKPDM&S
 [http://emailsignature/stationary/Icon_Development_Solutions.gif]

ICON plc made the following annotations.
------------------------------------------------------------------------------
This e-mail transmission may contain confidential or legally privileged information
that is intended only for the individual or entity named in the e-mail address. If you
are not the intended recipient, you are hereby notified that any disclosure, copying,
distribution, or reliance upon the contents of this e-mail is strictly prohibited. If
you have received this e-mail transmission in error, please reply to the sender, so that
ICON plc can arrange for proper delivery, and then please delete the message.
Thank You,
ICON plc
South County Business Park
Leopardstown
Dublin 18
Ireland
Registered number: 145835



From: owner-nmusers_at_globomaxnm.com [mailto:owner-nmusers_at_globomaxnm.com] On Behalf Of Anyue Yin
Sent: 22 December 2015 06:36
To: nmusers_at_globomaxnm.com
Subject: [NMusers] Problem of STS in NONMEM

Dear all,

I got a problem when using NONMEM with STS method. My aim is to get individual parameter estimates. I would like to let NONMEM estimate individual parameter one subject by one subject. My control stream is partially listed below, which is edited according to this thread http://www.cognigencorp.com/nonmem/nmo/topic035.html<http://www.cognigencorp.com/nonmem/nmo/topic035.html>

$SUBROUTINES ADVAN1 TRANS2

$PK
CL = THETA(1) * EXP(ETA(1))
V = THETA(2) * EXP(ETA(2))
S1 = V

$ERROR
IPRED = F
    W = F
    Y = IPRED*(1+EXP(ETA(3))*EPS(1))
 IRES = DV-IPRED
IWRES = IRES/W

$THETA
(10,20,30) ; CL
(10,80,100) ; V

$OMEGA
100 ; IIV CL
100 ; IIV V
100 ; IIV SIGMA

$SIGMA
1 FIXED ; PRO
$EST METHOD=1 INTER MAXEVAL=0 NOABORT SIG=3 PRINT=1 POSTHOC
$COV
$TABLE ID TIME AMT CMT DV PRED IPRED IWRES CL V ETA(1) ETA(2) ETA(3) ONEHEADER NOPRINT FILE=sdtab101.tab

The model can run successfully, and I got individual parameter estimates from output file. I found that for each subject the individual parameter estimate changed if I change the initial value of THETA. It seems that NONMEM regarded the initial value fixed (just like a fixed typical value). But I expected NONMEM to estimate individual parameters based on individual data rather than the initial value. Is there anything wrong with the code? Thanks in advance.

Anyue




image001.jpg
(image/jpeg attachment: image001.jpg)

Received on Tue Dec 22 2015 - 02:54:18 EST

The NONMEM Users Network is maintained by ICON plc. Requests to subscribe to the network should be sent to: nmusers-request_at_iconplc.com. Once subscribed, you may contribute to the discussion by emailing: nmusers_at_globomaxnm.com.