mgaughan/24_deb_pkg_gov
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

log1p nbinom save

Browse Source
This commit is contained in:
Matthew Gaughan 2024-05-09 19:47:38 -05:00
parent 92dcf50bd6
commit 5a21dd3e88
2 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

BIN
R/0509_log1p_gmodel.rda Normal file
View File

Binary file not shown.

6
R/readmeRDDAnalysis.R
View File

@ -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), #all_gmodel <- glmer.nb(count ~ D * I(week_offset)+ scaled_project_age + (D * I(week_offset) | upstream_vcs_link),
# control=glmerControl(optimizer="bobyqa", # control=glmerControl(optimizer="bobyqa",
# optCtrl=list(maxfun=2e5)), data=all_actions_data) # 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))) optCtrl=list(maxfun=1e5)))
summary(all_gmodel) summary(all_log1p_gmodel)
saveRDS(all_gmodel, "0509_gmodel.rda") saveRDS(all_log1p_gmodel, "0509_log1p_gmodel.rda")
#readRDS(path) #readRDS(path)
library(broom.mixed) library(broom.mixed)
test_condvals <- broom.mixed::tidy(all_gmodel, effects = "ran_vals", conf.int = TRUE) test_condvals <- broom.mixed::tidy(all_gmodel, effects = "ran_vals", conf.int = TRUE)