calibration.R

rm(list=ls())

library(dplyr)

## 
## Attaching package: 'dplyr'

## The following objects are masked from 'package:stats':
## 
##     filter, lag

## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union

library(ggplot2)

n.instances <- 30

dt_list <- as.list(1:n.instances)

dt_columns <- c(
  "time",
  "nintros", #deaths
  "number.of.positive.bc.agents",
  "number.of.hpos.agents",
  "number.of.hneg.agents",
  "number.of.diagnosed.cases",
  "number.of.diagnostic.referrals",
  "number.of.screening.referrals",
  "number.of.screen.completed",
  "number.of.dt.completed",
  "number.of.symptomatic",
  "number.of.navigated.agents",
  "time.with.cancer",
  "time.until.diagnosis",
  "time.until.diagnosis.navigated",
  "time.until.diagnosis.unnavigated",
  "number.of.diagnostic.referrals.at.t",
  "number.of.screening.visits.at.t"
  )

for (i in 1:n.instances){
  dt_list[[i]] <- read.table(paste0("../06_14_2020_Burnin_30_runs/", i,".data"))
}

which(unlist(lapply(dt_list, nrow) != 360))

## [1]  3  4 10 13 24

df <- bind_rows(dt_list[14:19])
colnames(df) <- dt_columns
df$source <- rep(14:19, each=360)

ggplot(df, aes(x=time, y=number.of.screening.visits.at.t, 
               colour=as.factor(source)))+
  geom_line()+
  theme_bw()

ggplot(df, aes(x=time, y=number.of.diagnostic.referrals.at.t, 
               colour=as.factor(source)))+
  geom_line()+
  theme_bw()

ggplot(df, aes(x=time, y=number.of.diagnostic.referrals.at.t/number.of.screening.visits.at.t, 
               colour=as.factor(source)))+
  geom_line()+
  theme_bw()

ggplot(df, aes(x=time, y=number.of.dt.completed/number.of.screen.completed, 
               colour=as.factor(source)))+
  geom_line()+
  theme_bw()

ggplot(df, aes(x=time, y=number.of.dt.completed, 
               colour=as.factor(source)))+
  geom_line()+
  theme_bw()

ggplot(df, aes(x=time, y=number.of.hpos.agents, 
               colour=as.factor(source)))+
  geom_line()+
  theme_bw()

ggplot(df, aes(x=time, y=number.of.hneg.agents, 
               colour=as.factor(source)))+
  geom_line()+
  theme_bw()

ggplot(df, aes(x=time, y=number.of.hneg.agents/number.of.hpos.agents, 
               colour=as.factor(source)))+
  geom_line()+
  theme_bw()

ggplot(df, aes(x=time, y=number.of.diagnosed.cases, 
               colour=as.factor(source)))+
  geom_line()+
  theme_bw()

ggplot(df, aes(x=time, y=number.of.positive.bc.agents, 
               colour=as.factor(source)))+
  geom_line()+
  theme_bw()

ggplot(df, aes(x=time, y=nintros, 
               colour=as.factor(source)))+
  geom_line()+
  theme_bw()