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

posted by ElMaestro  – Denmark, 2016-01-22 01:59 (2989 d 11:34 ago) – Posting: # 15842
Views: 21,907

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:

UA Flag
Activity
 Admin contact
22,957 posts in 4,819 threads, 1,636 registered users;
96 visitors (0 registered, 96 guests [including 7 identified bots]).
Forum time: 13:33 CET (Europe/Vienna)

With four parameters I can fit an elephant,
and with five I can make him wiggle his trunk.    John von Neumann

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