FDA loves replicate bears? [R for BE/BA]

posted by yjlee168 Homepage – Kaohsiung, Taiwan, 2009-04-23 08:08 (4609 d 18:15 ago) – Posting: # 3587
Views: 64,490

Dear D. Labes,

Thank you for your message. I find that currently bear's model dealing with replicate ABE data is equivalent to SAS Proc MIXED as the follows:
PROC MIXED;
CLASSES SEQ SUBJ PER TRT;
MODEL Y = SEQ PER TRT/ DDFM=SATTERTH;
RANDOM SUB=SUBJ G;
REPEATED/GRP=TRT SUB=SUBJ;
ESTIMATE 'T vs. R' TRT 1 -1/CL ALPHA=0.1;

(Ref.: Jones, B and Kenward MG. Design and Analysis of Cross-Over Trials, 2nd ed., Chapman & Hall/CRC, 2003). If using TYPE setting, either FA0(2) or CSH as recommended by FDA, we will not get the exactly same 90%CIs from running bear using lme(). I guess we need to find out how to comply with here.

All the best,
-- Yung-jin Lee
bear v2.9.0:- created by Hsin-ya Lee & Yung-jin Lee
Kaohsiung, Taiwan https://www.pkpd168.com/bear
Download link (updated) -> here

Complete thread:

Activity
 Admin contact
21,786 posts in 4,557 threads, 1,548 registered users;
online 16 (1 registered, 15 guests [including 5 identified bots]).
Forum time: Monday 01:23 CET (Europe/Vienna)

The history of statistics is like a telephone directory:
the plot is boring, full of numbers and the cast is endless.    Stephen Senn

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