diff --git a/R/contribRDDAnalysis.R b/R/contribRDDAnalysis.R
index 8c43fe2..3b1e0f4 100644
--- a/R/contribRDDAnalysis.R
+++ b/R/contribRDDAnalysis.R
@@ -58,9 +58,10 @@ library(optimx)
 library(lattice)
 #model
 print("fitting model")
-all_gmodel <- glmer.nb(log1p_count ~ D * week_offset + scaled_project_age + scaled_event_gap + (D * week_offset | upstream_vcs_link),
-              control=glmerControl(optimizer="bobyqa",
-                                           optCtrl=list(maxfun=2e5)), nAGQ=0, data=all_actions_data)
+#all_gmodel <- glmer.nb(log1p_count ~ D * week_offset + scaled_project_age + scaled_event_gap + (D * week_offset | upstream_vcs_link),
+#              control=glmerControl(optimizer="bobyqa",
+#                                           optCtrl=list(maxfun=2e5)), nAGQ=0, data=all_actions_data)
+all_gmodel <- readRDS("0512_contrib_all.rda")
 summary(all_gmodel)
 saveRDS(all_gmodel, "0512_contrib_all.rda")
 all_residuals <- residuals(all_gmodel)
diff --git a/R/final_models/0623_pop_rm_collab.rda b/R/final_models/0623_pop_rm_collab.rda
deleted file mode 100644
index 2f87015..0000000
Binary files a/R/final_models/0623_pop_rm_collab.rda and /dev/null differ
diff --git a/R/final_models/0623_pop_rm_collab_better.rda b/R/final_models/0623_pop_rm_collab_better.rda
deleted file mode 100644
index 3fa0f96..0000000
Binary files a/R/final_models/0623_pop_rm_collab_better.rda and /dev/null differ
diff --git a/R/final_models/0623_pop_rm_contrib.rda b/R/final_models/0623_pop_rm_contrib.rda
deleted file mode 100644
index ae62c72..0000000
Binary files a/R/final_models/0623_pop_rm_contrib.rda and /dev/null differ
diff --git a/R/readmeRDDAnalysis.R b/R/readmeRDDAnalysis.R
index 9ba7f49..34c2c62 100644
--- a/R/readmeRDDAnalysis.R
+++ b/R/readmeRDDAnalysis.R
@@ -38,6 +38,7 @@ windowed_data <- expanded_data |>
   mutate(D = ifelse(week > 27, 1, 0))
 #scale the age numbers 
 windowed_data$scaled_project_age <- scale(windowed_data$age_in_days)
+windowed_data$scaled_event_gap <- scale(event_gap)
 windowed_data$week_offset <- windowed_data$week - 27
 #break out the different types of commit actions that are studied
 all_actions_data <- windowed_data[which(windowed_data$observation_type == "all"),]
@@ -61,11 +62,11 @@ var(all_actions_data$count) # 268.4449
 mean (all_actions_data$count) # 3.757298
 median(all_actions_data$count) # 0
 print("fitting model")
-all_log1p_gmodel <- glmer.nb(log1p_count ~ D * week_offset+ scaled_project_age + (D * week_offset | upstream_vcs_link), data=all_actions_data, nAGQ=1,  control=glmerControl(optimizer="bobyqa",
-                           optCtrl=list(maxfun=1e5)))
-#all_log1p_gmodel <- readRDS("final_models/0510_rm_all.rda")
+all_log1p_gmodel <- glmer.nb(log1p_count ~ D * week_offset+ scaled_project_age + scaled_event_gap + (D * week_offset | upstream_vcs_link), data=all_actions_data, nAGQ=1,  control=glmerControl(optimizer="bobyqa",
+                         optCtrl=list(maxfun=1e5)))
+#all_log1p_gmodel <- readRDS("0624_log1p_nagq_gmodel_backup.rda")
 summary(all_log1p_gmodel)
-saveRDS(all_log1p_gmodel, "0624_log1p_nagq_gmodel_backup.rda")
+saveRDS(all_log1p_gmodel, "0624_eventgap_rm_rdd.rda")
 print("model fit")
 #I grouped the ranef D effects on 0624
 all_residuals <- residuals(all_log1p_gmodel)
@@ -84,7 +85,6 @@ g <- test_glmer_ranef_D |>
   geom_linerange(aes(ymin= conf.low, ymax= conf.high)) +
   theme_bw()
 g
-write.csv(test_glmer_ranef_D, "062424_readme_grouped.csv")
+write.csv(test_glmer_ranef_D, "062424_readme_grouped_1.csv")
 ggsave("0624caterpillar.png", g)
-# NOTE: below is the merge model for the same analysis, but it won't converge
 print("all pau")
\ No newline at end of file