melt_finemapping_results checks is alreadu melted now

RajLabMSSM · Oct 20, 2023 · aa2e176 · aa2e176
1 parent 59e4f95
commit aa2e176
Show file tree

Hide file tree

Showing 2 changed files with 9 additions and 2 deletions.
diff --git a/R/is_melted.R b/R/is_melted.R
@@ -0,0 +1,3 @@
+is_melted <- function(dat){
+    all(c("Method","PP","CS") %in% names(dat))
+}
diff --git a/R/melt_finemapping_results.R b/R/melt_finemapping_results.R
@@ -6,15 +6,19 @@
 #' @param dat Fine-mapping results data.
 #' @param verbose Print messages.
 #' 
-#' @return Melted \link[data.table]{data.table}
+#' @returns Modified \link[data.table]{data.table}
+#'  with the new columns "Method","PP", and "CS". 
 #' @export
 #' @importFrom data.table melt.data.table data.table
 #' @importFrom stats setNames
 #' @examples 
 #' finemap_melt <- echodata::melt_finemapping_results(dat = echodata::BST1)
 melt_finemapping_results <- function(dat,
                                      verbose=TRUE){
-
+    if(is_melted(dat)){
+        messager("Data is already melted. Skipping.",v=verbose)
+        return(dat)
+    }
     PP_cols <- grep("\\.PP$",colnames(dat), value = TRUE)
     CS_cols <- grep("\\.CS$",colnames(dat), value = TRUE)
     messager("Melting PP and CS from",length(CS_cols),