Advertisement
singmann

lsmeans custom contrast simple

Sep 9th, 2014
245
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
R 0.44 KB | None | 0 0
  1. require(lsmeans)
  2. data(md_12.1, package = "afex")  # data from Maxwell & Delaney (2004)
  3. options(contrasts = c("contr.sum", "contr.poly"))
  4.  
  5. str(md_12.1)
  6.  
  7. a1 <- aov(rt~angle*noise+Error(id/(angle*noise)), md_12.1)
  8. summary(a1)  # main effects and interaction match table 12.5 (p. 578)
  9.  
  10. (c1 <- lsmeans(a1, ~angle+noise))
  11.  
  12. contrast(c1, method = list(
  13.   abs_04_vs_8 = c(1, 1, -2, rep (0, 3))/2,
  14.   pres_04_vs_8 = c(rep (0, 3), 1, 1, -2)/2
  15.   ))
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement