confint and adjust in R [General Statistics]
Dear Anders!
Correct concerning four different results.
This behavior is described in the man page of lsmeans/pairs or summary.emmGrid.
But
This results in calling
❝
❝ library("lsmeans")
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
Regards,
Detlew
Complete thread:
- confint and adjust in R ElMaestro 2019-01-07 14:52 [General Statistics]
- confint and adjust in R d_labes 2019-01-07 17:02
- confint and adjust in R ElMaestro 2019-01-07 18:12
- confint and adjust in Rd_labes 2019-01-07 19:19
- confint and adjust in R ElMaestro 2019-01-07 22:22
- confint and adjust in Rd_labes 2019-01-07 19:19
- confint and adjust in R ElMaestro 2019-01-07 18:12
- confint and adjust in R d_labes 2019-01-07 17:02