Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- require(lsmeans)
- require(multcomp)
- afex::set_sum_contrasts()
- data(oats, package = "MASS")
- oats$Nf <- ordered(oats$N, levels = sort(levels(oats$N)))
- oats.aov <- aov(Y ~ Nf*V + Error(B/V), data = oats, qr = TRUE)
- (x <- lsmeans(oats.aov, ~Nf))
- (y <- pairs(x))
- summary(as.glht(lsmobj(bhat = x@bhat, V = x@V,levels = list(test = 1:10),linfct = rbind(x@linfct, y@linfct), df=x@dffun, dfargs = x@dfargs)), test = adjusted("free"))
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement