-
Notifications
You must be signed in to change notification settings - Fork 1
/
serology_data.R
33 lines (28 loc) · 1.63 KB
/
serology_data.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
###### DATA serology ######
story_data <- data.frame(
age_group = factor(c("0-4", "5-9","10-14", "15-19", "20-24", "25-29", "30-34", "35-39",
"40-44", "45-49", "50-54", "55-59", "60-64" , "65-69", "70-74", "75+"),
levels = c("0-4", "5-9","10-14", "15-19", "20-24", "25-29", "30-34", "35-39",
"40-44", "45-49", "50-54", "55-59", "60-64", "65-69", "70-74", "75+")),
data = c(0.007,0.038,0.027,0.03,0.077, rep(NA,11)),
type = "STORY",
lower = c(0.000,0.002,0.00,0.001,0.024,rep(NA,11)),
upper = c(0.058,0.101,0.077,0.084,0.159,rep(NA,11))
)
# story_data$data <- story_data$data /0.077
# story_data$lower <- story_data$lower /0.077
# story_data$upper <- story_data$upper /0.077
blood_data <- data.frame(
age_group = factor(c("0-4", "5-9","10-14", "15-19", "20-24", "25-29", "30-34", "35-39",
"40-44", "45-49", "50-54", "55-59", "60-64" , "65-69", "70-74", "75+"),
levels = c("0-4", "5-9","10-14", "15-19", "20-24", "25-29", "30-34", "35-39",
"40-44", "45-49", "50-54", "55-59", "60-64", "65-69", "70-74", "75+")),
data = c(rep(NA,4),0.102,0.102,0.093,0.093,0.079,0.079,0.078,0.078,0.063,0.063, rep(NA,2)),
type = "BT",
lower = c(rep(NA,4),0.08,0.08,0.072,0.072,0.059,0.059,0.058,0.058,0.043,0.043, rep(NA,2)),
upper = c(rep(NA,4),0.126,0.126,0.113,0.113,0.099,0.099,0.097,0.097,0.082,0.082, rep(NA,2))
)
# blood_data$data <- blood_data$data /0.102
# blood_data$lower <- blood_data$lower /0.102
# blood_data$upper <- blood_data$upper /0.102
sero_data <- rbind(story_data, blood_data)