Degrees of freedom 3-period full replicate [RSABE / ABEL]

posted by zizou – Plzeň, Czech Republic, 2015-07-06 19:25 (3633 d 15:10 ago) – Posting: # 15045
Views: 31,255

Dear Detlew,

I didn't try the code, so it was only what I expect in result (my mistake - I didn't write that for 2x2x3 replicate design I thought to filter data to only those used in GLM).

For 3-period 2-sequence replicate design with R replicated only in one sequence I would use code:
data var;
set replicate;
if (formulation='R' & sequence='TRR');
run;

proc glm data=var;
class subject period;
model logDATA= suject period;
run;
(Instead of the code which I copied from EMA Q&A in my previous post. I think there is no reason to have sequence in the GLModel for this case. And I am used to see degrees of freedom in ANOVA table which are giving us with real number of cases in evaluation: (number of subjects in evaluation-1) and (number of values in evaluation-1) in total)
In the example, due to we have 35 subjects in sequence with replicated R, I am expecting:
   source    df
   period     1
   subject   34
   Error     34
   Total     69
However, it does not affect anything important.

But if 72+1+34 results in total 108 (similar as 2+2=?, I think 5 is a very very good estimate, ... in accordance with the principles of doublethink "Sometimes they are five. Sometimes they are three. Sometimes they are all of them at once"), I have no power to know (same as Winston Smith has).

Btw. you are right about I did not include balance/unbalance option to table of dfs in my previous post (so it's true only for cases with no dropouts there).

Complete thread:

UA Flag
Activity
 Admin contact
23,424 posts in 4,927 threads, 1,672 registered users;
30 visitors (0 registered, 30 guests [including 19 identified bots]).
Forum time: 10:36 CEST (Europe/Vienna)

In matters of style, swim with the current;
in matters of principle, stand like a rock.    Thomas Jefferson

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