Description Usage Arguments Value Examples
If you run glmmboot on e.g. a grid of computers,
set return_coefs_instead = TRUE
for each.
Then enter them all here. Either just list them out,
or put them into one list and enter them.
1 | combine_resampled_lists(..., return_combined_list = FALSE)
|
... |
List of outputs to be combined, or just a bunch of output entries as separate unnamed arguments. |
return_combined_list |
Logical, default |
Returns the same output as bootstrap_ci
by default,
or the combined list (as if you had just run bootstrap_model once with
all resamples) if return_combined_list = TRUE
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 | data(test_data)
library(glmmTMB)
## where subj is some RE
test_model <- glmmTMB(y ~ x_var1 + (1 | subj),
data = test_data,
family = binomial)
output_list1 <- bootstrap_model(
test_model, base_data = test_data, 99, return_coefs_instead = TRUE)
output_list2 <- bootstrap_model(
test_model, base_data = test_data, 100, return_coefs_instead = TRUE)
output_list3 <- bootstrap_model(
test_model, base_data = test_data, 100, return_coefs_instead = TRUE)
combine_resampled_lists(output_list1, output_list2, output_list3)
num_blocks = 10
num_total_resamples = 299
reg_list <- list()
for(i in 1:num_blocks){
if(i < num_blocks){
block_resamples = floor((num_total_resamples + 1)/num_blocks)
} else {
block_resamples = floor((num_total_resamples + 1)/num_blocks - 1)
}
reg_list[[i]] = bootstrap_model(test_model,
base_data = test_data,
resamples = block_resamples,
return_coefs_instead = TRUE,
num_cores = 1) ## increase for parallel
}
boot_ci1 <- combine_resampled_lists(reg_list)
full_list <- combine_resampled_lists(reg_list, return_combined_list = TRUE)
boot_ci2 <- bootstrap_ci(full_list$base_coef_se,
full_list$resampled_coef_se)
identical(boot_ci1, boot_ci2)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.