Reading NONMEM output into S-Plus
From: R.Port@DKFZ-Heidelberg.de
Subject: Reading NONMEM output into S-Plus
Date: Mon, 5 Feb 2001 19:37:55 +0100
Hi Ganesh:
Re: Reading NONMEM output into S-Plus
I suggest you collecting the desired information in the NONMEM output *table*, then reading it into S-Plus using the read.table command. E.g. if the name of the NONMEM output table is 'NONMEM.output.table' and its first few lines look like this:
TABLE NO. 1 DATA RECS. 1 THROUGH 900
ID PS PAT MALE AGE0 BW
ADM TIME DRCH WDR ETA1 ETA2 ETA3
beta Emax ED50 tau mep DV PRED
RES WRES
1.0000E+00 1.0000E+00 1.0000E+00 1.0000E+00 1.0000E+01 2.5900E+01
1.0000E+00 0.0000E+00 2.0000E+00 7.5000E+03 -5.5092E-01 -2.0011E-01
5.0293E-01 6.6288E-06 0.0000E+00 1.0000E+00 8.4320E+01 6.9725E+00
7.4000E+00 7.4000E+00 0.0000E+00 3.6832E-01
1.0000E+00 1.0000E+00 1.0000E+00 1.0000E+00 1.0000E+01 2.6100E+01
1.0000E+00 2.0000E+00 0.0000E+00 7.5000E+03 -5.5092E-01 -2.0011E-01
5.0293E-01 6.6288E-06 0.0000E+00 1.0000E+00 8.4320E+01 7.0719E+00
8.6000E+00 7.5725E+00 1.0275E+00 1.6684E+00
1.0000E+00 1.0000E+00 1.0000E+00 1.0000E+00 1.0000E+01 2.6000E+01
1.0000E+00 4.0000E+00 0.0000E+00 7.5000E+03 -5.5092E-01 -2.0011E-01
5.0293E-01 6.6288E-06 0.0000E+00 1.0000E+00 8.4320E+01 7.1714E+00
6.5000E+00 7.7450E+00 -1.2450E+00 -9.1380E-01
then, in S-Plus, you would use:
> NONMEM.output.table_read.table('NONMEM.output.table',skip=1,header=T)
Don't forget to use the ONEHEADER option in $TABLE .
Good luck! Ruedi
----------------------------------------------------------------------------
Dr. R.E. Port, German Cancer Research Center, D-0200
P.O. Box 10 19 49, D-69009 Heidelberg, Germany
phone: +49-6221 42-3385
-3347
fax: -3382
e-mail: r.port@dkfz.de