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

posted by yjlee168 Homepage – Kaohsiung, Taiwan, 2009-05-21 15:35 (6227 d 19:54 ago) – Posting: # 3721
Views: 14,050

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,653 posts in 4,991 threads, 1,570 registered users;
270 visitors (0 registered, 270 guests [including 20 identified bots]).
Forum time: 11:30 CEST (Europe/Vienna)

The epistemological value of probability theory is based on the fact
that chance phenomena, considered collectively and on a grand scale,
create non-random regularity.    Andrey Kolmogorov

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