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

posted by d_labes  – Berlin, Germany, 2019-01-07 19:19 (992 d 09:06 ago) – Posting: # 19756
Views: 3,182

Dear Anders!

» » library("lsmeans")
I see that my suspicion was true. adjust= is an argument only available with package lsmeans (or another package overriding the behavior of confint() for other classes then "lm".
» ...
» confint(pairs(lsmeans(M, "Foo")), adjust="none")
» confint(pairs(lsmeans(M, "Foo")), adjust="bonferroni")
» confint(pairs(lsmeans(M, "Foo")), adjust="eat.my.shorts")
» confint(pairs(lsmeans(M, "Foo")))
» # HAPPY NEW YEAR [/code]
»
» It gives four different results here. You can see adjust is passed as an argument to the confint function so can't see how it would be part of lsmeans or pairs args.

Correct concerning four different results.
adjust="eat.my.shorts" is outstanding because this "method" isn't implemented but transferred via "..." and thus doesn't underlay any error checks. To react on this undesirable feature the programmer decided to set adjust= in such cases to adjust="bonferroni". Reasonable. Or?
This behavior is described in the man page of lsmeans/pairs or summary.emmGrid.
But adjust= is part of package lsmeans. Reason: class(pairs(lsmeans(M, "Foo"))).
This results in calling confint.emmGrid() with the possibility of reacting on multiplicity adjusted CI's instead of confint.lm() without the possibilty of multiplicity adjustments if you call confint().

Regards,

Detlew

Complete thread:

Activity
 Admin contact
21,699 posts in 4,537 threads, 1,542 registered users;
online 12 (0 registered, 12 guests [including 3 identified bots]).
Forum time: Sunday 05:26 CEST (Europe/Vienna)

A central lesson of science is that to understand complex issues
(or even simple ones), we must try to free our minds of dogma and
to guarantee the freedom to publish, to contradict, and to experiment.
Arguments from authority are unacceptable.    Carl Sagan

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