check code for Type III in R [🇷 for BE/BA]

posted by ElMaestro  – Denmark, 2015-02-24 10:59 (3731 d 17:42 ago) – Posting: # 14498
Views: 5,303

Hi WHM,

❝ if i want to get the SS of SEQ and residuals, should i use other function? or can get from drop1()?


Specifically for
a: lm.test1<-lm(log(Cmax)~TRT+PRD+SEQ+SID:SEQ, data=df)

you will see a zero value for Sequence.
If you want to mimic SAS' output, then I think you could try
b: lm.test1<-lm(log(Cmax)~TRT+PRD+SEQ, data=df)
and then you grab the F-value and P-value you get for Sequence in "b" and insert it in the anova table that you got with "a".

I hope this revolves it.

Pass or fail!
ElMaestro

Complete thread:

UA Flag
Activity
 Admin contact
23,424 posts in 4,927 threads, 1,669 registered users;
137 visitors (0 registered, 137 guests [including 7 identified bots]).
Forum time: 05:42 CEST (Europe/Vienna)

Only dead fish go with the current.    Scuba divers' proverb

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