RE : [NMusers] Probit function in NONMEM

Hello Yaning,

Sample code of using PHI (I did not use it yet in the current release of =
NONMEM 6 but previously we had to use it with verbatim code):

$PROB SINGLE-SUBJECT DICH DV SIM
$DATA data1
$INPUT DV
$PRED
IF (ICALL.EQ.4) THEN
" P=PHI(THETA(1))
  CALL RANDOM (2,R)
  IF (R.LE.P) THEN
    DV=1
  ELSE
    DV=0
  ENDIF
  RETURN
ENDIF
...
$THETA 1
$OMEGA 1 FIX
$SIM (1) (54543 UNIFORM)

$TABLE FILE=tablemode ONEHEADER

Hope this helps,

Samer Mouksassi
-------- Message d'origine--------
De: owner-nmusers_at_globomaxnm.com de la part de Bill Gillespie
Date: sam. 22/09/2007 13:16
À: Wang, Yaning
Cc: nmusers_at_globomaxnm.com
Objet : Re: [NMusers] Probit function in NONMEM
 
Hi Yaning,

See the PHI function in the NONMEM help:

> ___________________________________________________________________
> | |
> | PHI |
> |_________________________________________________________________|
>
> MEANING: PHI function
> CONTEXT: NONMEM utility routine
>
> USAGE:
> DOUBLE PRECISION PHI,X,CUMDIS
> C (If single precision is to be used):
> C REAL PHI,X
> CUMDIS= PHI(X)
>
> DISCUSSION:
> Gives the value of the cumulative distribution function.
>
> For X < 0, this is the integral from -INF to X of N(0,1).
> For X > 0, this is the integral from X TO INF of N(0,1).
> For X=0, this is .5.
>
> Input argument:
>
> X
>
> REFERENCES: None.
>
Regards,
Bill

Wang, Yaning wrote:
> Dear all:
> Is there a Probit function in NONMEM we can directly use?
> Thanks
> Yaning
>
>
> Yaning Wang, Ph.D.
> Team Leader, Pharmacometrics
> Office of Clinical Pharmacology
> Office of Translational Science
> Center for Drug Evaluation and Research
> U.S. Food and Drug Administration
> Phone: 301-796-1624
> Email: yaning.wang_at_fda.hhs.gov
>
> "The contents of this message are mine personally and do not =
necessarily
> reflect any position of the Government or the Food and Drug
> Administration."
>
>
>

Received on Sat Sep 22 2007 - 17:24:40 EDT