S×F vari­­ance: 3rd opinion [General Sta­tis­tics]

posted by AngusMcLean – USA, 2015-01-30 00:38 (3813 d 12:19 ago) – Posting: # 14330
Views: 41,461

❝ @Angus: I get the same complete 71 subjects (69 df) for T, 73 (71 df) for R, and 69 (67 df) as John reported in this post. From that I get exactly the variances he reported in SAS. The only difference is the Cinv for 67 df. I would get 49.16227 for p 0.05 instead of 0.95. Maybe I got sumfink wrong?

          dfd  Var_wt       Cinv

❝ dlat(T)    69  0.116539674  89.39120787

❝ dlat(R)    71  0.199313551  91.67023918

❝ ilat(T–R)  67  0.165897781  87.10807220


❝ Can you check your code a fourth time, please? When you copied the R-workflow and modified it for T, please check whether the fixed effect in intermediate dlat is indeed Sequence (and not empty!). Sometimes PHX “forgets” the model specification during copy/pasting.


I am able to resume this work. John sent me the list of subjects he used in SAS (total 71 subjects). I looked at the list of subjects I have and it is 69 subjects. The extra subjects he has over me are numbers 24 and 31.

I looked at the object within the "Copy of Prepare data sets for analysis….." entitled Dij complete rows Data Wizard. This is the point where incomplete data is excluded by the transformation. I see the commands for the transformation and reproduce them below.

Exclude where [LOG DATA R1] is NULL entire ROW
Exclude where [LOG DATA R2] is NULL entire ROW


For the transformation I changed over to LOG DATA T1 for x column and LOG DATA T2 for the Y column

John sent me his data set for Wt calculation and I compared with the subjects in mine.

This transformation excludes subjects 24 and 31 retained in John's data set so I get 69 instead of 71 subjects. Is it the correct exclusion criteria to apply?

The fixed effect is set to sequence in intermediate file mentioned above.

Angus

Complete thread:

UA Flag
Activity
 Admin contact
23,428 posts in 4,929 threads, 1,687 registered users;
52 visitors (0 registered, 52 guests [including 9 identified bots]).
Forum time: 13:57 CEST (Europe/Vienna)

If I’d observed all the rules,
I’d never have got anywhere.    Marilyn Monroe

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