More obvious to discuss the model.... for me... [Software]

posted by ElMaestro  – Denmark, 2022-11-10 01:34 (908 d 11:05 ago) – Posting: # 23359
Views: 3,135

Hi all,

The mixed model is for our purposes
Y= Xb + Zu + e

Z is the design matrix for the random effects expressed in u, typically the betweens.
Now, in BE if we decide to fit the model with an intercept in Z then obviously we will get a column of ones followed by a column with indices for either treatment A or treatment B, but not both. Under this model specification we are thus estimating the intercept (whose interpretation is not so straightforward if you ask me, but see Pharm Cat's example above) along with the betweens for A or B (implementation-dependent), but not both.
One the other hand, if we want to have separate variance estimates for A and B extracted directly from u once optimised, then we would make sure that the random effects are fit without intercept, implying Z without a column of ones, and Bob's your uncle.

I do not see any good arguments for fitting random effects with an intercept column in Z, but I do see arguments for doing it without. If the software takes care of everything from optimisation to generation of a CI behind the curtains, then there is no preference and software defaults will suffice whatever they are since the CI will be the same.

Pass or fail!
ElMaestro

Complete thread:

UA Flag
Activity
 Admin contact
23,424 posts in 4,927 threads, 1,672 registered users;
92 visitors (0 registered, 92 guests [including 11 identified bots]).
Forum time: 13:39 CEST (Europe/Vienna)

It’s always fun to have your models validated,
but is way more fun to have them trashed.
Finding out you are completely wrong
is a great part of science.    G. Randall Gladstone

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