API tools faq
paste
Advertisement
Guest User

Untitled

a guest
Jun 18th, 2019
74
0
Never
Add comment
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 2.36 KB | None | 0 0
raw download clone embed print report
  1. mod3 <- glmmTMB(sum.50 ~ location*trtmt + julian2 + min_aft_sunrise2 + mintemp2 + (1|week) + (1|site:day), data = wbnu_sub, ziformula=~1, family=poisson)
  2.  
  3. mod3.emm <- emmeans(mod3, ~ location|trtmt, adjust = "none")
  4. locdiffs <- contrast(mod3.emm, "revpairwise", adjust = "mvt")
  5. locdiffs
  6.  
  7. trtmt = pulsed_nofood:
  8. contrast estimate SE df t.ratio p.value
  9. near - far -0.273 0.586 124 -0.467 0.6416
  10.  
  11. trtmt = constant:
  12. contrast estimate SE df t.ratio p.value
  13. near - far 2.024 0.376 124 5.381 <.0001
  14.  
  15. trtmt = pulsed_food:
  16. contrast estimate SE df t.ratio p.value
  17. near - far 2.737 0.728 124 3.758 0.0003
  18.  
  19. Results are given on the log (not the response) scale.
  20.  
  21. c <- confint(contrast(locdiffs, "revpairwise", by = NULL, adjust = "mvt"))
  22. c
  23.  
  24. contrast estimate SE df lower.CL upper.CL
  25. near - far,constant - near - far,pulsed_nofood 2.297 0.696 124 0.652 3.94
  26. near - far,pulsed_food - near - far,pulsed_nofood 3.011 0.935 124 0.801 5.22
  27. near - far,pulsed_food - near - far,constant 0.714 0.820 124 -1.223 2.65
  28.  
  29. Results are given on the log (not the response) scale.
  30. Confidence level used: 0.95
  31. Conf-level adjustment: mvt method for 3 estimates
  32.  
  33. c.response <- confint(contrast(locdiffs, "revpairwise", by = NULL, adjust = "mvt", type = "response"))
  34.  
  35. contrast ratio SE df lower.CL upper.CL
  36. near / far,constant / near / far,pulsed_nofood 9.95 6.93 124 1.919 51.6
  37. near / far,pulsed_food / near / far,pulsed_nofood 20.31 18.99 124 2.228 185.1
  38. near / far,pulsed_food / near / far,constant 2.04 1.67 124 0.294 14.2
  39.  
  40. Confidence level used: 0.95
  41. Conf-level adjustment: mvt method for 3 estimates
  42. Intervals are back-transformed from the log scale
  43.  
  44. locdiffs.response <- contrast(mod3.emm, "revpairwise", adjust = "mvt", type = "response")
  45.  
  46. trtmt = pulsed_nofood:
  47. contrast ratio SE df lower.CL upper.CL
  48. near / far 0.761 0.446 124 0.238 2.43
  49.  
  50. trtmt = constant:
  51. contrast ratio SE df lower.CL upper.CL
  52. near / far 7.566 2.845 124 3.594 15.93
  53.  
  54. trtmt = pulsed_food:
  55. contrast ratio SE df lower.CL upper.CL
  56. near / far 15.447 11.252 124 3.653 65.32
  57.  
  58. Confidence level used: 0.95
  59. Intervals are back-transformed from the log scale
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement
Public Pastes
Advertisement
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy.  OK, I Understand
Not a member of Pastebin yet?
Sign Up, it unlocks many cool features!