X-Git-Url: https://git.auder.net/?p=morpheus.git;a=blobdiff_plain;f=reports%2FprintResults.R;h=abc2e9dcb982ffbea556cafaf02d7f7011af7e83;hp=b61084f4e19b93f8724535136c66cee6f97ec79c;hb=2989133a738da99b78594cb00afdd3d87916b8bd;hpb=2591fa8343c69ddb94dec5e55871d34c55eff9a3 diff --git a/reports/printResults.R b/reports/printResults.R index b61084f..abc2e9d 100644 --- a/reports/printResults.R +++ b/reports/printResults.R @@ -7,11 +7,11 @@ prms <- function(name, idx) mr[[idx]] = mr[[1]] p <- colMeans(do.call(rbind, lapply(mr[[idx]], function(m) m[1,]))) bVects <- lapply(mr[[idx]], function(m) m[2+d,]) - q98 <- quantile(sapply(bVects, function(bv) sum(abs(bv))), 0.98) + q98 <- Inf #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) + q98 <- Inf #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) list(p, beta, b, mr_params)