Intercept in mixed effects model [Software]

posted by PharmCat  – Russia, 2022-11-03 16:02 (979 d 21:55 ago) – Posting: # 23356
Views: 3,968

(edited on 2022-11-03 16:34)

Hello!

As Helmut said - the intercept term is practically always included in the model.

upd: HaHah)) I was mistaken... corrected interpretation is:

For random effect - you can include intercept, but what will it mean?

For example, if you have repeated effect model with factor 'trt' in sequence AABB for each subject and use random effect model ~ trt (correctly: ~ 0 + trt) with covariance structure "DIAG" you will get Z1 matrix (only 2 obs for subject, but you can make 4):

A B
1 0
0 1
...


With intercept your random effect model (~ 0 + trt) you will get Z2 matrix:

I B
1 0
1 1
...


G matrix same:

a  b
σ₁ 0
0  σ₂


and
σ₁ = 0.2
σ₂ = 0.3

V' = ZGZ'


So in the first case (no intercept), you will get variance estimate vector [a, b] where a - variance for A, b - variance for B.

Z1*G*Z1' =
 0.2  0.0
 0.0  0.3


So in the second case (with intercept), you will get variance estimate vector [a, b] where a - variance for intercept term (actually base level A), b - a "difference" of variance between A and B.

Z2*G*Z2' =
 0.2  0.2
 0.2  0.5


upd

Complete thread:

UA Flag
Activity
 Admin contact
23,428 posts in 4,929 threads, 1,682 registered users;
58 visitors (0 registered, 58 guests [including 15 identified bots]).
Forum time: 14:57 CEST (Europe/Vienna)

To know that we know what we know,
and to know that we do not know what we do not know,
that is true knowledge.    Nicolaus Copernicus

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