- b <- colMeans(do.call(rbind, lapply(mr[[idx]], function(m) m[2+d,])))
- L <- length(mr[[1]])
- beta <- (1/L) * Reduce("+", lapply(mr[[idx]], function(m) m[2:(d+1),]))
+ bVects <- lapply(mr[[idx]], function(m) m[2+d,])
+ q98 <- quantile(sapply(bVects, function(bv) sum(abs(bv))), 0.98)
+ bFiltered <- Filter(function(bv) sum(abs(bv)) < q98, bVects)
+ b <- colMeans(do.call(rbind, bFiltered))
+ betaMatrices <- lapply(mr[[idx]], function(m) m[2:(d+1),])
+ q98 <- quantile(sapply(betaMatrices, function(bm) sum(abs(bm))), 0.98)
+ bmFiltered <- Filter(function(bm) sum(abs(bm)) < q98, betaMatrices)
+ beta <- (1/length(bmFiltered)) * Reduce("+", bmFiltered)