eligible.R

# library(dplyr)

### Eligible or not 

combineInfo <- function(test1, test2){
	if(!is.na(test1)) return(test1)
	return(test2)
}

year <- 365.25

c_visits <- (c_visits
	%>% mutate(ageDays=as.numeric(visitdate - dateofbirth))
)

# We can do one-way age comparisons because criteria are nested
# (we are always more likely to treat younger children)
eligible <- (c_visits 
	%>%  rowwise()
	%>% mutate(eligible =
		visitdate < as.Date("2012-01-01") &(
			whostage >=3
			| ageDays <= 365
			| (ageDays <= round(1.5*year) & cd4percent<20)
			| (ageDays <= round(3*year) &
				combineInfo(cd4percent<20, cd4<750))
			| (ageDays <= round(5*year) &
				combineInfo(cd4percent<20, cd4<350))
			| combineInfo(cd4percent<15, cd4<200)
		)
		| (visitdate >= as.Date("2012-01-01")
			& (visitdate < as.Date("2015-04-01")) & (
				whostage >= 3
				| ageDays <= 365
				| (ageDays <= round(5*year) &
						 combineInfo(cd4percent<25,cd4<750))
				| cd4<350
			)        
		),
		e3 = ifelse((elig=="Eligible"),TRUE,eligible) ,
		e2 = (e3|eligible)
	)
	%>% ungroup 
)