Bioequivalence and Bioavailability Forum

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