Bioequivalence and Bioavailability Forum • Slightly off topic, but related :-)

Slightly off topic, but related :-) [Design Issues]

posted by Shuanghe – Spain, 2019-12-09 13:00 (1769 d 15:45 ago) – Posting: # 20964
Views: 6,898

Hi ElMaestro,

❝ Note e.g. that the two fits have different residuals and residual df's, which to me means incomplete subjects are not deleted (R does not know and is not being told something is incomplete; the full rank design matrix is still invertible and so on).

What about the following?


X3 <- Xm[-7,]


M3 <- lm(logCmax ~ factor(Seq)+factor(Subj)+factor(Trt)+factor(Per), data=X3)

anova(M3) 

lsmeans(M3, "Trt")

confint(pairs(lsmeans(M3, "Trt"), reverse =F), level=0.9)

Xm has missing period (1) for subject 4, X3 has no subject 4. compare anova(M2) and anova(M3), residual and df of residual are same. 90% CI also same. So wouldn't it mean that R deleted the extra period (2) of subject 4 in Xm automatically when doing BE evaluation? Lsmeans are different, so subject 4 period 2 was kept for that calculation. I woulds say that this behaviour is the same as SAS.

—
All the best,
Shuanghe

Complete thread:

Should those subjects have only one period data be included in BE analysis? ssussu 2019-12-06 10:58 [Design Issues]
- No way! But... Beholder 2019-12-06 12:36
  - No way! But... ElMaestro 2019-12-06 13:14
    - No way! But... Helmut 2019-12-06 14:47
      - No way! But... PharmCat 2019-12-06 23:46
        
        No way! But... wienui 2019-12-07 04:01
        
        EMA guideline: no way...: mittyri 2019-12-07 20:58
        
        EMA guideline: no way...: wienui 2019-12-09 05:25
        
        Slightly off topic, but related :-) ElMaestro 2019-12-08 01:31
        
        Slightly off topic, but related :-)Shuanghe 2019-12-09 12:00
        Slightly off topic, but related :-) PharmCat 2019-12-09 15:13
        
        Slightly off topic, but related :-) ElMaestro 2019-12-21 15:02
        
        Slightly off topic, but related :-) PharmCat 2019-12-22 01:12
        2.220446e-16 ≈ 0 Helmut 2019-12-22 10:37
        
        The optional tolerance argument ElMaestro 2019-12-23 14:37
        2.220446e-16 ≈ 0 PharmCat 2019-12-24 14:18
        
        Sum of residuals ~ ε Helmut 2019-12-24 14:54
        
        Sum of residuals ~ ε ElMaestro 2019-12-24 15:10
        
        Wrong terminology Helmut 2019-12-28 13:55
        
        Sum of residuals ~ ε PharmCat 2019-12-24 18:40