Hi R users,

I sent a message yesterday about NA in model estimates (
http://r.789695.n4.nabble.com/How-set-lm-to-don-t-return-NA-in-summary-td3722587.html).
If I use aov() instead of lm() I get no NA in model estimates and I use
gmodels::estimable() without problems. Ok!
Now I'm performing a lot of contrasts and I need correcting for
multiplicity. So, I can use multcomp::glht() for this. However, glht()
return an error message that is not compatible with my expectations. Someone
know I or has a suggestion for? Below some reproducible code.

# toy data
adi <- expand.grid(cult=gl(1,3,la=LETTERS[1]), fert=101)
fat <- expand.grid(cult=gl(2,3,la=LETTERS[2:3]), fert=seq(50,150,50))
da <- rbind(adi, fat)
da$y <- da$fert+rnorm(nrow(da),0,10)

# plot
require(lattice)
xyplot(y~fert|cult, da)

# adjust complete factorial
m0 <- aov(y~cult*fert, subset(da, cult!="A"))
summary(m0)
coef(m0)

# adjust incomplete factorial
m1 <- aov(y~cult*fert, da)
summary(m1)
coef(m1)

require(multcomp)
glht(m0, linfct=matrix(c(1,1,10,0), nrow=1))   # work
glht(m1, linfct=matrix(c(1,1,0,10,0), nrow=1)) # don't work

Thank you.
Walmes.

==========================================================================
Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
Departamento de Estatística - Universidade Federal do Paraná
fone: (+55) 41 3361 3573
VoIP: (3361 3600) 1053 1173
e-mail: walmes@ufpr.br
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================

	[[alternative HTML version deleted]]

