Re: Reading NONMEM output into S-Plus

From: "Niclas Jonsson" <Niclas.Jonsson@farmbio.uu.se> Subject: Re: Reading NONMEM output into S-Plus Date: Wed, 7 Feb 2001 08:37:58 +0000 (GMT) If you want to get hold of the parameter estimates in a NONMEM output file you could use teh read.lst function that comes with Xpose, e.g: > tmp <- read.lst("run13.lst") > tmp $term: [1] "MINIMIZATION SUCCESSFUL" [2] "NO. OF FUNCTION EVALUATIONS USED: 208" [3] "NO. OF SIG. DIGITS IN FINAL EST.: 3.9" $ofv: [1] 2116.521 $thetas: [1] 20.30000 84.90000 1.77000 0.00566 $omega: $omega[[1]]: [1] 0.136 $omega[[2]]: [1] 0.000 0.182 $omega[[3]]: [1] 0.00 0.00 1.39 $sigma: $sigma[[1]]: [1] 0.0238 $sethetas: [1] 1.07000 4.51000 0.39000 0.00317 $seomegas: $seomegas[[1]]: [1] 0.0297 $seomegas[[2]]: [1] 0.0000 0.0299 $seomegas[[3]]: [1] 0.000 0.000 0.663 $sesigmas: $sesigmas[[1]]: [1] 0.00469 In other words, you will get back a list with various components like ofv, thetas and so on. If you want to access, for example, the thetas from the tmp object created above, you would type: > tmp$thetas and that would give you a vector of the theta estimates, i.e. [1] 20.30000 84.90000 1.77000 0.00566 Best, Niclas -- Department of Pharmacy Uppsala University Box 580 SE-751 23 Uppsala Sweden Phone: +46 18 471 43 85 Fax: +46 18 471 40 03 Mobile: +46 70 485 61 98 E-mail: niclas.jonsson@biof.uu.se