Logistic Mixed models with NONMEM

From: lewis@c255.ucsf.EDU (LSheiner) Subject: Logistic Mixed models with NONMEM Date: 23 Jan 1998 19:44:58 -0500 To reply to Doug's comments - I'm not saying there wouldn't be problems, but Doug's concern is not entirely justified ... THETA(3) multiplies ETA(1), so that, equivalently, GAMMA=THETA(1)*EXP(ETA(1))**THETA(3), or GAMMA = a*RV**b, Where RV = log-normally distributed random variable. The two constants, a and b, are not, in principle, unidentifiable. In fact, the model I used arose in a particular problem because the more general model (below) showed that ETA(1) and ETA(2) had a correlation of 1 in the particular data set I was analysing. Here is the more general case (which I will now put in the repository to avoid further confusion): $PRED ; NEXT LINE IS VERBATIM CODE " CALL GETETA(ETA) GAMMA=THETA(1)*EXP(ETA(1)) IF(GAMMA.GT.10) EXIT D50 = THETA(2)*EXP(ETA(2)) IF(D50.LT..00001) EXIT LDOS= LOG(DOSE) A = GAMMA*LDOS - GAMMA*LOG(D50) B = EXP(A) P = B/(1+B) Y = P**DV*(1-P)**(1-DV) ..... $OMEGA BLOCK (2) .2 .1 .2 > > 2. I've created a very simple PRED record (also below) > for calculation of a global mean in a mixed logistic model > (ie, the covariate DOSE is no longer present). Should it work? > $SUBROUTINES CONTR=contr CCONTR=ccontr > $PRED > ; NEXT LINE IS VERBATIM CODE > " CALL GETETA(ETA) > A = THETA(1)+ETA(1) > B = EXP(A) > P = B/(1+B) > Y = P**DV*(1-P)**(1-DV) > Looks fine to me. > > 3. The documentation for the PRED record from the NONMEM > user's archives (copied below) says that it > implements a 'sigmoid emax model' . > I know little pharmacology so I don't know > what that model is. How can I find out (Online > resources, if any, would be especially helpful)? > See the README.009 file in the same directory. Also, you can read about the various types of PD models, for example (and among many other places) in: Black JW, Leff P: Operational models of pharmacological antagonism. Proc Roy Soc B, 220:141-162, 1983, or Holford NHG, Sheiner LB: (a) Pharmacokinetic and Pharmacodynamic modelling in vivo. CRC-Critical Reviews in Bioengineering, 5:273, 1981; (b) Understanding the Dose-Effect Relationship: Clinical application of pharmacokinetic-pharmacodynamic models. Clinical Pharmacokinetics 6: 429, 1981; (c) Kinetics of pharmacologic Response. Pharmac. Ther. 16:143, 1982, and in Rowland, M & Tucker, G (eds), Pharmacokinetics: Theory and Methodology, Pergammon Press, Oxford, 1986.