Guest User

Untitled

a guest
Feb 25th, 2018
95
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.55 KB | None | 0 0
  1. # 2 levels, treatment contrasts (default)
  2. # here, by luck, our baseline level happens to be alphabetically first, so we get a sensible
  3. # assignment of treatment=1.
  4. foo <- factor(c("treatment", "control"))
  5. contrasts(foo)
  6. # treatment
  7. # control 0
  8. # treatment 1
  9.  
  10. # 2 levels, sum contrasts (better)
  11. contrasts(foo) <- contr.sum
  12. contrasts(foo)
  13. # [,1]
  14. # control 1
  15. # treatment -1
  16.  
  17. # 2 levels, deviation contrasts (preferred way)
  18. contrasts(foo) <- matrix(c(-0.5, 0.5), nrow=2, dimnames=list(c("control", "treatment"), "predictor"))
  19. contrasts(foo)
  20. # predictor
  21. # control -0.5
  22. # treatment 0.5
  23.  
  24. # 3 levels, dummy coding (default)
  25. # here, by sheer luck, a sensible choice for baseline category happens to also be the alphabetically
  26. # first level of the factor, so you get a reasonable dummy coding.
  27. foo <- factor(c("disease+placebo", "disease+treatment", "controls+placebo"))
  28. contrasts(foo)
  29. # disease+placebo disease+treatment
  30. # controls+placebo 0 0
  31. # disease+placebo 1 0
  32. # disease+treatment 0 1
  33.  
  34. # 3 levels, planned contrasts
  35. foo <- factor(c("disease+placebo", "disease+treatment", "controls+placebo"))
  36. contrasts(foo) <- matrix(c(0, -1/2, 1/2, -2/3, 1/3, 1/3), nrow=3,
  37. dimnames=list(levels(foo), c("treatment", "disease")))
  38. contrasts(foo)
  39. # treatment disease
  40. # controls+placebo 0.0 -0.6666667
  41. # disease+placebo -0.5 0.3333333
  42. # disease+treatment 0.5 0.3333333
Add Comment
Please, Sign In to add comment