回复: Block versus diagonal omega
Dear Al Berg
Take what you cited as an example.If -0.16(correlation between ETA2 and
ETA3) was removed,the number
0.31
0.91 0.63 representing meaning would not exist.Because the ETA1 was
correlated to ETA2,
there is impossible to be correlated to ETA3 when
-0.16 was removed.That is why 0.31 does exist.And so forth.
Here is the block:
$OMEGA BLOCK(2)
1.00
0.60 1.00
$OMEGA BLOCK(2)
1.00
0.58 1.00
The covariance being included could only reduce the variability of
parameters in
the diagonal elements,in my opinion.So regarding your last question,my reply is
no.
This issue have been discussed several times before,maybe you could index
on "covariance" to find it out in our mail list.I am not sure that I am
misleading you,
for I am a nmuser for only 2 years.But indeed,there is no consensus on the off-
diagonal element.The only consensus is that the model should fit the data as
far
as possible.
hongbo ye
from nanjing city.
2010-08-27
yhb5442387
发件人: "Berg, Alexander K., Pharm.D., Ph.D." <[email protected]>
发送时间: 2010-08-26 22:04
主 题: RE: [NMusers] Block versus diagonal omega
收件人: "yhb5442387" <[email protected]>
Hongbo -
Thank you for your reply to my question. Regarding your suggestions, I
understand what you mean but I'm unsure how to apply them in the following
situation. (I'm not really clear how to work with the off-diagonal elements as
I'm new to NONMEM, so this may be a trivial question.) For example, I have a
2-comp model (Cl, V1, Q, V2) with a correlation matrix as follows:
ETA1 ETA2 ETA3 ETA4
1.00
0.60 1.00
0.31 -0.16 1.00
0.91 0.63 0.58 1.00
How do I simplify this? If I decide that I want to remove only the correlation
between ETA2 and ETA3 (which is the smallest), how can I do that without
running into an error in NONMEM? Also, is it even necessary to look at the
off-diagonal elements at this point since I'm in the structural model building
portion of my analysis?
Thank you for your help and feedback -
Al Berg
Quoted reply history
From: [email protected] [mailto:[email protected]] On
Behalf Of yhb5442387
Sent: Thursday, August 26, 2010 4:38 AM
To: nmusers
Subject: RE: [NMusers] Block versus diagonal omega
Hi Al,
Serge's suggestion is available in practise,however when we are considering
to add one covariate such as body weight to the parameter-CL,e.g.,the number of
off-diagonal elements retained in the base model may be different from the one
in the covariate model .As I have noticed,one or above off-diagonal elements
could cross the zero cutoff again and should be excluded from the OMEGA block
structure.So it is hardly to keep the constructure of OMEGA block same during
the model improvment .
In my opinion,if all the diagonal elements fall in an acceptable
interval,such as the CV of parameter is within 50%,there is no need to insert
the off-diagonal element.The off-diagonal element which means covariance
between the diagonal paramete,represents the correlation between them.So
another way is to refer to the scatter plots between ETAs estimated by the
model with diagonal elements .Which off-diagonal element is included depends
on the correlation between two ETAs in the scaterr plots.
Most frequently,the diagonal elements are enough.Do not be worried about that.
By the way, when we discussed the off-diagonal issue,we should not forget the
basic purpose of model building-to make the model predictive performance to be
in accordance with the observed values as far as possible.
Jeroen,
Do you mean the off-diagonal elements instead of diagonal elements when you
mentioned in the second paragraph,because I would like to believe the the
off-digonal elements are more difficult to estimate
hongbo ye
from nanjing city.
2010-08-26
yhb5442387
发件人: "Serge Guzy" <[email protected]>
发送时间: 2010-08-26 04:46
主 题: RE: [NMusers] Block versus diagonal omega
收件人: "Berg, Alexander K., Pharm.D., Ph.D." <[email protected]>,
<[email protected]>
I am not sure there is one single statistical test you can use like we do with
covariate selection (forward followed by backward deletion method).
The easiest way to deal with this problem would be first to use a stable method
like importance sampling assisted by MAP estimation (IMPMAP in NONMEM7) and
getting the full variance covariance matrix and correlation matrix. NONMEM7
will give you also like SADAPT the standard errors associated with each
correlation coefficient. A way to categorize these correlation coefficients
would be to look at each correlation mean +- 2 standard errors and see if it
crosses the zero cutoff. If so, you would assume this correlation not to be
statistically significant. Once all the not statistically significant
correlations are deleted, you have your new blocks to be considered (I guess
you have sometimes to change the order of your parameters to define this new
block in NONMEM7) and you refit your model with this new blocks. Of course,
this is an approximation but at least it allows you ranking the most important
correlations based on both their mean but also their corresponding standard
errors.
A pure diagonal variance covariance matrix will affect the outcome of your
subsequent simulations and usually would inflate the response variability
across the population as important correlations are may be missing.
Serge Guzy; Ph.D
President, CEO; POP_PHARM; INC;
www.poppharm.com
[email protected]
510 684 87 40
From: [email protected] [mailto:[email protected]] On
Behalf Of Berg, Alexander K., Pharm.D., Ph.D.
Sent: Wednesday, August 25, 2010 12:20 PM
To: [email protected]
Subject: [NMusers] Block versus diagonal omega
Hello -
I was curious if someone from the group could perhaps describe the basis for
deciding whether to use a block (variance and covariance) versus diagonal
(variance only) form of omega. Specifically, what tests if any can be
performed to decide between the two forms and are there certain situations
where one is preferred over the other as I often see only the diagonal form
used. Any help would be much appreciated -
Al Berg, PhD/PharmD
Clinical Pharmacology Fellow
Mayo Clinic - Rochester
[email protected]
The information contained in this email message may contain confidential or
legally privileged information and is intended solely for the use of the named
recipient(s). No confidentiality or privilege is waived or lost by any
transmission error. If the reader of this message is not the intended
recipient, please immediately delete the e-mail and all copies of it from your
system, destroy any hard copies of it and notify the sender either by telephone
or return e-mail. Any direct or indirect use, disclosure, distribution,
printing, or copying of any part of this message is prohibited. Any views
expressed in this message are those of the individual sender, except where the
message states otherwise and the sender is authorized to state them to be the
views of XOMA.