GLM in R and the power to know [🇷 for BE/BA]

posted by d_labes  – Berlin, Germany, 2008-09-25 18:45 (6060 d 03:27 ago) – Posting: # 2411
Views: 60,145

Dear ElMaestro, dear HS,

I cannot follow your discussion exactly because I have no experience with R up to now.

Without entering combat in the war of religions ;-) about type1, type3 and so on 'Sum of squares' (Ssq) and and their appropriateness:

SAS Proc GLM behaves like R, I think.

If the model is specified as
(1) log(PKparm) = sequence treatment period subject
you will get DF=0 and type3 Ssq = missing or zero (I do not remember my first attempt, in which I tried this model also.

Only if the model is specified as
(2) log(PKparm) = sequence treatment period subject(sequence)
where subject(sequence) means subject nested within sequence, an appropriate DF and a value for the type 3 Ssq is obtained. This is because the Ssq(sequence) and Ssq(subject(sequence)) sum up to the Ssq(subject) which means model (1) is over-specified.

So I presume that R gives analog results if the model is specified according to formula (2).
I'm not sure (as already said, ≤zero experience with R) but think I found nesting operator in R is | (pipe).
Thus do me the favor and try it in R.
—
Regards,

Detlew
Thread locked

Complete thread:

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

The whole purpose of education is
to turn mirrors into windows.    Sydney J. Harris

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