scratchpad/sapply_map_situation.R

In meghapsimatrix/wildmeta: Cluster Wild Bootstrapping for Meta-Analysis

library(clubSandwich)
library(robumeta)
library(metafor)
library(tidyverse)

devtools::load_all()

source("R/helpers.R")
source("R/plot_wildmeta.R")


load("scratchpad/boots.Rdata")

boots

full <- rma.mv(yi = d ~ 0 + study_type + hrs + test,
                     V = V,
                     random = ~ study_type| study,
                     data = SATcoaching)

constraints <- constrain_equal(1:3, coefs = coef(full))

sapply(boots,
       FUN = get_boot_F,
       full_model = full,
       C_mat = constraints)


sapply(boots,
       FUN = get_boot_F.rma.mv,
       full_model = full,
       C_mat = constraints)

boots_1 <- boots[[5]]


get_boot_F(full,
           boots_1,
           constraints)

y_boot <- boots_1
full_model <- full
C_mat <- constraints

y_new <- rep(NA, length = nrow(full_model$X.f))
y_new[full_model$not.na] <- y_boot


boot_mod <- tryCatch(update(full_model, formula = y_new ~ .),
                     error = function(e) NA)

if(inherits(boot_mod, "rma.mv")){

  cov_mat <- clubSandwich::vcovCR(boot_mod, type = "CR1")

  res <- clubSandwich::Wald_test(boot_mod,
                                 constraints = C_mat,
                                 vcov = cov_mat,
                                 test = "Naive-F")

  res <- res$Fstat


} else{

  res <- NA
}

res