diff --git a/R/0509_log1p_gmodel.rda b/R/0509_log1p_gmodel.rda new file mode 100644 index 0000000..183b560 Binary files /dev/null and b/R/0509_log1p_gmodel.rda differ diff --git a/R/readmeRDDAnalysis.R b/R/readmeRDDAnalysis.R index dd32eb1..5d9ad0a 100644 --- a/R/readmeRDDAnalysis.R +++ b/R/readmeRDDAnalysis.R @@ -82,10 +82,10 @@ mean (all_actions_data$log1p_count) # 0.6426873 #all_gmodel <- glmer.nb(count ~ D * I(week_offset)+ scaled_project_age + (D * I(week_offset) | upstream_vcs_link), # control=glmerControl(optimizer="bobyqa", # optCtrl=list(maxfun=2e5)), data=all_actions_data) -all_gmodel <- glmer.nb(count ~ D * I(week_offset)+ scaled_project_age + (D * I(week_offset) | upstream_vcs_link), data=all_actions_data, nAGQ=0, control=glmerControl(optimizer="bobyqa", +all_log1p_gmodel <- glmer.nb(log1p_count ~ D * I(week_offset)+ scaled_project_age + (D * I(week_offset) | upstream_vcs_link), data=all_actions_data, nAGQ=0, control=glmerControl(optimizer="bobyqa", optCtrl=list(maxfun=1e5))) -summary(all_gmodel) -saveRDS(all_gmodel, "0509_gmodel.rda") +summary(all_log1p_gmodel) +saveRDS(all_log1p_gmodel, "0509_log1p_gmodel.rda") #readRDS(path) library(broom.mixed) test_condvals <- broom.mixed::tidy(all_gmodel, effects = "ran_vals", conf.int = TRUE)