Epsilon shrinkage and IWRES

Dear all, Can someone enlighten me on a NONMEM implementation? Since NONMEM 7, epsilon shrinkage is reported in the output, and the manual states that it is calculated as: 100%*[1-SD(IWRES)] in accordance with "Karlsson MO and Savic RM. Diagnosing Model Diagnostics. Clinical Pharmacology and Therapeutics, 2007; 82(1): 17-20" as one would expect. However, as far as I know, the user needs to manually code the IWRES bit using the "Uppsala implementation": IPRED = F IRES = DV - IPRED W = THETA(3) IWRES = IRES/W Y = IPRED+W*EPS(1) SIGMA 1 FIXED So how does NONMEM arrive at epsilon shrinkage? And is there by any chance an (undocumented?) IWRES item that can be exported in a table file as well? Kind regards, Rik Rik Schoemaker, PhD Exprimo NV Tel: +31 (0)20 4416410 E-mail: [email protected] Web: www.exprimo.com This e-mail is confidential. It is also privileged or otherwise protected by work product immunity or other legal rules. The information is intended to be for use of the individual or entity named above. If you are not the intended recipient, please be aware that any disclosure, copying, distribution or use of the contents of this information is prohibited. You should therefore delete this message from your computer system. If you have received the message in error, please notify us by reply e-mail. The integrity and security of this message cannot be guaranteed on the Internet. Thank you for your co-operation. ---------------------------- This e-mail message has been scanned for Viruses by Norman Virus Control and it's Content by MailMarshal