Bioequivalence and Bioavailability Forum • check code for Type III in R

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

posted by ElMaestro – Denmark, 2015-02-24 10:59 (4135 d 21:07 ago) – Posting: # 14498
Views: 6,618

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:

check code for Type III in R WHM 2015-02-16 07:54
- check code for Type III in R ElMaestro 2015-02-16 12:17
  - check code for Type III in R WHM 2015-02-24 01:56
    - check code for Type III in RElMaestro 2015-02-24 09:59