CLS method

1 messages 1 people Latest: Jul 13, 2004

CLS method

From: John T. Mondick Date: July 13, 2004 technical
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 _______________________________________________________