In [1]:
agestr <- read.csv('agestr.csv')
yr <- 2000
offset <- c(rep(5, 17), 15)
In [2]:
agestr
In [3]:
as1 <- agestr[agestr[1]==yr, -(1:4)][1,]/offset
as1 <- splinefun(c(cumsum(offset) - offset/2, 100), c(as1, 0), method = "natural")
as1 <- round(as1(0:99))
as1 <- as1 * agestr[agestr[1]==yr, 4][1]/sum(as1)
as2 <- agestr[agestr[1]==yr, -(1:4)][2,]/offset
as2 <- splinefun(c(cumsum(offset) - offset/2, 100), c(as2, 0), method = "natural")
as2 <- round(as2(0:99))
as2 <- as2 * agestr[agestr[1]==yr, 4][2]/sum(as2)
par(mfrow=c(2,1))
plot(as1)
lines(cumsum(offset) - offset/2, agestr[agestr[1]==yr, -(1:4)][1,]/offset)
plot(as2)
lines(cumsum(offset) - offset/2, agestr[agestr[1]==yr, -(1:4)][2,]/offset)
In [4]:
sum(as1*0:99)/sum(as1)
sum(as2*0:99)/sum(as2)
In [11]:
asnew <- rbind(data.frame(Sex="Male", Age=0:99, Weight=as1), data.frame(Sex="Female", Age=0:99, Weight=as2))
In [12]:
write.csv(asnew, paste0("AgeStr", yr, ".csv"), row.names = F)