major folder reorganisation, R pkg is now epclust/ at first level. Experimental usage...

[epclust.git] / old_C_code / stage2 / src / unused / aux2.r

diff --git a/old_C_code/stage2/src/unused/aux2.r b/old_C_code/stage2/src/unused/aux2.r
new file mode 100644 (file)
index 0000000..d7bdb79
--- /dev/null
+++ b/old_C_code/stage2/src/unused/aux2.r
@@ -0,0 +1,25 @@
+# File: aux2.r
+
+dayimpact <- function(vec) {
+  num <- tapply(vec, rep(1:365, each = 48), mean, na.rm = TRUE) 
+  den <- tapply(vec, rep(1:365, each = 48), max, na.rm = TRUE) 
+  return(num / den)
+}
+
+nightimpact <- function(vec) { 
+  idnuits <- logical(48)
+  idnuits[c(1:14, 47:48)] <- TRUE
+  nuits <-  vec[rep(idnuits, 365)] 
+  num <- tapply(nuits, rep(1:365, each = sum(idnuits)), mean, na.rm = TRUE)
+  den <- tapply(vec, rep(1:365, each = 48), mean, na.rm = TRUE) 
+  return(num / den / 3)
+}
+
+lunchimpact <- function(vec) { 
+  idlunchs <- logical(48)
+  idlunchs[25:30] <- TRUE
+  lunchs <- vec[rep(idlunchs, 365)]
+  num <- tapply(lunchs, rep(1:365, each = sum(idlunchs)), mean, na.rm = TRUE)
+  den <- tapply(vec, rep(1:365, each = 48), mean, na.rm = TRUE) 
+  return(num / den / 8)
+}