Bioequivalence and Bioavailability Forum • R code for analyzing classical 2X2 crossover designed bioequivalence data

R code for analyzing classical 2X2 crossover designed bioequivalence data [🇷 for BE/BA]

posted by ElMaestro – Denmark, 2016-01-22 01:59 (3010 d 18:44 ago) – Posting: # 15842
Views: 22,068

Hi Lizhao,

it is not too difficult.
You can use lm in the form of something like:

Mod = lm(lnCmax ~Trt+Seq+Subj+Per)

An easy option for the direct treatment effects is to use
Mod = lm(lnCmax ~0+Trt+Seq+Subj+Per)
and then you have the treatment effects directly out of the fixed effects or model summaries for Mod.

To generate a type III anova you can use drop1(Mod, test="F")
To generate the confidence interval, you just need the critical value at your chosen alpha and df, the fixed treatment effects (estimates thereof), and the SE which is calculated from the model's MSE and from the sequence sample sizes (which determine df too). You can find some code via this link.

Good luck. If you get stuck, send an email.

—
Pass or fail!
ElMaestro

Complete thread:

R code for analyzing classical 2X2 crossover designed bioequivalence data lizhao 2016-01-22 00:13
- R code for analyzing classical 2X2 crossover designed bioequivalence dataElMaestro 2016-01-22 00:59
  - R code for analyzing classical 2X2 crossover designed bioequivalence data lizhao 2016-01-22 01:26
    - On random effects and bogus statements ElMaestro 2016-01-22 12:13
      - FDA 'own' SAS code for 2x2 d_labes 2016-01-27 09:57
        
        FDA 'own' SAS code for 2x2 nobody 2016-01-27 10:09
        
        OT: PM to ElMaestro Helmut 2016-01-27 13:39
        On coffee and such... ElMaestro 2016-01-27 13:47
      - Geometric means ± SD Helmut 2016-01-27 13:51
        
        Geometric means ± SD nobody 2016-01-27 15:20
        
        Money makes the world go ’round Helmut 2016-01-27 16:11
        
        Money makes the world go ’round nobody 2016-01-27 16:31
        
        Fixed or random effect kumarnaidu 2016-11-04 14:05
        
        PROC GLM fixes Subject mittyri 2016-11-04 23:00
        
        PROC GLM fixes Subject kumarnaidu 2016-11-05 05:22
        
        PROC GLM fixes Subject mittyri 2016-11-05 10:13
        
        PROC GLM fixes Subject kumarnaidu 2016-11-07 06:24
        
        PROC GLM fixes Subject mittyri 2016-11-07 11:57
        
        PROC GLM fixes Subject kumarnaidu 2016-11-07 12:42
        
        Geometric means ± SD Ben 2016-11-08 21:20
        
        Misunderstanding Helmut 2016-11-09 15:29
        
        Misunderstanding Ben 2016-11-14 19:40
        
        RTFM Helmut 2016-11-14 21:50