RE: OMEGA HAS A NONZERO BLOCK
From:Pascal Girard
Subject:RE: [NMusers] OMEGA HAS A NONZERO BLOCK
Date:Wed, 2 Oct 2002 00:25:42 -0700
Hi Steve,
I often encountered the problem when I want to resample OMEGAs from MVN distribution using
a non-diagonal final OMEGA. I fixed it definitively by writing my own INFN subroutine which outputs
final OMEGA in an ASCII file, with more digits than the one provided by NONMEM output.
By the way this is really not rocket science. It would be much more efficient to have an
option in NONMEM that outputs all NONMEM estimates, with sufficient precisions, (Objective
functions, THETAs, OMEGAs and SIGMAs, SE, COV matrix ..., indexed by the PROB #) in an ASCII,
EXCEL compatible, file. Software like SAS statistical procedures or Splus functions offer this
functionality. This would be a kind of sophisticated Model Specification File, but in ASCII,
rather than binary coding. (Hope Stuart is listening -:)
All my best,
Pascal Girard
Tel/Fax +33 (0)437 37 80 29