From 91efcfe9c2209f267468f33fed63d5224d81927b Mon Sep 17 00:00:00 2001 From: Alexander Fischer Date: Sun, 25 Jun 2023 19:42:58 +0200 Subject: [PATCH] more on cluster fixed effects #124 --- R/boot_algo_fastnreliable.R | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/R/boot_algo_fastnreliable.R b/R/boot_algo_fastnreliable.R index 4b855b55..7a704d23 100644 --- a/R/boot_algo_fastnreliable.R +++ b/R/boot_algo_fastnreliable.R @@ -100,6 +100,8 @@ boot_algo_fastnreliable <- function( G <- N_G_bootcluster <- length(unique(bootcluster[[1]])) k <- length(R) + + bootstrap_type_x <- paste0(substr(bootstrap_type, 1, 1), "x") if(impose_null){ # WCR1x or WCR3x @@ -132,7 +134,7 @@ boot_algo_fastnreliable <- function( function(g) Matrix::t(X_list[[g]]) %*% y_list[[g]] ) - + tXX <- Reduce("+", tXgXg) # crossprod(X) tXy <- Reduce("+", tXgyg) # t(X) %*% y tXXinv <- solve(tXX) @@ -288,7 +290,7 @@ boot_algo_fastnreliable <- function( score_diff <- scores_boot - scores_g_boot[g,] delta_diff[g,] <- - ( + as.vector( (inv_tXX_tXgXg[[g]] %*% score_diff) - delta_b_star )^2 }