bmFiltered <- Filter(function(bm) sum(abs(bm)) < q98, betaMatrices)
beta <- (1/length(bmFiltered)) * Reduce("+", bmFiltered)
list(p, beta, b, mr_params)
bmFiltered <- Filter(function(bm) sum(abs(bm)) < q98, betaMatrices)
beta <- (1/length(bmFiltered)) * Reduce("+", bmFiltered)
list(p, beta, b, mr_params)