From c6478783ce058f9c2bcbe1a85ec31d9640645989 Mon Sep 17 00:00:00 2001 From: Alex Date: Fri, 10 Oct 2025 14:59:42 -0600 Subject: [PATCH] Working on death trend simulations --- Scripts/Data_Load.r | 23 ++++++++++++++++++++++- 1 file changed, 22 insertions(+), 1 deletion(-) diff --git a/Scripts/Data_Load.r b/Scripts/Data_Load.r index 86e262d..35d6716 100644 --- a/Scripts/Data_Load.r +++ b/Scripts/Data_Load.r @@ -269,5 +269,26 @@ GET_MORTALITY_DATA("Data/Mortality_Rates/Male/I_85+.csv","Male",85,Inf) ) ####NEXT STEP JOIN MORTALITY TO POPULATION!!!!!!!!!!!!!!! LIN_MORTALITY <- MORTALITY_DATA_ALL %>% filter(County=="Lincoln") -DEM_2020 %>% filter(County=='Lincoln') +LIN_DEM <- DEM_2020 %>% filter(County=='Lincoln') +CURRENT <- LIN_DEM[1,] +VALUES <- LIN_MORTALITY %>% filter(Sex=="Male",Min_Age<=0,Max_Age>=0) %>% select(Death_Rate,Rate_SD,Trend,Trend_SD) +VALUES +VALUES$RATE_SD +TREND_SIM <- function(NUM_YEARS_FORWARD,VALUES){ + RES <- c() + C_YEAR_VALUE <- VALUES$Death_Rate + for(i in 1:100){ + RES[length(RES)+1] <- C_YEAR_VALUE + C_YEAR_VALUE <- rnorm(1,mean=VALUES$Trend,sd=VALUES$Trend_SD)+C_YEAR_VALUE + } + return(RES) +} +TREND_SIM(100,VALUES) + + +rnorm(50,mean=VALUES$Death_Rate,sd=VALUES$Rate_SD) +YEARS_AHEAD<-10 +for(i in 0:YEARS_AHEAD) +rnorm(50,mean=VALUES$Death_Rate,sd=VALUES$Rate_SD) +