NONMEM Users Network Archive

Hosted by Cognigen

AW: Meta-analysis with Nonmem

From: Garmann, Dirk <Dirk.Garmann>
Date: Thu, 25 Feb 2010 11:10:25 +0100

Dear Leonid, Andreas,
I would suggest to simulate each study with N subjects based on the known=
 distribution, mean value (DV), SD,N. This would consider the number of s=
ubjects in each study as well as the study specific variability.
The datasets can be merged and used for a combined NONMEM analysis

If the number of subjects in a study is small, it might be useful to repe=
at the simulation and produce different combined datasets

The NONMEM results based on all combined datasets can be compared/analyze=

Best regards

Dirk Garmann, PhD
Clinical Scientific Expert /Pharmacokineticist
Merz Pharmaceuticals
Eckenheimer Landstrasse 100
60318 Frankfurt
Phone +49 (69) 1503 720

-----Ursprüngliche Nachricht-----
Von: owner-nmusers
m Auftrag von Leonid Gibiansky
Gesendet: Wednesday, February 24, 2010 4:45 PM
An: alindauer-research
Cc: nmusers
Betreff: Re: [NMusers] Meta-analysis with Nonmem

I think if you have only study-level data, you can treat each study as a
"subject", build meta-population dose-PD model in the usual way,

IPRED=Emax (or some other) function of DOSE

and modify the error part as follows: For each data point, you know the
mean value (DV), SD, and N. From SD and N you can get an estimate of
standard error of the mean as SE=SD/sqrt(N) (computed for each data
point and included into the data file). Then your error model would inclu=


I am not sure whether you need to fix SIGMA
1 FIXED ; for EPS(1)

thus assuming that all error in your model comes from the "assay", or
estimate it thus allowing for unexplained model misspecification and
extra error. I would try both ways to see the difference.

This is the simplest version that can be further improved (? or at
least, made more complicated) by adding the study effect on error (thus
accounting for possible differences in study populations;
Y=IPRED+SE*EXP(ETA(1))*EPS(1)), etc.


Leonid Gibiansky, Ph.D.
President, QuantPharm LLC
e-mail: LGibiansky at
tel: (301) 767 5566

> I wanted to investigate the dose-response relationship (Emax model) of =
> drug with NONMEM, based on data from literature (i.e. a meta-analysis).=
> However, I am not quite sure how to deal with the different levels of
> random effects. Suppose I have 10 studies of different size where
> different doses were given and the response is presented as average
> change of a biomarker +/- standard deviation for each dose level. How
> would I incorporate the standard deviation of the biomarker measurement=
> reported in each study for each dose level and how would I account for
> the different number of patients in the study?
> I would greatly appreciate your help, maybe with a NM-code snippet or
> reference to a paper where something similar has been done.
> Thanks in advance, Andreas.
> Ferrer
> Andreas Lindauer
> Pharmacokineticist
> Pharmacokinetics and Metabolism
> R&D Center. Ferrer Internacional S.A.
> Juan de Sada 32, 08028 Barcelona
> alindauer-research
> Recicla ¿Necesita imprimir este mensaje? Protejamos el medio
> ambiente. Li cal imprimir aquest missatge? Protegim el medi ambient.=
> Do you need to print this message? Let's protect the environment.
> Este mensaje, y en su caso, cualquier fichero anexo al mismo, puede
> contener información confidencial, siendo para uso exclusivo del
> destinatario, quedando prohibida su divulgación, copia o distribuci=
ón a
> terceros sin la autorización expresa del remitente. Si Vd. ha recibid=
> este mensaje erróneamente, se ruega lo notifique al remitente y proce=
> a su borrado. Gracias por su colaboración.
> This message and its annexed files may contain confidential information=
> which is exclusively for the use of the addressee. It is strictly
> forbidden to distribute copies to third parties without the explicit
> permission of the sender. If you receive this message by mistake, pleas=
> notify it to the sender and make sure to delete it. Thank you for your
> kind cooperation.

Merz Pharmaceuticals GmbH, Frankfurt am Main
Amtsgericht Frankfurt am Main, HRB 53808
Geschäftsführung: Dr. Martin Zügel (Vors.), Dr. Alexander Gebauer,
Dr. Karsten Schlemm, Dr. Eugen Wilbert

Die vorgenannten Angaben der E-Mail haben grundsätzlich nur informative=
n Charakter. Dies ist kein Anerkenntnis, dass es sich beim Inhalt dieser =
E-Mail um eine rechtsverbindliche Erklärung der entsprechenden Gesellsc=
haft der Merz Gruppe handelt, es sei denn dies ist ausdrücklich als sol=
ches formuliert. Erklärungen, die eine Gesellschaft der Merz Gruppe ver=
pflichten sollen, bedürfen jeweils der Unterschrift durch zwei zeichnun=
gsberechtigte Personen dieser Gesellschaft.
Diese E-Mail enthält vertrauliche und/oder rechtlich geschützte Info=
rmationen. Wenn Sie nicht der richtige Adressat sind oder diese E-Mail ir=
rtümlich erhalten haben, informieren Sie bitte sofort den Absender und=
 vernichten Sie diese E-Mail. Das unerlaubte Kopieren sowie die unbefugt=
e Weitergabe dieser E-Mail ist nicht gestattet.

This e-mail may contain confidential and/or privileged information. If y=
ou are not the intended recipient (or have received this e-mail in error=
) please notify the sender immediately and destroy this email. Any unaut=
horised copying, disclosure or distribution of the material in this e-ma=
il is strictly forbidden. (MRZ/2010)
Received on Thu Feb 25 2010 - 05:10:25 EST

The NONMEM Users Network is maintained by ICON plc. Requests to subscribe to the network should be sent to:

Once subscribed, you may contribute to the discussion by emailing: