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 firstname.lastname@example.org ____________________________ Clinical Pharmacology - Pharmacometrics Novo Nordisk A/S, Krogshoejvej 53A 2880 Bagsvaerd, Denmark ******* 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:email@example.com 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!