Advertisement
singmann

combining multiple contrasts in lsmeans

Nov 5th, 2014
196
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
R 0.43 KB | None | 0 0
  1. require(lsmeans)
  2. require(multcomp)
  3.  
  4. afex::set_sum_contrasts()
  5. data(oats, package = "MASS")
  6.  
  7. oats$Nf <- ordered(oats$N, levels = sort(levels(oats$N)))
  8. oats.aov <- aov(Y ~ Nf*V + Error(B/V), data = oats, qr = TRUE)
  9.  
  10. (x <- lsmeans(oats.aov, ~Nf))
  11. (y <- pairs(x))
  12.  
  13. 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