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 2.37 KB | None | 0 0
  1. library(tidyverse)
  2. rounded.data = data %>%
  3. mutate(MeanTemp = round(MeanTemp*5)/5) %>%
  4. group_by(MeanTemp) %>%
  5. summarise(Snow = sum(SnowBinary), NoSnow = n() - Snow)
  6.  
  7.  
  8. rounded.model = glm(cbind(Snow, NoSnow) ~ MeanTemp, data = rounded.data, family = binomial())
  9.  
  10. # Deviance GOF test
  11. pchisq(rounded.model$deviance, rounded.model$df.residual, lower.tail = F)
  12.  
  13. cut.data = data %>%
  14. mutate(MeanTemp = cut(MeanTemp,seq(-25,25,5))) %>%
  15. group_by(MeanTemp) %>%
  16. summarise(Snow = sum(SnowBinary), NoSnow = n() - Snow)
  17.  
  18. cut.model = glm(cbind(Snow, NoSnow) ~ MeanTemp, data = cut.data, family = binomial())
  19.  
  20. snow = tibble::tribble(
  21. ~MeanTemp, ~SnowBinary,
  22. 3.5, 0,
  23. 3.7, 0,
  24. 2.5, 0,
  25. 4.9, 0,
  26. 5, 0,
  27. -0.3, 0,
  28. -2.7, 0,
  29. -2.8, 0,
  30. -1.7, 0,
  31. -4, 1,
  32. -5.5, 0,
  33. -6.2, 1,
  34. -9.8, 1,
  35. -11.9, 0,
  36. -9.6, 0,
  37. -5, 1,
  38. -6, 0,
  39. -0.4, 0,
  40. 3.4, 0,
  41. -0.9, 0,
  42. -2, 0,
  43. -2.6, 1,
  44. -1.7, 0,
  45. -6.9, 1,
  46. -9.9, 1,
  47. -15.9, 0,
  48. -16, 0,
  49. -14.9, 0,
  50. -12.6, 1,
  51. -15.8, 1,
  52. -18.4, 1,
  53. -15.6, 1,
  54. -13.4, 0,
  55. -13.4, 0,
  56. -14.2, 1,
  57. -18.3, 0,
  58. -19.9, 1,
  59. -13.8, 0,
  60. -1.1, 1,
  61. -7.4, 1,
  62. -3.6, 0,
  63. 8.1, 0,
  64. 0.4, 0,
  65. -12, 0,
  66. -10.8, 0,
  67. -9.6, 1,
  68. -9.2, 1,
  69. -10.8, 1,
  70. -7.9, 0,
  71. -2.6, 0,
  72. 0.8, 0,
  73. -1.1, 0,
  74. 4.5, 0,
  75. 1.5, 0,
  76. -6.3, 0,
  77. -6.9, 1,
  78. 0.1, 0,
  79. 3.4, 0,
  80. 0.7, 0,
  81. -5.9, 0,
  82. -10.6, 1,
  83. -4.9, 0,
  84. -4.6, 1,
  85. -11.1, 1,
  86. -6.6, 1,
  87. -4.8, 1
  88. )
Add Comment
Please, Sign In to add comment