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

posted by zizou – Plzeň, Czech Republic, 2015-07-06 19:25 (3446 d 21:24 ago) – Posting: # 15045
Views: 29,863

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,336 posts in 4,902 threads, 1,669 registered users;
23 visitors (0 registered, 23 guests [including 11 identified bots]).
Forum time: 15:50 CET (Europe/Vienna)

Biostatistician. One who has neither the intellect for mathematics
nor the commitment for medicine but likes to dabble in both.    Stephen Senn

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