Splines

3 messages 2 people Latest: Jan 17, 2019

Splines

From: Andre Jackson Date: January 16, 2019 technical
I looked at a readme text posted on line (attached) that originated from Icon some time ago: https://www.google.com/url?sa=t&rct=j&q=&esrc=s&source=web&cd=6&cad=rja&uac t=8&ved=2ahUKEwiQ4ozXuvLfAhWCZd8KHVW2DRIQFjAFegQIARAB&url=https%3A%2F%2Fnonm em.iconplc.com%2Fnonmem%2Fconstrained%2520longitudinal%2520spline%2Freadme.t xt&usg=AOvVaw07-r6tKLZsz_5GMOTYBMuM readme.txt ... https://www.google.com/url?sa=t&rct=j&q=&esrc=s&source=web&cd=6&cad=rja&uac t=8&ved=2ahUKEwiQ4ozXuvLfAhWCZd8KHVW2DRIQFjAFegQIARAB&url=https%3A%2F%2Fnonm em.iconplc.com%2Fnonmem%2Fconstrained%2520longitudinal%2520spline%2Freadme.t xt&usg=AOvVaw07-r6tKLZsz_5GMOTYBMuM https://www.google.com/url?sa=t&rct=j&q=&esrc=s&source=web&cd=6&cad=rja&uac t=8&ved=2ahUKEwiQ4ozXuvLfAhWCZd8KHVW2DRIQFjAFegQIARAB&url=https%3A%2F%2Fnonm em.iconplc.com%2Fnonmem%2Fconstrained%2520longitudinal%2520spline%2Freadme.t xt&usg=AOvVaw07-r6tKLZsz_5GMOTYBMuM https://nonmem.iconplc.com/nonmem/constrained%20longitudinal%20spline/readme .txt 1. https://webcache.googleusercontent.com/search?q=cache:xWJ5Iy0lqgsJ: https:// nonmem.iconplc.com/nonmem/constrained%2520longitudinal%2520spline/readme.txt +&cd=6&hl=en&ct=clnk&gl=us Cached Nov 18, 1997 - The longitudinal splines are constrained to exhibit (i) decreasing tail (ii) ... random effects denoted ETAj in NONMEM, where ETAj is distributed ... You've visited this page 5 times. Last visit: 1/15/19 When I tried to run the code for a single dose I got an error message. The control and data (1 subject) are attached. The error message was that the characters in error were: Pack $268 SUBROUTINES EEROR WHILE OPENING OR COPYING USED SUPPLIED SUBROUTINE Control $PROBLEM Fits the single-dose PK data with the CLS method $DATA ..\SDSPLINE_DAT.CSV IGNORE=I $INPUT ID TIME DV AMT MDV $SUB OTHER=PACK OTHER=single ; on SUN/UNIX $SUB OTHER=PACK.f OTHER=single.f ; on PC $PRED "FIRST " double precision ggg(10),col(4),desp(20) " character fname(4)*12 NB2=3 ORDER=4 NDESP=5 FRAC=.05 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) ; any error model ; Enter nonzero values of H(i,j+1), j=1,nb2-1, when the INTERACTION ; option is used in $EST. ; ;" do j=1,nb2-1 ;" h(1,j+1)=ggg(j) ;" enddo JUNK=ETA(1)+EPS(1) ; NMTRAN requirement $OMEGA .0625 .0625 ; diagonal ; dimension = NB2-1 due to the end-point constraint $SIGMA .0625 .09 $EST MAX=9999 POSTHOC $THETA (0 FIXED) ; THETA(1) = 0 due to the end-point constraint (0 2 INF) ; THETA(2) (0 2 INF) ; THETA(3) (0 2 INF) ; THETA(4) (0 2 INF) ; THETA(5) (0 2 INF) ; THETA(6) ;$SIMULATION(2792293 NEW) ONLYSIM SUBPROBLEMS=100 ;$ESTIMATION MAX=5000,PRINT=5 ;$COVARIANCE $TABLE TIME ID Y DV ONEHEADER NOPRINT FILE=MODEL1.fit ID TIME DV AMT MDV 1 0 . 50 1 1 0.5 6.82 . 0 1 2.26 8.66 . 0 1 3.58 7.25 . 0 1 5.91 4.32 . 0 1 10 5.95 . 0

