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

posted by ElMaestro  – Belgium?, 2019-01-07 18:12  – Posting: # 19754
Views: 1,374

Dear d_labes,

» 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.

Thanks a lot for commenting :-)

Can you try this:

library("lsmeans")
set.seed(12341)
v1=c(rep("A",10),rep( "B",10), rep( "C", 10))
v2=c(rep("Q",10),rep( "W",10), rep( "E",10))
Foo=sample(v1)
Bar=sample(v2)
y=runif(30)
M=lm(y~factor(Foo)+factor(Bar))

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


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.

Le tits now.

Best regards,
ElMaestro

Complete thread:

Activity
 Admin contact
20,246 posts in 4,261 threads, 1,397 registered users;
online 5 (0 registered, 5 guests [including 1 identified bots]).
Forum time (Europe/Vienna): 01:25 CET

To call the statistician after the experiment is done
may be no more than asking him to perform a postmortem examination:
he may be able to say what the experiment died of.    R.A. Fisher

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