CLS method
From: "John T. Mondick" jmond001@umaryland.edu
Subject:[NMusers] CLS method
Date: Tue, July 13, 2004 3:15 pm
Hi everyone,
I am attempting to run the Constrained Longitudinal Splines method submitted
by Park, Veroota and Sheiner in the NONMEM repository. When I try to run the
single dose example (example 1), I get some peculiar errors and the output
files contain no predictions or parameters. I have attached the error
message and the control stream for reference. The control stream is taken
directly from the example.
I am running NONMEM under Windows XP, Visual Fortran v6.6. Does anyone have
any suggestions on how to remedy this?
Thanks,
John Mondick
WARNINGS AND ERRORS (IF ANY) FOR PROBLEM 1
(WARNING 2) NM-TRAN INFERS THAT THE DATA ARE POPULATION.
(WARNING 31) $OMEGA INCLUDES A NON-FIXED INITIAL ESTIMATE CORRESPONDING TO
AN ETA THAT IS NOT USED IN ABBREVIATED CODE.
(WARNING 40) $THETA INCLUDES A NON-FIXED INITIAL ESTIMATE CORRESPONDING TO
A THETA THAT IS NOT USED IN ABBREVIATED CODE.
1 file(s) copied.
Starting nonmem execution ...
MONITORING OF SEARCH:
0ITERATION NO.: 0 OBJECTIVE VALUE: 0.1392E+05 NO. OF FUNC.
EVALS.: 9
CUMULATIVE NO. OF FUNC. EVALS.: 9
PARAMETER: 0.1000E+00 0.1000E+00 0.1000E+00 0.1000E+00 0.1000E+00
0.1000E+00 0.1000E+00 0.1000E+00 0.1000E+00
GRADIENT: 0.1962E+05 0.3795E+05 -0.2637E+06 0.3315E+05 0.1750E+05 -
0.7580E+05 -0.9156E+03 -0.3095E+05 -0.1722E+06
0ITERATION NO.: 51 OBJECTIVE VALUE: 0.2352E+03 NO. OF FUNC.
EVALS.: 0
CUMULATIVE NO. OF FUNC. EVALS.: 554
PARAMETER: -0.2095E-01 0.2578E-01 0.2042E+00 0.3728E-01 0.3911E-05
0.1286E+00 -0.8732E-01 0.1248E-04 0.5427E+00
GRADIENT: 0.5720E-01 0.1689E-02 0.3732E-02 0.5587E-01 0.2354E-01
0.2148E-02 0.897 4E-03 0.3023E-02 -0.1217E-01
forrtl: error (65): floating invalid
Image PC Routine Line Source
nonmem.exe 00403BC2 Unknown Unknown Unknown
nonmem.exe 00405118 Unknown Unknown Unknown
nonmem.exe 004182F3 Unknown Unknown Unknown
nonmem.exe 0040570D Unknown Unknown Unknown
nonmem.exe 004B7189 Unknown Unknown Unknown
nonmem.exe 004A84F4 Unknown Unknown Unknown
kernel32.dll 77E814C7 Unknown Unknown Unknown
1 file(s) copied.
-------------------------------------------------------------------------
$PROBLEM single dose
$INPUT ID TIME DV AMT MDV
$DATA ex1.dat
$SUB OTHER=pack.F OTHER=single.F
$PRED
"FIRST
" double precision ggg(10),col(4),desp(20)
" character fname(4)*12
NB2=3
ORDER=4
FRAC=.05
NDESP=5
STEP=.1
TMAX=36
" DESP(1)=.5
" DESP(2)=2.26
" DESP(3)=3.58
" DESP(4)=5.91
" DESP(5)=10
" COL(1)=1
" COL(2)=2
" COL(3)=3
" COL(4)=5
" FNAME(1)='pred.rep'
" FNAME(2)='ipred.rep'
" FNAME(3)='mipk.rep'
" FNAME(4)='para.rep'
" call single(icall,newind,theta,datrec,frac,nb2,order,time,
" x eta,col,desp,ndesp,step,tmax,fname,cls,ggg,id,0,icro)
" do j=1,nb2-1
" g(j,1)=ggg(j)
" enddo
Y=CLS*(1+EPS(1))+EPS(2)
JUNK=ETA(1)+EPS(1)
$OMEGA .0625 .0625
$SIGMA .0625 .09
$EST MAX=9999 POSTHOC
$THETA
(0 FIX)
(0 2 INF)
(0 2 INF)
(0 2 INF)
(0 2 INF)
(0 2 INF)
$TABLE ID TIME DV PRED
ONEHE NOPRINT FILE=examp.TXT
_______________________________________________________