Commit | Line | Data |
---|---|---|
ad642dc6 BA |
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 | } |