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

posted by yjlee168 Homepage – Kaohsiung, Taiwan, 2010-06-03 15:04 (5443 d 05:43 ago) – Posting: # 5425
Views: 24,075

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.2:- 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,424 posts in 4,927 threads, 1,673 registered users;
90 visitors (0 registered, 90 guests [including 0 identified bots]).
Forum time: 20:47 CEST (Europe/Vienna)

There are two possible outcomes: if the result confirms the
hypothesis, then you’ve made a measurement. If the result is
contrary to the hypothesis, then you’ve made a discovery.    Enrico Fermi

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