[R] pairwise comparisons with multcomp package
gaiarrido
gaiarrido at usal.es
Mon Feb 13 18:16:07 CET 2012
Hi,
I've got this model and following Hothorn et al advices, I used glht for a
post hoc comparison
> modezqM<-glm(rojos~estacion*zona3,quasipoisson,subset=(edadysexo=="M"))
> anova(modezqM,test="F")
Df Deviance Resid. Df Resid.
Dev F Pr(>F)
NULL 293 41148
estacion 1 3647.0 292 37501
23.6250 1.927e-06 ***
zona3 2 931.4 290 36569
3.0167 0.050510.
estacion:zona3 2 2421.7 288 34148
7.8438 0.000482 ***
> summary(glht(modezqM,linfct=mcp(zona3="Tukey",interaction_average=T)))
Estimate Std. Error z-value
Pr(>|z|)
norte - turist == 0 -0.05182 0.23701 -0.219 0.9727
occid - turist == 0 -1.05970 0.41618 -2.546 0.0271 *
occid - norte == 0 -1.00788 0.45830 -2.199 0.0665 .
I'm trying to understanding the way this command made pairwise comparisons,
but, where does this z-value comes from?
I'm suspected this not come from Tukey's test, but i don't know it's Wald or
something i've never use
I know it is not an habitual doubt.
Sorry and thanks in advance
I suspect it
-----
Mario Garrido Escudero
PhD student
Dpto. de Biología Animal, Ecología, Parasitología, Edafología y Qca. Agrícola
Universidad de Salamanca
--
View this message in context: http://r.789695.n4.nabble.com/pairwise-comparisons-with-multcomp-package-tp4384499p4384499.html
Sent from the R help mailing list archive at Nabble.com.
More information about the R-help
mailing list