Splines

From: Andre Jackson Date: January 17, 2019 technical
There is a readme text posted on line (link below) that originated from Icon some time ago: https://www.google.com/url?sa=t&rct=j&q=&esrc=s&source=web&cd=6&cad=rja&uac t=8&ved=2ahUKEwiQ4ozXuvLfAhWCZd8KHVW2DRIQFjAFegQIARAB&url=https%3A%2F%2Fnonm em.iconplc.com%2Fnonmem%2Fconstrained%2520longitudinal%2520spline%2Freadme.t xt&usg=AOvVaw07-r6tKLZsz_5GMOTYBMuM readme.txt ... https://www.google.com/url?sa=t&rct=j&q=&esrc=s&source=web&cd=6&cad=rja&uac t=8&ved=2ahUKEwiQ4ozXuvLfAhWCZd8KHVW2DRIQFjAFegQIARAB&url=https%3A%2F%2Fnonm em.iconplc.com%2Fnonmem%2Fconstrained%2520longitudinal%2520spline%2Freadme.t xt&usg=AOvVaw07-r6tKLZsz_5GMOTYBMuM https://www.google.com/url?sa=t&rct=j&q=&esrc=s&source=web&cd=6&cad=rja&uac t=8&ved=2ahUKEwiQ4ozXuvLfAhWCZd8KHVW2DRIQFjAFegQIARAB&url=https%3A%2F%2Fnonm em.iconplc.com%2Fnonmem%2Fconstrained%2520longitudinal%2520spline%2Freadme.t xt&usg=AOvVaw07-r6tKLZsz_5GMOTYBMuM https://nonmem.iconplc.com/nonmem/constrained%20longitudinal%20spline/readme .txt 1. https://webcache.googleusercontent.com/search?q=cache:xWJ5Iy0lqgsJ: https:// nonmem.iconplc.com/nonmem/constrained%2520longitudinal%2520spline/readme.txt +&cd=6&hl=en&ct=clnk&gl=us Cached Nov 18, 1997 - The longitudinal splines are constrained to exhibit (i) decreasing tail (ii) ... random effects denoted ETAj in NONMEM, where ETAj is distributed ... You've visited this page 5 times. Last visit: 1/15/19 When I tried to run the code for a single dose I got an error message and since the code referenced Icon, I thought that I would ask you before sending it to NMusers. The control and data are attached. The error message was that the characters in error were: Pack $268 SUBROUTINES EEROR WHILE OPENING OR COPYING USED SUPPLIED SUBROUTINE SDSPLINE_DAT.csv Description: MS-Excel spreadsheet TEST.ctl Description: Binary data

Splines

From: Robert Bauer Date: January 17, 2019 technical
Regarding Andre Jackson's query about using the Spline example in the NONMEM repository, the reason for the error message is that the appropriate OTHER routines, such as pack.f and single.f, that are referenced in the control stream file, need to be available. When using an example, it is good idea to copy all of the pertinent files from that directory. But this leads to another matter, that these OTHER files that some examples use are coded in Fortran 77. These would need to be re-coded into F90 files in order to be used in NONMEM 7. I will soon be putting out another message to nmusers about a utility that will do that. Robert J. Bauer, Ph.D. Senior Director Pharmacometrics R&D ICON Early Phase 820 W. Diamond Avenue Suite 100 Gaithersburg, MD 20878 Office: (215) 616-6428 Mobile: (925) 286-0769 [email protected]<mailto:[email protected]> http://www.iconplc.com/