Variance=0 [RSABE / ABEL]

posted by Helmut Homepage – Vienna, Austria, 2013-05-03 18:07 (4008 d 06:45 ago) – Posting: # 10535
Views: 21,837

Dear Detlew!

❝ Please note the 'Variance' parameter in case of AUCinf.


… and the negative one obtained in PHX. Seems that SAS’ RLME-engine forces negative values to zero.

❝ See this post to notice that we really need 'Variance'=0 in our model (Type CS is only an approximation to that end, hoping the 'Variance' parameter is fitted with a value near zero).


Ooh – that one. :-D

❝ So we eventually have the optimizer to tell that for AUCt.


I’m not very optimistic whether this is possible in PHX at all; I will ask Pharsight. In PHX for linear mixed effects models the initial estimates are derived by the method of moments:
log AUCt

Starting estimates of variance parameters:
csDiag_11                           0.00881650
csBlock_11                          0.188267
Var(Period*Formulation*Subject)_21  0.0519883
Var(Period*Formulation*Subject)_22  0.0703479

log AUCi

Starting estimates of variance parameters:
csDiag_11                          -0.00961080
csBlock_11                          0.148549
Var(Period*Formulation*Subject)_21  0.0620715
Var(Period*Formulation*Subject)_22  0.0894799


Alternatively one can state initial variances. If I set csDiag_11 (PHX’ terminology) to zero (whilst keeping the others), the optimizer iterates happily around (four iterations instead of two) – only to end up with the same final estimates…

Dif-tor heh smusma 🖖🏼 Довге життя Україна! [image]
Helmut Schütz
[image]

The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes

Complete thread:

UA Flag
Activity
 Admin contact
22,993 posts in 4,828 threads, 1,656 registered users;
57 visitors (0 registered, 57 guests [including 5 identified bots]).
Forum time: 00:53 CEST (Europe/Vienna)

If you don’t like something change it;
if you can’t change it, change the way you think about it.    Mary Engelbreit

The Bioequivalence and Bioavailability Forum is hosted by
BEBAC Ing. Helmut Schütz
HTML5