-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathtable1.R
52 lines (35 loc) · 1.63 KB
/
table1.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
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
library(dplyr)
db2020_res <- research %>% filter(year_firstpub == 2020)
db2021_res <- research %>% filter(year_firstpub == 2021)
db2022_res <- research %>% filter(year_firstpub == 2022)
db2020_rct <- rct %>% filter(year_firstpub == 2020)
db2021_rct <- rct %>% filter(year_firstpub == 2021)
db2022_rct <- rct %>% filter(year_firstpub == 2022)
db2020_rev <- review %>% filter(year_firstpub == 2020)
db2021_rev <- review %>% filter(year_firstpub == 2021)
db2022_rev <- review %>% filter(year_firstpub == 2022)
kable(round(epi.prev(pos = length(db2021_res$is_coi_pred[db2021_res$is_coi_pred == TRUE]),
tested = nrow(db2021_res),
se = 0.992,
sp = 0.995)$ap,
1))
kable(round(epi.prev(pos = length(db2021_res$is_fund_pred[db2021_res$is_fund_pred == TRUE]),
tested = nrow(db2020_res),
se = 0.997,
sp = 0.981)$ap,
1))
kable(round(epi.prev(pos = length(db2020_res$is_register_pred[db2020_res$is_register_pred == TRUE]),
tested = nrow(db2020_res),
se = 0.955,
sp = 0.997)$ap,
1))
kable(round(epi.prev(pos = length(db2020_res$is_open_data[db2020_res$is_open_data == TRUE]),
tested = nrow(db2020_res),
se = 0.758,
sp = 0.986)$ap,
1))
kable(round(epi.prev(pos = length(db2020_res$is_open_code[db2020_res$is_open_code == TRUE]),
tested = nrow(db2020_res),
se = 0.587,
sp = 0.997)$ap,
1))