inst/doc/conversions.R

In incidence: Compute, Handle, Plot and Model Incidence of Dated Events

## ---- echo = FALSE------------------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>", 
  fig.width=7, 
  fig.height=5 
)

## ----example------------------------------------------------------------------
library(outbreaks)
library(incidence)
dat <- ebola_sim$linelist$date_of_onset
i_14 <- incidence(dat, interval = "2 epiweeks", groups = ebola_sim$linelist$gender)
i_14
plot(i_14, border = "white")

## -----------------------------------------------------------------------------
as.data.frame(i_14)

## ---- long--------------------------------------------------------------------
df <- as.data.frame(i_14, long = TRUE)
head(df)
tail(df)

## example of custom plot using steps:
library(ggplot2)
ggplot(df, aes(x = dates, y = counts)) + geom_step(aes(color = groups))

## ---- iso---------------------------------------------------------------------
i_7 <- incidence(dat, interval = "week")
i_7
plot(i_7, border = "white")
head(as.data.frame(i_7))
tail(as.data.frame(i_7))


## ---- conversions-------------------------------------------------------------
args(incidence:::as.incidence.matrix)

## -----------------------------------------------------------------------------
vec <- c(1,2,3,0,3,2,4,1,2,1)

i <- as.incidence(vec)
i

plot(vec, type = "s")
plot(i, border = "white")


## -----------------------------------------------------------------------------
i <- as.incidence(vec, interval = 7)
i
plot(i, border = "white")


## -----------------------------------------------------------------------------
i$dates

## ---- round_trip--------------------------------------------------------------
## convertion: incidence --> data.frame:
i_14
df <- as.data.frame(i_14)
head(df)
tail(df)

## conversion: data.frame --> incidence
new_i <- as.incidence(df[group_names(i_14)], df$dates, interval = "2 epiweeks")
new_i

## check round trip
identical(new_i, i_14)