diff --git a/R/0510_pop_contrib_collab.rda b/R/0510_pop_contrib_collab.rda
new file mode 100644
index 0000000..e980ed3
Binary files /dev/null and b/R/0510_pop_contrib_collab.rda differ
diff --git a/R/0510_pop_contrib_contrib.rda b/R/0510_pop_contrib_contrib.rda
new file mode 100644
index 0000000..dddb034
Binary files /dev/null and b/R/0510_pop_contrib_contrib.rda differ
diff --git a/R/0510_pop_rm_collab.rda b/R/0510_pop_rm_collab.rda
new file mode 100644
index 0000000..1bb0da8
Binary files /dev/null and b/R/0510_pop_rm_collab.rda differ
diff --git a/R/0510_pop_rm_contrib.rda b/R/0510_pop_rm_contrib.rda
new file mode 100644
index 0000000..8fdd753
Binary files /dev/null and b/R/0510_pop_rm_contrib.rda differ
diff --git a/R/popRDDAnalyssis.R b/R/popRDDAnalyssis.R
index 951fd78..2dd3e74 100644
--- a/R/popRDDAnalyssis.R
+++ b/R/popRDDAnalyssis.R
@@ -36,14 +36,22 @@ contrib_pop_contrib <- expanded_contrib_data[which(expanded_contrib_data$is_coll
 #import models 
 library(lme4)
 library(optimx)
-collab_readme_model <- lmer(log1pcount ~ after_doc + (1| upstream_vcs_link), data=collab_pop_readme, REML=FALSE)
+collab_readme_model <- glmer.nb(log1pcount ~ after_doc + (after_doc| upstream_vcs_link), data=collab_pop_readme)
 summary(collab_readme_model)
-contrib_readme_model <- lmer(log1pcount ~ after_doc + ( 1| upstream_vcs_link), data=contrib_pop_readme, REML=FALSE)
+saveRDS(collab_readme_model, "0510_pop_rm_collab.rda")
+crm_residuals <- residuals(collab_readme_model)
+qqnorm(crm_residuals)
+contrib_readme_model <- glmer.nb(log1pcount ~ after_doc + (after_doc| upstream_vcs_link), data=contrib_pop_readme)
 summary(contrib_readme_model)
-collab_contrib_model <- lmer(log1pcount ~ after_doc + ( 1 | upstream_vcs_link), data=collab_pop_contrib, REML=FALSE)
+saveRDS(contrib_readme_model, "0510_pop_rm_contrib.rda")
+conrm_residuals <- residuals(contrib_readme_model)
+qqnorm(conrm_residuals)
+collab_contrib_model <- glmer.nb(log1pcount ~ after_doc + (after_doc| upstream_vcs_link), data=collab_pop_contrib)
 summary(collab_contrib_model)
-contrib_contrib_model <- lmer(log1pcount ~ after_doc + ( 1 | upstream_vcs_link), data=contrib_pop_contrib, REML=FALSE)
+saveRDS(collab_contrib_model, "0510_pop_contrib_collab.rda")
+contrib_contrib_model <- glmer.nb(log1pcount ~ after_doc + (after_doc| upstream_vcs_link), data=contrib_pop_contrib)
 summary(contrib_contrib_model)
+saveRDS(contrib_contrib_model, "0510_pop_contrib_contrib.rda")
 
 library(ggplot2)
 expanded_readme_data |>