confint and adjust in R [General Sta­tis­tics]

posted by d_labes  – Berlin, Germany, 2019-01-07 16:02 (573 d 19:22 ago) – Posting: # 19749
Views: 2,548

Dear ElMaestro,

» I am exploring the confint function in R and the oppotunities for multiplicity adjustment.
» In fact, the confint function allows an adjust argument :
» In my current models I can apply adjust= e.g. "bonferroni", "none", "sidak", "mvt" and "tukey" (which is the default).

Sorry, but AFAIK the default confint method or that for lm() objects has no adjust argument, or more precisely spoken accepts any string for adjust= but without any effect. Example:
> fit <- lm(100/mpg ~ disp + hp + wt + am, data = mtcars)
> confint(fit)
                   2.5 %      97.5 %
(Intercept) -0.774822875 2.256118188
disp        -0.002867999 0.008273849
hp          -0.001400580 0.011949674
wt           0.380088737 1.622517536
am          -0.614677730 0.926307310

> confint(fit, adjust="")
                   2.5 %      97.5 %
(Intercept) -0.774822875 2.256118188
disp        -0.002867999 0.008273849
hp          -0.001400580 0.011949674
wt           0.380088737 1.622517536
am          -0.614677730 0.926307310

There is no error thrown because confint() has an ... argument, for which there is no error check available.

What you describe must be an effect of using an add-on package with comes with the possibility to consider multiplicity in calculation of confidence intervals.
If you are lucky, the documentation of that package will contain details of what the adjust= argument will do. Look into the ./doc sub-directory of the package and look for vignettes or other documentation material.
But here I must say that being lucky is not so often the case:no: in R. Many authors of packages assume that you have sufficient knowledge of the statistical methods underlying a package they create. Eventually you may figure out whats going on if you look at the code of the confint.class() function, where class is the model class, defined in the package you uses.



Complete thread:

 Admin contact
20,961 posts in 4,373 threads, 1,459 registered users;
online 7 (0 registered, 7 guests [including 4 identified bots]).
Forum time: Monday 11:24 UTC (Europe/Vienna)

Garbage in, garbage out.
It’s very very simple.    Anders Fuglsang

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