R/showtable.daily.R

In cyhsuTN/COVID19: Prediction of Cases Volumne and Hospitalization by SIR Model

#' Tables to show daily cases volumne and daily hospitalization.
#' @description Tables to show daily cases volumne and daily hospitalization.
#' @import stats


#' @param showdays Input. Number of days to project.
#' @param windowsize Input. Display by every \emph{windowsize} days.
#' @param casevolumne Input. Object from function \emph{fitSIR}.
#' @param dailyhosp Input. Object from function \emph{Prehos.daily}.


#' @return \item{short.dailycv}{Short table of daily cases volumne.}
#' @return \item{short.dailyhos}{short table of predicted daily hospitalization.}

#' @examples ## To predicte 100 days from today (dayFT=100).
#' @examples casevolumne <- fitSIR(susceptible=4119405, Infected=3733, inihos=14,
#' @examples      hosrate=2.5, hms=15, inidbt=4, mrt=14, sdis=30, dayFT=100)
#' @examples dailyhosp <- Prehos.daily(casevolumne, hosrate=2.5, icurate=0.75, venrate=0.5, hms=15)
#' @examples ## show the first 60 days/ every 2 days
#' @examples showtable.daily(showdays=60, windowsize=7, casevolumne, dailyhosp)
#'

#' @export

showtable.daily <- function(showdays=60, windowsize=7, casevolumne, dailyhosp) {
  maxint <- floor((showdays+1)/ windowsize) + 1
  daily.cv <- casevolumne$result
  daily.cv[,3] <- c(0,-diff(casevolumne$result[,2]))
  daily.cv[,4] <- c(0, diff(casevolumne$result[,4]))
  short.dailycv <- daily.cv[1+ windowsize*(1:maxint - 1),]
  short.dailyhos <- dailyhosp[1+ windowsize*(1:maxint - 1),]
  list(short.dailycv=short.dailycv, short.dailyhos=short.dailyhos)
}