Revert to previous x_init settings in optimParams (keeping the initial one)
[morpheus.git] / reports / printResults.R
index b61084f..abc2e9d 100644 (file)
@@ -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)