random vs. fixed [Software]

posted by Helmut Homepage – Vienna, Austria, 2011-10-07 23:35 (5365 d 09:42 ago) – Posting: # 7447
Views: 8,133

Hi Yicaoting & ElMaestro!

❝ ❝ For dataset 3, results from WNL:

❝ ❝ LSM_R: 82.5594 (WNL)  vs  82.5594 (SAS)

❝ ❝ LSM_T: 79.6926 (WNL)  vs  79.2074 (SAS)


❝ Hrmmmmmmmfff...


Yessir.

❝ However, there is -at least in theory- an alternative when one value is missing for a period in one (or more) subject(s) and that is to try a maximum likelihood approach where you specify subject as random in the mixed model and trt+seq+per all fixed. When I do that in R, I actually can reproduce your values from WNL (but I do not have WNL on my machine so cannot play around). It could thus be that WNL actually uses a mixed model to obtain the estimates?


Right guess. Phoenix/WinNonlin’s default in BE is:
fixed:  Sequence+Formulation+Period
Random: Subject(Sequence)


Therefore we get REML estimates. If we delete the random effect and specify the model as ‘all fixed’
Sequence+Formulation+Period+Subject(Sequence)
WNL will spit out exactly SAS’ results (LSM_T, SE, Diff, and CI) for dataset 3.

❝ Someone, go read the manual?!


Wasn’t necessary. ;-)

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,654 posts in 4,992 threads, 1,571 registered users;
151 visitors (0 registered, 151 guests [including 22 identified bots]).
Forum time: 09:18 CEST (Europe/Vienna)

Always listen to experts.
They’ll tell you what can’t be done and why.
Then do it.    Robert A. Heinlein

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