confint() for difference of effect levels [🇷 for BE/BA]

posted by d_labes  – Berlin, Germany, 2017-08-29 14:23 (2762 d 06:43 ago) – Posting: # 17749
Views: 8,698

Dear Öberster Größter Meister,

❝ let us say we have an lm object for a BE analysis in a study where we have n treatments.

Ok. I assume we have.

❝ Is there a way to use confint() to derive CI's for the difference of the i'th and j'th treatments (where obviously i!=j, and i is within 1..n and j is within 1..n) ?

No. confint() gives you confidence intervals for the model parameters only.

❝ If not, is there another built-in function that achieves this goal?

Built in I know of none, if you mean a base R installation. Only from add-on packages.
Have a look at package multcomp.
Or if you prefer using such beasts like "least square means" have a look at package lsmeans.

BTW: All-at-once is not the preferred method of a big regulatory body.
Pairwise lm objects (throwing away the data under those treatments not under consideration) will give you the opportunity to use confint() to derive CI's.

Regards,

Detlew

Complete thread:

UA Flag
Activity
 Admin contact
23,409 posts in 4,921 threads, 1,669 registered users;
22 visitors (0 registered, 22 guests [including 5 identified bots]).
Forum time: 20:06 CET (Europe/Vienna)

Facts are meaningless. You could use facts to prove
anything that’s even remotely true!    Homer Simpson

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