I am in the process of building a function for drug efficacy that depends both on age of administration and years since administration. The idea is to figure out the efficacy of a drug at the time you get it and then the decrease in efficacy as time goes on; for example, get drug at age 50, and then follow each year and see how efficacy changes. A modified version of the function is below.
What I would like to have at the end of this is a data.frame for drug efficacy with each column being the age at which the drug was administered (e.g., 50, 51, 52, ... , 85) and each row being the efficacy at each year after dosage (e.g., 1, 2, ... , 10).
I have tried a variety of things but have gotten basically nowhere. This code below is the closest I have got, but is still way off. If need be, the function can be restructured to two parts if that's easier.
Also if there is an easy way to name the columns that would be great, but that's the least of my concerns right now.
Any help would be greatly appreciated. Thank you
b2 <- 0.28
a1 <- 0.02
a2 <- 0.04
e <- exp(1)
p_tot <- 19132 + 18827 + 14505
p1 <- 19132
p2 <- 18827
p3 <- 14505
years_i <- seq(1, 10, by = 1)
age_j <- seq(50, 85, by = 1)
VEij <- function(age_j, years_i){
b1j <- ((a1 - a2*age_j)*p_tot + b2*(p1 + 2*p2 + 3*p3)) / p_tot
Efficacy <- b1j - b2*years_i
return(Efficacy)
}
EffDF <- data.frame(years_i)
for (x in age_j){
for (y in years_i){
EffDF <- cbind(EffDF, VEij(x, y))
}
}
