Variance=0 [RSABE / ABEL]

posted by Helmut Homepage – Vienna, Austria, 2013-05-03 18:07 (4431 d 09:30 ago) – Posting: # 10535
Views: 25,490

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
23,424 posts in 4,927 threads, 1,672 registered users;
34 visitors (0 registered, 34 guests [including 10 identified bots]).
Forum time: 03:38 CEST (Europe/Vienna)

Medical researches can be divided into two sorts:
those who think that meta is better and those
who believe that pooling is fooling.    Stephen Senn

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