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

# 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)
}
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	}