From:"TKT (Thomas Klitgaard)"
Subject:[NMusers] Dumping of internally used parameter values to file
Date:Thu, February 28, 2002 4:49 pm
Title: Dumping of internally used parameter values to file
Dear nmusers,
Could any of you inform me, if it is possible to have the parameter
estimates produced by NONMEM written to a file, with all decimals
used (internally) by the program, instead of "just" the rounded
off values in the output-file?
Thanks in advance
Thomas Klitgaard
Pharmacometrician
Direct +45 44424960
Fax +45 44421940
tkt@novonordisk.com
Dumping of internally used parameter values to file
From:Nick Holford
Subject:Re: [NMusers] Dumping of internally used parameter values to file
Date:Thu, February 28, 2002 6:44 pm
I have tried to follow the logic from OUTPT, PRNT3, USWLSM and finally PRNT. My best
guess is that you change this statement in PRNT.for:
DATA VAR/'1PE9','.2',' 9(','1H.)'/
This statement determines the format when the values of THETA etc are printed in the
output listing. However, my experiments have all failed. Any changes I have tried
e.g. '1PE12' instead of '1PE9' give me a string of *s in the output. Some advice is
needed from a FORTRAN guru.
Nick
--
Nick Holford, Divn Pharmacology & Clinical Pharmacology
University of Auckland, 85 Park Rd, Private Bag 92019, Auckland, New Zealand
email:n.holford@auckland.ac.nz tel:+64(9)373-7599x6730 fax:373-7556
http://www.health.auckland.ac.nz/pharmacology/staff/nholford/
From:Pascal Girard
Subject:RE: [NMusers] Dumping of internally used parameter values to file
Date:Fri, March 1, 2002 2:19 am
Dear Thomas,
Iwould not change the NONMEM FORTRAN code. The easiest way to have
parameters with higher precisions is to write your own INFN.f fortran
subroutine, which will replace the dummy one that comes with PREDDPP. Within
INFN you will have access to all final parameter estimates, placed in
various COMMON, which in turn you can write, with the desired format, to an
ASCII file.
Way to call INFN, depends whether you use PREDPP or not.
1) If using PREDPP just change the dummy INFN in $SUBROUTINE by:
$SUBROUTINE ADVANxx TRANSxx INFN=MYINFN.f
Then MYINFN subroutine will be automatically call at the end of each
problem.
2) If not using PREDPP, you will have to control the call to INFN. For
exemple:
$SUBROUTINE OTHER=MYINFN.f
IF (ICALL.EQ.0) CALL MYINFN( ...) ; e.g. initialization for a multiple
problem run
...
IF (ICALL.EQ.3) CALL MYINFN( ...) ; output parameters with desired
precisions at the end of a problem
You have lots of help in the nmhelp for this.
Best regards,
Pascal Girard, PHARSIGHT Corp.
Tel/Fax +33 (0)437 37 80 29
From:"Rombout, Dr. Ferdinand"
Subject:AW: [NMusers] Dumping of internally used parameter values to file
Date:Fri, March 1, 2002 2:46 am
The answer is simple but not easy to perform. You have to look for the
Fortran formatting statements in the file NMD and change those. I used that
in the past to change for example the output of the objective function
during the iteration process to see if it changes (which is difficult to see
if the OF is larger than 10000)and to change the digits outputted for THETA
in the final results. It is a bit trial and error to find the right
statements for the numbers you need.
But remember, according the license agreement from Globomax you have to pay
a lot if you change the code without permission!