Hello:
I am running a model using Perl Speaks Nonmem which runs ok and has good
shrinkage. However, when I run the command with (-CWRES) I get the
following error.
AN ERROR WAS FOUND IN THE CONTROL STATEMENTS.
300 VERBATIM "LAST MUST BE AFTER LAST ABBREVIATED CODE
I searched the archives but couldn't find an answer. I have my error block
code listed. My suspicion is that may be where the problem may exist since
the model does run ok.
$ERROR
IF(CMT.EQ.3) IPRED=A(3)/S3
W=SQRT((THETA(13)*F)**2+THETA(14)**2); equation var 2 proost1
CP=IPRED
IRES=DV-IPRED
IWRES=IRES/W
Y=CP +CP*ERR(1)*W + ERR(2)
Verbatim Code
Nowadays you can just specify CWRES in the $TABLE statement of NONMEM to get those residuals, circumventing the issue you ran into. See the nonmem75 pdf for further residuals that are available such as NPDE. (if you really want to investigate what goes wrong, go into the NM_run1 folder and look at the code that PsN has generated for you in psn.mod)
For all on the nm list: best wishes for 2024!
Best regards,
Jeroen
http://pd-value.com
[email protected]
@PD_value
+31 6 23118438
-- More value out of your data!
Quoted reply history
On 08-01-2024 14:12, [email protected] wrote:
> Hello:
>
> I am running a model using Perl Speaks Nonmem which runs ok and has good shrinkage. However, when I run the command with (-CWRES) I get the following error.
>
> AN ERROR WAS FOUND IN THE CONTROL STATEMENTS.
>
> 300 VERBATIM "LAST MUST BE AFTER LAST ABBREVIATED CODE
>
> I searched the archives but couldn't find an answer. I have my error block code listed. My suspicion is that may be where the problem may exist since the model does run ok.
>
> $ERROR
>
> IF(CMT.EQ.3) IPRED=A(3)/S3
>
> W=SQRT((THETA(13)*F)**2+THETA(14)**2); equation var 2 proost1
>
> CP=IPRED
>
> IRES=DV-IPRED
>
> IWRES=IRES/W
>
> Y=CP +CP*ERR(1)*W + ERR(2)
Dear Andre,
You may want to look at how you code your error model.
However, for the question at hand: The need for this option in PsN is legacy, I
believe.
With modern nonmem (versions), you can just request CWRES on $TABLE, and you do
not have to specify how this should be calculated (neither does PsN need to do
this).
Best wishes
Jakob
Jakob Ribbing, Ph.D.
Principal Consultant, Pharmetheus AB
Cell/Mobile: +46 (0)70 514 33 77
[email protected]
www.pharmetheus.com
Phone, Office: +46 (0)18 513 328
Uppsala Science Park, Dag Hammarskjölds väg 36B
SE-752 37 Uppsala, Sweden
Quoted reply history
> On 8 Jan 2024, at 14:12, <[email protected]> <[email protected]>
> wrote:
>
> Hello:
>
> I am running a model using Perl Speaks Nonmem which runs ok and has good
> shrinkage. However, when I run the command with (-CWRES) I get the
> following error.
>
> AN ERROR WAS FOUND IN THE CONTROL STATEMENTS.
> 300 VERBATIM "LAST MUST BE AFTER LAST ABBREVIATED CODE
>
> I searched the archives but couldn't find an answer. I have my error block
> code listed. My suspicion is that may be where the problem may exist since
> the model does run ok.
>
> $ERROR
>
> IF(CMT.EQ.3) IPRED=A(3)/S3
> W=SQRT((THETA(13)*F)**2+THETA(14)**2); equation var 2 proost1
> CP=IPRED
> IRES=DV-IPRED
> IWRES=IRES/W
> Y=CP +CP*ERR(1)*W + ERR(2)
--
*This communication is confidential and is only intended for the use of the
individual or entity to which it is directed. It may contain information
that is privileged and exempt from disclosure under applicable law. If you
are not the intended recipient please notify us immediately. Please do not
copy it or disclose its contents to any other person.*
*Any personal data
will be processed in accordance with Pharmetheus' privacy notice, available
here https://pharmetheus.com/privacy-policy/.**
*