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

posted by zizou – Plzeň, Czech Republic, 2015-07-06 19:25 (3269 d 23:10 ago) – Posting: # 15045
Views: 28,397

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,057 posts in 4,840 threads, 1,641 registered users;
112 visitors (0 registered, 112 guests [including 7 identified bots]).
Forum time: 18:36 CEST (Europe/Vienna)

You should treat as many patients as possible with the new drugs
while they still have the power to heal.    Armand Trousseau

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