Joachim, Yes, I included dummy observations before each new dosing time, therefore, it is maybe the easiest way to read my table in R and let R determine the minimum biomarker conc for each ID
Thank you very much, Friederike Von: Joachim Grevel [mailto:jgre...@o2.co.uk] Gesendet: Mittwoch, 6. April 2011 15:39 An: Friederike Kanefendt; nmusers@globomaxnm.com Betreff: RE: [NMusers] Modeling Cmin and Tmin Friederike, You can just as well enrich the individual records of each ID with additional times using EVID=2. Then you get predictions for all compartments at these additional times for each ID with its optimal set of ETA estimates without altering the OBFV or the estimates. You can then list these predictions in the output table and play with them as you like. Hope this helps you furher, Joachim From: Friederike Kanefendt [mailto:friederike.kanefe...@uni-bonn.de] Sent: 06 April 2011 14:05 To: Joachim Grevel; nmusers@globomaxnm.com Subject: AW: [NMusers] Modeling Cmin and Tmin Dear Joachim, thank you for your explanations. My plan was not to determine the typical value or distribution for this parameters (as THETA or ETA) but to determine the minimum biomarker concentration for each ID with the appropriate time for further correlation studies with outcome etc. I tried to use NONMEM to determine Cmin while individual parameters/profiles are estimated. I read something in NMusers about Cmax and thought that this is maybe also possible for Cmin and Tmin <http://www.cognigencorp.com/nonmem/current/2008-January/0733.html> http://www.cognigencorp.com/nonmem/current/2008-January/0733.html http://www.cognigencorp.com/nonmem/current/2007-December/0716.html If this is not possible I would try this using simulations and R Thanks King regards, Friederike Von: owner-nmus...@globomaxnm.com [mailto:owner-nmus...@globomaxnm.com] Im Auftrag von Joachim Grevel Gesendet: Mittwoch, 6. April 2011 11:49 An: Friederike Kanefendt; nmusers@globomaxnm.com Betreff: RE: [NMusers] Modeling Cmin and Tmin Dear Friederike, Please, step back for one minute and try to understand what you want to find out: You are applying a nonlinear mixed-effects modelling program to obtain parameter values for your model. All parameters in your model are random variables. Some have a mean of zero (ETAs, EPSs) some should have a mean <> 0 (THETAs). Your results are not point estimates but distributions. Therefore you should approach the question of Cmin and Tmin through simulation and report the confidence interval (or other stats) for your parameters of interest. Thus construct a simulation data set with dense time points around the expected Tmin, and run $SIMULATION (240311) ONLYSIMULATION SUBPROBLEMS=1000 $TABLE ID TIME ...... NOHEADER NOAPPEND NOPRINT FILE=xxxx.tab with all your final parameter estimates as initial values in your control file. Then analyse your 100 (or 1000) simulated data files with R or SAS and report the distribution of the Cmin and Tmin you found. Please, come back with more questions in case I have confused you, Joachim Joachim Grevel, PhD Scientific Director BAST Inc Limited BioCity Nottingham Pennyfoot Street Nottingham, NG1 1GF Tel: +44 (0)115 8120497 From: owner-nmus...@globomaxnm.com [mailto:owner-nmus...@globomaxnm.com] On Behalf Of Friederike Kanefendt Sent: 06 April 2011 09:56 To: nmusers@globomaxnm.com Subject: [NMusers] Modeling Cmin and Tmin Dear NMUser, I am a beginner in NONMEM and I have a question regarding modeling Cmin and Tmin. I have data of a biomarker which decreases after drug intake and I want to determine the minimum change from baseline (Cmin) as well as the corresponding time (Tmin). I tried different models but the estimated values for these parameters are not plausible or always zero. Has someone experience with this or any idea to solve this problem? Thanks a lot Best regards Friederike Attached I send you a part of my code: $SUBROUTINE ADVAN6 TOL=3 ... $DES ... DADT(5) = KIN*INH-KOUT*A(5) ; Biomarker CMT IF(TIME.EQ.0) THEN CMIN = BASE ; Baseline conc also estimated by NM TMIN=0 ENDIF CC=A(5) IF(CC.LT.CMIN) THEN CMIN= CC TMIN= T ENDIF REL = CMIN/BASE ; change relative to baseline ... $ESTIMATION SIG=2 PRINT=1 METHOD=1 INTER MAXEVAL=0 NOABORT Friederike Kanefendt Pharmacist, PhD-Student University of Bonn -Clinical Pharmacy- An der Immenburg 4 D-53121 Bonn Phone: +49 (0)228 73-5781 Fax: +49(0) 228 73-9757 friederike.kanefe...@uni-bonn.de
