-
Notifications
You must be signed in to change notification settings - Fork 2
/
README.Rmd
126 lines (90 loc) · 5.4 KB
/
README.Rmd
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
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
---
output: github_document
bibliography: vignettes/references.bib
link-citations: true
editor_options:
markdown:
wrap: 72
---
```{r, include = FALSE}
knitr::opts_chunk$set(
collapse = TRUE,
comment = "#>",
fig.path = "man/figures/README-",
out.width = "100%"
)
```
## *{{ packagename }}*: Tools for estimating vaccine effectiveness and vaccine related metrics <img src="man/figures/logo.png" align="right" width="130"/>
<!-- badges: start -->
[![License: MIT](https://img.shields.io/badge/License-MIT-blue.svg)](https://opensource.org/license/mit)
[![R-CMD-check](https://github.com/{{ gh_repo }}/actions/workflows/R-CMD-check.yaml/badge.svg)](https://github.com/{{ gh_repo }}/actions/workflows/R-CMD-check.yaml)
[![Codecov test coverage](https://codecov.io/gh/{{ gh_repo }}/branch/main/graph/badge.svg)](https://app.codecov.io/gh/{{ gh_repo }}?branch=main)
[![Project Status: WIP – Initial development is in progress, but there has not yet been a stable, usable release suitable for the public.](https://www.repostatus.org/badges/latest/wip.svg)](https://www.repostatus.org/#wip)
[![lifecycle-maturing](https://raw.githubusercontent.com/reconverse/reconverse.github.io/master/images/badge-maturing.svg)](https://www.reconverse.org/lifecycle.html#concept)
<!-- badges: end -->
*{{ packagename }}* is developed at [Pontificia Universidad Javeriana](https://www.javeriana.edu.co/inicio) as part of the
[Epiverse-TRACE initiative](https://data.org/initiatives/epiverse/).
*{{ packagename }}* is an R package that offers tools for estimating vaccine effectiveness (VE), using a series of epidemiological designs including cohort studies, test-negative case-control, and screening methods [@bookvaccine]. The current version of the package provides a set of features for preparing, visualizing, and managing cohort data, estimating vaccine effectiveness, and assessing the performance of the models. Test-negative design and screening method will be included in future versions.
## Installation
Our stable versions are released on CRAN, and can be installed using:
```r
install.packages("vaccineff", build_vignettes = TRUE)
```
The current development version of _{{ packagename }}_ can be installed from [GitHub](https://github.com/) using the `pak` package.
```{r, eval = FALSE}
if(!require("pak")) install.packages("pak")
pak::pak("{{ gh_repo }}")
```
Or using the `remotes` package
```{r, eval = FALSE }
if(!require("remotes")) install.packages("remotes")
remotes::install_github("{{ gh_repo }}")
```
## Quick start
```{r cleaning, include = FALSE, echo = TRUE}
library(vaccineff)
```
*{{ packagename }}* provides a minimal cohort dataset that can be used to test out the models.
```{r example, include = TRUE, echo = TRUE}
# Load example `cohortdata` included in the package
data("cohortdata")
head(cohortdata, 5)
```
The function `make_vaccineff_data` allows defining different aspects of the study design—such as vaccination dates, immunization delays, and potential confounding factors—and creates an object of class `vaccineff_data`. Its output is used to estimate VE using a Cox model regression invoked by the function `estimate_vaccineff`, which returns the object `vaccineff`.
The proportional hazard assumption can be tested both through the Schoenfeld test and visually using the `plot` method by setting `type = "loglog"`.
```{r}
# Create `vaccineff_data`
vaccineff_data <- make_vaccineff_data(
data_set = cohortdata,
outcome_date_col = "death_date",
censoring_date_col = "death_other_causes",
vacc_date_col = "vaccine_date_2",
vaccinated_status = "v",
unvaccinated_status = "u",
immunization_delay = 15,
end_cohort = as.Date("2021-12-31"),
match = TRUE,
exact = "sex",
nearest = c(age = 1)
)
# Estimate VE
ve <- estimate_vaccineff(vaccineff_data, at = 180)
# Print summary of VE
summary(ve)
# Generate loglog plot to check proportional hazards
plot(ve, type = "loglog")
```
## Package vignettes
More details on how to use _{{ packagename }}_ can be found in the [online documentation as package vignettes](https://epiverse-trace.github.io/{{ packagename }}/), and in the articles "Get Started with vaccineff" and "Introduction to cohort design with vaccineff".
## Help
To report a bug or to request a new feature please open an [issue](https://github.com/{{ gh_repo }}/issues/new/choose).
## Contribute
Contributions to _{{ packagename }}_ are welcomed. Contributions are welcome via [pull
requests](https://github.com/{{ gh_repo }}/pulls).
Contributors to the project include:
**Authors**: [David Santiago Quevedo](https://github.com/davidsantiagoquevedo) and [Zulma M. Cucunubá](https://github.com/zmcucunuba) (maintainer)
**Contributors**:
[Geraldine Gómez](https://github.com/GeraldineGomez), [Pratik Gupte](https://github.com/pratikunterwegs), [Érika J. Cantor](https://github.com/ErikaCantor), [Santiago Loaiza](https://github.com/santilo9513), [Jaime A. Pavlich-Mariscal](https://github.com/jpavlich), [Hugo Gruson](https://github.com/Bisaloo), [Chris Hartgerink](https://github.com/chartgerink), [Felipe Segundo Abril-Bermúdez](https://github.com/fsabrilb), [Joshua W. Lambert](https://github.com/joshwlambert), [Julian Otero](https://github.com/jd-otero)
## Code of conduct
Please note that the _{{ packagename }}_ project is released with a [Contributor Code of Conduct](https://github.com/epiverse-trace/.github/blob/main/CODE_OF_CONDUCT.md). By contributing to this project, you agree to abide by its terms.
## References