Handle alpha matrix with only 1 row

AnqiWang2021 · Apr 16, 2024 · 266fe48 · 266fe48
1 parent 7d8d07a
commit 266fe48
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/R/susie_utils.R b/R/susie_utils.R
@@ -662,7 +662,7 @@ susie_get_cs_attainable <- function(res, coverage = 0.95, ethres = 20, ...) {
     return(H)
   }
   # Get attainable coverage
-  alpha_attainable <- apply(res$alpha, 2, function(x) ifelse(x == max(x), x, 0))
+  alpha_attainable <- apply(res$alpha, 2, function(x) ifelse(x == max(x), x, 0), simplify = FALSE)
 
   # Calculate the coverage attained
   coverage_attained <- rowSums(alpha_attainable)