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--------
Quoted reply history
De: [EMAIL PROTECTED] de la part de Bill Gillespie
Date: sam. 22/09/2007 13:16
À: Wang, Yaning
Cc: [email protected]
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: [EMAIL PROTECTED]
>
> "The contents of this message are mine personally and do not necessarily
> reflect any position of the Government or the Food and Drug
> Administration."
>
>
>