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
_______________________________________________________