X-Git-Url: https://git.auder.net/?a=blobdiff_plain;f=reports%2FprintResults.R;h=abc2e9dcb982ffbea556cafaf02d7f7011af7e83;hb=5af71d43f3f2dba21c6667939fcff88923af3b7b;hp=b61084f4e19b93f8724535136c66cee6f97ec79c;hpb=07f2d0453144f5b639dcc28f1df40a141db7da79;p=morpheus.git 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)