SAS-like Type III SS sequence with bear? [🇷 for BE/BA]

posted by yjlee168 Homepage – Kaohsiung, Taiwan, 2010-06-03 15:04 (5854 d 00:08 ago) – Posting: # 5425
Views: 29,719

Dear ElMaestro and all,

I made some code modifications and got the output like as follows.

Statistical analysis (ANOVA(lm))   
-----------------------------------------------------------
   Dependent Variable: log(Cmax)                           

Analysis of Variance Table

Response: log(Cmax)
          Df   Sum Sq  Mean Sq F value Pr(>F)
prd        1 0.027294 0.027294  1.7187 0.2144
drug       1 0.025583 0.025583  1.6110 0.2284
subj(seq) 12 0.255270 0.021272  1.3395 0.3103
Residuals 12 0.190565 0.015880               

Analysis of Variance Table

Response: log(Cmax)
           DF    Type I SS  Mean Square  F Value  Pr > F
prd         1     0.027294     0.027294   1.7187 0.21439
drug        1     0.025583     0.025583   1.6110 0.22842
subj(seq)  12     0.255270     0.021272   1.3395 0.31028

Analysis of Variance Table

Response: log(Cmax)
           DF  Type III SS  Mean Square  F Value  Pr > F
prd         1     0.027294     0.027294   1.7187 0.21439
drug        1     0.025583     0.025583   1.6110 0.22842
subj(seq)  12     0.255270     0.021272   1.3395 0.31028
-----------------------------------------------------------

Tests of Hypothesis using the Type I MS for
SUBJECT(SEQUENCE) as an error term
            Df  Sum Sq   Mean Sq F value Pr(>F)
seq          1 0.00002 0.0000172   9e-04 0.9764
Residuals   26 0.49871 0.0191813               

Tests of Hypothesis using the Type III MS for
SUBJECT(SEQUENCE) as an error term
            Df  Sum Sq   Mean Sq F value Pr(>F)
seq          1 0.00002 0.0000172   9e-04 0.9764
Residuals   26 0.49871 0.0191813


❝ ❝ seq 0 0.000000 0.21364 -104.52


❝ Since you ask I'd say this line is not informative; Seq will always have a

❝ zero (well, effectively zero) SS with type III when you are using R and

❝ drop1. SAS gives you another value for reasons explained above. [...]



What do you think about this revision? Thanks.

All the best,
-- Yung-jin Lee
bear v2.9.6:- created by Hsin-ya Lee & Yung-jin Lee
Kaohsiung, Taiwan https://www.pkpd168.com/bear
Download link (updated) -> here

Complete thread:

UA Flag
Activity
 Admin contact
23,653 posts in 4,991 threads, 1,570 registered users;
117 visitors (0 registered, 117 guests [including 27 identified bots]).
Forum time: 15:12 CEST (Europe/Vienna)

To propose that poor design can be corrected by subtle analysis techniques
is contrary to good scientific thinking.    Stuart J. Pocock

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