[NMusers] Numerical integreaion and change in clearance over time

2017-12-19 Thread Pavel Belo
Hello NONMEM Users,   When clearance changes over time, the classic analytical solution for the 2-compartment linear model is not correct.  Nevertheless, if the change is only ~30%, one can say it is still acceptable.  Time required to run an analytical and numerical models is drastically diff

RE: [NMusers] Simultaneous pk model of 2 drugs

2016-09-05 Thread Pavel Belo
    Hello Chris,   What is the point of modeling 2 drugs in one NONMEM code?  Do the drugs interact?  For example, you can have 2 monoclonal antibodies competing for or binding to the same target and/or concentration of one drug changes elimination rate or some other parameters of th

RE: [NMusers] Calculation of half life

2016-03-30 Thread Pavel Belo
 You can add formulas to NONMEM code.  All individual noncompartmental parameters can be calculated within NONMEM.  An example of the NONMEM code is below.  tH     DADT(6) = A(8) ;calculates AUC if A is concentration   DATD(8) = [put here the right part of the differential equation which calc

[NMusers] model of error, linearisation, additive error, and VPC

2016-01-29 Thread Pavel Belo
Hello NONMEM Users, When I tried to print log-scaled VPC, there was an error message about negative values.  It can be caused by an additive error and/or linearization of the error model when NONMEM transforms  Y  = F*DEXP(ERR(1)*SD1)  into Y  = F+ F*ERR(1)*SD1.  The error is inflated at

[NMusers] setup of parallel processing and supporting software - help wanted

2015-12-08 Thread Pavel Belo
Hello The Team,   We hear different opinions about effectiveness of parallel processing with NONMEM from very helpful to less helpful.  It can be task dependent.   How useful is it in phase 3 for basic and covariate models, as well as for bootstrapping?    We reached a non-exploratory (produc

[NMusers] message "A VALUE IN THE ESTIMATION RECORD IS INAPPROPRIATE" whene there is no $sigma record

2015-11-26 Thread Pavel Belo
 Hello The NONMEM Team,   I removed $sigma record because ERR(1) and ERR(2) were not used (home-made likelihood was used instead; SD was parameterized as THETA) and got "A VALUE IN THE ESTIMATION RECORD IS INAPPROPRIATE" (see below).  The message sounds rather emotional than mathematical ("in

Re: [NMusers] PRED for BLQ-like observations

2015-11-26 Thread Pavel Belo
- Leonid Gibiansky, Ph.D. President, QuantPharm LLC web:www.quantpharm.com e-mail: LGibiansky at quantpharm.com tel:(301) 767 5566 On 11/20/2015 5:15 PM, Pavel Belo wrote: That is perfect! P. On Fri, Nov 20, 2015 at 04:04 PM, Bauer, Robert wrote: Unfortunately adding records to est

Re: [NMusers] PRED for BLQ-like observations

2015-11-20 Thread Pavel Belo
be largely ignored in the estimation step and not contribute much to run time. Nick On 21-Nov-15 08:59, Pavel Belo wrote: Thank you Bill, In my case it exactly doubles the number of records... The records are daily measures and the code is running slow enough. I'll split the code into es

RE: [NMusers] PRED for BLQ-like observations

2015-11-20 Thread Pavel Belo
expect that the MDV=1 records would be largely ignored in the estimation step and not contribute much to run time. Nick On 21-Nov-15 08:59, Pavel Belo wrote: Thank you Bill, In my case it exactly doubles the number of records... The records are daily measures and the code is running s

[NMusers] mu coding for variability in SD

2015-11-20 Thread Pavel Belo
Hello The NONMEM Professionals,   The model requires SAEM.  Is mu coding useful when each subject has to have different SD (SDi=SD0*exp(etai))?   It may look something like   MU_8  =THETA(8)  SD  = DEXP(MU_8 +ETA(8))   $THETA (0.1,0.0413,50)  ;SD1 $THETA (0.001 FIX) 

RE: [NMusers] PRED for BLQ-like observations

2015-11-20 Thread Pavel Belo
IPRED in the rows with MDV=1.   Thanks,   Bill   From: owner-nmus...@globomaxnm.com [mailto:owner-nmus...@globomaxnm.com] On Behalf Of Pavel Belo Sent: Friday, November 20, 2015 11:47 AM To: nmusers@globomaxnm.com Subject: [NMusers] PRED for BLQ-like observations   Hello The NONMEM

[NMusers] PRED for BLQ-like observations

2015-11-20 Thread Pavel Belo
Hello The NONMEM Users,   When we use M3-like approach, the outputs has PRED for non-missing observations and something else for BLQ (is that PRED=CUMD?).  As in the diagnostic figures PRED for BLQs looks like noise, I remove them.  It is not always perfect, but OK in for most frequent cases. 

Re: [NMusers] eigenvalues

2015-11-16 Thread Pavel Belo
> @PD_value +31 6 23118438 -- More value out of your data! Op 06-11-15 om 17:05 schreef Pavel Belo: NONMEM demonstrated very large differences in objective function when variability or correlations were added or removed.  Monolix demonstrated close-to-insignificant

Re: [NMusers] eigenvalues

2015-11-06 Thread Pavel Belo
matrix.  Monolix teachers recommend to be careful and try a diagonal matrix first.    Thanks, Pavel    On Fri, Nov 06, 2015 at 08:42 AM, Pavel Belo wrote:        Hello Jeroen,   Thank you for your response.  It was a practical question.  I understand the theory.  What is the reason different

Re: [NMusers] eigenvalues

2015-11-06 Thread Pavel Belo
-value.com <http://pd-value.com> jer...@pd-value.com <mailto:jer...@pd-value.com> @PD_value +31 6 23118438 -- More value out of your data! On 11/04/2015 05:55 PM, Pavel Belo wrote: Hello NONMEM Users,   I try to make sense of the results and one of the ways

RE: [NMusers] Job Posting

2015-11-04 Thread Pavel Belo
Dear Ashley, What is the location of the position? Kind regards, Pavel Kovalenko 914-4622515 On Tue, Aug 11, 2015 at 11:30 AM, Ashley Weavers wrote: Dear all,   Hope this email finds you well. Currently I have a senior level Pharmacometrician position from a top, global pharmaceutic

[NMusers] eigenvalues

2015-11-04 Thread Pavel Belo
Hello NONMEM Users,   I try to make sense of the results and one of the ways to do it is to compare the same or similar models across software packages.  5x5 full omega matrix is used because it was prohibitive to remove some insignificant correlations from the matrix without removing significa

Re: [NMusers] OMEGA matrix

2014-09-30 Thread Pavel Belo
er-nmus...@globomaxnm.com> ] On Behalf Of Standing Joseph (GREAT ORMOND STREET HOSPITAL FOR CHILDREN NHS FOUNDATION TRUST) Sent: Friday, September 26, 2014 09:15 To: Kowalski, Ken; 'Eleveld, DJ'; 'Pavel Belo'; nmusers@globomaxnm.com <mailto:nmusers@globomaxnm.com> Su

Re: [NMusers] OMEGA matrix

2014-09-30 Thread Pavel Belo
h (GREAT ORMOND STREET HOSPITAL FOR CHILDREN NHS FOUNDATION TRUST) Sent: Friday, September 26, 2014 09:15 To: Kowalski, Ken; 'Eleveld, DJ'; 'Pavel Belo'; nmusers@globomaxnm.com <mailto:nmusers@globomaxnm.com> Subject: RE: [NMusers] OMEGA matrix Dear Pavel, To answer your

[NMusers] OMEGA matrix

2014-09-25 Thread Pavel Belo
Hello Nonmem Community,   It seems like NONMEM developers may advise to start with full OMEGA matrix at the beginning of model development.  Monolix developers may advise to start with a diagonal matrix.  Is there something different in NONMEM SAEM algorithms that makes model stable when a

RE: [NMusers] SAEM and IMP

2014-05-15 Thread Pavel Belo
ivariate normal, but then use a 1-d pseudo-random sequence to generate the chi-square random variable.    From: owner-nmus...@globomaxnm.com [mailto:owner-nmus...@globomaxnm.com] On Behalf Of Emmanuel Chigutsa Sent: Thursday, May 15, 2014 1:03 PM To: Pavel Belo; nmusers@globomaxnm.com Sub

[NMusers] SAEM and IMP

2014-05-15 Thread Pavel Belo
Hello NONMEM Users, As SAEM does not provide a useful objective function, the manuals recommend using IMP after SAEM. It works well in many cases when IMP works well. When IMP works well, SAEM is not always needed. SAEM is really needed when the other methods do not work well. The issue i

Re: [NMusers] call random: 372 "R" IN CALL RANDOM MAY NOT BE USED GLOBALLY IN ANOTHER BLOCK.

2014-03-27 Thread Pavel Belo
to the error message. Can you send me your complete control stream?   On Wed, Mar 26, 2014, at 05:23 PM, Pavel Belo wrote: Thank you Nick.   It does not work.  The error message is the same.  It is strange.  It almost looks like a bag or an installation issue.  I tried it on

Re: [NMusers] call random: 372 "R" IN CALL RANDOM MAY NOT BE USED GLOBALLY IN ANOTHER BLOCK.

2014-03-26 Thread Pavel Belo
gical test is required      HM=1    ELSE      HM=0    ENDIF ENDIF  On 27/03/2014 11:23 a.m., Pavel Belo wrote:  Hello NONMEM73 Users,   I try to use call random: $PK IF (ICALL.EQ.4.AND.NEWIND.NE.2) CALL RANDOM(2,R) IF (ICALL.EQ.4.AND.R.LE.0.1855) HM=1 IF (ICALL.EQ.4.AND.R.GT.0.1855)

[NMusers] call random: 372 "R" IN CALL RANDOM MAY NOT BE USED GLOBALLY IN ANOTHER BLOCK.

2014-03-26 Thread Pavel Belo
 Hello NONMEM73 Users,   I try to use call random: $PK IF (ICALL.EQ.4.AND.NEWIND.NE.2) CALL RANDOM(2,R) IF (ICALL.EQ.4.AND.R.LE.0.1855) HM=1 IF (ICALL.EQ.4.AND.R.GT.0.1855) HM=0    It does not matter how I try to code it, I am receiving the same error message:   372  "R" IN CA

RE: [NMusers] backward integration from t-a to t

2014-01-17 Thread Pavel Belo
=K0*A(2)-K0*A(3) ... DADT(X)=K0*A(X-1)-K0*A(X) AUCapprox=A(2)+...+A(X) This will prolong the shape of AUCapprox(t). It could be a bit simpler and smoother than tlag implementation Leonid Original email: - From: Pavel Belo non...@optonline.net Date: Thu, 16 Jan 2014 13:

Re: [NMusers] backward integration from t-a to t

2014-01-17 Thread Pavel Belo
r than tlag implementation Leonid Original email: - From: Pavel Belo non...@optonline.net <mailto:non...@optonline.net> Date: Thu, 16 Jan 2014 13:05:54 -0500 (EST) To: lgibian...@quantpharm.com <mailto:lgibian...@quantpharm.com> , nmusers@globomaxnm.c

RE: [NMusers] backward integration from t-a to t

2014-01-16 Thread Pavel Belo
Leonid Original email: ----- From: Pavel Belo non...@optonline.net Date: Tue, 14 Jan 2014 13:45:18 -0500 (EST) To: robert.ba...@iconplc.com, nmusers@globomaxnm.com Subject: [NMusers] backward integration from t-a to t Dear Robert,  Efficacy is frequently considered a function of AUC. (A

RE: [NMusers] backward integration from t-a to t

2014-01-16 Thread Pavel Belo
t start of treatment)   Are you aware of the standard models for effect delay that one would commonly consider and why did you dismiss these?   Best regards   Jakob   From: owner-nmus...@globomaxnm.com [mailto:owner-nmus...@globomaxnm.com] On Behalf Of Pavel Belo Sent: 14 January 2014 18:45

[NMusers] Sr. Statistical Programmer\Statistical Programmer I\II - PK/PD - position can be adjusted for knowledge and expirience

2014-01-15 Thread Pavel Belo
Sr. Statistical Programmer\Statistical Programmer I\II Job Location: Tarrytown, NY Auto req ID: 1208BR Job Posting Title: Sr. Statistical Programmer/Statistical Programmer I/II Person Type: Regular Posting Category: Research & Discovery Location: Tarrytown Job Summary and Essential Functi

RE: [NMusers] backward integration from t-a to t

2014-01-15 Thread Pavel Belo
you aware of the standard models for effect delay that one would commonly consider and why did you dismiss these?   Best regards   Jakob   From: owner-nmus...@globomaxnm.com [mailto:owner-nmus...@globomaxnm.com] On Behalf Of Pavel Belo Sent: 14 January 2014 18:45 To: Bauer, Robert C

[NMusers] backward integration from t-a to t

2014-01-15 Thread Pavel Belo
Please never mind. Someone suggested tose tlag, which is a simple and efficient enough way to do it.   Regards, Pavel        Dear Robert,   Efficacy is frequently considered a function of AUC.  (AUC is just an integral. It is obvious how to calculate AUC any software whic

[NMusers] backward integration from t-a to t

2014-01-14 Thread Pavel Belo
Dear Robert,   Efficacy is frequently considered a function of AUC.  (AUC is just an integral. It is obvious how to calculate AUC any software which can solve ODE.)  A disadvantage of this model of efficacy is that the effect is irreversable because AUC of concentration can only increa

RE: [NMusers] PsN 3.7.6 released

2013-12-16 Thread Pavel Belo
Hello Kajsa, Do you have an example of PsN configuration file for NONMEM7.3/windows? Whan I try to install PsN, it asks me for NONMEM directory for the psn.config file. When I provide it, PsN states it cannot find the nmfe file in run or other directories. Pavel On Wed, Dec 04, 2013 at 09:

RE: [NMusers] Getting rid of correlation issues between CL and volume parameters

2013-12-05 Thread Pavel Belo
Hello Nele,   You can try Ke, V2, K23 and K32 instead of Cl, Q, V2 and V3. It  may mostly remove the correlation between V2 and Cl because  Cl=V2*Ke. If F1 is not a parameter in your model and you give  SC doses, it is the same. You have (Cl/F1) = (V2/F1)*Ke, i.e.  your parameters are "app