EpiModel / ARTnet

Network and Epidemic Parameterization with the ARTnet Study Data
10 stars 5 forks source link

Age pyramid inputs for ARTnet #45

Closed smjenness closed 1 year ago

smjenness commented 1 year ago

@karina-d-wallrafen : opening up this separate issue here as you are continuing to work on this.

I already created the age.pyramid argument for build_netstats, so we will just need to plug in your updated data and test the mechanism to update the parameters.

But one thing arose related to the age distributions I mentioned today at our team meeting: it will be necessary to reweight the upper age group (that group above the age of sexual cessation) or else all of the target stats get inconsistent.

With a uniform age distribution, we can reweight the overall size of that upper age group with any ratio like so:

age1 <- rep(1/length(15:64) * 0.995, length(15:64))
age2 <- rep(1/length(65:99) * 0.005, length(65:99))

age.pyramid <- c(age1, age2)

where 0.005 is the fraction of the total population that should be within the ages of 65 and 99. I'd like you to create a generalized weighting function for build_netstats that specifies this proportion for any age distribution (a uniform like the above, a proper age pyramid as you're working on, or anything else).