Mean as intercept; model matrices [🇷 for BE/BA]

posted by yjlee168 Homepage – Kaohsiung, Taiwan, 2009-05-21 15:35 (6239 d 18:58 ago) – Posting: # 3721
Views: 14,107

dear Elmaestro,

So far we still cannot figure out why the intercept is not the mean of lnY. However, you put Subj as one of regressors in your lm() functions. In Chow & Liu's model, the S(i,k) is the intersubject residuals of Subj i on Sequence k, not Subj itself. Could this be the reason that these two models might be different? What we use in bear with your example is like the following:
Snot <- lm(lnY~Trt+Subj:Seq+Per+Seq) for a 2x2x2 ABE. It works o.k..

❝ Snot = lm(lnY~Trt+Subj+Seq+Per) ## seems to be standard syntax in

❝ 1. Chow & Liu's model specification is

❝ lnY = Mean + S(i,k) + P(j) + F(j,k) + C(j-1,k) + E(i,j,k)


All the best,
-- Yung-jin Lee
bear v2.9.6:- created by Hsin-ya Lee & Yung-jin Lee
Kaohsiung, Taiwan https://www.pkpd168.com/bear
Download link (updated) -> here

Complete thread:

UA Flag
Activity
 Admin contact
23,655 posts in 4,993 threads, 1,570 registered users;
104 visitors (0 registered, 104 guests [including 12 identified bots]).
Forum time: 10:34 CEST (Europe/Vienna)

Science is simply common sense at its best that is,
rigidly accurate in observation, and
merciless to fallacy in logic.    Thomas Henry Huxley

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