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

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

posted by ElMaestro – Denmark, 2019-12-08 02:31 (1950 d 22:47 ago) – Posting: # 20961
Views: 9,113

Hi all,

❝ ❝ Correct. If we would be allowed (pun!) to use a mixed-effects model. Patterson and Jones argued against this doubtful practice and the ethical implications of discarding data.

❝

❝ Really, there is no problem in mixed model approach to use all data if model constructed correctly.

But would we really need a mixed model ??
I seem to recall a thread about it some time ago. I believe the normal linear model fits without trouble to a 222BE dataset with a missing period, so why would it be necessary to use a mixed model for that situation at all?
As I recall it mixed model done with ML are proven to give unbiased estimates, but it has never been proven that REML fits are unbiased (and REML is generally used in BE, whenever a discussion of the mixed model comes into play), hence I can't imagine that bias would be a reason for picking the mixed model over a linear model when we have a simple missing period situation?!? Am curious to learn more about this.

Here's an example in R:

rm(list=ls(all=TRUE))

library("emmeans")

set.seed(123456)

logCmax=rnorm(10,6, 2)

Subj=c(1,1,2,2,3,3,4,4,5,5)

Seq=c("AB", "AB", "BA", "BA", "AB", "AB", "BA","BA", "BA", "BA")

Per=c(rep(c(1,2), 5))

Trt=substr(Seq, Per, Per)

X=data.frame(Subj, Per, Trt, Seq, logCmax)

X

M1=lm(logCmax ~ factor(Seq)+factor(Subj)+factor(Trt)+factor(Per), data=X)

anova(M1) 

lsmeans(M1, "Trt")

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



Xm=X[-7,] ##let's delete a period

Xm

M2=lm(logCmax ~ factor(Seq)+factor(Subj)+factor(Trt)+factor(Per), data=Xm)

anova(M2) 

lsmeans(M2, "Trt")

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

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).
Many thanks for any light you can shed onto this.

Edit: I added set.seed(123456) to your code in order to get reproducible results. [Helmut]

—
Pass or fail!
ElMaestro

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