API tools faq
paste
Guest User

Untitled

a guest
Jan 19th, 2018
98
0
Never
Add comment
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 4.78 KB | None | 0 0
raw download clone embed print report
  1. structure(list(hosptg = structure(c(3L, 3L, 1L, 3L, 3L, 3L, 3L,
  2. 3L, 3L, 3L, 3L, 3L, 3L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
  3. 3L, 3L, 3L, 1L, 3L, 3L, 3L, 3L, 2L, 3L, 3L, 3L, 3L, 1L, 1L, 3L,
  4. 3L, 3L, 3L, 1L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
  5. 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 3L, 3L,
  6. 3L, 3L, 2L, 3L, 1L, 3L, 1L, 3L, 3L, 1L, 3L, 3L, 3L, 3L, 3L, 1L,
  7. 3L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L), .Label = c("1",
  8. "2", "3"), class = "factor"), quarter.adm = structure(c(4L, 11L,
  9. 3L, 12L, 7L, 8L, 12L, 9L, 1L, 11L, 7L, 1L, 2L, 2L, 10L, 10L,
  10. 8L, 11L, 6L, 1L, 4L, 6L, 10L, 10L, 6L, 11L, 11L, 7L, 3L, 6L,
  11. 10L, 12L, 7L, 6L, 6L, 3L, 6L, 12L, 4L, 4L, 12L, 1L, 6L, 5L, 11L,
  12. 9L, 4L, 4L, 3L, 10L, 4L, 8L, 10L, 3L, 7L, 1L, 12L, 5L, 4L, 6L,
  13. 6L, 3L, 9L, 7L, 8L, 3L, 7L, 8L, 7L, 6L, 5L, 11L, 9L, 11L, 1L,
  14. 4L, 6L, 5L, 5L, 6L, 5L, 5L, 11L, 3L, 4L, 12L, 12L, 1L, 9L, 9L,
  15. 6L, 9L, 1L, 4L, 8L, 1L, 5L, 2L, 9L, 11L), .Label = c("2011Q1",
  16. "2011Q2", "2011Q3", "2011Q4", "2012Q1", "2012Q2", "2012Q3", "2012Q4",
  17. "2013Q1", "2013Q2", "2013Q3", "2013Q4"), class = "factor"), g.mdc = c("08",
  18. "05", "09", "08", "14", "15", "15", "11", "09", "01", "08", "11",
  19. "16", "14", "08", "06", "08", "06", "06", "08", "15", "14", "14",
  20. "08", "11", "09", "08", "08", "06", "06", "06", "08", "03", "05",
  21. "05", "15", "02", "05", "08", "04", "04", "10", "06", "01", "08",
  22. "05", "03", "06", "01", "01", "06", "08", "08", "04", "12", "05",
  23. "01", "15", "08", "01", "08", "01", "05", "15", "15", "01", "06",
  24. "15", "01", "08", "01", "05", "08", "02", "15", "03", "06", "05",
  25. "05", "03", "09", "08", "11", "12", "06", "04", "08", "01", "06",
  26. "01", "08", "06", "15", "05", "08", "07", "08", "13", "08", "08"
  27. ), sex = structure(c(2L, 2L, 2L, 2L, 2L, 1L, 2L, 1L, 2L, 2L,
  28. 2L, 1L, 2L, 2L, 1L, 2L, 1L, 1L, 2L, 2L, 1L, 2L, 2L, 2L, 1L, 2L,
  29. 2L, 2L, 1L, 1L, 2L, 1L, 1L, 2L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 1L,
  30. 1L, 1L, 1L, 2L, 1L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 1L, 1L,
  31. 2L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 1L, 2L, 2L, 1L,
  32. 1L, 1L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 2L, 1L, 2L,
  33. 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L), .Label = c("0", "1"), class = "factor"),
  34. age = c(23L, 83L, 51L, 54L, 37L, 0L, 0L, 82L, 45L, 88L, 84L,
  35. 58L, 41L, 33L, 71L, 79L, 67L, 42L, 73L, 66L, 0L, 26L, 38L,
  36. 65L, 31L, 87L, 38L, 38L, 77L, 44L, 54L, 74L, 38L, 70L, 44L,
  37. 0L, 78L, 65L, 56L, 85L, 70L, 83L, 89L, 46L, 39L, 34L, 5L,
  38. 85L, 18L, 5L, 41L, 73L, 18L, 41L, 75L, 77L, 36L, 0L, 84L,
  39. 83L, 58L, 93L, 83L, 0L, 0L, 2L, 49L, 0L, 55L, 46L, 40L, 81L,
  40. 60L, 51L, 0L, 22L, 78L, 69L, 75L, 65L, 31L, 15L, 79L, 87L,
  41. 72L, 78L, 48L, 16L, 81L, 63L, 84L, 17L, 0L, 60L, 60L, 74L,
  42. 44L, 44L, 53L, 71L), deceased = structure(c(1L, 1L, 1L, 1L,
  43. 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
  44. 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L,
  45. 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
  46. 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
  47. 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
  48. 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L,
  49. 1L, 1L, 1L, 1L, 1L, 1L), .Label = c("0", "1"), class = "factor"),
  50. n.group = c(3L, 2L, 1L, 1L, 1L, 3L, 1L, 1L, 1L, 1L, 2L, 1L,
  51. 1L, 1L, 3L, 2L, 3L, 1L, 3L, 1L, 1L, 1L, 1L, 3L, 1L, 1L, 1L,
  52. 2L, 1L, 3L, 2L, 1L, 1L, 2L, 2L, 1L, 1L, 1L, 3L, 1L, 2L, 1L,
  53. 1L, 2L, 2L, 2L, 1L, 3L, 3L, 1L, 3L, 3L, 3L, 1L, 1L, 1L, 1L,
  54. 1L, 3L, 1L, 3L, 3L, 2L, 1L, 3L, 3L, 1L, 3L, 1L, 3L, 2L, 2L,
  55. 2L, 1L, 2L, 1L, 3L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 3L, 2L, 1L,
  56. 1L, 1L, 1L, 3L, 1L, 2L, 1L, 3L, 1L, 2L, 1L, 2L, 2L)), class = c("tbl_df",
  57. "tbl", "data.frame"), .Names = c("hosptg", "quarter.adm", "g.mdc",
  58. "sex", "age", "deceased", "n.group"), row.names = c(NA, -100L
  59. ))
  60.  
  61. red111 <- red11 %>% group_by(hosptg, quarter.adm, g.mdc) %>%
  62. mutate(n= n()
  63. , female = mean(sex == '1', na.rm=T)
  64. , age = mean(age, na.rm=T)
  65. , m.mortf = mean(deceased == '1', na.rm=T)
  66. , w.mortf = weighted.mean(deceased == '1', n.group,
  67. na.rm=T))
  68.  
  69. red211 <- red11 %>% group_by(hosptg, quarter.adm, g.mdc) %>%
  70. summarize(n= n()
  71. , female = mean(sex == '1', na.rm=T)
  72. , age = mean(age, na.rm=T)
  73. , m.mortf = mean(deceased == '1', na.rm=T)
  74. , w.mortf = weighted.mean(deceased == '1', n.group, na.rm=T))
  75.  
  76. mean(red11$deceased == 1, na.rm=T) [1] 0.02
  77.  
  78. sum(red211$n) [1] 170
  79.  
  80. > mean(red111$female) [1] 0.52
  81.  
  82. > mean(red111$w.mortf) [1] 0.02
  83.  
  84. > mean(red111$m.mortf) [1] 0.02
  85.  
  86. sum(red211$n) [1] 100
  87.  
  88. mean(red211$female) [1] 0.4977169
  89.  
  90. mean(red211$w.mortf) [1] 0.02739726
  91.  
  92. mean(red211$m.mortf) [1] 0.02739726
Add Comment
Please, Sign In to add comment
Public Pastes
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!