old_C_code/stage2_UNFINISHED/src/unused/aux2.r

   1 # File: aux2.r
   2
   3 dayimpact <- function(vec) {
   4   num <- tapply(vec, rep(1:365, each = 48), mean, na.rm = TRUE)
   5   den <- tapply(vec, rep(1:365, each = 48), max, na.rm = TRUE)
   6   return(num / den)
   7 }
   8
   9 nightimpact <- function(vec) {
  10   idnuits <- logical(48)
  11   idnuits[c(1:14, 47:48)] <- TRUE
  12   nuits <-  vec[rep(idnuits, 365)]
  13   num <- tapply(nuits, rep(1:365, each = sum(idnuits)), mean, na.rm = TRUE)
  14   den <- tapply(vec, rep(1:365, each = 48), mean, na.rm = TRUE)
  15   return(num / den / 3)
  16 }
  17
  18 lunchimpact <- function(vec) {
  19   idlunchs <- logical(48)
  20   idlunchs[25:30] <- TRUE
  21   lunchs <- vec[rep(idlunchs, 365)]
  22   num <- tapply(lunchs, rep(1:365, each = sum(idlunchs)), mean, na.rm = TRUE)
  23   den <- tapply(vec, rep(1:365, each = 48), mean, na.rm = TRUE)
  24   return(num / den / 8)
  25 }